Changeset 12995


Ignore:
Timestamp:
2020-05-30T08:17:37+02:00 (4 months ago)
Author:
dancopsey
Message:

Do not use isnow in forming zkappa_i but still use it for the coupler

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/UKMO/NEMO_4.0.1_remove_0.1m_snow_test/src/ICE/icethd_zdf_bl99.F90

    r12834 r12995  
    124124      REAL(wp), DIMENSION(jpij,0:nlay_s)   ::   zkappa_s    ! Kappa factor in the snow 
    125125      REAL(wp), DIMENSION(jpij,0:nlay_s)   ::   zeta_s      ! Eta factor in the snow 
     126      REAL(wp), DIMENSION(jpij)            ::   zkappa_combi ! Combined snow and ice surface conductivity 
    126127      REAL(wp), DIMENSION(jpij,nlay_i+3)   ::   zindterm    ! 'Ind'ependent term 
    127128      REAL(wp), DIMENSION(jpij,nlay_i+3)   ::   zindtbis    ! Temporary 'ind'ependent term 
     
    317318         END DO 
    318319         DO ji = 1, npti   ! Snow-ice interface 
    319             IF ( .NOT. l_T_converged(ji) ) & 
    320                zkappa_i(ji,0) = zkappa_s(ji,nlay_s) * isnow(ji) + zkappa_i(ji,0) * ( 1._wp - isnow(ji) ) 
     320            IF ( .NOT. l_T_converged(ji) ) THEN 
     321               ! Calculate combined surface snow and ice conductivity to pass through the coupler 
     322               zkappa_combi(ji) = zkappa_s(ji,0) * isnow(ji) + zkappa_i(ji,0) * ( 1._wp - isnow(ji) ) 
     323 
     324               ! If there is snow then use the same snow-ice interface conductivity for the top layer of ice 
     325               IF ( h_s_1d(ji) > 0._wp ) THEN 
     326                  zkappa_i(ji,0) = zkappa_s(ji,nlay_s) 
     327               END IF 
     328 
    321329         END DO 
    322330         ! 
     
    847855      DO ji = 1, npti 
    848856         IF( h_i_1d(ji) > 0.1_wp ) THEN 
    849             cnd_ice_1d(ji) = 2._wp * zkappa_i(ji,0) 
     857            cnd_ice_1d(ji) = 2._wp * zkappa_combi(ji) 
    850858         ELSE 
    851859            cnd_ice_1d(ji) = 2._wp * ztcond_i(ji,0) * 10._wp 
Note: See TracChangeset for help on using the changeset viewer.