Changeset 3274 for palm/trunk/SOURCE/surface_mod.f90
- Timestamp:
- Sep 24, 2018 3:42:55 PM (6 years ago)
- File:
-
- 1 edited
-
palm/trunk/SOURCE/surface_mod.f90 (modified) (21 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/surface_mod.f90
r3254 r3274 26 26 ! ----------------- 27 27 ! $Id$ 28 ! Modularization of all bulk cloud physics code components 29 ! 30 ! 3254 2018-09-17 10:53:57Z suehring 28 31 ! Remove redundant subroutine argument 29 32 ! … … 207 210 USE chem_modules 208 211 209 USE control_parameters 212 USE control_parameters 210 213 211 214 USE indices, & … … 505 508 !< in the domain (required to activiate RTM) 506 509 510 LOGICAL :: surf_bulk_cloud_model = .FALSE. !< use cloud microphysics 511 LOGICAL :: surf_microphysics_morrison = .FALSE. !< use 2-moment Morrison (add. prog. eq. for nc and qc) 512 LOGICAL :: surf_microphysics_seifert = .FALSE. !< use 2-moment Seifert and Beheng scheme 513 507 514 508 515 SAVE … … 543 550 PUBLIC bc_h, ind_pav_green, ind_veg_wall, ind_wat_win, ns_h_on_file, & 544 551 ns_v_on_file, surf_def_h, surf_def_v, surf_lsm_h, surf_lsm_v, & 545 surf_usm_h, surf_usm_v, surf_type, vertical_surfaces_exist 552 surf_usm_h, surf_usm_v, surf_type, vertical_surfaces_exist, & 553 surf_bulk_cloud_model, surf_microphysics_morrison, surf_microphysics_seifert 546 554 ! 547 555 !-- Public subroutines and functions … … 1090 1098 ! 1091 1099 !-- 1092 IF ( cloud_physics .AND.microphysics_morrison) THEN1100 IF ( surf_bulk_cloud_model .AND. surf_microphysics_morrison) THEN 1093 1101 DEALLOCATE ( surfaces%qcs ) 1094 1102 DEALLOCATE ( surfaces%ncs ) … … 1098 1106 ! 1099 1107 !-- 1100 IF ( cloud_physics .AND.microphysics_seifert) THEN1108 IF ( surf_bulk_cloud_model .AND. surf_microphysics_seifert) THEN 1101 1109 DEALLOCATE ( surfaces%qrs ) 1102 1110 DEALLOCATE ( surfaces%nrs ) … … 1217 1225 ! 1218 1226 !-- 1219 IF ( cloud_physics .AND.microphysics_morrison) THEN1227 IF ( surf_bulk_cloud_model .AND. surf_microphysics_morrison) THEN 1220 1228 ALLOCATE ( surfaces%qcs(1:surfaces%ns) ) 1221 1229 ALLOCATE ( surfaces%ncs(1:surfaces%ns) ) … … 1225 1233 ! 1226 1234 !-- 1227 IF ( cloud_physics .AND.microphysics_seifert) THEN1235 IF ( surf_bulk_cloud_model .AND. surf_microphysics_seifert) THEN 1228 1236 ALLOCATE ( surfaces%qrs(1:surfaces%ns) ) 1229 1237 ALLOCATE ( surfaces%nrs(1:surfaces%ns) ) … … 1286 1294 ! 1287 1295 !-- 1288 IF ( cloud_physics .AND.microphysics_morrison) THEN1296 IF ( surf_bulk_cloud_model .AND. surf_microphysics_morrison) THEN 1289 1297 DEALLOCATE ( surfaces%qcsws ) 1290 1298 DEALLOCATE ( surfaces%ncsws ) … … 1292 1300 ! 1293 1301 !-- 1294 IF ( cloud_physics .AND.microphysics_seifert) THEN1302 IF ( surf_bulk_cloud_model .AND. surf_microphysics_seifert) THEN 1295 1303 DEALLOCATE ( surfaces%qrsws ) 1296 1304 DEALLOCATE ( surfaces%nrsws ) … … 1358 1366 ! 1359 1367 !-- 1360 IF ( cloud_physics .AND.microphysics_morrison) THEN1368 IF ( surf_bulk_cloud_model .AND. surf_microphysics_morrison) THEN 1361 1369 ALLOCATE ( surfaces%qcsws(1:surfaces%ns) ) 1362 1370 ALLOCATE ( surfaces%ncsws(1:surfaces%ns) ) … … 1364 1372 ! 1365 1373 !-- 1366 IF ( cloud_physics .AND.microphysics_seifert) THEN1374 IF ( surf_bulk_cloud_model .AND. surf_microphysics_seifert) THEN 1367 1375 ALLOCATE ( surfaces%qrsws(1:surfaces%ns) ) 1368 1376 ALLOCATE ( surfaces%nrsws(1:surfaces%ns) ) … … 1473 1481 DEALLOCATE ( surfaces%vpt1 ) 1474 1482 1475 IF ( cloud_physics .AND.microphysics_seifert) THEN1483 IF ( surf_bulk_cloud_model .AND. surf_microphysics_seifert) THEN 1476 1484 DEALLOCATE ( surfaces%qcs ) 1477 1485 DEALLOCATE ( surfaces%ncs ) … … 1480 1488 ENDIF 1481 1489 1482 IF ( cloud_physics .AND.microphysics_seifert) THEN1490 IF ( surf_bulk_cloud_model .AND. surf_microphysics_seifert) THEN 1483 1491 DEALLOCATE ( surfaces%qrs ) 1484 1492 DEALLOCATE ( surfaces%nrs ) … … 1598 1606 ALLOCATE ( surfaces%vpt1(1:surfaces%ns) ) 1599 1607 1600 IF ( cloud_physics .AND.microphysics_seifert) THEN1608 IF ( surf_bulk_cloud_model .AND. surf_microphysics_seifert) THEN 1601 1609 ALLOCATE ( surfaces%qcs(1:surfaces%ns) ) 1602 1610 ALLOCATE ( surfaces%ncs(1:surfaces%ns) ) … … 1605 1613 ENDIF 1606 1614 1607 IF ( cloud_physics .AND.microphysics_seifert) THEN1615 IF ( surf_bulk_cloud_model .AND. surf_microphysics_seifert) THEN 1608 1616 ALLOCATE ( surfaces%qrs(1:surfaces%ns) ) 1609 1617 ALLOCATE ( surfaces%nrs(1:surfaces%ns) ) … … 2133 2141 IF ( humidity ) THEN 2134 2142 surf%qs(num_h) = 0.0_wp 2135 IF ( cloud_physics .AND.microphysics_morrison) THEN2143 IF ( surf_bulk_cloud_model .AND. surf_microphysics_morrison) THEN 2136 2144 surf%qcs(num_h) = 0.0_wp 2137 2145 surf%ncs(num_h) = 0.0_wp … … 2141 2149 2142 2150 ENDIF 2143 IF ( cloud_physics .AND.microphysics_seifert) THEN2151 IF ( surf_bulk_cloud_model .AND. surf_microphysics_seifert) THEN 2144 2152 surf%qrs(num_h) = 0.0_wp 2145 2153 surf%nrs(num_h) = 0.0_wp … … 2314 2322 IF ( humidity ) THEN 2315 2323 surf%qsws(num_h) = 0.0_wp 2316 IF ( cloud_physics .AND.microphysics_morrison ) THEN2324 IF ( surf_bulk_cloud_model .AND. surf_microphysics_morrison ) THEN 2317 2325 surf%ncsws(num_h) = 0.0_wp 2318 2326 surf%qcsws(num_h) = 0.0_wp 2319 2327 ENDIF 2320 IF ( cloud_physics .AND.microphysics_seifert ) THEN2328 IF ( surf_bulk_cloud_model .AND. surf_microphysics_seifert ) THEN 2321 2329 surf%nrsws(num_h) = 0.0_wp 2322 2330 surf%qrsws(num_h) = 0.0_wp … … 2449 2457 ! 2450 2458 !-- Following wall fluxes are assumed to be zero 2451 IF ( cloud_physics .AND.microphysics_morrison) THEN2459 IF ( surf_bulk_cloud_model .AND. surf_microphysics_morrison) THEN 2452 2460 surf%qcs(num_v) = 0.0_wp 2453 2461 surf%ncs(num_v) = 0.0_wp … … 2456 2464 surf%ncsws(num_v) = 0.0_wp 2457 2465 ENDIF 2458 IF ( cloud_physics .AND.microphysics_seifert) THEN2466 IF ( surf_bulk_cloud_model .AND. surf_microphysics_seifert) THEN 2459 2467 surf%qrs(num_v) = 0.0_wp 2460 2468 surf%nrs(num_v) = 0.0_wp
Note: See TracChangeset
for help on using the changeset viewer.
