Changeset 2375 for palm/trunk/SOURCE/microphysics_mod.f90
- Timestamp:
- Aug 29, 2017 2:10:28 PM (7 years ago)
- File:
-
- 1 edited
-
palm/trunk/SOURCE/microphysics_mod.f90 (modified) (19 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/microphysics_mod.f90
r2318 r2375 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Improved aerosol initilization and some minor bugfixes 28 ! for droplet sedimenation 29 ! 30 ! 2318 2017-07-20 17:27:44Z suehring 27 31 ! Get topography top index via Function call 28 32 ! … … 188 192 REAL(wp) :: xrmax = 5.0E-6_wp !< maximum rain drop site (kg) 189 193 190 REAL(wp) :: c_sedimentation = 2.0_wp !< Courant number of sedimentation process 191 REAL(wp) :: dpirho_l !< 6.0 / ( pi * rho_l ) 192 REAL(wp) :: dt_micro !< microphysics time step 193 REAL(wp) :: na_init = 100.0E6_wp !< Total particle/aerosol concentration (cm-3) 194 REAL(wp) :: nc_const = 70.0E6_wp !< cloud droplet concentration 195 REAL(wp) :: dt_precipitation = 100.0_wp !< timestep precipitation (s) 196 REAL(wp) :: sed_qc_const !< const. for sedimentation of cloud water 197 REAL(wp) :: pirho_l !< pi * rho_l / 6.0; 194 REAL(wp) :: c_sedimentation = 2.0_wp !< Courant number of sedimentation process 195 REAL(wp) :: dpirho_l !< 6.0 / ( pi * rho_l ) 196 REAL(wp) :: dry_aerosol_radius = 0.05E-6_wp !< dry aerosol radius 197 REAL(wp) :: dt_micro !< microphysics time step 198 REAL(wp) :: sigma_bulk = 2.0_wp !< width of aerosol spectrum 199 REAL(wp) :: na_init = 100.0E6_wp !< Total particle/aerosol concentration (cm-3) 200 REAL(wp) :: nc_const = 70.0E6_wp !< cloud droplet concentration 201 REAL(wp) :: dt_precipitation = 100.0_wp !< timestep precipitation (s) 202 REAL(wp) :: sed_qc_const !< const. for sedimentation of cloud water 203 REAL(wp) :: pirho_l !< pi * rho_l / 6.0; 198 204 199 205 REAL(wp), DIMENSION(:), ALLOCATABLE :: nc_1d !< … … 210 216 211 217 PUBLIC cloud_water_sedimentation, collision_turbulence, & 212 curvature_solution_effects_bulk, c_sedimentation, dt_precipitation, & 213 limiter_sedimentation, na_init, nc_const, sigma_gc, & 218 curvature_solution_effects_bulk, c_sedimentation, & 219 dry_aerosol_radius, dt_precipitation, & 220 limiter_sedimentation, na_init, nc_const, sigma_gc, sigma_bulk, & 214 221 ventilation_effect 215 222 … … 290 297 291 298 USE cloud_parameters, & 292 ONLY: rho_l299 ONLY: molecular_weight_of_solute, rho_l, rho_s, vanthoff 293 300 294 301 USE control_parameters, & 295 ONLY: microphysics_morrison, microphysics_seifert 302 ONLY: aerosol_nacl, aerosol_c3h4o4 , aerosol_nh4no3, & 303 microphysics_morrison, microphysics_seifert 296 304 297 305 USE indices, & … … 313 321 ENDIF 314 322 323 ! 324 !-- Set constants for certain aerosol type 325 IF ( microphysics_morrison ) THEN 326 IF ( aerosol_nacl ) THEN 327 molecular_weight_of_solute = 0.05844_wp 328 rho_s = 2165.0_wp 329 vanthoff = 2.0_wp 330 ELSEIF ( aerosol_c3h4o4 ) THEN 331 molecular_weight_of_solute = 0.10406_wp 332 rho_s = 1600.0_wp 333 vanthoff = 1.37_wp 334 ELSEIF ( aerosol_nh4no3 ) THEN 335 molecular_weight_of_solute = 0.08004_wp 336 rho_s = 1720.0_wp 337 vanthoff = 2.31_wp 338 ENDIF 339 ENDIF 340 315 341 ! 316 342 !-- Pre-calculate frequently calculated fractions of pi and rho_l … … 509 535 510 536 USE cloud_parameters, & 511 ONLY: hyrho, l_d_cp, l_d_r, l_v, rho_l, r_v, t_d_pt 537 ONLY: hyrho, l_d_cp, l_d_r, l_v, molecular_weight_of_solute, & 538 molecular_weight_of_water, rho_l, rho_s, r_v, t_d_pt, & 539 vanthoff 512 540 513 541 USE constants, & … … 524 552 USE control_parameters, & 525 553 ONLY: simulated_time 526 527 USE particle_attributes, &528 ONLY: molecular_weight_of_solute, molecular_weight_of_water, rho_s, &529 log_sigma, vanthoff530 554 531 555 IMPLICIT NONE … … 545 569 REAL(wp) :: n_ccn !< 546 570 REAL(wp) :: q_s !< 547 REAL(wp) :: rd0 !<548 571 REAL(wp) :: s_0 !< 549 572 REAL(wp) :: sat !< 550 573 REAL(wp) :: sat_max !< 551 574 REAL(wp) :: sigma !< 575 REAL(wp) :: sigma_act !< 552 576 REAL(wp) :: t_int !< 553 577 REAL(wp) :: t_l !< … … 616 640 ! 617 641 !-- Prescribe power index that describes the soluble fraction 618 !-- of an aerosol particle (beta) and mean geometric radius of 619 !-- dry aerosol spectrum 642 !-- of an aerosol particle (beta) 620 643 !-- (see: Morrison + Grabowski, 2007, JAS, 64) 621 beta_act = 0.5_wp622 rd0 = 0.05E-6_wp644 beta_act = 0.5_wp 645 sigma_act = sigma_bulk**( 1.0_wp + beta_act ) 623 646 ! 624 647 !-- Calculate mean geometric supersaturation (s_0) with 625 648 !-- parameterization by Khvorostyanov and Curry (2006) 626 s_0 = rd0 **(- ( 1.0_wp + beta_act ) ) *&649 s_0 = dry_aerosol_radius **(- ( 1.0_wp + beta_act ) ) * & 627 650 ( 4.0_wp * afactor**3 / ( 27.0_wp * bfactor ) )**0.5_wp 628 651 … … 632 655 !-- (see: Khvorostyanov + Curry, 2006, J. Geo. Res., 111) 633 656 n_ccn = ( na_init / 2.0_wp ) * ( 1.0_wp - ERF( & 634 LOG( s_0 / sat ) / ( SQRT(2.0_wp) * log_sigma(1) ) ) )657 LOG( s_0 / sat ) / ( SQRT(2.0_wp) * LOG(sigma_act) ) ) ) 635 658 activ = MAX( ( n_ccn - nc(k,j,i) ) / dt_micro, 0.0_wp ) 636 659 ENDIF … … 1357 1380 1358 1381 sed_qc(nzt+1) = 0.0_wp 1382 sed_nc(nzt+1) = 0.0_wp 1359 1383 1360 1384 DO i = nxlg, nxrg … … 1386 1410 ENDIF 1387 1411 1388 IF ( qc(k,j,i) > eps_sb ) THEN1412 IF ( qc(k,j,i) > eps_sb .AND. nc_sedi > eps_mr ) THEN 1389 1413 sed_qc(k) = sed_qc_const * nc_sedi**( -2.0_wp / 3.0_wp ) * & 1390 1414 ( qc(k,j,i) * hyrho(k) )**( 5.0_wp / 3.0_wp ) * & … … 1948 1972 1949 1973 USE cloud_parameters, & 1950 ONLY: hyrho, l_d_cp, l_d_r, l_v, rho_l, r_v, t_d_pt 1974 ONLY: hyrho, l_d_cp, l_d_r, l_v, molecular_weight_of_solute, & 1975 molecular_weight_of_water, rho_l, rho_s, r_v, t_d_pt, & 1976 vanthoff 1951 1977 1952 1978 USE constants, & … … 1963 1989 USE control_parameters, & 1964 1990 ONLY: simulated_time 1965 1966 USE particle_attributes, &1967 ONLY: molecular_weight_of_solute, molecular_weight_of_water, rho_s, &1968 log_sigma, vanthoff1969 1991 1970 1992 IMPLICIT NONE … … 1984 2006 REAL(wp) :: n_ccn !< 1985 2007 REAL(wp) :: q_s !< 1986 REAL(wp) :: rd0 !<1987 2008 REAL(wp) :: s_0 !< 1988 2009 REAL(wp) :: sat !< 1989 2010 REAL(wp) :: sat_max !< 1990 2011 REAL(wp) :: sigma !< 2012 REAL(wp) :: sigma_act !< 1991 2013 REAL(wp) :: t_int !< 1992 2014 REAL(wp) :: t_l !< … … 2052 2074 ! 2053 2075 !-- Prescribe power index that describes the soluble fraction 2054 !-- of an aerosol particle (beta) and mean geometric radius of 2055 !-- dry aerosol spectrum 2076 !-- of an aerosol particle (beta). 2056 2077 !-- (see: Morrison + Grabowski, 2007, JAS, 64) 2057 beta_act = 0.5_wp2058 rd0 = 0.05E-6_wp2078 beta_act = 0.5_wp 2079 sigma_act = sigma_bulk**( 1.0_wp + beta_act ) 2059 2080 ! 2060 2081 !-- Calculate mean geometric supersaturation (s_0) with 2061 2082 !-- parameterization by Khvorostyanov and Curry (2006) 2062 s_0 = rd0 **(- ( 1.0_wp + beta_act ) ) *&2083 s_0 = dry_aerosol_radius **(- ( 1.0_wp + beta_act ) ) * & 2063 2084 ( 4.0_wp * afactor**3 / ( 27.0_wp * bfactor ) )**0.5_wp 2064 2085 … … 2067 2088 !-- supersaturation and taking Koehler theory into account 2068 2089 !-- (see: Khvorostyanov + Curry, 2006, J. Geo. Res., 111) 2069 n_ccn = ( na_init / 2.0_wp ) * ( 1.0_wp - ERF( 2070 LOG( s_0 / sat ) / ( SQRT(2.0_wp) * log_sigma(1) ) ) )2090 n_ccn = ( na_init / 2.0_wp ) * ( 1.0_wp - ERF( & 2091 LOG( s_0 / sat ) / ( SQRT(2.0_wp) * LOG(sigma_act) ) ) ) 2071 2092 activ = MAX( ( n_ccn ) / dt_micro, 0.0_wp ) 2072 2093 … … 2709 2730 2710 2731 sed_qc(nzt+1) = 0.0_wp 2732 sed_nc(nzt+1) = 0.0_wp 2733 2711 2734 2712 2735 DO k = nzt, nzb+1, -1 … … 2735 2758 ENDIF 2736 2759 2737 IF ( qc_1d(k) > eps_sb ) THEN2760 IF ( qc_1d(k) > eps_sb .AND. nc_sedi > eps_mr ) THEN 2738 2761 sed_qc(k) = sed_qc_const * nc_sedi**( -2.0_wp / 3.0_wp ) * & 2739 2762 ( qc_1d(k) * hyrho(k) )**( 5.0_wp / 3.0_wp ) * flag
Note: See TracChangeset
for help on using the changeset viewer.
