Changeset 5589 for LMDZ6/branches/contrails/libf/phylmd/lmdz_lscp_ini.f90
- Timestamp:
- Mar 26, 2025, 6:05:40 PM (3 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
TabularUnified LMDZ6/branches/contrails/libf/phylmd/lmdz_lscp_ini.f90 ¶
r5580 r5589 162 162 !$OMP THREADPRIVATE(ok_weibull_warm_clouds) 163 163 164 INTEGER, SAVE, PROTECTED :: iflag_cloud_sublim_pdf=4 ! iflag for the distribution of water inside ice clouds165 !$OMP THREADPRIVATE(iflag_cloud_sublim_pdf)166 167 164 REAL, SAVE, PROTECTED :: depo_coef_cirrus=.7 ! [-] deposition coefficient for growth of ice crystals in cirrus clouds 168 165 !$OMP THREADPRIVATE(depo_coef_cirrus) … … 173 170 REAL, SAVE, PROTECTED :: std_subl_pdf_lscp=2. ! [%] standard deviation of the gaussian distribution of water inside ice clouds 174 171 !$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 clouds177 !$OMP THREADPRIVATE(mu_subl_pdf_lscp)178 172 179 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 … … 260 254 !$OMP THREADPRIVATE(ok_corr_vap_evasub) 261 255 256 LOGICAL, SAVE, PROTECTED :: ok_growth_precip_deposition=.FALSE. 257 !$OMP THREADPRIVATE(ok_growth_precip_deposition) 258 262 259 REAL, SAVE, PROTECTED :: cld_lc_lsc_snow=2.e-5 ! snow autoconversion coefficient, stratiform. default from Chaboureau and PInty 2006 263 260 !$OMP THREADPRIVATE(cld_lc_lsc_snow) … … 311 308 !$OMP THREADPRIVATE(r_snow) 312 309 310 REAL, SAVE, PROTECTED :: expo_tau_auto_snow=0.1 311 !$OMP THREADPRIVATE(expo_tau_auto_snow) 312 313 313 REAL, SAVE, PROTECTED :: tau_auto_snow_min=100. ! Snow autoconversion minimal timescale (when liquid) [s] 314 314 !$OMP THREADPRIVATE(tau_auto_snow_min) … … 344 344 !$OMP THREADPRIVATE(snow_fallspeed_cld) 345 345 346 LOGICAL, SAVE, PROTECTED :: ok_ precip_resuspension=.FALSE.! Flag to activate the resuspension of ice crystals347 !$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 resuspended350 !$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) 351 351 !--End of the parameters for poprecip 352 352 … … 448 448 CALL getin_p('ok_poprecip',ok_poprecip) 449 449 CALL getin_p('ok_corr_vap_evasub',ok_corr_vap_evasub) 450 CALL getin_p('ok_growth_precip_deposition',ok_growth_precip_deposition) 450 451 CALL getin_p('rain_int_min',rain_int_min) 451 452 CALL getin_p('gamma_agg',gamma_agg) … … 467 468 CALL getin_p('snow_fallspeed_clr',snow_fallspeed_clr) 468 469 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) 471 472 ! for condensation and ice supersaturation 472 473 CALL getin_p('ok_unadjusted_clouds',ok_unadjusted_clouds) 473 474 CALL getin_p('ok_weibull_warm_clouds',ok_weibull_warm_clouds) 474 CALL getin_p('iflag_cloud_sublim_pdf',iflag_cloud_sublim_pdf)475 475 CALL getin_p('depo_coef_cirrus',depo_coef_cirrus) 476 476 CALL getin_p('capa_cond_cirrus',capa_cond_cirrus) 477 477 CALL getin_p('std_subl_pdf_lscp',std_subl_pdf_lscp) 478 CALL getin_p('mu_subl_pdf_lscp',mu_subl_pdf_lscp)479 478 CALL getin_p('beta_pdf_lscp',beta_pdf_lscp) 480 479 CALL getin_p('temp_thresh_pdf_lscp',temp_thresh_pdf_lscp) … … 553 552 WRITE(lunout,*) 'lscp_ini, ok_poprecip', ok_poprecip 554 553 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 555 555 WRITE(lunout,*) 'lscp_ini, rain_int_min:', rain_int_min 556 556 WRITE(lunout,*) 'lscp_ini, gamma_agg:', gamma_agg … … 566 566 WRITE(lunout,*) 'lscp_ini, snow_fallspeed_clr:', snow_fallspeed_clr 567 567 WRITE(lunout,*) 'lscp_ini, snow_fallspeed_cld:', snow_fallspeed_cld 568 WRITE(lunout,*) 'lscp_ini, ok_ precip_resuspension:', ok_precip_resuspension569 WRITE(lunout,*) 'lscp_ini, snow_thresh_resuspension:', snow_thresh_resuspension568 WRITE(lunout,*) 'lscp_ini, ok_ice_sedim:', ok_ice_sedim 569 WRITE(lunout,*) 'lscp_ini, ice_fallspeed_sedim:', ice_fallspeed_sedim 570 570 ! for condensation and ice supersaturation 571 571 WRITE(lunout,*) 'lscp_ini, ok_ice_supersat:', ok_ice_supersat 572 572 WRITE(lunout,*) 'lscp_ini, ok_unadjusted_clouds:', ok_unadjusted_clouds 573 573 WRITE(lunout,*) 'lscp_ini, ok_weibull_warm_clouds:', ok_weibull_warm_clouds 574 WRITE(lunout,*) 'lscp_ini, iflag_cloud_sublim_pdf:', iflag_cloud_sublim_pdf575 574 WRITE(lunout,*) 'lscp_ini, depo_coef_cirrus:', depo_coef_cirrus 576 575 WRITE(lunout,*) 'lscp_ini, capa_cond_cirrus:', capa_cond_cirrus 577 576 WRITE(lunout,*) 'lscp_ini, std_subl_pdf_lscp:', std_subl_pdf_lscp 578 WRITE(lunout,*) 'lscp_ini, mu_subl_pdf_lscp:', mu_subl_pdf_lscp579 577 WRITE(lunout,*) 'lscp_ini, beta_pdf_lscp:', beta_pdf_lscp 580 578 WRITE(lunout,*) 'lscp_ini, temp_thresh_pdf_lscp:', temp_thresh_pdf_lscp
Note: See TracChangeset
for help on using the changeset viewer.