Changeset 11488 for NEMO/branches/2019/ENHANCE-02_ISF_nemo/src/OCE/ISF
- Timestamp:
- 2019-08-30T19:59:38+02:00 (5 years ago)
- Location:
- NEMO/branches/2019/ENHANCE-02_ISF_nemo/src/OCE/ISF
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/ENHANCE-02_ISF_nemo/src/OCE/ISF/isfcavgam.F90
r11486 r11488 111 111 ! 112 112 ! compute ustar (AD15 eq. 27) 113 zustar(:,:) = SQRT( pCd(:,:) * ( putbl(:,:) * putbl(:,:) + pvtbl(:,:) * pvtbl(:,:) + pke2 ) ) 113 zustar(:,:) = SQRT( pCd(:,:) * ( putbl(:,:) * putbl(:,:) + pvtbl(:,:) * pvtbl(:,:) + pke2 ) ) * mskisf_cav(:,:) 114 114 ! 115 115 ! Compute gammats -
NEMO/branches/2019/ENHANCE-02_ISF_nemo/src/OCE/ISF/isfcavmlt.F90
r11425 r11488 137 137 ! 138 138 ! thermal driving 139 zthd (:,:) = pttbl(:,:) - ztfrz(:,:)139 zthd (:,:) = ( pttbl(:,:) - ztfrz(:,:) ) * mskisf_cav(:,:) 140 140 ! 141 141 ! compute ocean-ice heat flux and then derive fwf assuming that ocean heat flux equal latent heat … … 211 211 ! 212 212 ! thermal driving 213 zthd(ji,jj) = pttbl(ji,jj) - ztfrz(ji,jj)213 zthd(ji,jj) = ( pttbl(ji,jj) - ztfrz(ji,jj) ) 214 214 ! 215 215 ! compute the upward water and heat flux (eq. 24 and eq. 26) … … 218 218 pqhc (ji,jj) = rcp * pqfwf(ji,jj) * ztfrz(ji,jj) ! heat content flux (> 0 out) 219 219 ! 220 zqcon(ji,jj) = zeps3 * ( ztfrz(ji,jj) - rtsurf ) ! conductive flux through the ice (> 0 out)220 zqcon(ji,jj) = zeps3 * ( ztfrz(ji,jj) - rtsurf ) 221 221 ! 222 222 END DO … … 224 224 ! 225 225 ! output conductive heat flux through the ice 226 CALL iom_put('qconisf', zqcon )226 CALL iom_put('qconisf', zqcon(:,:) * mskisf_cav(:,:) ) 227 227 ! 228 228 ! output thermal driving 229 CALL iom_put('isfthermald_cav', zthd )229 CALL iom_put('isfthermald_cav', zthd(:,:) * mskisf_cav(:,:) ) 230 230 ! 231 231 END SUBROUTINE isfcav_mlt_3eq
Note: See TracChangeset
for help on using the changeset viewer.