Changeset 5202 for LMDZ6/branches/cirrus/libf/phylmd/lmdz_lscp_ini.F90
- Timestamp:
- Sep 20, 2024, 12:32:04 PM (31 hours ago)
- Location:
- LMDZ6/branches/cirrus
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/cirrus
- Property svn:mergeinfo changed
-
LMDZ6/branches/cirrus/libf/phylmd/lmdz_lscp_ini.F90
r5165 r5202 67 67 !$OMP THREADPRIVATE(iflag_t_glace) 68 68 69 INTEGER, SAVE, PROTECTED :: iflag_cloudth_vert=0 ! option for determining cloud fraction and content in convective boundary layers69 INTEGER, SAVE, PROTECTED :: iflag_cloudth_vert=0 ! option for determining cloud fraction and content in convective boundary layers 70 70 !$OMP THREADPRIVATE(iflag_cloudth_vert) 71 71 72 INTEGER, SAVE, PROTECTED :: iflag_gammasat=0 ! which threshold for homogeneous nucleation below -40oC72 INTEGER, SAVE, PROTECTED :: iflag_gammasat=0 ! which threshold for homogeneous nucleation below -40oC 73 73 !$OMP THREADPRIVATE(iflag_gammasat) 74 74 75 INTEGER, SAVE, PROTECTED :: iflag_rain_incloud_vol=0 ! use of volume cloud fraction for rain autoconversion75 INTEGER, SAVE, PROTECTED :: iflag_rain_incloud_vol=0 ! use of volume cloud fraction for rain autoconversion 76 76 !$OMP THREADPRIVATE(iflag_rain_incloud_vol) 77 77 78 INTEGER, SAVE, PROTECTED :: iflag_bergeron=0 ! bergeron effect for liquid precipitation treatment78 INTEGER, SAVE, PROTECTED :: iflag_bergeron=0 ! bergeron effect for liquid precipitation treatment 79 79 !$OMP THREADPRIVATE(iflag_bergeron) 80 80 81 INTEGER, SAVE, PROTECTED :: iflag_fisrtilp_qsat=0 ! qsat adjustment (iterative) during autoconversion81 INTEGER, SAVE, PROTECTED :: iflag_fisrtilp_qsat=0 ! qsat adjustment (iterative) during autoconversion 82 82 !$OMP THREADPRIVATE(iflag_fisrtilp_qsat) 83 83 84 INTEGER, SAVE, PROTECTED :: iflag_pdf=0 ! type of subgrid scale qtot pdf84 INTEGER, SAVE, PROTECTED :: iflag_pdf=0 ! type of subgrid scale qtot pdf 85 85 !$OMP THREADPRIVATE(iflag_pdf) 86 86 87 INTEGER, SAVE, PROTECTED :: iflag_autoconversion=0 ! autoconversion option 87 INTEGER, SAVE, PROTECTED :: iflag_icefrac=0 ! which phase partitioning function to use 88 !$OMP THREADPRIVATE(iflag_icefrac) 89 90 INTEGER, SAVE, PROTECTED :: iflag_autoconversion=0 ! autoconversion option 88 91 !$OMP THREADPRIVATE(iflag_autoconversion) 89 92 90 LOGICAL, SAVE, PROTECTED :: reevap_ice=.false. ! no liquid precip for T< threshold 93 94 LOGICAL, SAVE, PROTECTED :: reevap_ice=.false. ! no liquid precip for T< threshold 91 95 !$OMP THREADPRIVATE(reevap_ice) 92 96 93 REAL, SAVE, PROTECTED :: cld_lc_lsc=2.6e-4 ! liquid autoconversion coefficient, stratiform rain97 REAL, SAVE, PROTECTED :: cld_lc_lsc=2.6e-4 ! liquid autoconversion coefficient, stratiform rain 94 98 !$OMP THREADPRIVATE(cld_lc_lsc) 95 99 … … 118 122 !$OMP THREADPRIVATE(coef_eva) 119 123 120 REAL, SAVE, PROTECTED :: coef_sub ! tuning coefficient ice precip sublimation124 REAL, SAVE, PROTECTED :: coef_sub ! tuning coefficient ice precip sublimation 121 125 !$OMP THREADPRIVATE(coef_sub) 122 126 … … 124 128 !$OMP THREADPRIVATE(expo_eva) 125 129 126 REAL, SAVE, PROTECTED :: expo_sub ! tuning coefficient ice precip sublimation130 REAL, SAVE, PROTECTED :: expo_sub ! tuning coefficient ice precip sublimation 127 131 !$OMP THREADPRIVATE(expo_sub) 128 132 … … 226 230 !$OMP THREADPRIVATE(thresh_precip_frac) 227 231 232 REAL, SAVE, PROTECTED :: tau_mixenv=100000 ! Homogeneization time of mixed phase clouds [s] 233 !$OMP THREADPRIVATE(tau_mixenv) 234 235 REAL, SAVE, PROTECTED :: capa_crystal=1. ! Sursaturation of ice part in mixed phase clouds [-] 236 !$OMP THREADPRIVATE(capa_crystal) 237 238 REAL, SAVE, PROTECTED :: lmix_mpc=1000 ! Length of turbulent zones in Mixed Phase Clouds [m] 239 !$OMP THREADPRIVATE(lmix_mpc) 240 241 REAL, SAVE, PROTECTED :: naero5=0.5 ! Number concentration of aerosol larger than 0.5 microns [scm-3] 242 !$OMP THREADPRIVATE(naero5) 243 244 REAL, SAVE, PROTECTED :: gamma_snwretro = 0. ! Proportion of snow taken into account in ice retroaction in icefrac_turb [-] 245 !$OMP THREADPRIVATE(gamma_snwretro) 246 247 REAL, SAVE, PROTECTED :: gamma_taud = 1. ! Tuning coeff for tau_dissipturb [-] 248 !$OMP THREADPRIVATE(gamma_taud) 249 228 250 REAL, SAVE, PROTECTED :: gamma_col=1. ! A COMMENTER TODO [-] 229 251 !$OMP THREADPRIVATE(gamma_col) … … 235 257 !$OMP THREADPRIVATE(gamma_rim) 236 258 237 REAL, SAVE, PROTECTED :: rho_rain=1000. ! A COMMENTER TODO[kg/m3]259 REAL, SAVE, PROTECTED :: rho_rain=1000. ! Rain density [kg/m3] 238 260 !$OMP THREADPRIVATE(rho_rain) 239 261 240 REAL, SAVE, PROTECTED :: rho_ice=920. ! A COMMENTER TODO[kg/m3]262 REAL, SAVE, PROTECTED :: rho_ice=920. ! Ice density [kg/m3] 241 263 !$OMP THREADPRIVATE(rho_ice) 242 264 243 REAL, SAVE, PROTECTED :: r_rain=500.E-6 ! A COMMENTER TODO[m]265 REAL, SAVE, PROTECTED :: r_rain=500.E-6 ! Rain droplets radius for POPRECIP [m] 244 266 !$OMP THREADPRIVATE(r_rain) 245 267 246 REAL, SAVE, PROTECTED :: r_snow=1.E-3 ! A COMMENTER TODO[m]268 REAL, SAVE, PROTECTED :: r_snow=1.E-3 ! Ice crystals radius for POPRECIP [m] 247 269 !$OMP THREADPRIVATE(r_snow) 248 270 249 REAL, SAVE, PROTECTED :: tau_auto_snow_min=100. ! A COMMENTER TODO [s]271 REAL, SAVE, PROTECTED :: tau_auto_snow_min=100. ! A COMMENTER TODO [s] 250 272 !$OMP THREADPRIVATE(tau_auto_snow_min) 251 273 … … 256 278 !$OMP THREADPRIVATE(eps) 257 279 258 REAL, SAVE, PROTECTED :: gamma_melt=1. ! A COMMENTER TODO [-]280 REAL, SAVE, PROTECTED :: gamma_melt=1. ! A COMMENTER TODO [-] 259 281 !$OMP THREADPRIVATE(gamma_melt) 260 282 261 REAL, SAVE, PROTECTED :: alpha_freez=4. ! A COMMENTER TODO [-]283 REAL, SAVE, PROTECTED :: alpha_freez=4. ! A COMMENTER TODO [-] 262 284 !$OMP THREADPRIVATE(alpha_freez) 263 285 264 REAL, SAVE, PROTECTED :: beta_freez=0.1 ! A COMMENTER TODO [m-3.s-1]286 REAL, SAVE, PROTECTED :: beta_freez=0.1 ! A COMMENTER TODO [m-3.s-1] 265 287 !$OMP THREADPRIVATE(beta_freez) 266 288 267 REAL, SAVE, PROTECTED :: gamma_freez=1. ! A COMMENTER TODO [-]289 REAL, SAVE, PROTECTED :: gamma_freez=1. ! A COMMENTER TODO [-] 268 290 !$OMP THREADPRIVATE(gamma_freez) 269 291 270 REAL, SAVE, PROTECTED :: rain_fallspeed=4. ! A COMMENTER TODO [m/s]292 REAL, SAVE, PROTECTED :: rain_fallspeed=4. ! A COMMENTER TODO [m/s] 271 293 !$OMP THREADPRIVATE(rain_fallspeed) 272 294 273 REAL, SAVE, PROTECTED :: rain_fallspeed_clr ! A COMMENTER TODO [m/s]295 REAL, SAVE, PROTECTED :: rain_fallspeed_clr ! A COMMENTER TODO [m/s] 274 296 !$OMP THREADPRIVATE(rain_fallspeed_clr) 275 297 276 REAL, SAVE, PROTECTED :: rain_fallspeed_cld ! A COMMENTER TODO [m/s]298 REAL, SAVE, PROTECTED :: rain_fallspeed_cld ! A COMMENTER TODO [m/s] 277 299 !$OMP THREADPRIVATE(rain_fallspeed_cld) 278 300 279 REAL, SAVE, PROTECTED :: snow_fallspeed=1. ! A COMMENTER TODO [m/s]301 REAL, SAVE, PROTECTED :: snow_fallspeed=1. ! A COMMENTER TODO [m/s] 280 302 !$OMP THREADPRIVATE(snow_fallspeed) 281 303 282 REAL, SAVE, PROTECTED :: snow_fallspeed_clr ! A COMMENTER TODO [m/s]304 REAL, SAVE, PROTECTED :: snow_fallspeed_clr ! A COMMENTER TODO [m/s] 283 305 !$OMP THREADPRIVATE(snow_fallspeed_clr) 284 306 285 REAL, SAVE, PROTECTED :: snow_fallspeed_cld ! A COMMENTER TODO [m/s]307 REAL, SAVE, PROTECTED :: snow_fallspeed_cld ! A COMMENTER TODO [m/s] 286 308 !$OMP THREADPRIVATE(snow_fallspeed_cld) 287 309 !--End of the parameters for poprecip … … 325 347 RLMLT=RLMLT_in 326 348 RTT=RTT_in 327 R G=RG_in349 RV=RV_in 328 350 RVTMP2=RVTMP2_in 329 351 RPI=RPI_in … … 347 369 CALL getin_p('iflag_fisrtilp_qsat',iflag_fisrtilp_qsat) 348 370 CALL getin_p('iflag_pdf',iflag_pdf) 371 CALL getin_p('iflag_icefrac',iflag_icefrac) 349 372 CALL getin_p('reevap_ice',reevap_ice) 350 373 CALL getin_p('cld_lc_lsc',cld_lc_lsc) … … 368 391 CALL getin_p('dist_liq',dist_liq) 369 392 CALL getin_p('tresh_cl',tresh_cl) 393 CALL getin_p('tau_mixenv',tau_mixenv) 394 CALL getin_p('capa_crystal',capa_crystal) 395 CALL getin_p('lmix_mpc',lmix_mpc) 396 CALL getin_p('naero5',naero5) 397 CALL getin_p('gamma_snwretro',gamma_snwretro) 398 CALL getin_p('gamma_taud',gamma_taud) 370 399 CALL getin_p('iflag_oldbug_fisrtilp',iflag_oldbug_fisrtilp) 371 400 CALL getin_p('temp_nowater',temp_nowater) … … 430 459 WRITE(lunout,*) 'lscp_ini, iflag_fisrtilp_qsat:', iflag_fisrtilp_qsat 431 460 WRITE(lunout,*) 'lscp_ini, iflag_pdf', iflag_pdf 461 WRITE(lunout,*) 'lscp_ini, iflag_icefrac', iflag_icefrac 432 462 WRITE(lunout,*) 'lscp_ini, reevap_ice', reevap_ice 433 463 WRITE(lunout,*) 'lscp_ini, cld_lc_lsc', cld_lc_lsc … … 448 478 WRITE(lunout,*) 'lscp_ini, dist_liq', dist_liq 449 479 WRITE(lunout,*) 'lscp_ini, tresh_cl', tresh_cl 480 WRITE(lunout,*) 'lscp_ini, tau_mixenv', tau_mixenv 481 WRITE(lunout,*) 'lscp_ini, capa_crystal', capa_crystal 482 WRITE(lunout,*) 'lscp_ini, lmix_mpc', lmix_mpc 483 WRITE(lunout,*) 'lscp_ini, naero5', naero5 484 WRITE(lunout,*) 'lscp_ini, gamma_snwretro', gamma_snwretro 485 WRITE(lunout,*) 'lscp_ini, gamma_taud', gamma_taud 450 486 WRITE(lunout,*) 'lscp_ini, iflag_oldbug_fisrtilp', iflag_oldbug_fisrtilp 451 487 WRITE(lunout,*) 'lscp_ini, fl_cor_ebil', fl_cor_ebil
Note: See TracChangeset
for help on using the changeset viewer.