Changeset 13295 for NEMO/trunk/src/OCE/DYN/dynldf_iso.F90
- Timestamp:
- 2020-07-10T20:24:21+02:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/OCE/DYN/dynldf_iso.F90
r13237 r13295 128 128 IF( ln_dynldf_hor .AND. ln_traldf_iso ) THEN 129 129 ! 130 DO_3D _00_00(1, jpk )130 DO_3D( 0, 0, 0, 0, 1, jpk ) 131 131 uslp (ji,jj,jk) = - ( gdept(ji+1,jj,jk,Kbb) - gdept(ji ,jj ,jk,Kbb) ) * r1_e1u(ji,jj) * umask(ji,jj,jk) 132 132 vslp (ji,jj,jk) = - ( gdept(ji,jj+1,jk,Kbb) - gdept(ji ,jj ,jk,Kbb) ) * r1_e2v(ji,jj) * vmask(ji,jj,jk) … … 168 168 169 169 IF( ln_zps ) THEN ! z-coordinate - partial steps : min(e3u) 170 DO_2D _00_01170 DO_2D( 0, 0, 0, 1 ) 171 171 zabe1 = ( ahmt(ji,jj,jk)+rn_ahm_b ) * e2t(ji,jj) & 172 172 & * MIN( e3u(ji ,jj,jk,Kmm), & … … 183 183 END_2D 184 184 ELSE ! other coordinate system (zco or sco) : e3t 185 DO_2D _00_01185 DO_2D( 0, 0, 0, 1 ) 186 186 zabe1 = ( ahmt(ji,jj,jk)+rn_ahm_b ) & 187 187 & * e2t(ji,jj) * e3t(ji,jj,jk,Kmm) * r1_e1t(ji,jj) … … 199 199 200 200 ! j-flux at f-point 201 DO_2D _10_10201 DO_2D( 1, 0, 1, 0 ) 202 202 zabe2 = ( ahmf(ji,jj,jk) + rn_ahm_b ) & 203 203 & * e1f(ji,jj) * e3f(ji,jj,jk) * r1_e2f(ji,jj) … … 219 219 ! i-flux at f-point | t | 220 220 221 DO_2D _00_10221 DO_2D( 0, 0, 1, 0 ) 222 222 zabe1 = ( ahmf(ji,jj,jk) + rn_ahm_b ) & 223 223 & * e2f(ji,jj) * e3f(ji,jj,jk) * r1_e1f(ji,jj) … … 235 235 ! j-flux at t-point 236 236 IF( ln_zps ) THEN ! z-coordinate - partial steps : min(e3u) 237 DO_2D _01_10237 DO_2D( 0, 1, 1, 0 ) 238 238 zabe2 = ( ahmt(ji,jj,jk)+rn_ahm_b ) * e1t(ji,jj) & 239 239 & * MIN( e3v(ji,jj ,jk,Kmm), & … … 250 250 END_2D 251 251 ELSE ! other coordinate system (zco or sco) : e3t 252 DO_2D _01_10252 DO_2D( 0, 1, 1, 0 ) 253 253 zabe2 = ( ahmt(ji,jj,jk)+rn_ahm_b ) & 254 254 & * e1t(ji,jj) * e3t(ji,jj,jk,Kmm) * r1_e2t(ji,jj) … … 268 268 ! Second derivative (divergence) and add to the general trend 269 269 ! ----------------------------------------------------------- 270 DO_2D _00_00270 DO_2D( 0, 0, 0, 0 ) 271 271 puu(ji,jj,jk,Krhs) = puu(ji,jj,jk,Krhs) + ( ziut(ji+1,jj) - ziut(ji,jj ) & 272 272 & + zjuf(ji ,jj) - zjuf(ji,jj-1) ) * r1_e1e2u(ji,jj) &
Note: See TracChangeset
for help on using the changeset viewer.