- Timestamp:
- 2013-08-14T10:54:56+02:00 (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90
r3998 r4006 1253 1253 DO jj = 1, jpj 1254 1254 DO ji = 1, jpi 1255 ztaper = EXP( -(gphit(ji,jj)/8._wp)**2 )1255 ztaper = EXP( -(gphit(ji,jj)/8._wp)**2._wp ) 1256 1256 hbatt(ji,jj) = rn_sbot_max * ztaper + hbatt(ji,jj) * ( 1._wp - ztaper ) 1257 1257 END DO … … 1368 1368 fsde3w(:,:,:) = gdep3w(:,:,:) 1369 1369 ! 1370 where (e3t (:,:,:).eq.0.0) e3t(:,:,:) = 1. 01371 where (e3u (:,:,:).eq.0.0) e3u(:,:,:) = 1. 01372 where (e3v (:,:,:).eq.0.0) e3v(:,:,:) = 1. 01373 where (e3f (:,:,:).eq.0.0) e3f(:,:,:) = 1. 01374 where (e3w (:,:,:).eq.0.0) e3w(:,:,:) = 1. 01375 where (e3uw (:,:,:).eq.0.0) e3uw(:,:,:) = 1. 01376 where (e3vw (:,:,:).eq.0.0) e3vw(:,:,:) = 1. 01370 where (e3t (:,:,:).eq.0.0) e3t(:,:,:) = 1._wp 1371 where (e3u (:,:,:).eq.0.0) e3u(:,:,:) = 1._wp 1372 where (e3v (:,:,:).eq.0.0) e3v(:,:,:) = 1._wp 1373 where (e3f (:,:,:).eq.0.0) e3f(:,:,:) = 1._wp 1374 where (e3w (:,:,:).eq.0.0) e3w(:,:,:) = 1._wp 1375 where (e3uw (:,:,:).eq.0.0) e3uw(:,:,:) = 1._wp 1376 where (e3vw (:,:,:).eq.0.0) e3vw(:,:,:) = 1._wp 1377 1377 1378 1378 #if defined key_agrif … … 1752 1752 ENDDO 1753 1753 ! 1754 CALL lbc_lnk(e3t ,'T',1.) ; CALL lbc_lnk(e3u ,'T',1.)1755 CALL lbc_lnk(e3v ,'T',1.) ; CALL lbc_lnk(e3f ,'T',1.)1756 CALL lbc_lnk(e3w ,'T',1.)1757 CALL lbc_lnk(e3uw,'T',1.) ; CALL lbc_lnk(e3vw,'T',1.)1758 !1759 1754 ! ! ============= 1760 1755 … … 1853 1848 !!---------------------------------------------------------------------- 1854 1849 ! 1855 pf = ( TANH( rn_theta * ( -(pk-0.5_wp) / REAL(jpkm1 ) + rn_thetb ) ) &1850 pf = ( TANH( rn_theta * ( -(pk-0.5_wp) / REAL(jpkm1,wp) + rn_thetb ) ) & 1856 1851 & - TANH( rn_thetb * rn_theta ) ) & 1857 1852 & * ( COSH( rn_theta ) & … … 1879 1874 ! 1880 1875 IF ( rn_theta == 0 ) then ! uniform sigma 1881 pf1 = - ( pk1 - 0.5_wp ) / REAL( jpkm1 )1876 pf1 = - ( pk1 - 0.5_wp ) / REAL( jpkm1,wp ) 1882 1877 ELSE ! stretched sigma 1883 pf1 = ( 1._wp - pbb ) * ( SINH( rn_theta*(-(pk1-0.5_wp)/REAL(jpkm1 )) ) ) / SINH( rn_theta ) &1884 & + pbb * ( (TANH( rn_theta*( (-(pk1-0.5_wp)/REAL(jpkm1 )) + 0.5_wp) ) - TANH( 0.5_wp * rn_theta ) ) &1878 pf1 = ( 1._wp - pbb ) * ( SINH( rn_theta*(-(pk1-0.5_wp)/REAL(jpkm1,wp)) ) ) / SINH( rn_theta ) & 1879 & + pbb * ( (TANH( rn_theta*( (-(pk1-0.5_wp)/REAL(jpkm1,wp)) + 0.5_wp) ) - TANH( 0.5_wp * rn_theta ) ) & 1885 1880 & / ( 2._wp * TANH( 0.5_wp * rn_theta ) ) ) 1886 1881 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.