Changeset 13169


Ignore:
Timestamp:
2020-06-27T19:10:19+02:00 (2 weeks ago)
Author:
clem
Message:

adapt conductivities for met-office. Here a minimum of 1mm of snow is used, but we may need to increase it to 1cm (only for met-office specific fields)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/r4.0-HEAD_r12713_clem_dan_fixcpl/src/ICE/icethd_zdf_bl99.F90

    r13148 r13169  
    133133      REAL(wp), DIMENSION(jpij)            ::   za_s_fra    ! ice fraction covered by snow  
    134134      REAL(wp), DIMENSION(jpij)            ::   isnow       ! snow presence (1) or not (0)  
     135      REAL(wp), DIMENSION(jpij)            ::   isnow_comb  ! snow presence for met-office  
    135136      ! 
    136137      ! Mono-category 
     
    182183            isnow (ji) = 0._wp 
    183184         ENDIF 
     185         ! for Met-Office 
     186         IF( h_s_1d(ji) < zh_min ) THEN 
     187            isnow_comb(ji) = h_s_1d(ji) / zh_min 
     188         ELSE 
     189            isnow_comb(ji) = 1._wp 
     190         ENDIF 
    184191      END DO 
    185192      ! clem: we should apply correction on snow thickness to take into account snow fraction 
     
    331338         DO ji = 1, npti   ! Snow-ice interface 
    332339            IF ( .NOT. l_T_converged(ji) ) THEN 
    333                ! Calculate combined surface snow and ice conductivity to pass through the coupler 
    334                zkappa_comb(ji) = isnow(ji) * zkappa_s(ji,0) + ( 1._wp - isnow(ji) ) * zkappa_i(ji,0) 
     340               ! Calculate combined surface snow and ice conductivity to pass through the coupler (met-office) 
     341               zkappa_comb(ji) = isnow_comb(ji) * zkappa_s(ji,0) + ( 1._wp - isnow_comb(ji) ) * zkappa_i(ji,0) 
    335342               ! If there is snow then use the same snow-ice interface conductivity for the top layer of ice 
    336343               IF( h_s_1d(ji) > 0._wp )   zkappa_i(ji,0) = zkappa_s(ji,nlay_s) 
     
    896903      DO ji = 1, npti 
    897904         IF( h_i_1d(ji) >= zhi_ssl ) THEN 
    898 !!clem metO            cnd_ice_1d(ji) = 2._wp * zkappa_comb(ji) 
    899             cnd_ice_1d(ji) = 2._wp * zkappa_i(ji,0) 
     905            cnd_ice_1d(ji) = 2._wp * zkappa_comb(ji) 
     906            !!cnd_ice_1d(ji) = 2._wp * zkappa_i(ji,0) 
    900907         ELSE 
    901908            cnd_ice_1d(ji) = 2._wp * ztcond_i(ji,0) / zhi_ssl ! cnd_ice is capped by: cond_i/zhi_ssl 
Note: See TracChangeset for help on using the changeset viewer.