Changeset 4533 for palm/trunk/SOURCE/bulk_cloud_model_mod.f90
- Timestamp:
- May 14, 2020 2:46:46 PM (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/bulk_cloud_model_mod.f90
r4532 r4533 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Reformat intrinsic function 28 ! 29 ! 4532 2020-05-14 13:41:35Z schwenkel 27 30 ! Calculate mean droplet radius assuming gamma distibution for condensation 28 31 ! … … 4473 4476 !-- following (Seifert and Beheng, 2006, Kogan and Khairoutdinov, 2000, 4474 4477 !-- Seifert and Stevens, 2010) 4475 rc = MAX( alpha_rc * gamma(nu + 1.33_wp) / gamma(nu + 1.0_wp) * &4478 rc = MAX( alpha_rc * GAMMA(nu + 1.33_wp) / GAMMA(nu + 1.0_wp) * & 4476 4479 ( 3.0_wp * qc(k,j,i) / & 4477 4480 ( 4.0_wp * pi * rho_l * ( nu + 2.0_wp ) * nc(k,j,i) ) & … … 4562 4565 !-- in Seifert and Stevens, 2010 ) is switched off. Minimum radius is set to 1µm following 4563 4566 !-- (Seifert and Beheng, 2006, Kogan and Khairoutdinov, 2000, Seifert and Stevens, 2010) 4564 rc = MAX( alpha_rc * gamma(nu + 1.33_wp) / gamma(nu + 1.0_wp) * &4567 rc = MAX( alpha_rc * GAMMA(nu + 1.33_wp) / GAMMA(nu + 1.0_wp) * & 4565 4568 ( 3.0_wp * qc(k,j,i) / ( 4.0_wp * pi * rho_l * ( nu + 2.0_wp ) * nc(k,j,i) ) & 4566 4569 )**0.33_wp, 1.0E-6_wp )
Note: See TracChangeset
for help on using the changeset viewer.