Changeset 14475 for NEMO/trunk/src/OCE/DYN
- Timestamp:
- 2021-02-17T17:19:28+01:00 (3 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/OCE/DYN/dynatf_qco.F90
r14224 r14475 43 43 USE isf_oce , ONLY: ln_isf ! ice shelf 44 44 USE isfdynatf , ONLY: isf_dynatf ! ice shelf volume filter correction subroutine 45 USE zdfdrg , ONLY: ln_drgice_imp, rCdU_top 45 46 ! 46 47 USE in_out_manager ! I/O manager … … 101 102 REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: zue, zve 102 103 REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) :: zua, zva 104 REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: zutau, zvtau 103 105 !!---------------------------------------------------------------------- 104 106 ! … … 239 241 ENDIF 240 242 ! 243 IF ( iom_use("utau") ) THEN 244 IF ( ln_drgice_imp.OR.ln_isfcav ) THEN 245 ALLOCATE(zutau(jpi,jpj)) 246 DO_2D( 0, 0, 0, 0 ) 247 jk = miku(ji,jj) 248 zutau(ji,jj) = utau(ji,jj) + 0.5_wp * rho0 * ( rCdU_top(ji+1,jj)+rCdU_top(ji,jj) ) * puu(ji,jj,jk,Kaa) 249 END_2D 250 CALL iom_put( "utau", zutau(:,:) ) 251 DEALLOCATE(zutau) 252 ELSE 253 CALL iom_put( "utau", utau(:,:) ) 254 ENDIF 255 ENDIF 256 ! 257 IF ( iom_use("vtau") ) THEN 258 IF ( ln_drgice_imp.OR.ln_isfcav ) THEN 259 ALLOCATE(zvtau(jpi,jpj)) 260 DO_2D( 0, 0, 0, 0 ) 261 jk = mikv(ji,jj) 262 zvtau(ji,jj) = vtau(ji,jj) + 0.5_wp * rho0 * ( rCdU_top(ji,jj+1)+rCdU_top(ji,jj) ) * pvv(ji,jj,jk,Kaa) 263 END_2D 264 CALL iom_put( "vtau", zvtau(:,:) ) 265 DEALLOCATE(zvtau) 266 ELSE 267 CALL iom_put( "vtau", vtau(:,:) ) 268 ENDIF 269 ENDIF 270 ! 241 271 IF(sn_cfctl%l_prtctl) CALL prt_ctl( tab3d_1=puu(:,:,:,Kaa), clinfo1=' nxt - puu(:,:,:,Kaa): ', mask1=umask, & 242 272 & tab3d_2=pvv(:,:,:,Kaa), clinfo2=' pvv(:,:,:,Kaa): ' , mask2=vmask )
Note: See TracChangeset
for help on using the changeset viewer.