Changeset 13295 for NEMO/trunk/tests/CANAL/MY_SRC/trazdf.F90
- Timestamp:
- 2020-07-10T20:24:21+02:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/tests/CANAL/MY_SRC/trazdf.F90
r12740 r13295 156 156 IF( l_ldfslp ) THEN ! isoneutral diffusion: add the contribution 157 157 IF( ln_traldf_msc ) THEN ! MSC iso-neutral operator 158 DO_3D _00_00(2, jpkm1 )158 DO_3D( 0, 0, 0, 0, 2, jpkm1 ) 159 159 zwt(ji,jj,jk) = zwt(ji,jj,jk) + akz(ji,jj,jk) 160 160 END_3D 161 161 ELSE ! standard or triad iso-neutral operator 162 DO_3D _00_00(2, jpkm1 )162 DO_3D( 0, 0, 0, 0, 2, jpkm1 ) 163 163 zwt(ji,jj,jk) = zwt(ji,jj,jk) + ah_wslp2(ji,jj,jk) 164 164 END_3D … … 168 168 ! Diagonal, lower (i), upper (s) (including the bottom boundary condition since avt is masked) 169 169 IF( ln_zad_Aimp ) THEN ! Adaptive implicit vertical advection 170 DO_3D _00_00(1, jpkm1 )170 DO_3D( 0, 0, 0, 0, 1, jpkm1 ) 171 171 zzwi = - p2dt * zwt(ji,jj,jk ) / e3w(ji,jj,jk ,Kmm) 172 172 zzws = - p2dt * zwt(ji,jj,jk+1) / e3w(ji,jj,jk+1,Kmm) … … 177 177 END_3D 178 178 ELSE 179 DO_3D _00_00(1, jpkm1 )179 DO_3D( 0, 0, 0, 0, 1, jpkm1 ) 180 180 zwi(ji,jj,jk) = - p2dt * zwt(ji,jj,jk ) / e3w(ji,jj,jk,Kmm) 181 181 zws(ji,jj,jk) = - p2dt * zwt(ji,jj,jk+1) / e3w(ji,jj,jk+1,Kmm) … … 203 203 ! used as a work space array: its value is modified. 204 204 ! 205 DO_2D _00_00205 DO_2D( 0, 0, 0, 0 ) 206 206 zwt(ji,jj,1) = zwd(ji,jj,1) 207 207 END_2D 208 DO_3D _00_00(2, jpkm1 )208 DO_3D( 0, 0, 0, 0, 2, jpkm1 ) 209 209 zwt(ji,jj,jk) = zwd(ji,jj,jk) - zwi(ji,jj,jk) * zws(ji,jj,jk-1) / zwt(ji,jj,jk-1) 210 210 END_3D … … 212 212 ENDIF 213 213 ! 214 DO_2D _00_00214 DO_2D( 0, 0, 0, 0 ) 215 215 pt(ji,jj,1,jn,Kaa) = e3t(ji,jj,1,Kbb) * pt(ji,jj,1,jn,Kbb) + p2dt * e3t(ji,jj,1,Kmm) * pt(ji,jj,1,jn,Krhs) 216 216 END_2D 217 DO_3D _00_00(2, jpkm1 )217 DO_3D( 0, 0, 0, 0, 2, jpkm1 ) 218 218 zrhs = e3t(ji,jj,jk,Kbb) * pt(ji,jj,jk,jn,Kbb) + p2dt * e3t(ji,jj,jk,Kmm) * pt(ji,jj,jk,jn,Krhs) ! zrhs=right hand side 219 219 pt(ji,jj,jk,jn,Kaa) = zrhs - zwi(ji,jj,jk) / zwt(ji,jj,jk-1) * pt(ji,jj,jk-1,jn,Kaa) 220 220 END_3D 221 221 ! 222 DO_2D _00_00222 DO_2D( 0, 0, 0, 0 ) 223 223 pt(ji,jj,jpkm1,jn,Kaa) = pt(ji,jj,jpkm1,jn,Kaa) / zwt(ji,jj,jpkm1) * tmask(ji,jj,jpkm1) 224 224 END_2D 225 DO_3DS _00_00(jpk-2, 1, -1 )225 DO_3DS( 0, 0, 0, 0, jpk-2, 1, -1 ) 226 226 pt(ji,jj,jk,jn,Kaa) = ( pt(ji,jj,jk,jn,Kaa) - zws(ji,jj,jk) * pt(ji,jj,jk+1,jn,Kaa) ) & 227 227 & / zwt(ji,jj,jk) * tmask(ji,jj,jk)
Note: See TracChangeset
for help on using the changeset viewer.