Ignore:
Timestamp:
2020-07-09T15:19:01+02:00 (3 months ago)
Author:
ayoung
Message:

Merging changes for ticket #2407. Patch for iceberg melting stability. Fully SETTE tested: this merge introduces a change to the results in ORCA2_ICE_PISCES.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/src/OCE/ICB/icbthm.F90

    r13226 r13281  
    2020   USE phycst         ! NEMO physical constants 
    2121   USE sbc_oce 
     22   USE eosbn2         ! equation of state 
    2223   USE lib_fortran, ONLY : DDPDD 
    2324 
     
    5051      INTEGER  ::   ii, ij 
    5152      REAL(wp) ::   zM, zT, zW, zL, zSST, zVol, zLn, zWn, zTn, znVol, zIC, zDn 
     53      REAL(wp) ::   zSSS, zfzpt 
    5254      REAL(wp) ::   zMv, zMe, zMb, zmelt, zdvo, zdva, zdM, zSs, zdMe, zdMb, zdMv 
    5355      REAL(wp) ::   zMnew, zMnew1, zMnew2, zheat_hcflux, zheat_latent, z1_12 
     
    8587         CALL icb_utl_interp( pt%xi, pt%e1, pt%uo, pt%ui, pt%ua, pt%ssh_x,   & 
    8688            &                 pt%yj, pt%e2, pt%vo, pt%vi, pt%va, pt%ssh_y,   & 
    87             &                 pt%sst, pt%cn, pt%hi, zff ) 
     89            &                 pt%sst, pt%cn, pt%hi, zff, pt%sss ) 
    8890         ! 
    8991         zSST = pt%sst 
     92         zSSS = pt%sss 
     93         CALL eos_fzp(zSSS,zfzpt)                       ! freezing point 
    9094         zIC  = MIN( 1._wp, pt%cn + rn_sicn_shift )     ! Shift sea-ice concentration       !!gm ??? 
    9195         zM   = pt%mass 
     
    109113 
    110114         ! Melt rates in m/s (i.e. division by rday) 
    111          zMv = MAX( 7.62d-3*zSST+1.29d-3*(zSST**2)                    , 0._wp ) * z1_rday   ! Buoyant convection at sides (eqn M.A10) 
    112          zMb = MAX( 0.58_wp*(zdvo**0.8_wp)*(zSST+4.0_wp)/(zL**0.2_wp) , 0._wp ) * z1_rday   ! Basal turbulent melting     (eqn M.A7 ) 
    113          zMe = MAX( z1_12*(zSST+2.)*zSs*(1._wp+COS(rpi*(zIC**3)))     , 0._wp ) * z1_rday   ! Wave erosion                (eqn M.A8 ) 
     115         zMv = MAX( 7.62d-3*zSST+1.29d-3*(zSST**2)                    , 0._wp ) * z1_rday      ! Buoyant convection at sides (eqn M.A10) 
     116         IF ( zSST > zfzpt ) THEN                                                              ! Calculate basal melting only if SST above freezing point   
     117            zMb = MAX( 0.58_wp*(zdvo**0.8_wp)*(zSST+4.0_wp)/(zL**0.2_wp) , 0._wp ) * z1_rday   ! Basal turbulent melting     (eqn M.A7 ) 
     118         ELSE 
     119            zMb = 0._wp                                                                        ! No basal melting if SST below freezing point      
     120         ENDIF 
     121         zMe = MAX( z1_12*(zSST+2.)*zSs*(1._wp+COS(rpi*(zIC**3)))     , 0._wp ) * z1_rday      ! Wave erosion                (eqn M.A8 ) 
    114122 
    115123         IF( ln_operator_splitting ) THEN      ! Operator split update of volume/mass 
Note: See TracChangeset for help on using the changeset viewer.