- Timestamp:
- 2021-02-08T21:55:49+01:00 (3 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2021/dev_r14318_RK3_stage1/src/OCE/DIA/diawri.F90
r14239 r14418 250 250 CALL iom_put( "sbv", z2d ) ! bottom j-current 251 251 ENDIF 252 253 IF( ln_zad_Aimp ) ww = ww + wi ! Recombine explicit and implicit parts of vertical velocity for diagnostic output 254 CALL iom_put( "woce", ww ) ! vertical velocity 255 252 253 IF( ln_zad_Aimp ) THEN ! Recombine explicit and implicit parts of vertical velocity for diagnostic output 254 CALL iom_put( "woce", ww + wi ) ! vertical velocity 255 ELSE 256 CALL iom_put( "woce", ww ) ! vertical velocity 257 ENDIF 258 256 259 IF( iom_use('w_masstr') .OR. iom_use('w_masstr2') ) THEN ! vertical mass transport & its square value 257 260 ! Caution: in the VVL case, it only correponds to the baroclinic mass transport. 258 261 z2d(:,:) = rho0 * e1e2t(:,:) 259 DO jk = 1, jpk 260 z3d(:,:,jk) = ww(:,:,jk) * z2d(:,:) 261 END DO 262 IF( ln_zad_Aimp ) THEN 263 DO jk = 1, jpk 264 z3d(:,:,jk) = ( ww(:,:,jk) + wi(:,:,jk) ) * z2d(:,:) 265 END DO 266 ELSE 267 DO jk = 1, jpk 268 z3d(:,:,jk) = ww(:,:,jk) * z2d(:,:) 269 END DO 270 ENDIF 262 271 CALL iom_put( "w_masstr" , z3d ) 263 272 IF( iom_use('w_masstr2') ) CALL iom_put( "w_masstr2", z3d(:,:,:) * z3d(:,:,:) ) 264 273 ENDIF 265 !266 IF( ln_zad_Aimp ) ww = ww - wi ! Remove implicit part of vertical velocity that was added for diagnostic output267 274 268 275 CALL iom_put( "avt" , avt ) ! T vert. eddy diff. coef.
Note: See TracChangeset
for help on using the changeset viewer.