Changeset 1822 for palm/trunk/SOURCE/prognostic_equations.f90
- Timestamp:
- Apr 7, 2016 7:49:42 AM (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/prognostic_equations.f90
r1818 r1822 19 19 ! Current revisions: 20 20 ! ------------------ 21 ! 21 ! Kessler microphysics scheme moved to microphysics. 22 22 ! 23 23 ! Former revisions: … … 190 190 USE control_parameters, & 191 191 ONLY: call_microphysics_at_all_substeps, cloud_physics, & 192 cloud_top_radiation, constant_diffusion, dp_external, &192 cloud_top_radiation, constant_diffusion, dp_external, & 193 193 dp_level_ind_b, dp_smooth_factor, dpdxy, dt_3d, humidity, & 194 i cloud_scheme, inflow_l, intermediate_timestep_count,&194 inflow_l, intermediate_timestep_count, & 195 195 intermediate_timestep_count_max, large_scale_forcing, & 196 large_scale_subsidence, neutral, nudging, ocean, outflow_l,&197 outflow_s, passive_scalar, precipitation,&198 prho_reference, prho_reference, prho_reference, pt_reference,&199 p t_reference, pt_reference, scalar_advec, scalar_advec,&200 s imulated_time, sloping_surface,&196 large_scale_subsidence, microphysics_seifert, & 197 microphysics_sat_adjust, neutral, nudging, ocean, outflow_l, & 198 outflow_s, passive_scalar, prho_reference, prho_reference, & 199 prho_reference, pt_reference, pt_reference, pt_reference, & 200 scalar_advec, scalar_advec, simulated_time, sloping_surface, & 201 201 timestep_scheme, tsc, use_subsidence_tendencies, & 202 202 use_upstream_for_tke, wall_heatflux, & … … 247 247 USE buoyancy_mod, & 248 248 ONLY: buoyancy, buoyancy_acc 249 250 USE calc_precipitation_mod, &251 ONLY: calc_precipitation252 249 253 250 USE calc_radiation_mod, & … … 271 268 USE diffusion_w_mod, & 272 269 ONLY: diffusion_w, diffusion_w_acc 273 274 USE impact_of_latent_heat_mod, &275 ONLY: impact_of_latent_heat276 270 277 271 USE kinds … … 377 371 DO j = nys, nyn 378 372 ! 379 !-- If required, calculate cloud microphysic al impacts (two-moment scheme)380 IF ( cloud_physics .AND. icloud_scheme == 0 .AND.&373 !-- If required, calculate cloud microphysics 374 IF ( cloud_physics .AND. .NOT. microphysics_sat_adjust .AND. & 381 375 ( intermediate_timestep_count == 1 .OR. & 382 376 call_microphysics_at_all_substeps ) & … … 584 578 IF ( cloud_top_radiation ) THEN 585 579 CALL calc_radiation( i, j ) 586 ENDIF587 588 !589 !-- If required compute impact of latent heat due to precipitation590 IF ( cloud_physics .AND. icloud_scheme == 1 .AND. &591 precipitation ) THEN592 CALL impact_of_latent_heat( i, j )593 580 ENDIF 594 581 … … 729 716 730 717 ! 731 !-- If required compute decrease of total water content due to732 !-- precipitation733 IF ( cloud_physics .AND. icloud_scheme == 1 .AND. &734 precipitation ) THEN735 CALL calc_precipitation( i, j )736 ENDIF737 738 !739 718 !-- Sink or source of scalar concentration due to canopy elements 740 719 IF ( plant_canopy ) CALL plant_canopy_model( i, j, 5 ) … … 788 767 !-- If required, calculate prognostic equations for rain water content 789 768 !-- and rain drop concentration 790 IF ( cloud_physics .AND. icloud_scheme == 0 .AND. & 791 precipitation ) THEN 769 IF ( cloud_physics .AND. microphysics_seifert ) THEN 792 770 ! 793 771 !-- Calculate prognostic equation for rain water content … … 807 785 CALL diffusion_s( i, j, qr, qrsws, qrswst, wall_qrflux ) 808 786 809 CALL user_actions( i, j, 'qr-tendency' )810 787 ! 811 788 !-- Prognostic equation for rain water content … … 848 825 CALL diffusion_s( i, j, nr, nrsws, nrswst, wall_nrflux ) 849 826 850 CALL user_actions( i, j, 'nr-tendency' )851 827 ! 852 828 !-- Prognostic equation for rain drop concentration … … 971 947 972 948 ! 973 !-- If required, calculate cloud microphysical impacts (two-moment scheme)974 IF ( cloud_physics .AND. icloud_scheme == 0 .AND.&949 !-- If required, calculate cloud microphysical impacts 950 IF ( cloud_physics .AND. .NOT. microphysics_sat_adjust .AND. & 975 951 ( intermediate_timestep_count == 1 .OR. & 976 952 call_microphysics_at_all_substeps ) & … … 1254 1230 1255 1231 ! 1256 !-- If required compute impact of latent heat due to precipitation1257 IF ( cloud_physics .AND. icloud_scheme == 1 .AND. precipitation ) THEN1258 CALL impact_of_latent_heat1259 ENDIF1260 1261 !1262 1232 !-- Consideration of heat sources within the plant canopy 1263 1233 IF ( plant_canopy .AND. ( cthf /= 0.0_wp ) ) THEN … … 1459 1429 1460 1430 ! 1461 !-- If required compute decrease of total water content due to1462 !-- precipitation1463 IF ( cloud_physics .AND. icloud_scheme == 1 .AND. precipitation ) THEN1464 CALL calc_precipitation1465 ENDIF1466 1467 !1468 1431 !-- Sink or source of scalar concentration due to canopy elements 1469 1432 IF ( plant_canopy ) CALL plant_canopy_model( 5 ) … … 1530 1493 !-- If required, calculate prognostic equations for rain water content 1531 1494 !-- and rain drop concentration 1532 IF ( cloud_physics .AND. icloud_scheme == 0 .AND. precipitation) THEN1495 IF ( cloud_physics .AND. microphysics_seifert ) THEN 1533 1496 1534 1497 CALL cpu_log( log_point(52), 'qr-equation', 'start' ) … … 1565 1528 1566 1529 CALL diffusion_s( qr, qrsws, qrswst, wall_qrflux ) 1567 1568 CALL user_actions( 'qr-tendency' )1569 1530 1570 1531 ! … … 1639 1600 1640 1601 CALL diffusion_s( nr, nrsws, nrswst, wall_nrflux ) 1641 1642 CALL user_actions( 'nr-tendency' )1643 1602 1644 1603 ! … … 1821 1780 ! 1822 1781 !-- If required, calculate cloud microphysical impacts (two-moment scheme) 1823 IF ( cloud_physics .AND. icloud_scheme == 0 .AND.&1782 IF ( cloud_physics .AND. .NOT. microphysics_sat_adjust .AND. & 1824 1783 ( intermediate_timestep_count == 1 .OR. & 1825 1784 call_microphysics_at_all_substeps ) & … … 2079 2038 2080 2039 ! 2081 !-- If required compute impact of latent heat due to precipitation2082 IF ( cloud_physics .AND. icloud_scheme == 1 .AND. precipitation ) THEN2083 CALL impact_of_latent_heat2084 ENDIF2085 2086 !2087 2040 !-- Consideration of heat sources within the plant canopy 2088 2041 IF ( plant_canopy .AND. ( cthf /= 0.0_wp ) ) THEN … … 2256 2209 2257 2210 ! 2258 !-- If required compute decrease of total water content due to2259 !-- precipitation2260 IF ( cloud_physics .AND. icloud_scheme == 1 .AND. precipitation ) THEN2261 CALL calc_precipitation2262 ENDIF2263 2264 !2265 2211 !-- Sink or source of scalar concentration due to canopy elements 2266 2212 IF ( plant_canopy ) CALL plant_canopy_model( 5 ) … … 2311 2257 !-- If required, calculate prognostic equations for rain water content 2312 2258 !-- and rain drop concentration 2313 IF ( cloud_physics .AND. icloud_scheme == 0 .AND. precipitation) THEN2259 IF ( cloud_physics .AND. microphysics_seifert ) THEN 2314 2260 2315 2261 CALL cpu_log( log_point(52), 'qr-equation', 'start' ) … … 2346 2292 CALL diffusion_s( qr, qrsws, qrswst, wall_qrflux ) 2347 2293 2348 CALL user_actions( 'qr-tendency' )2349 2350 2294 ! 2351 2295 !-- Prognostic equation for rain water content … … 2403 2347 2404 2348 CALL diffusion_s( nr, nrsws, nrswst, wall_nrflux ) 2405 2406 CALL user_actions( 'nr-tendency' )2407 2349 2408 2350 !
Note: See TracChangeset
for help on using the changeset viewer.