Changeset 13295 for NEMO/trunk/src/OCE/DIA/diawri.F90
- Timestamp:
- 2020-07-10T20:24:21+02:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/OCE/DIA/diawri.F90
r13286 r13295 175 175 CALL iom_put( "sst", ts(:,:,1,jp_tem,Kmm) ) ! surface temperature 176 176 IF ( iom_use("sbt") ) THEN 177 DO_2D _11_11177 DO_2D( 1, 1, 1, 1 ) 178 178 ikbot = mbkt(ji,jj) 179 179 z2d(ji,jj) = ts(ji,jj,ikbot,jp_tem,Kmm) … … 185 185 CALL iom_put( "sss", ts(:,:,1,jp_sal,Kmm) ) ! surface salinity 186 186 IF ( iom_use("sbs") ) THEN 187 DO_2D _11_11187 DO_2D( 1, 1, 1, 1 ) 188 188 ikbot = mbkt(ji,jj) 189 189 z2d(ji,jj) = ts(ji,jj,ikbot,jp_sal,Kmm) … … 199 199 zztmp = rho0 * 0.25 200 200 z2d(:,:) = 0._wp 201 DO_2D _00_00201 DO_2D( 0, 0, 0, 0 ) 202 202 zztmp2 = ( ( rCdU_bot(ji+1,jj)+rCdU_bot(ji ,jj) ) * uu(ji ,jj,mbku(ji ,jj),Kmm) )**2 & 203 203 & + ( ( rCdU_bot(ji ,jj)+rCdU_bot(ji-1,jj) ) * uu(ji-1,jj,mbku(ji-1,jj),Kmm) )**2 & … … 214 214 CALL iom_put( "ssu", uu(:,:,1,Kmm) ) ! surface i-current 215 215 IF ( iom_use("sbu") ) THEN 216 DO_2D _11_11216 DO_2D( 1, 1, 1, 1 ) 217 217 ikbot = mbku(ji,jj) 218 218 z2d(ji,jj) = uu(ji,jj,ikbot,Kmm) … … 224 224 CALL iom_put( "ssv", vv(:,:,1,Kmm) ) ! surface j-current 225 225 IF ( iom_use("sbv") ) THEN 226 DO_2D _11_11226 DO_2D( 1, 1, 1, 1 ) 227 227 ikbot = mbkv(ji,jj) 228 228 z2d(ji,jj) = vv(ji,jj,ikbot,Kmm) … … 254 254 255 255 IF ( iom_use("sstgrad") .OR. iom_use("sstgrad2") ) THEN 256 DO_2D _00_00256 DO_2D( 0, 0, 0, 0 ) 257 257 zztmp = ts(ji,jj,1,jp_tem,Kmm) 258 258 zztmpx = ( ts(ji+1,jj,1,jp_tem,Kmm) - zztmp ) * r1_e1u(ji,jj) + ( zztmp - ts(ji-1,jj ,1,jp_tem,Kmm) ) * r1_e1u(ji-1,jj) … … 270 270 IF( iom_use("heatc") ) THEN 271 271 z2d(:,:) = 0._wp 272 DO_3D _11_11(1, jpkm1 )272 DO_3D( 1, 1, 1, 1, 1, jpkm1 ) 273 273 z2d(ji,jj) = z2d(ji,jj) + e3t(ji,jj,jk,Kmm) * ts(ji,jj,jk,jp_tem,Kmm) * tmask(ji,jj,jk) 274 274 END_3D … … 278 278 IF( iom_use("saltc") ) THEN 279 279 z2d(:,:) = 0._wp 280 DO_3D _11_11(1, jpkm1 )280 DO_3D( 1, 1, 1, 1, 1, jpkm1 ) 281 281 z2d(ji,jj) = z2d(ji,jj) + e3t(ji,jj,jk,Kmm) * ts(ji,jj,jk,jp_sal,Kmm) * tmask(ji,jj,jk) 282 282 END_3D … … 286 286 IF ( iom_use("eken") ) THEN 287 287 z3d(:,:,jpk) = 0._wp 288 DO_3D _00_00(1, jpkm1 )288 DO_3D( 0, 0, 0, 0, 1, jpkm1 ) 289 289 zztmp = 0.25_wp * r1_e1e2t(ji,jj) / e3t(ji,jj,jk,Kmm) 290 290 z3d(ji,jj,jk) = zztmp * ( uu(ji-1,jj,jk,Kmm)**2 * e2u(ji-1,jj) * e3u(ji-1,jj,jk,Kmm) & … … 312 312 IF( iom_use("u_heattr") ) THEN 313 313 z2d(:,:) = 0._wp 314 DO_3D _00_00(1, jpkm1 )314 DO_3D( 0, 0, 0, 0, 1, jpkm1 ) 315 315 z2d(ji,jj) = z2d(ji,jj) + z3d(ji,jj,jk) * ( ts(ji,jj,jk,jp_tem,Kmm) + ts(ji+1,jj,jk,jp_tem,Kmm) ) 316 316 END_3D … … 321 321 IF( iom_use("u_salttr") ) THEN 322 322 z2d(:,:) = 0.e0 323 DO_3D _00_00(1, jpkm1 )323 DO_3D( 0, 0, 0, 0, 1, jpkm1 ) 324 324 z2d(ji,jj) = z2d(ji,jj) + z3d(ji,jj,jk) * ( ts(ji,jj,jk,jp_sal,Kmm) + ts(ji+1,jj,jk,jp_sal,Kmm) ) 325 325 END_3D … … 339 339 IF( iom_use("v_heattr") ) THEN 340 340 z2d(:,:) = 0.e0 341 DO_3D _00_00(1, jpkm1 )341 DO_3D( 0, 0, 0, 0, 1, jpkm1 ) 342 342 z2d(ji,jj) = z2d(ji,jj) + z3d(ji,jj,jk) * ( ts(ji,jj,jk,jp_tem,Kmm) + ts(ji,jj+1,jk,jp_tem,Kmm) ) 343 343 END_3D … … 348 348 IF( iom_use("v_salttr") ) THEN 349 349 z2d(:,:) = 0._wp 350 DO_3D _00_00(1, jpkm1 )350 DO_3D( 0, 0, 0, 0, 1, jpkm1 ) 351 351 z2d(ji,jj) = z2d(ji,jj) + z3d(ji,jj,jk) * ( ts(ji,jj,jk,jp_sal,Kmm) + ts(ji,jj+1,jk,jp_sal,Kmm) ) 352 352 END_3D … … 357 357 IF( iom_use("tosmint") ) THEN 358 358 z2d(:,:) = 0._wp 359 DO_3D _00_00(1, jpkm1 )359 DO_3D( 0, 0, 0, 0, 1, jpkm1 ) 360 360 z2d(ji,jj) = z2d(ji,jj) + e3t(ji,jj,jk,Kmm) * ts(ji,jj,jk,jp_tem,Kmm) 361 361 END_3D … … 365 365 IF( iom_use("somint") ) THEN 366 366 z2d(:,:)=0._wp 367 DO_3D _00_00(1, jpkm1 )367 DO_3D( 0, 0, 0, 0, 1, jpkm1 ) 368 368 z2d(ji,jj) = z2d(ji,jj) + e3t(ji,jj,jk,Kmm) * ts(ji,jj,jk,jp_sal,Kmm) 369 369 END_3D
Note: See TracChangeset
for help on using the changeset viewer.