Changeset 1831 for palm/trunk/SOURCE/microphysics.f90
- Timestamp:
- Apr 7, 2016 1:15:51 PM (8 years ago)
- File:
-
- 1 edited
-
palm/trunk/SOURCE/microphysics.f90 (modified) (15 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/microphysics.f90
r1823 r1831 19 19 ! Current revisions: 20 20 ! ------------------ 21 ! 21 ! turbulence renamed collision_turbulence, 22 ! drizzle renamed cloud_water_sedimentation. cloud_water_sedimentation also 23 ! avaialble for microphysics_kessler. 22 24 ! 23 25 ! Former revisions: … … 163 165 164 166 USE cloud_parameters, & 165 ONLY: c p, hyrho, prr, pt_d_t, r_d, t_d_pt167 ONLY: cloud_water_sedimentation, cp, hyrho, prr, pt_d_t, r_d, t_d_pt 166 168 167 169 USE control_parameters, & 168 ONLY: call_microphysics_at_all_substeps, d rizzle, dt_3d, dt_micro,&170 ONLY: call_microphysics_at_all_substeps, dt_3d, dt_micro, & 169 171 g, intermediate_timestep_count, large_scale_forcing, & 170 172 lsf_surf, microphysics_kessler, microphysics_seifert, & … … 223 225 224 226 CALL autoconversion_kessler 227 IF ( cloud_water_sedimentation ) CALL sedimentation_cloud 225 228 226 229 ELSEIF ( microphysics_seifert ) THEN … … 232 235 CALL evaporation_rain 233 236 CALL sedimentation_rain 234 235 IF ( drizzle ) CALL sedimentation_cloud 237 IF ( cloud_water_sedimentation ) CALL sedimentation_cloud 236 238 237 239 ENDIF … … 306 308 USE cloud_parameters, & 307 309 ONLY: a_1, a_2, a_3, b_1, b_2, b_3, beta_cc, c_1, c_2, c_3, & 308 c_const, dpirho_l, eps_sb, hyrho, k_cc, kin_vis_air,&309 nc_const, x0310 c_const, collision_turbulence, dpirho_l, eps_sb, hyrho, & 311 k_cc, kin_vis_air, nc_const, x0 310 312 311 313 USE control_parameters, & 312 ONLY: dt_micro, rho_surface , turbulence314 ONLY: dt_micro, rho_surface 313 315 314 316 USE cpulog, & … … 371 373 !-- Parameterized turbulence effects on autoconversion (Seifert, 372 374 !-- Nuijens and Stevens, 2010) 373 IF ( turbulence ) THEN375 IF ( collision_turbulence ) THEN 374 376 ! 375 377 !-- Weight averaged radius of cloud droplets: … … 495 497 496 498 USE cloud_parameters, & 497 ONLY: eps_sb, hyrho, k_cr0499 ONLY: collision_turbulence, eps_sb, hyrho, k_cr0 498 500 499 501 USE control_parameters, & 500 ONLY: dt_micro, rho_surface , turbulence502 ONLY: dt_micro, rho_surface 501 503 502 504 USE cpulog, & … … 537 539 !-- Nuijens and Stevens, 2010). The factor of 1.0E4 is needed to 538 540 !-- convert the dissipation rate (diss) from m2 s-3 to cm2 s-3. 539 IF ( turbulence ) THEN541 IF ( collision_turbulence ) THEN 540 542 k_cr = k_cr0 * ( 1.0_wp + 0.05_wp * & 541 543 MIN( 600.0_wp, & … … 1165 1167 1166 1168 USE cloud_parameters, & 1167 ONLY: cp, hyrho, nc_const, prr, pt_d_t, r_d, t_d_pt 1169 ONLY: cloud_water_sedimentation, cp, hyrho, nc_const, prr, pt_d_t, & 1170 r_d, t_d_pt 1168 1171 1169 1172 USE control_parameters, & 1170 ONLY: call_microphysics_at_all_substeps, d rizzle, dt_3d, dt_micro,&1173 ONLY: call_microphysics_at_all_substeps, dt_3d, dt_micro, & 1171 1174 g, intermediate_timestep_count, large_scale_forcing, & 1172 1175 lsf_surf, microphysics_seifert, microphysics_kessler, & … … 1236 1239 1237 1240 CALL autoconversion_kessler( i,j ) 1241 IF ( cloud_water_sedimentation ) CALL sedimentation_cloud( i,j ) 1238 1242 1239 1243 ELSEIF ( microphysics_seifert ) THEN … … 1245 1249 CALL evaporation_rain( i,j ) 1246 1250 CALL sedimentation_rain( i,j ) 1247 1248 IF ( drizzle ) CALL sedimentation_cloud( i,j ) 1251 IF ( cloud_water_sedimentation ) CALL sedimentation_cloud( i,j ) 1249 1252 1250 1253 ENDIF … … 1326 1329 USE cloud_parameters, & 1327 1330 ONLY: a_1, a_2, a_3, b_1, b_2, b_3, beta_cc, c_1, c_2, c_3, & 1328 c_const, dpirho_l, eps_sb, hyrho, kin_vis_air, k_cc, x0 1331 c_const, collision_turbulence, dpirho_l, eps_sb, hyrho, & 1332 kin_vis_air, k_cc, x0 1329 1333 1330 1334 USE control_parameters, & 1331 ONLY: dt_micro, rho_surface , turbulence1335 ONLY: dt_micro, rho_surface 1332 1336 1333 1337 USE grid_variables, & … … 1382 1386 !-- Parameterized turbulence effects on autoconversion (Seifert, 1383 1387 !-- Nuijens and Stevens, 2010) 1384 IF ( turbulence ) THEN1388 IF ( collision_turbulence ) THEN 1385 1389 ! 1386 1390 !-- Weight averaged radius of cloud droplets: … … 1495 1499 1496 1500 USE cloud_parameters, & 1497 ONLY: eps_sb, hyrho, k_cr01501 ONLY: collision_turbulence, eps_sb, hyrho, k_cr0 1498 1502 1499 1503 USE control_parameters, & 1500 ONLY: dt_micro, rho_surface , turbulence1504 ONLY: dt_micro, rho_surface 1501 1505 1502 1506 USE indices, & … … 1529 1533 !-- Nuijens and Stevens, 2010). The factor of 1.0E4 is needed to 1530 1534 !-- convert the dissipation rate (diss) from m2 s-3 to cm2 s-3. 1531 IF ( turbulence ) THEN1535 IF ( collision_turbulence ) THEN 1532 1536 k_cr = k_cr0 * ( 1.0_wp + 0.05_wp * & 1533 1537 MIN( 600.0_wp, &
Note: See TracChangeset
for help on using the changeset viewer.
