Changeset 1849 for palm/trunk/SOURCE/modules.f90
 Timestamp:
 Apr 8, 2016 11:33:18 AM (5 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

palm/trunk/SOURCE/modules.f90
r1846 r1849 19 19 ! Current revisions: 20 20 !  21 ! 22 ! 21 ! bfactor, mass_of_solute, molecular_weight_of_solute, molecular_weight_of_water, 22 ! vanthoff moved to mod_particle_attributes. 23 ! dt_micro and several cloud_parameters moved to microphysics_mod. 24 ! 1dmicrophysics profiles moved to microphysics_mod. 25 ! 23 26 ! Former revisions: 24 27 !  … … 394 397 c_u_m, c_u_m_l, c_v_m, c_v_m_l, c_w_m, c_w_m_l, ddzu, ddzu_pres, & 395 398 dd2zu, dzu, ddzw, dzw, hyp, inflow_damping_factor, l_grid, & 396 nc_1d, nr_1d, ptdf_x, ptdf_y, p_surf, pt_surf, pt_1d, pt_init, & 397 qsws_surf, q_1d, q_init, q_surf, qc_1d, qr_1d, rdf, rdf_sc, & 398 ref_state, sa_init, shf_surf, timenudge, time_surf, time_vert, & 399 tmp_tnudge, ug, u_init, u_nzb_p1_for_vfc, vg, v_init, & 399 ptdf_x, ptdf_y, p_surf, pt_surf, pt_init, qsws_surf, q_init, q_surf, & 400 rdf, rdf_sc, ref_state, sa_init, shf_surf, timenudge, time_surf, & 401 time_vert, tmp_tnudge, ug, u_init, u_nzb_p1_for_vfc, vg, v_init, & 400 402 v_nzb_p1_for_vfc, w_subs, zu, zw 401 403 … … 407 409 mean_inflow_profiles, nrs, nrsws, nrswst, & 408 410 ol, & !< Obukhov length 409 p tnudge, pt_slope_ref,&411 precipitation_amount, precipitation_rate, ptnudge, pt_slope_ref, & 410 412 qnudge, qs, qsws, qswst, qswst_remote, qrs, qrsws, qrswst, & 411 413 saswsb, saswst, shf, tnudge, td_lsa_lpt, td_lsa_q, td_sub_lpt, & … … 421 423 diss_l_v, diss_l_w, flux_l_e, flux_l_nr, flux_l_pt, flux_l_q, & 422 424 flux_l_qr, flux_l_sa, flux_l_u, flux_l_v, flux_l_w, kh, km, & 423 l_wall, p _loc, tend, tric,&425 l_wall, prr, p_loc, tend, tric, & 424 426 u_m_l, u_m_n, u_m_r, u_m_s, v_m_l, v_m_n, v_m_r, v_m_s, w_m_l, & 425 427 w_m_n, w_m_r, w_m_s … … 490 492 USE kinds 491 493 492 LOGICAL :: cloud_water_sedimentation = .FALSE., & !< cloud water sedimentation in bulk cloud physics 493 limiter_sedimentation = .TRUE., & !< sedimentation limiter in bulk cloud physics 494 collision_turbulence = .FALSE., & !< turbulence effects in bulk cloud physics 495 ventilation_effect = .TRUE. !< ventilation effect in bulk cloud physics 496 497 498 REAL(wp) :: a_1 = 8.69E4_wp, & !< coef. in turb. parametrization (cm2 s3) 499 a_2 = 7.38E5_wp, & !< coef. in turb. parametrization (cm2 s3) 500 a_3 = 1.40E2_wp, & !< coef. in turb. parametrization 501 a_term = 9.65_wp, & !< coef. for terminal velocity (m s1) 502 a_vent = 0.78_wp, & !< coef. for ventilation effect 503 b_1 = 11.45E6_wp, & !< coef. in turb. parametrization (m) 504 b_2 = 9.68E6_wp, & !< coef. in turb. parametrization (m) 505 b_3 = 0.62_wp, & !< coef. in turb. parametrization 506 b_term = 9.8_wp, & !< coef. for terminal velocity (m s1) 507 b_vent = 0.308_wp, & !< coef. for ventilation effect 508 beta_cc = 3.09E4_wp, & !< coef. in turb. parametrization (cm2 s3) 509 bfactor, & !< solution effect on diffusional growth 510 c_1 = 4.82E6_wp, & !< coef. in turb. parametrization (m) 511 c_2 = 4.8E6_wp, & !< coef. in turb. parametrization (m) 512 c_3 = 0.76_wp, & !< coef. in turb. parametrization 513 c_const = 0.93_wp, & !< const. in Taylormicroscale Reynolds number 514 c_evap = 0.7_wp, & !< constant in evaporation 515 c_sedimentation = 2.0_wp, & !< Courant number of sedimentation process 516 c_term = 600.0_wp, & !< coef. for terminal velocity (m1) 517 cof(6) = (/ 76.18009172947146_wp, & !< coefficients in the 518 86.50532032941677_wp, & !< numerical 519 24.01409824083091_wp, & !< calculation of the 520 1.231739572450155_wp, & !< gamma function 521 0.1208650973866179E2_wp, & 522 0.5395239384953E5_wp /), & 523 cp = 1005.0_wp, & !< heat capacity of dry air (J kg1 K1) 524 diff_coeff_l = 0.23E4_wp, & !< diffusivity of water vapor (m2 s1) 525 effective_coll_efficiency, & !< effective collision efficiency 526 eps_ros = 1.0E4_wp, & !< accuracy of Rosenbrock method 527 eps_sb = 1.0E20_wp, & !< threshold in twomoments scheme 528 k_cc = 9.44E09_wp, & !< const. cloudcloud kernel (m3 kg2 s1) 529 k_cr0 = 4.33_wp, & !< const. cloudrain kernel (m3 kg1 s1) 530 k_rr = 7.12_wp, & !< const. rainrain kernel (m3 kg1 s1) 531 k_br = 1000._wp, & !< const. in breakup parametrization (m1) 532 k_st = 1.2E8_wp, & !< const. in drizzle parametrization (m1 s1) 533 kappa_rr = 60.7_wp, & !< const. in collision kernel (kg1/3) 534 kin_vis_air = 1.4086E5_wp, & !< kin. viscosity of air (m2 s1) 494 REAL(wp) :: cp = 1005.0_wp, & !< heat capacity of dry air (J kg1 K1) 535 495 l_v = 2.5E+06_wp, & !< latent heat of vaporization (J kg1) 536 496 l_d_cp, l_d_r, l_d_rv, & !< l_v / cp, l_v / r_d, l_v / r_v 537 mass_of_solute = 1.0E17_wp, & !< soluted NaCl (kg)538 molecular_weight_of_solute = 0.05844_wp, & !< mol. m. NaCl (kg mol1)539 molecular_weight_of_water = 0.01801528_wp, & !< mol. m. H2O (kg mol1)540 nc_const = 70.0E6_wp, & !< cloud droplet concentration541 prec_time_const = 0.001_wp, & !< coef. in Kessler scheme542 pirho_l, dpirho_l, & !< pi * rho_l / 6.0; 6.0 / ( pi * rho_l )543 497 rho_l = 1.0E3_wp, & !< density of water (kg m3) 544 ql_crit = 0.0005_wp, & !< coef. in Kessler scheme545 498 r_d = 287.0_wp, & !< sp. gas const. dry air (J kg1 K1) 546 r_v = 461.51_wp, & !< sp. gas const. water vapor (J kg1 K1) 547 schmidt = 0.71_wp, & !< Schmidt number 548 schmidt_p_1d3, & !< schmidt**( 1.0 / 3.0 ) 549 sed_qc_const, & !< const. for sedimentation of cloud water 550 sigma_gc = 1.3_wp, & !< lognormal geometric standard deviation 551 stp = 2.5066282746310005_wp, & !< parameter in gamma function 552 thermal_conductivity_l = 2.43E2_wp, & !< therm. cond. air (J m1 s1 K1) 553 vanthoff = 2.0_wp, & !< van't Hoff factor for NaCl 554 x0 = 2.6E10_wp, & !< separating drop mass (kg) 555 xrmin = 2.6E10_wp, & !< minimum rain drop size (kg) 556 xrmax = 5.0E6_wp, & !< maximum rain drop site (kg) 557 dt_precipitation = 100.0_wp, & !< timestep precipitation (s) 558 w_precipitation = 9.65_wp !< maximum terminal velocity (m s1) 499 r_v = 461.51_wp !< sp. gas const. water vapor (J kg1 K1) 500 559 501 560 502 REAL(wp), DIMENSION(:), ALLOCATABLE :: hyrho, pt_d_t, t_d_pt 561 562 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: precipitation_amount, &563 precipitation_rate564 !565 ! 3D array of precipitation rate566 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: prr567 503 568 504 SAVE … … 807 743 dt_do2d_xy = 9999999.9_wp, dt_do2d_xz = 9999999.9_wp, & 808 744 dt_do2d_yz = 9999999.9_wp, dt_do3d = 9999999.9_wp, dt_dvrp = 9999999.9_wp, & 809 dt_max = 20.0_wp, dt_micro = 1.0_wp,&745 dt_max = 20.0_wp, & 810 746 dt_restart = 9999999.9_wp, & 811 747 dt_run_control = 60.0_wp, dt_3d = 1.0_wp, dz = 1.0_wp, &
Note: See TracChangeset
for help on using the changeset viewer.