Changeset 5383 for LMDZ6/trunk/libf/phylmd/lmdz_lscp_ini.f90
- Timestamp:
- Dec 5, 2024, 11:09:43 AM (8 days ago)
- File:
-
- 1 edited
-
LMDZ6/trunk/libf/phylmd/lmdz_lscp_ini.f90 (modified) (6 diffs)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmd/lmdz_lscp_ini.f90
r5268 r5383 222 222 !$OMP THREADPRIVATE(rain_int_min) 223 223 224 REAL, SAVE, PROTECTED :: thresh_precip_frac=1.E-6 ! precipitation fraction threshold TODO[-]224 REAL, SAVE, PROTECTED :: thresh_precip_frac=1.E-6 ! precipitation fraction threshold [-] 225 225 !$OMP THREADPRIVATE(thresh_precip_frac) 226 226 227 REAL, SAVE, PROTECTED :: tau_mixenv=100000 ! Homogeneization time of mixed phase clouds [s] 228 !$OMP THREADPRIVATE(tau_mixenv) 229 230 REAL, SAVE, PROTECTED :: capa_crystal=1. ! Sursaturation of ice part in mixed phase clouds [-] 227 REAL, SAVE, PROTECTED :: capa_crystal=1. ! Crystal capacitance (shape factor) for lscp_icefrac_turb [-] 231 228 !$OMP THREADPRIVATE(capa_crystal) 232 233 REAL, SAVE, PROTECTED :: lmix_mpc=1000 ! Length of turbulent zones in Mixed Phase Clouds [m]234 !$OMP THREADPRIVATE(lmix_mpc)235 229 236 230 REAL, SAVE, PROTECTED :: naero5=0.5 ! Number concentration of aerosol larger than 0.5 microns [scm-3] 237 231 !$OMP THREADPRIVATE(naero5) 238 232 239 REAL, SAVE, PROTECTED :: gamma_snwretro = 0. ! Proportion of snow taken into account in ice retroaction in icefrac_turb [-]233 REAL, SAVE, PROTECTED :: gamma_snwretro = 0. ! Proportion of snow taken into account in ice retroaction in lscp_icefrac_turb [-] 240 234 !$OMP THREADPRIVATE(gamma_snwretro) 241 235 242 REAL, SAVE, PROTECTED :: gamma_taud = 1. ! Tuning coeff for tau_dissipturb [-]236 REAL, SAVE, PROTECTED :: gamma_taud = 1. ! Tuning coeff for Lagrangian decorrelation timescale in lscp_icefrac_turb [-] 243 237 !$OMP THREADPRIVATE(gamma_taud) 244 238 245 REAL, SAVE, PROTECTED :: gamma_col=1. ! A COMMENTER TODO[-]239 REAL, SAVE, PROTECTED :: gamma_col=1. ! Tuning coefficient for rain collection efficiency (poprecip) [-] 246 240 !$OMP THREADPRIVATE(gamma_col) 247 241 248 REAL, SAVE, PROTECTED :: gamma_agg=1. ! A COMMENTER TODO[-]242 REAL, SAVE, PROTECTED :: gamma_agg=1. ! Tuning coefficient for snow aggregation efficiency (poprecip) [-] 249 243 !$OMP THREADPRIVATE(gamma_agg) 250 244 251 REAL, SAVE, PROTECTED :: gamma_rim=1. ! A COMMENTER TODO[-]245 REAL, SAVE, PROTECTED :: gamma_rim=1. ! Tuning coefficient for riming efficiency (poprecip) [-] 252 246 !$OMP THREADPRIVATE(gamma_rim) 247 248 REAL, SAVE, PROTECTED :: gamma_melt=1. ! Tuning coefficient for snow melting efficiency (poprecip) [-] 249 !$OMP THREADPRIVATE(gamma_melt) 250 251 REAL, SAVE, PROTECTED :: gamma_freez=1. ! Tuning coefficient for rain collision freezing efficiency (poprecip) [-] 252 !$OMP THREADPRIVATE(gamma_freez) 253 253 254 254 REAL, SAVE, PROTECTED :: rho_rain=1000. ! Rain density [kg/m3] … … 258 258 !$OMP THREADPRIVATE(rho_ice) 259 259 260 REAL, SAVE, PROTECTED :: r_rain=500.E-6 ! Rain droplets radius for POPRECIP[m]260 REAL, SAVE, PROTECTED :: r_rain=500.E-6 ! Rain droplets radius (poprecip) [m] 261 261 !$OMP THREADPRIVATE(r_rain) 262 262 263 REAL, SAVE, PROTECTED :: r_snow=1.E-3 ! Ice crystals radius for POPRECIP[m]263 REAL, SAVE, PROTECTED :: r_snow=1.E-3 ! Ice crystals radius (poprecip) [m] 264 264 !$OMP THREADPRIVATE(r_snow) 265 265 266 REAL, SAVE, PROTECTED :: tau_auto_snow_min=100. ! A COMMENTER TODO[s]266 REAL, SAVE, PROTECTED :: tau_auto_snow_min=100. ! Snow autoconversion minimal timescale (when liquid) [s] 267 267 !$OMP THREADPRIVATE(tau_auto_snow_min) 268 268 269 REAL, SAVE, PROTECTED :: tau_auto_snow_max=1000. ! A COMMENTER TODO[s]269 REAL, SAVE, PROTECTED :: tau_auto_snow_max=1000. ! Snow autoconversion minimal timescale (when only ice) [s] 270 270 !$OMP THREADPRIVATE(tau_auto_snow_max) 271 271 272 REAL, SAVE, PROTECTED :: eps=1.E-10 ! A COMMENTER TODO[-]272 REAL, SAVE, PROTECTED :: eps=1.E-10 ! Treshold 0 [-] 273 273 !$OMP THREADPRIVATE(eps) 274 274 275 REAL, SAVE, PROTECTED :: gamma_melt=1. ! A COMMENTER TODO [-] 276 !$OMP THREADPRIVATE(gamma_melt) 277 278 REAL, SAVE, PROTECTED :: alpha_freez=4. ! A COMMENTER TODO [-] 275 REAL, SAVE, PROTECTED :: alpha_freez=4. ! Slope of exponential for immersion freezing timescale [-] 279 276 !$OMP THREADPRIVATE(alpha_freez) 280 277 281 REAL, SAVE, PROTECTED :: beta_freez=0.1 ! A COMMENTER TODO [m-3.s-1]278 REAL, SAVE, PROTECTED :: beta_freez=0.1 ! Inv.time immersion freezing [s-1] 282 279 !$OMP THREADPRIVATE(beta_freez) 283 280 284 REAL, SAVE, PROTECTED :: gamma_freez=1. ! A COMMENTER TODO [-] 285 !$OMP THREADPRIVATE(gamma_freez) 286 287 REAL, SAVE, PROTECTED :: rain_fallspeed=4. ! A COMMENTER TODO [m/s] 281 REAL, SAVE, PROTECTED :: rain_fallspeed=4. ! Rain fall velocity [m/s] 288 282 !$OMP THREADPRIVATE(rain_fallspeed) 289 283 290 REAL, SAVE, PROTECTED :: rain_fallspeed_clr ! A COMMENTER TODO[m/s]284 REAL, SAVE, PROTECTED :: rain_fallspeed_clr ! Rain fall velocity in clear sky [m/s] 291 285 !$OMP THREADPRIVATE(rain_fallspeed_clr) 292 286 293 REAL, SAVE, PROTECTED :: rain_fallspeed_cld ! A COMMENTER TODO[m/s]287 REAL, SAVE, PROTECTED :: rain_fallspeed_cld ! Rain fall velocity in cloudy sky [m/s] 294 288 !$OMP THREADPRIVATE(rain_fallspeed_cld) 295 289 296 REAL, SAVE, PROTECTED :: snow_fallspeed=1. ! A COMMENTER TODO[m/s]290 REAL, SAVE, PROTECTED :: snow_fallspeed=1. ! Snow fall velocity [m/s] 297 291 !$OMP THREADPRIVATE(snow_fallspeed) 298 292 299 REAL, SAVE, PROTECTED :: snow_fallspeed_clr ! A COMMENTER TODO[m/s]293 REAL, SAVE, PROTECTED :: snow_fallspeed_clr ! Snow fall velocity in clear sky [m/s] 300 294 !$OMP THREADPRIVATE(snow_fallspeed_clr) 301 295 302 REAL, SAVE, PROTECTED :: snow_fallspeed_cld ! A COMMENTER TODO[m/s]296 REAL, SAVE, PROTECTED :: snow_fallspeed_cld ! Snow fall velocity in cloudy sky [m/s] 303 297 !$OMP THREADPRIVATE(snow_fallspeed_cld) 304 298 !--End of the parameters for poprecip … … 386 380 CALL getin_p('dist_liq',dist_liq) 387 381 CALL getin_p('tresh_cl',tresh_cl) 388 CALL getin_p('tau_mixenv',tau_mixenv)389 382 CALL getin_p('capa_crystal',capa_crystal) 390 CALL getin_p('lmix_mpc',lmix_mpc)391 383 CALL getin_p('naero5',naero5) 392 384 CALL getin_p('gamma_snwretro',gamma_snwretro) … … 403 395 CALL getin_p('gamma_freez',gamma_freez) 404 396 CALL getin_p('gamma_melt',gamma_melt) 397 CALL getin_p('tau_auto_snow_max',tau_auto_snow_max) 398 CALL getin_p('tau_auto_snow_min',tau_auto_snow_min) 405 399 CALL getin_p('r_snow',r_snow) 406 400 CALL getin_p('rain_fallspeed',rain_fallspeed) … … 471 465 WRITE(lunout,*) 'lscp_ini, dist_liq', dist_liq 472 466 WRITE(lunout,*) 'lscp_ini, tresh_cl', tresh_cl 473 WRITE(lunout,*) 'lscp_ini, tau_mixenv', tau_mixenv474 467 WRITE(lunout,*) 'lscp_ini, capa_crystal', capa_crystal 475 WRITE(lunout,*) 'lscp_ini, lmix_mpc', lmix_mpc476 468 WRITE(lunout,*) 'lscp_ini, naero5', naero5 477 469 WRITE(lunout,*) 'lscp_ini, gamma_snwretro', gamma_snwretro … … 489 481 WRITE(lunout,*) 'lscp_ini, gamma_freez:', gamma_freez 490 482 WRITE(lunout,*) 'lscp_ini, gamma_melt:', gamma_melt 483 WRITE(lunout,*) 'lscp_ini, tau_auto_snow_max:',tau_auto_snow_max 484 WRITE(lunout,*) 'lscp_ini, tau_auto_snow_min:',tau_auto_snow_min 491 485 WRITE(lunout,*) 'lscp_ini, r_snow:', r_snow 492 486 WRITE(lunout,*) 'lscp_ini, rain_fallspeed_clr:', rain_fallspeed_clr
Note: See TracChangeset
for help on using the changeset viewer.
