Ignore:
Timestamp:
Mar 26, 2025, 6:05:40 PM (3 months ago)
Author:
aborella
Message:

Multiple changes:

  • added new radiative diagnostics for contrails
  • added ok_rad_contrail option to allow for a double call of RRTM (w/ and w/o contrails)
  • transformed resuspension of snow into ice sedimentation (poprecip)
  • some modifications in poprecip in line with the ones from EV
  • cleaned sublimation of ice clouds in lmdz_lscp_condensation, option ok_ice_supersat
  • aviation emissions can now be read with IOIPSL (in lon/lat mode)
File:
1 edited

Legend:

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

    r5580 r5589  
    162162  !$OMP THREADPRIVATE(ok_weibull_warm_clouds)
    163163
    164   INTEGER, SAVE, PROTECTED :: iflag_cloud_sublim_pdf=4       ! iflag for the distribution of water inside ice clouds
    165   !$OMP THREADPRIVATE(iflag_cloud_sublim_pdf)
    166 
    167164  REAL, SAVE, PROTECTED :: depo_coef_cirrus=.7               ! [-] deposition coefficient for growth of ice crystals in cirrus clouds
    168165  !$OMP THREADPRIVATE(depo_coef_cirrus)
     
    173170  REAL, SAVE, PROTECTED :: std_subl_pdf_lscp=2.              ! [%] standard deviation of the gaussian distribution of water inside ice clouds
    174171  !$OMP THREADPRIVATE(std_subl_pdf_lscp)
    175 
    176   REAL, SAVE, PROTECTED :: mu_subl_pdf_lscp=1./3.            ! [-] shape factor of the gamma distribution of water inside ice clouds
    177   !$OMP THREADPRIVATE(mu_subl_pdf_lscp)
    178172 
    179173  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
     
    260254  !$OMP THREADPRIVATE(ok_corr_vap_evasub)
    261255
     256  LOGICAL, SAVE, PROTECTED :: ok_growth_precip_deposition=.FALSE.
     257  !$OMP THREADPRIVATE(ok_growth_precip_deposition)
     258
    262259  REAL, SAVE, PROTECTED :: cld_lc_lsc_snow=2.e-5            ! snow autoconversion coefficient, stratiform. default from  Chaboureau and PInty 2006
    263260  !$OMP THREADPRIVATE(cld_lc_lsc_snow)
     
    311308  !$OMP THREADPRIVATE(r_snow)
    312309
     310  REAL, SAVE, PROTECTED :: expo_tau_auto_snow=0.1
     311  !$OMP THREADPRIVATE(expo_tau_auto_snow)
     312
    313313  REAL, SAVE, PROTECTED :: tau_auto_snow_min=100.           ! Snow autoconversion minimal timescale (when liquid) [s]
    314314  !$OMP THREADPRIVATE(tau_auto_snow_min)
     
    344344  !$OMP THREADPRIVATE(snow_fallspeed_cld)
    345345
    346   LOGICAL, SAVE, PROTECTED :: ok_precip_resuspension=.FALSE.! Flag to activate the resuspension of ice crystals
    347   !$OMP THREADPRIVATE(ok_precip_resuspension)
    348 
    349   REAL, SAVE, PROTECTED :: snow_thresh_resuspension=1.e-5   ! [kg/m2/s] Threshold flux below which part of the snow flux will be resuspended
    350   !$OMP THREADPRIVATE(snow_thresh_resuspension)
     346  LOGICAL, SAVE, PROTECTED :: ok_ice_sedim=.FALSE.          ! Flag to activate the sedimentation of ice crystals
     347  !$OMP THREADPRIVATE(ok_ice_sedim)
     348
     349  REAL, SAVE, PROTECTED :: ice_fallspeed_sedim=0.1          ! Ice fallspeed velocity for sedimentation [m/s]
     350  !$OMP THREADPRIVATE(ice_fallspeed_sedim)
    351351  !--End of the parameters for poprecip
    352352
     
    448448    CALL getin_p('ok_poprecip',ok_poprecip)
    449449    CALL getin_p('ok_corr_vap_evasub',ok_corr_vap_evasub)
     450    CALL getin_p('ok_growth_precip_deposition',ok_growth_precip_deposition)
    450451    CALL getin_p('rain_int_min',rain_int_min)
    451452    CALL getin_p('gamma_agg',gamma_agg)
     
    467468    CALL getin_p('snow_fallspeed_clr',snow_fallspeed_clr)
    468469    CALL getin_p('snow_fallspeed_cld',snow_fallspeed_cld)
    469     CALL getin_p('ok_precip_resuspension',ok_precip_resuspension)
    470     CALL getin_p('snow_thresh_resuspension',snow_thresh_resuspension)
     470    CALL getin_p('ok_ice_sedim',ok_ice_sedim)
     471    CALL getin_p('ice_fallspeed_sedim',ice_fallspeed_sedim)
    471472    ! for condensation and ice supersaturation
    472473    CALL getin_p('ok_unadjusted_clouds',ok_unadjusted_clouds)
    473474    CALL getin_p('ok_weibull_warm_clouds',ok_weibull_warm_clouds)
    474     CALL getin_p('iflag_cloud_sublim_pdf',iflag_cloud_sublim_pdf)
    475475    CALL getin_p('depo_coef_cirrus',depo_coef_cirrus)
    476476    CALL getin_p('capa_cond_cirrus',capa_cond_cirrus)
    477477    CALL getin_p('std_subl_pdf_lscp',std_subl_pdf_lscp)
    478     CALL getin_p('mu_subl_pdf_lscp',mu_subl_pdf_lscp)
    479478    CALL getin_p('beta_pdf_lscp',beta_pdf_lscp)
    480479    CALL getin_p('temp_thresh_pdf_lscp',temp_thresh_pdf_lscp)
     
    553552    WRITE(lunout,*) 'lscp_ini, ok_poprecip', ok_poprecip
    554553    WRITE(lunout,*) 'lscp_ini, ok_corr_vap_evasub', ok_corr_vap_evasub
     554    WRITE(lunout,*) 'lscp_ini, ok_growth_precip_deposition', ok_growth_precip_deposition
    555555    WRITE(lunout,*) 'lscp_ini, rain_int_min:', rain_int_min
    556556    WRITE(lunout,*) 'lscp_ini, gamma_agg:', gamma_agg
     
    566566    WRITE(lunout,*) 'lscp_ini, snow_fallspeed_clr:', snow_fallspeed_clr
    567567    WRITE(lunout,*) 'lscp_ini, snow_fallspeed_cld:', snow_fallspeed_cld
    568     WRITE(lunout,*) 'lscp_ini, ok_precip_resuspension:', ok_precip_resuspension
    569     WRITE(lunout,*) 'lscp_ini, snow_thresh_resuspension:', snow_thresh_resuspension
     568    WRITE(lunout,*) 'lscp_ini, ok_ice_sedim:', ok_ice_sedim
     569    WRITE(lunout,*) 'lscp_ini, ice_fallspeed_sedim:', ice_fallspeed_sedim
    570570    ! for condensation and ice supersaturation
    571571    WRITE(lunout,*) 'lscp_ini, ok_ice_supersat:', ok_ice_supersat
    572572    WRITE(lunout,*) 'lscp_ini, ok_unadjusted_clouds:', ok_unadjusted_clouds
    573573    WRITE(lunout,*) 'lscp_ini, ok_weibull_warm_clouds:', ok_weibull_warm_clouds
    574     WRITE(lunout,*) 'lscp_ini, iflag_cloud_sublim_pdf:', iflag_cloud_sublim_pdf
    575574    WRITE(lunout,*) 'lscp_ini, depo_coef_cirrus:', depo_coef_cirrus
    576575    WRITE(lunout,*) 'lscp_ini, capa_cond_cirrus:', capa_cond_cirrus
    577576    WRITE(lunout,*) 'lscp_ini, std_subl_pdf_lscp:', std_subl_pdf_lscp
    578     WRITE(lunout,*) 'lscp_ini, mu_subl_pdf_lscp:', mu_subl_pdf_lscp
    579577    WRITE(lunout,*) 'lscp_ini, beta_pdf_lscp:', beta_pdf_lscp
    580578    WRITE(lunout,*) 'lscp_ini, temp_thresh_pdf_lscp:', temp_thresh_pdf_lscp
Note: See TracChangeset for help on using the changeset viewer.