Changeset 3241 for palm/trunk/SOURCE/microphysics_mod.f90
- Timestamp:
- Sep 12, 2018 3:02:00 PM (6 years ago)
- File:
-
- 1 edited
-
palm/trunk/SOURCE/microphysics_mod.f90 (modified) (19 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/microphysics_mod.f90
r3026 r3241 25 25 ! ----------------- 26 26 ! $Id$ 27 ! unused variables removed 28 ! 29 ! 3026 2018-05-22 10:30:53Z schwenkel 27 30 ! Changed the name specific humidity to mixing ratio, since we are computing 28 31 ! mixing ratios. … … 201 204 REAL(wp) :: k_br = 1000.0_wp !< const. in breakup parametrization (m-1) 202 205 REAL(wp) :: k_st = 1.2E8_wp !< const. in drizzle parametrization (m-1 s-1) 203 REAL(wp) :: kappa_rr = 60.7_wp !< const. in collision kernel (kg-1/3)204 206 REAL(wp) :: kin_vis_air = 1.4086E-5_wp !< kin. viscosity of air (m2 s-1) 205 207 REAL(wp) :: prec_time_const = 0.001_wp !< coef. in Kessler scheme (s-1) … … 212 214 REAL(wp) :: xamin = 5.24E-19_wp !< average aerosol mass (kg) (~ 0.05µm) 213 215 REAL(wp) :: xcmin = 4.18E-15_wp !< minimum cloud drop size (kg) (~ 1µm) 214 REAL(wp) :: xcmax = 2.6E-10_wp !< maximum cloud drop size (kg) (~ 40µm)215 216 REAL(wp) :: xrmin = 2.6E-10_wp !< minimum rain drop size (kg) 216 217 REAL(wp) :: xrmax = 5.0E-6_wp !< maximum rain drop site (kg) … … 556 557 557 558 USE arrays_3d, & 558 ONLY: hyp, nc, nr, pt, q, qc, qr559 ONLY: nc 559 560 560 561 USE cloud_parameters, & 561 ONLY: hyrho, l_d_cp, l_d_r, l_v, molecular_weight_of_solute, & 562 molecular_weight_of_water, rho_l, rho_s, r_v, t_d_pt, & 563 vanthoff 564 565 USE constants, & 566 ONLY: pi 562 ONLY: molecular_weight_of_solute, & 563 molecular_weight_of_water, rho_l, rho_s, r_v, vanthoff 567 564 568 565 USE cpulog, & … … 570 567 571 568 USE diagnostic_quantities_mod, & 572 ONLY: e_s, magnus, q_s, sat, supersaturation, t_l569 ONLY: magnus, sat, supersaturation, t_l 573 570 574 571 USE indices, & … … 576 573 577 574 USE kinds 578 579 USE control_parameters, &580 ONLY: simulated_time581 575 582 576 IMPLICIT NONE … … 684 678 685 679 USE arrays_3d, & 686 ONLY: hyp, nr, pt,q, qc, qr, nc680 ONLY: q, qc, qr, nc 687 681 688 682 USE cloud_parameters, & 689 ONLY: hyrho, l_d_cp, l_ d_r, l_v, r_v, t_d_pt683 ONLY: hyrho, l_d_cp, l_v, r_v 690 684 691 685 USE constants, & … … 702 696 703 697 USE kinds 704 705 USE control_parameters, &706 ONLY: simulated_time707 698 708 699 IMPLICIT NONE … … 716 707 REAL(wp) :: dc !< 717 708 REAL(wp) :: evap !< 718 REAL(wp) :: evap_nc !<719 709 REAL(wp) :: g_fac !< 720 710 REAL(wp) :: nc_0 !< … … 1173 1163 1174 1164 USE arrays_3d, & 1175 ONLY: hyp, nr, pt, q,qc, qr1165 ONLY: nr, qc, qr 1176 1166 1177 1167 USE cloud_parameters, & 1178 ONLY: hyrho, l_d_cp, l_ d_r, l_v, r_v, t_d_pt1168 ONLY: hyrho, l_d_cp, l_v, r_v 1179 1169 1180 1170 USE constants, & … … 1185 1175 1186 1176 USE diagnostic_quantities_mod, & 1187 ONLY: e_s, magnus, q_s,sat, supersaturation, t_l1177 ONLY: e_s, magnus, sat, supersaturation, t_l 1188 1178 1189 1179 USE indices, & … … 1461 1451 INTEGER(iwp) :: k !< running index z direction 1462 1452 INTEGER(iwp) :: k_run !< 1463 INTEGER(iwp) :: l !< running index of surface type1464 1453 INTEGER(iwp) :: m !< running index surface elements 1465 1454 INTEGER(iwp) :: surf_e !< End index of surface elements at (j,i)-gridpoint … … 1713 1702 precipitation_amount_interval, time_do2d_xy 1714 1703 1715 USE indices, &1716 ONLY: nxl, nxr, nys, nyn, nzb, nzt, wall_flags_01717 1718 1704 USE kinds 1719 1705 … … 1727 1713 INTEGER(iwp) :: k !< running index y direction 1728 1714 INTEGER(iwp) :: m !< running index surface elements 1729 INTEGER(iwp) :: surf_e !< End index of surface elements at (j,i)-gridpoint1730 INTEGER(iwp) :: surf_s !< Start index of surface elements at (j,i)-gridpoint1731 1715 1732 1716 IF ( ( dt_do2d_xy - time_do2d_xy ) < precipitation_amount_interval .AND.& … … 1945 1929 1946 1930 USE arrays_3d, & 1947 ONLY: hyp , nr, pt, q, qc, qr, nc1931 ONLY: hyp 1948 1932 1949 1933 USE cloud_parameters, & 1950 ONLY: hyrho, l_d_cp, l_d_r, l_v, molecular_weight_of_solute,&1934 ONLY: l_d_cp, l_d_r, molecular_weight_of_solute, & 1951 1935 molecular_weight_of_water, rho_l, rho_s, r_v, t_d_pt, & 1952 1936 vanthoff 1953 1937 1954 USE constants, &1955 ONLY: pi1956 1957 1938 USE cpulog, & 1958 ONLY: cpu_log , log_point_s1939 ONLY: cpu_log 1959 1940 1960 1941 USE indices, & 1961 ONLY: n xlg, nxrg, nysg, nyng, nzb, nzt1942 ONLY: nzb, nzt 1962 1943 1963 1944 USE kinds 1964 1965 USE control_parameters, &1966 ONLY: simulated_time1967 1945 1968 1946 IMPLICIT NONE … … 2085 2063 2086 2064 USE arrays_3d, & 2087 ONLY: hyp , nr, pt, q, qc, qr, nc2065 ONLY: hyp 2088 2066 2089 2067 USE cloud_parameters, & … … 2094 2072 2095 2073 USE cpulog, & 2096 ONLY: cpu_log , log_point_s2074 ONLY: cpu_log 2097 2075 2098 2076 USE indices, & 2099 ONLY: n xlg, nxrg, nysg, nyng, nzb, nzt2077 ONLY: nzb, nzt 2100 2078 2101 2079 USE kinds 2102 2103 USE control_parameters, &2104 ONLY: simulated_time2105 2080 2106 2081 IMPLICIT NONE … … 2116 2091 REAL(wp) :: e_s !< 2117 2092 REAL(wp) :: evap !< 2118 REAL(wp) :: evap_nc !<2119 2093 REAL(wp) :: g_fac !< 2120 2094 REAL(wp) :: nc_0 !< … … 3027 3001 precipitation_amount_interval, time_do2d_xy 3028 3002 3029 USE indices, &3030 ONLY: nzb, nzt, wall_flags_03031 3032 3003 USE kinds 3033 3004
Note: See TracChangeset
for help on using the changeset viewer.
