Ignore:
Timestamp:
Apr 13, 2025, 7:10:19 PM (3 months ago)
Author:
aborella
Message:
  • changed treatment of prognostic variables for prognostic clouds
  • adapted sedimentation and autoconversion for prognostic cirrus clouds
  • cloud mixing, ice sedimentation and ISSR diagnosis are now consistent with the water vapor PDF
  • simplified assumptions for ice crystals deposition / sublimation
  • first version of the coupling between prognostic cirrus clouds and deep convection
  • added persistent contrail cirrus clouds in radiative diagnostics
File:
1 edited

Legend:

Unmodified
Added
Removed
  • TabularUnified LMDZ6/branches/contrails/libf/phylmd/lmdz_lscp_ini.f90

    r5602 r5609  
    162162  !$OMP THREADPRIVATE(ok_weibull_warm_clouds)
    163163
     164  REAL, SAVE, PROTECTED :: ffallv_issr                       ! tuning coefficient crystal fall velocity, cirrus clouds (with ISSR)
     165  !$OMP THREADPRIVATE(ffallv_issr)
     166
    164167  REAL, SAVE, PROTECTED :: depo_coef_cirrus=.7               ! [-] deposition coefficient for growth of ice crystals in cirrus clouds
    165168  !$OMP THREADPRIVATE(depo_coef_cirrus)
     
    168171  !$OMP THREADPRIVATE(capa_cond_cirrus)
    169172
    170   REAL, SAVE, PROTECTED :: std_subl_pdf_lscp=2.              ! [%] standard deviation of the gaussian distribution of water inside ice clouds
    171   !$OMP THREADPRIVATE(std_subl_pdf_lscp)
    172  
    173   REAL, SAVE, PROTECTED :: beta_pdf_lscp=1.E-3               ! [SI] tuning coefficient for the standard deviation of the PDF of water vapor in the clear sky region
     173  REAL, SAVE, PROTECTED :: mu_subl_pdf_lscp=1./3.            ! [-] factor for the ice distribution inside cirrus clouds
     174  !$OMP THREADPRIVATE(mu_subl_pdf_lscp)
     175 
     176  REAL, SAVE, PROTECTED :: beta_pdf_lscp=8.75E-4             ! [SI] tuning coefficient for the standard deviation of the PDF of water vapor in the clear sky region
    174177  !$OMP THREADPRIVATE(beta_pdf_lscp)
    175178 
     
    192195  !$OMP THREADPRIVATE(b_homofreez)
    193196
    194   REAL, SAVE, PROTECTED :: delta_hetfreez=1.                 ! [-] value between 0 and 1 to simulate for heterogeneous freezing.
     197  REAL, SAVE, PROTECTED :: delta_hetfreez=0.85               ! [-] value between 0 and 1 to simulate for heterogeneous freezing.
    195198  !$OMP THREADPRIVATE(delta_hetfreez)
    196199 
     
    344347  !$OMP THREADPRIVATE(ok_ice_sedim)
    345348
    346   REAL, SAVE, PROTECTED :: ice_fallspeed_sedim=0.1          ! Ice fallspeed velocity for sedimentation [m/s]
    347   !$OMP THREADPRIVATE(ice_fallspeed_sedim)
     349  REAL, SAVE, PROTECTED :: fallice_sedim=1.                 ! Tuning factor for ice fallspeed velocity for sedimentation [-]
     350  !$OMP THREADPRIVATE(fallice_sedim)
    348351  !--End of the parameters for poprecip
    349352
     
    466469    CALL getin_p('snow_fallspeed_cld',snow_fallspeed_cld)
    467470    CALL getin_p('ok_ice_sedim',ok_ice_sedim)
    468     CALL getin_p('ice_fallspeed_sedim',ice_fallspeed_sedim)
     471    CALL getin_p('fallice_sedim',fallice_sedim)
    469472    ! for condensation and ice supersaturation
    470473    CALL getin_p('ok_unadjusted_clouds',ok_unadjusted_clouds)
    471474    CALL getin_p('ok_weibull_warm_clouds',ok_weibull_warm_clouds)
     475    ffallv_issr=ffallv_lsc
     476    CALL getin_p('ffallv_issr',ffallv_issr)
    472477    CALL getin_p('depo_coef_cirrus',depo_coef_cirrus)
    473478    CALL getin_p('capa_cond_cirrus',capa_cond_cirrus)
    474     CALL getin_p('std_subl_pdf_lscp',std_subl_pdf_lscp)
     479    CALL getin_p('mu_subl_pdf_lscp',mu_subl_pdf_lscp)
    475480    CALL getin_p('beta_pdf_lscp',beta_pdf_lscp)
    476481    CALL getin_p('temp_thresh_pdf_lscp',temp_thresh_pdf_lscp)
     
    562567    WRITE(lunout,*) 'lscp_ini, snow_fallspeed_cld:', snow_fallspeed_cld
    563568    WRITE(lunout,*) 'lscp_ini, ok_ice_sedim:', ok_ice_sedim
    564     WRITE(lunout,*) 'lscp_ini, ice_fallspeed_sedim:', ice_fallspeed_sedim
     569    WRITE(lunout,*) 'lscp_ini, fallice_sedim:', fallice_sedim
    565570    ! for condensation and ice supersaturation
    566571    WRITE(lunout,*) 'lscp_ini, ok_ice_supersat:', ok_ice_supersat
    567572    WRITE(lunout,*) 'lscp_ini, ok_unadjusted_clouds:', ok_unadjusted_clouds
    568573    WRITE(lunout,*) 'lscp_ini, ok_weibull_warm_clouds:', ok_weibull_warm_clouds
     574    WRITE(lunout,*) 'lscp_ini, ffallv_issr', ffallv_issr
    569575    WRITE(lunout,*) 'lscp_ini, depo_coef_cirrus:', depo_coef_cirrus
    570576    WRITE(lunout,*) 'lscp_ini, capa_cond_cirrus:', capa_cond_cirrus
    571     WRITE(lunout,*) 'lscp_ini, std_subl_pdf_lscp:', std_subl_pdf_lscp
     577    WRITE(lunout,*) 'lscp_ini, mu_subl_pdf_lscp:', mu_subl_pdf_lscp
    572578    WRITE(lunout,*) 'lscp_ini, beta_pdf_lscp:', beta_pdf_lscp
    573579    WRITE(lunout,*) 'lscp_ini, temp_thresh_pdf_lscp:', temp_thresh_pdf_lscp
Note: See TracChangeset for help on using the changeset viewer.