Changeset 3874
- Timestamp:
- Apr 8, 2019 4:53:48 PM (6 years ago)
- Location:
- palm/trunk/SOURCE
- Files:
-
- 4 edited
-
Makefile (modified) (1 diff)
-
bulk_cloud_model_mod.f90 (modified) (5 diffs)
-
module_interface.f90 (modified) (7 diffs)
-
prognostic_equations.f90 (modified) (8 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/Makefile
r3872 r3874 1510 1510 advec_ws.o \ 1511 1511 advec_w_up.o \ 1512 bulk_cloud_model_mod.o \1513 1512 buoyancy.o \ 1514 1513 chemistry_model_mod.o \ -
palm/trunk/SOURCE/bulk_cloud_model_mod.f90
r3870 r3874 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Implemented non_transport_physics module interfaces 28 ! 29 ! 3870 2019-04-08 13:44:34Z knoop 27 30 ! Moving prognostic equations of bcm into bulk_cloud_model_mod 28 31 ! … … 360 363 bcm_header, & 361 364 bcm_actions, & 362 bcm_ actions_micro, &365 bcm_non_transport_physics, & 363 366 bcm_prognostic_equations, & 364 367 bcm_3d_data_averaging, & … … 420 423 END INTERFACE bcm_actions 421 424 422 INTERFACE bcm_ actions_micro423 MODULE PROCEDURE bcm_ actions_micro424 MODULE PROCEDURE bcm_ actions_micro_ij425 END INTERFACE bcm_ actions_micro425 INTERFACE bcm_non_transport_physics 426 MODULE PROCEDURE bcm_non_transport_physics 427 MODULE PROCEDURE bcm_non_transport_physics_ij 428 END INTERFACE bcm_non_transport_physics 426 429 427 430 INTERFACE bcm_prognostic_equations … … 1174 1177 !> Control of microphysics for all grid points 1175 1178 !------------------------------------------------------------------------------! 1176 SUBROUTINE bcm_actions_micro 1177 1178 IMPLICIT NONE 1179 1180 IF ( large_scale_forcing .AND. lsf_surf ) THEN 1181 ! 1182 !-- Calculate vertical profile of the hydrostatic pressure (hyp) 1183 hyp = barometric_formula(zu, pt_surface * exner_function(surface_pressure * 100.0_wp), surface_pressure * 100.0_wp) 1184 d_exner = exner_function_invers(hyp) 1185 exner = 1.0_wp / exner_function_invers(hyp) 1186 hyrho = ideal_gas_law_rho_pt(hyp, pt_init) 1187 ! 1188 !-- Compute reference density 1189 rho_surface = ideal_gas_law_rho(surface_pressure * 100.0_wp, pt_surface * exner_function(surface_pressure * 100.0_wp)) 1179 SUBROUTINE bcm_non_transport_physics 1180 1181 1182 CALL cpu_log( log_point(51), 'microphysics', 'start' ) 1183 1184 IF ( .NOT. microphysics_sat_adjust .AND. & 1185 ( intermediate_timestep_count == 1 .OR. & 1186 call_microphysics_at_all_substeps ) ) & 1187 THEN 1188 1189 IF ( large_scale_forcing .AND. lsf_surf ) THEN 1190 ! 1191 !-- Calculate vertical profile of the hydrostatic pressure (hyp) 1192 hyp = barometric_formula(zu, pt_surface * exner_function(surface_pressure * 100.0_wp), surface_pressure * 100.0_wp) 1193 d_exner = exner_function_invers(hyp) 1194 exner = 1.0_wp / exner_function_invers(hyp) 1195 hyrho = ideal_gas_law_rho_pt(hyp, pt_init) 1196 ! 1197 !-- Compute reference density 1198 rho_surface = ideal_gas_law_rho(surface_pressure * 100.0_wp, pt_surface * exner_function(surface_pressure * 100.0_wp)) 1199 ENDIF 1200 1201 ! 1202 !-- Compute length of time step 1203 IF ( call_microphysics_at_all_substeps ) THEN 1204 dt_micro = dt_3d * weight_pres(intermediate_timestep_count) 1205 ELSE 1206 dt_micro = dt_3d 1207 ENDIF 1208 1209 ! 1210 !-- Reset precipitation rate 1211 IF ( intermediate_timestep_count == 1 ) prr = 0.0_wp 1212 1213 ! 1214 !-- Compute cloud physics 1215 IF ( microphysics_kessler ) THEN 1216 1217 CALL autoconversion_kessler 1218 IF ( cloud_water_sedimentation ) CALL sedimentation_cloud 1219 1220 ELSEIF ( microphysics_seifert ) THEN 1221 1222 CALL adjust_cloud 1223 IF ( microphysics_morrison ) CALL activation 1224 IF ( microphysics_morrison ) CALL condensation 1225 CALL autoconversion 1226 CALL accretion 1227 CALL selfcollection_breakup 1228 CALL evaporation_rain 1229 CALL sedimentation_rain 1230 IF ( cloud_water_sedimentation ) CALL sedimentation_cloud 1231 1232 ENDIF 1233 1234 CALL calc_precipitation_amount 1235 1190 1236 ENDIF 1191 1237 1192 ! 1193 !-- Compute length of time step 1194 IF ( call_microphysics_at_all_substeps ) THEN 1195 dt_micro = dt_3d * weight_pres(intermediate_timestep_count) 1196 ELSE 1197 dt_micro = dt_3d 1198 ENDIF 1199 1200 ! 1201 !-- Reset precipitation rate 1202 IF ( intermediate_timestep_count == 1 ) prr = 0.0_wp 1203 1204 ! 1205 !-- Compute cloud physics 1206 IF ( microphysics_kessler ) THEN 1207 1208 CALL autoconversion_kessler 1209 IF ( cloud_water_sedimentation ) CALL sedimentation_cloud 1210 1211 ELSEIF ( microphysics_seifert ) THEN 1212 1213 CALL adjust_cloud 1214 IF ( microphysics_morrison ) CALL activation 1215 IF ( microphysics_morrison ) CALL condensation 1216 CALL autoconversion 1217 CALL accretion 1218 CALL selfcollection_breakup 1219 CALL evaporation_rain 1220 CALL sedimentation_rain 1221 IF ( cloud_water_sedimentation ) CALL sedimentation_cloud 1222 1223 ENDIF 1224 1225 CALL calc_precipitation_amount 1226 1227 END SUBROUTINE bcm_actions_micro 1238 CALL cpu_log( log_point(51), 'microphysics', 'stop' ) 1239 1240 END SUBROUTINE bcm_non_transport_physics 1228 1241 1229 1242 … … 1234 1247 !------------------------------------------------------------------------------! 1235 1248 1236 SUBROUTINE bcm_actions_micro_ij( i, j ) 1237 1238 IMPLICIT NONE 1249 SUBROUTINE bcm_non_transport_physics_ij( i, j ) 1250 1239 1251 1240 1252 INTEGER(iwp) :: i !< 1241 1253 INTEGER(iwp) :: j !< 1242 1254 1243 IF ( large_scale_forcing .AND. lsf_surf ) THEN 1244 ! 1245 !-- Calculate vertical profile of the hydrostatic pressure (hyp) 1246 hyp = barometric_formula(zu, pt_surface * exner_function(surface_pressure * 100.0_wp), surface_pressure * 100.0_wp) 1247 d_exner = exner_function_invers(hyp) 1248 exner = 1.0_wp / exner_function_invers(hyp) 1249 hyrho = ideal_gas_law_rho_pt(hyp, pt_init) 1250 ! 1251 !-- Compute reference density 1252 rho_surface = ideal_gas_law_rho(surface_pressure * 100.0_wp, pt_surface * exner_function(surface_pressure * 100.0_wp)) 1255 IF ( .NOT. microphysics_sat_adjust .AND. & 1256 ( intermediate_timestep_count == 1 .OR. & 1257 call_microphysics_at_all_substeps ) ) & 1258 THEN 1259 1260 IF ( large_scale_forcing .AND. lsf_surf ) THEN 1261 ! 1262 !-- Calculate vertical profile of the hydrostatic pressure (hyp) 1263 hyp = barometric_formula(zu, pt_surface * exner_function(surface_pressure * 100.0_wp), surface_pressure * 100.0_wp) 1264 d_exner = exner_function_invers(hyp) 1265 exner = 1.0_wp / exner_function_invers(hyp) 1266 hyrho = ideal_gas_law_rho_pt(hyp, pt_init) 1267 ! 1268 !-- Compute reference density 1269 rho_surface = ideal_gas_law_rho(surface_pressure * 100.0_wp, pt_surface * exner_function(surface_pressure * 100.0_wp)) 1270 ENDIF 1271 1272 ! 1273 !-- Compute length of time step 1274 IF ( call_microphysics_at_all_substeps ) THEN 1275 dt_micro = dt_3d * weight_pres(intermediate_timestep_count) 1276 ELSE 1277 dt_micro = dt_3d 1278 ENDIF 1279 ! 1280 !-- Reset precipitation rate 1281 IF ( intermediate_timestep_count == 1 ) prr(:,j,i) = 0.0_wp 1282 1283 ! 1284 !-- Compute cloud physics 1285 IF( microphysics_kessler ) THEN 1286 1287 CALL autoconversion_kessler_ij( i,j ) 1288 IF ( cloud_water_sedimentation ) CALL sedimentation_cloud_ij( i,j ) 1289 1290 ELSEIF ( microphysics_seifert ) THEN 1291 1292 CALL adjust_cloud_ij( i,j ) 1293 IF ( microphysics_morrison ) CALL activation_ij( i,j ) 1294 IF ( microphysics_morrison ) CALL condensation_ij( i,j ) 1295 CALL autoconversion_ij( i,j ) 1296 CALL accretion_ij( i,j ) 1297 CALL selfcollection_breakup_ij( i,j ) 1298 CALL evaporation_rain_ij( i,j ) 1299 CALL sedimentation_rain_ij( i,j ) 1300 IF ( cloud_water_sedimentation ) CALL sedimentation_cloud_ij( i,j ) 1301 1302 ENDIF 1303 1304 CALL calc_precipitation_amount_ij( i,j ) 1305 1253 1306 ENDIF 1254 1307 1255 ! 1256 !-- Compute length of time step 1257 IF ( call_microphysics_at_all_substeps ) THEN 1258 dt_micro = dt_3d * weight_pres(intermediate_timestep_count) 1259 ELSE 1260 dt_micro = dt_3d 1261 ENDIF 1262 ! 1263 !-- Reset precipitation rate 1264 IF ( intermediate_timestep_count == 1 ) prr(:,j,i) = 0.0_wp 1265 1266 ! 1267 !-- Compute cloud physics 1268 IF( microphysics_kessler ) THEN 1269 1270 CALL autoconversion_kessler_ij( i,j ) 1271 IF ( cloud_water_sedimentation ) CALL sedimentation_cloud_ij( i,j ) 1272 1273 ELSEIF ( microphysics_seifert ) THEN 1274 1275 CALL adjust_cloud_ij( i,j ) 1276 IF ( microphysics_morrison ) CALL activation_ij( i,j ) 1277 IF ( microphysics_morrison ) CALL condensation_ij( i,j ) 1278 CALL autoconversion_ij( i,j ) 1279 CALL accretion_ij( i,j ) 1280 CALL selfcollection_breakup_ij( i,j ) 1281 CALL evaporation_rain_ij( i,j ) 1282 CALL sedimentation_rain_ij( i,j ) 1283 IF ( cloud_water_sedimentation ) CALL sedimentation_cloud_ij( i,j ) 1284 1285 ENDIF 1286 1287 CALL calc_precipitation_amount_ij( i,j ) 1288 1289 END SUBROUTINE bcm_actions_micro_ij 1308 END SUBROUTINE bcm_non_transport_physics_ij 1290 1309 1291 1310 -
palm/trunk/SOURCE/module_interface.f90
r3873 r3874 147 147 bcm_header, & 148 148 bcm_actions, & 149 bcm_non_transport_physics, & 149 150 bcm_prognostic_equations, & 150 151 bcm_swap_timelevel, & … … 382 383 module_interface_header, & 383 384 module_interface_actions, & 385 module_interface_non_transport_physics, & 384 386 module_interface_prognostic_equations, & 385 387 module_interface_swap_timelevel, & … … 443 445 MODULE PROCEDURE module_interface_actions_ij 444 446 END INTERFACE module_interface_actions 447 448 INTERFACE module_interface_non_transport_physics 449 MODULE PROCEDURE module_interface_non_transport_physics 450 MODULE PROCEDURE module_interface_non_transport_physics_ij 451 END INTERFACE module_interface_non_transport_physics 445 452 446 453 INTERFACE module_interface_prognostic_equations … … 872 879 ! Description: 873 880 ! ------------ 874 !> Perform module-specific actions while in time-integration ( cache-optimized)881 !> Perform module-specific actions while in time-integration (vector-optimized) 875 882 !------------------------------------------------------------------------------! 876 883 SUBROUTINE module_interface_actions( location ) … … 893 900 ! Description: 894 901 ! ------------ 895 !> Perform module-specific actions while in time-integration ( vector-optimized)902 !> Perform module-specific actions while in time-integration (cache-optimized) 896 903 !------------------------------------------------------------------------------! 897 904 SUBROUTINE module_interface_actions_ij( i, j, location ) … … 916 923 ! Description: 917 924 ! ------------ 925 !> Compute module-specific prognostic_equations (vector-optimized) 926 !------------------------------------------------------------------------------! 927 SUBROUTINE module_interface_non_transport_physics() 928 929 930 IF ( bulk_cloud_model ) CALL bcm_non_transport_physics() 931 932 933 END SUBROUTINE module_interface_non_transport_physics 934 935 936 !------------------------------------------------------------------------------! 937 ! Description: 938 ! ------------ 918 939 !> Compute module-specific prognostic_equations (cache-optimized) 940 !------------------------------------------------------------------------------! 941 SUBROUTINE module_interface_non_transport_physics_ij( i, j ) 942 943 944 INTEGER(iwp), INTENT(IN) :: i !< grid index in x-direction 945 INTEGER(iwp), INTENT(IN) :: j !< grid index in y-direction 946 947 948 IF ( bulk_cloud_model ) CALL bcm_non_transport_physics( i, j ) 949 950 951 END SUBROUTINE module_interface_non_transport_physics_ij 952 953 954 !------------------------------------------------------------------------------! 955 ! Description: 956 ! ------------ 957 !> Compute module-specific prognostic_equations (vector-optimized) 919 958 !------------------------------------------------------------------------------! 920 959 SUBROUTINE module_interface_prognostic_equations() … … 933 972 ! Description: 934 973 ! ------------ 935 !> Compute module-specific prognostic_equations ( vector-optimized)974 !> Compute module-specific prognostic_equations (cache-optimized) 936 975 !------------------------------------------------------------------------------! 937 976 SUBROUTINE module_interface_prognostic_equations_ij( i, j, i_omp_start, tn ) -
palm/trunk/SOURCE/prognostic_equations.f90
r3872 r3874 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Added non_transport_physics module interfaces and moved bcm code into it 28 ! 29 ! 3872 2019-04-08 15:03:06Z knoop 27 30 ! Moving prognostic equations of bcm into bulk_cloud_model_mod 28 31 ! … … 378 381 ug, u_init, u_p, v, vg, vpt, v_init, v_p, w, w_p 379 382 380 USE bulk_cloud_model_mod, &381 ONLY: call_microphysics_at_all_substeps, bulk_cloud_model, &382 bcm_actions_micro, microphysics_sat_adjust383 384 383 USE buoyancy_mod, & 385 384 ONLY: buoyancy … … 441 440 USE module_interface, & 442 441 ONLY: module_interface_actions, & 442 module_interface_non_transport_physics, & 443 443 module_interface_prognostic_equations 444 444 … … 478 478 ONLY: wtm_tendencies 479 479 480 IMPLICIT NONE 480 481 481 482 PRIVATE … … 510 511 511 512 512 IMPLICIT NONE513 514 513 INTEGER(iwp) :: i !< 515 514 INTEGER(iwp) :: i_omp_start !< … … 649 648 ENDIF 650 649 651 ENDIF 652 653 ! 654 !-- If required, calculate cloud microphysics 655 IF ( bulk_cloud_model .AND. .NOT. microphysics_sat_adjust .AND. & 656 ( intermediate_timestep_count == 1 .OR. & 657 call_microphysics_at_all_substeps ) ) & 658 THEN 659 !$OMP PARALLEL PRIVATE (i,j) 660 !$OMP DO 661 DO i = nxlg, nxrg 662 DO j = nysg, nyng 663 CALL bcm_actions_micro( i, j ) 664 ENDDO 650 ENDIF 651 652 !$OMP PARALLEL PRIVATE (i,j) 653 !$OMP DO 654 DO i = nxlg, nxrg 655 DO j = nysg, nyng 656 ! 657 !-- Calculate non transport physics for all other modules 658 CALL module_interface_non_transport_physics( i, j ) 665 659 ENDDO 666 !$OMP END PARALLEL667 ENDIF660 ENDDO 661 !$OMP END PARALLEL 668 662 669 663 ! … … 1238 1232 1239 1233 1240 IMPLICIT NONE1241 1242 1234 INTEGER(iwp) :: i !< 1243 1235 INTEGER(iwp) :: ib !< index for aerosol size bins (salsa) … … 1299 1291 ENDIF 1300 1292 ENDIF 1301 ENDIF1302 1303 !1304 !-- If required, calculate cloud microphysical impacts1305 IF ( bulk_cloud_model .AND. .NOT. microphysics_sat_adjust .AND. &1306 ( intermediate_timestep_count == 1 .OR. &1307 call_microphysics_at_all_substeps ) &1308 ) THEN1309 CALL cpu_log( log_point(51), 'microphysics', 'start' )1310 CALL bcm_actions_micro1311 CALL cpu_log( log_point(51), 'microphysics', 'stop' )1312 1293 ENDIF 1294 1295 ! 1296 !-- Calculate non transport physics for all other modules 1297 CALL module_interface_non_transport_physics 1313 1298 1314 1299 !
Note: See TracChangeset
for help on using the changeset viewer.
