Changeset 8768


Ignore:
Timestamp:
2017-11-21T15:12:09+01:00 (3 years ago)
Author:
clem
Message:

add effective conductivity in LIM

Location:
branches/UKMO/dev_r8183_ICEMODEL_svn_removed/NEMOGCM/NEMO/LIM_SRC_3
Files:
5 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r8183_ICEMODEL_svn_removed/NEMOGCM/NEMO/LIM_SRC_3/ice.F90

    r8752 r8768  
    313313   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)   ::   rn_amax_2d     !: maximum ice concentration 2d array 
    314314   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   ftr_ice        !: transmitted solar radiation under ice 
     315   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   cnd_ice        !: effective conductivity at the top of ice/snow (Jules coupling) [W.m-2.K-1] 
    315316 
    316317   !!---------------------------------------------------------------------- 
     
    449450      ! * Ice global state variables 
    450451      ii = ii + 1 
    451       ALLOCATE( ftr_ice(jpi,jpj,jpl) ,                                                 & 
     452      ALLOCATE( ftr_ice(jpi,jpj,jpl) ,cnd_ice(jpi,jpj,jpl),                           & 
    452453         &      h_i   (jpi,jpj,jpl) , a_i   (jpi,jpj,jpl) , v_i   (jpi,jpj,jpl) ,     & 
    453454         &      v_s    (jpi,jpj,jpl) , h_s  (jpi,jpj,jpl) , t_su  (jpi,jpj,jpl) ,     & 
  • branches/UKMO/dev_r8183_ICEMODEL_svn_removed/NEMOGCM/NEMO/LIM_SRC_3/ice1d.F90

    r8752 r8768  
    3838    
    3939   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) ::   qml_ice_1d     !: heat available for snow / ice surface melting [W/m2]  
    40    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) ::   qcn_ice_1d     !: heat available for snow / ice surface sublimation¬ [W/m2]  
    41    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) ::   qsr_ice_tr_1d  !: solar flux transmitted below the ice surface¬ [W/m2]  
     40   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) ::   qcn_ice_1d     !: heat available for snow / ice surface sublimation [W/m2]  
     41   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) ::   qsr_ice_tr_1d  !: solar flux transmitted below the ice surface [W/m2]  
     42   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) ::   cnd_ice_1d     !: conductivity at the top of ice/snow (Jules coupling) [W/K/m2] 
    4243 
    4344   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) ::   hfx_sum_1d 
     
    180181         &      qns_ice_1d(jpij)  ,   & 
    181182         &      qml_ice_1d(jpij), qcn_ice_1d(jpij) , qsr_ice_tr_1d(jpij) , & 
    182          &      t_bo_1d   (jpij) ,                                         & 
     183         &      cnd_ice_1d(jpij), t_bo_1d   (jpij) ,                       & 
    183184         &      hfx_sum_1d(jpij) , hfx_bom_1d(jpij) ,hfx_bog_1d(jpij) ,    &  
    184185         &      hfx_dif_1d(jpij) , hfx_opw_1d(jpij) ,hfx_dyn_1d(jpij) ,    & 
  • branches/UKMO/dev_r8183_ICEMODEL_svn_removed/NEMOGCM/NEMO/LIM_SRC_3/icestp.F90

    r8752 r8768  
    413413 
    414414      tau_icebfr(:,:) = 0._wp; ! landfast ice param only (clem: important to keep the init here) 
     415      cnd_ice(:,:,:) = 0._wp   ! effective conductivity at the top of ice/snow (Jules coupling) 
    415416       
    416417   END SUBROUTINE diag_set0 
  • branches/UKMO/dev_r8183_ICEMODEL_svn_removed/NEMOGCM/NEMO/LIM_SRC_3/icethd.F90

    r8752 r8768  
    527527         CALL tab_1d_2d( npti, nptidx(1:npti), qns_ice_1d(1:npti), qns_ice(:,:,kl) ) 
    528528         CALL tab_1d_2d( npti, nptidx(1:npti), ftr_ice_1d(1:npti), ftr_ice(:,:,kl) ) 
    529          ! 
     529         ! effective conductivity (for Jules coupling) 
     530         CALL tab_1d_2d( npti, nptidx(1:npti), cnd_ice_1d(1:npti), cnd_ice(:,:,kl) ) 
    530531         ! SIMIP diagnostics          
    531532         CALL tab_1d_2d( npti, nptidx(1:npti), t_si_1d      (1:npti), t_si(:,:,kl) ) 
  • branches/UKMO/dev_r8183_ICEMODEL_svn_removed/NEMOGCM/NEMO/LIM_SRC_3/icethd_zdf.F90

    r8752 r8768  
    388388            END DO 
    389389         END DO 
    390           
     390 
    391391         ! 
    392392         !----------------------------------------! 
     
    909909      ! 11) Jules coupling: reset inner snow and ice temperatures, update conduction fluxes 
    910910      !--------------------------------------------------------------------------------------- 
     911      ! effective conductivity (10cm treshold to avoid cnd_ice to be too big) 
     912      DO ji = 1, npti 
     913         IF( zh_s(ji) >= 0.1_wp ) THEN   ;   cnd_ice_1d(ji) = zkappa_s(ji,0) * 2._wp 
     914         ELSE                            ;   cnd_ice_1d(ji) = zkappa_i(ji,0) * 2._wp 
     915         ENDIF 
     916      END DO 
    911917      ! 
    912918      IF ( k_jules == np_zdf_jules_SND ) THEN   ! --- Jules coupling in "SND" mode 
Note: See TracChangeset for help on using the changeset viewer.