New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 11488 for NEMO/branches/2019/ENHANCE-02_ISF_nemo/src/OCE/ISF – NEMO

Ignore:
Timestamp:
2019-08-30T19:59:38+02:00 (5 years ago)
Author:
mathiot
Message:

ENHANCE-02_ISF: mask isf output in the open ocean + activation of bdy and tide in WED025 (ticket #2142)

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  
    111111      ! 
    112112      ! 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(:,:) 
    114114      ! 
    115115      ! Compute gammats 
  • NEMO/branches/2019/ENHANCE-02_ISF_nemo/src/OCE/ISF/isfcavmlt.F90

    r11425 r11488  
    137137      ! 
    138138      ! thermal driving 
    139       zthd (:,:) = pttbl(:,:) - ztfrz(:,:) 
     139      zthd (:,:) = ( pttbl(:,:) - ztfrz(:,:) ) * mskisf_cav(:,:) 
    140140      ! 
    141141      ! compute ocean-ice heat flux and then derive fwf assuming that ocean heat flux equal latent heat 
     
    211211            ! 
    212212            ! thermal driving 
    213             zthd(ji,jj) = pttbl(ji,jj) - ztfrz(ji,jj) 
     213            zthd(ji,jj) = ( pttbl(ji,jj) - ztfrz(ji,jj) ) 
    214214            ! 
    215215            ! compute the upward water and heat flux (eq. 24 and eq. 26) 
     
    218218            pqhc (ji,jj) = rcp      * pqfwf(ji,jj) * ztfrz(ji,jj)                             ! heat content   flux (> 0 out) 
    219219            ! 
    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 ) 
    221221            ! 
    222222         END DO 
     
    224224      ! 
    225225      ! output conductive heat flux through the ice 
    226       CALL iom_put('qconisf', zqcon) 
     226      CALL iom_put('qconisf', zqcon(:,:) * mskisf_cav(:,:) ) 
    227227      ! 
    228228      ! output thermal driving 
    229       CALL iom_put('isfthermald_cav', zthd) 
     229      CALL iom_put('isfthermald_cav', zthd(:,:) * mskisf_cav(:,:) ) 
    230230      ! 
    231231   END SUBROUTINE isfcav_mlt_3eq 
Note: See TracChangeset for help on using the changeset viewer.