Changeset 4881 for palm/trunk
- Timestamp:
- Feb 19, 2021 10:05:08 PM (3 years ago)
- Location:
- palm/trunk/SOURCE
- Files:
-
- 4 edited
-
chem_modules.f90 (modified) (3 diffs)
-
chem_photolysis_mod.f90 (modified) (2 diffs)
-
chemistry_model_mod.f90 (modified) (11 diffs)
-
surface_mod.f90 (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/chem_modules.f90
r4828 r4881 21 21 ! Current revisions: 22 22 ! ----------------- 23 ! 24 ! 23 ! 24 ! 25 25 ! Former revisions: 26 26 ! ----------------- 27 27 ! $Id$ 28 ! removed unnecessarty namelist parameters 29 ! 30 ! 31 ! 4828 2021-01-05 11:21:41Z Giersch 28 32 ! Enable output of vertical fluxes of chemical species. 29 33 ! … … 224 228 LOGICAL :: bc_radiation_cs_s = .FALSE. !< flag for indicating a radiation/neumann 225 229 !< condition at the south boundary 226 LOGICAL :: constant_top_csflux(99) = .TRUE. !< internal flag, set to .FALSE. if no227 !< top_csflux is prescribed228 230 LOGICAL :: constant_csflux(99) = .TRUE. !< internal flag, set to .FALSE. if no 229 231 !< surface_csflux is prescribed 230 LOGICAL :: call_chem_at_all_substeps = .FALSE. !< namelist parameter: ....??? 231 LOGICAL :: chem_debug0 = .FALSE. !< namelist parameter: flag for minimum print 232 !< output 233 LOGICAL :: chem_debug1 = .FALSE. !< namelist parameter: flag for print output 234 LOGICAL :: chem_debug2 = .FALSE. !< namelist parameter: flag for further print 235 !< output 232 LOGICAL :: call_chem_at_all_substeps = .FALSE. !< namelist parameter: Never set true 233 !< except for tests 236 234 LOGICAL :: chem_gasphase_on = .TRUE. !< namelist parameter: flag to switch off 237 235 !< chemical reactions … … 248 246 !< model 249 247 250 REAL(wp) :: cs_surface_initial_change(99) = 0.0_wp !< namelist parameter: ...??? 251 REAL(wp) :: cs_vertical_gradient(99,10) = 0.0_wp !< namelist parameter: ...??? 252 REAL(wp) :: cs_vertical_gradient_level(99,10) = -999999.9_wp !< namelist parameter: ...??? 253 REAL(wp) :: emiss_factor_main(99) = -9999.0_wp !< namelist parameter: ...??? 254 REAL(wp) :: emiss_factor_side(99) = -9999.0_wp !< namelist parameter: ...??? 255 REAL(wp) :: surface_csflux(99) = 0.0_wp !< namelist parameter: ...??? 256 REAL(wp) :: top_csflux(99) = 0.0_wp !< namelist parameter: ...??? 257 REAL(wp) :: wall_csflux(99,0:5) = 0.0_wp !< namelist parameter: ...??? 248 REAL(wp) :: emiss_factor_main(99) = -9999.0_wp !< namelist parameter: weighting factor 249 !< for LOD 0 parameterized emissions 250 REAL(wp) :: emiss_factor_side(99) = -9999.0_wp !< namelist parameter: weighting factor 251 !< for LOD 0 parameterized emissions 252 REAL(wp) :: surface_csflux(99) = 0.0_wp !< namelist parameter: surface emission 253 !< flux or LOD 0 parameterized emissions 254 REAL(wp) :: wall_csflux(99,0:5) = 0.0_wp !< namelist parameter: ...??? 258 255 259 256 REAL(wp), DIMENSION(99) :: cs_surface = 0.0_wp !< namelist parameter: chem species -
palm/trunk/SOURCE/chem_photolysis_mod.f90
r4828 r4881 20 20 ! Current revisions: 21 21 ! ----------------- 22 ! 23 ! 22 ! 23 ! 24 24 ! Former revisions: 25 25 ! ----------------- 26 26 ! $Id$ 27 ! removed chem_debug parameter 28 ! 29 ! 30 ! 4828 2021-01-05 11:21:41Z Giersch 27 31 ! further re-formatting to follow the PALM coding standard 28 32 ! … … 78 82 USE chem_modules, & 79 83 ONLY: phot_frequen, photolysis_scheme 80 81 USE chem_modules, &82 ONLY: chem_debug283 84 84 85 USE kinds -
palm/trunk/SOURCE/chemistry_model_mod.f90
r4860 r4881 26 26 ! ----------------- 27 27 ! $Id$ 28 ! removed unnecessarty namelist parameters and commented output statements 29 ! and cs_surface_initial_change 30 ! 31 ! 32 ! 4860 2021-02-01 08:10:59Z raasch 28 33 ! further re-numbering of message IDs 29 34 ! … … 297 302 ! ------------ 298 303 !> Chemistry model for PALM-4U 299 !> @todo Extend chem_species type by nspec and nvar as addititional elements (RF)300 !> @todo Check possibility to reduce dimension of chem_species%conc from nspec to nvar (RF)301 304 !> @todo Adjust chem_rrd_local to CASE structure of others modules. It is not allowed to use the 302 305 !> chemistry model in a precursor run and additionally not using it in a main run … … 308 311 !> @todo slight differences in passive scalar and chem spcs when chem reactions turned off. Need to 309 312 !> be fixed. bK 310 !> @todo test nesting for chem spcs, was implemented by suehring (kanani)311 313 !> @todo chemistry error messages 312 314 ! … … 1521 1523 (variable(char_len-2:) == '_xz') .OR. & 1522 1524 (variable(char_len-2:) == '_yz') ) ) THEN 1523 !1524 !-- todo: remove or replace by "CALL message" mechanism (kanani)1525 ! IF(myid == 0) WRITE(6,*) 'Output of species ' // TRIM( variable ) // &1526 ! TRIM( chem_species(lsp)%name )1527 1525 IF (av == 0) THEN 1528 1526 DO i = nxl, nxr … … 1638 1636 DO lsp = 1, nspec 1639 1637 IF ( TRIM( spec_name ) == TRIM( chem_species(lsp)%name) ) THEN 1640 !1641 !-- todo: remove or replace by "CALL message" mechanism (kanani)1642 ! IF(myid == 0 .AND. chem_debug0 ) WRITE(6,*) 'Output of species ' // TRIM( variable ) // &1643 ! TRIM( chem_species(lsp)%name )1644 1638 IF (av == 0) THEN 1645 1639 DO i = nxl, nxr … … 1717 1711 DO lsp=1,nspec 1718 1712 IF (TRIM( spec_name ) == TRIM( chem_species(lsp)%name) ) THEN 1719 !1720 !-- todo: remove or replace by "CALL message" mechanism (kanani)1721 ! IF(myid == 0 .AND. chem_debug0 ) WRITE(6,*) 'Output of species ' // TRIM( variable ) // &1722 ! TRIM( chem_species(lsp)%name )1723 1713 IF (av == 0) THEN 1724 1714 IF ( .NOT. mask_surface(mid) ) THEN … … 2272 2262 ENDIF 2273 2263 ! 2274 !-- If required, change the surface chem spcs at the start of the 3D run2275 IF ( cs_surface_initial_change(1) /= 0.0_wp ) THEN2276 DO lsp = 1, nspec2277 DO i = nxlg, nxrg2278 DO j = nysg, nyng2279 flag = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_total_0(nzb,j,i), 0 ) )2280 chem_species(lsp)%conc(nzb,j,i) = chem_species(lsp)%conc(nzb,j,i) + &2281 cs_surface_initial_change(lsp) * flag2282 ENDDO2283 ENDDO2284 ENDDO2285 ENDIF2286 2287 2264 ENDIF 2288 2265 ! … … 2295 2272 DO lsp = 1, nphot 2296 2273 phot_frequen(lsp)%name = phot_names(lsp) 2297 !2298 !-- todo: remove or replace by "CALL message" mechanism (kanani)2299 !-- IF( myid == 0 ) THEN2300 !-- WRITE(6,'(a,i4,3x,a)') 'Photolysis: ',lsp,TRIM( phot_names(lsp) )2301 !-- ENDIF2302 2274 phot_frequen(lsp)%freq(nzb:nzt+1,nysg:nyng,nxlg:nxrg) => freq_1(:,:,:,lsp) 2303 2275 ENDDO … … 2546 2518 bc_cs_t, & 2547 2519 call_chem_at_all_substeps, & 2548 chem_debug0, &2549 chem_debug1, &2550 chem_debug2, &2551 2520 chem_gasphase_on, & 2552 2521 chem_mechanism, & … … 2555 2524 cs_profile, & 2556 2525 cs_surface, & 2557 cs_surface_initial_change, &2558 cs_vertical_gradient_level, &2559 2526 daytype_mdh, & 2560 2527 deposition_dry, & … … 2575 2542 photolysis_scheme, & 2576 2543 wall_csflux, & 2577 cs_vertical_gradient, &2578 top_csflux, &2579 2544 surface_csflux, & 2580 2545 surface_csflux_name, & -
palm/trunk/SOURCE/surface_mod.f90
r4877 r4881 25 25 ! ----------------- 26 26 ! $Id$ 27 ! removed constant_top_csflux option 28 ! 29 ! 30 ! 4877 2021-02-17 16:17:35Z suehring 27 31 ! bugfix in init_single_surface_properties for vertical surfaces 28 32 ! … … 2928 2932 IF ( passive_scalar .AND. constant_top_scalarflux ) surf%ssws(num_h) = top_scalarflux * & 2929 2933 rho_air_zw(nzt+1) 2930 !2931 !-- Prescribe top chemical species' flux2932 DO lsp = 1, nvar2933 IF ( air_chemistry .AND. constant_top_csflux(lsp) ) THEN2934 surf%cssws(lsp,num_h) = top_csflux(lsp) * rho_air_zw(nzt+1)2935 ENDIF2936 ENDDO2937 2934 ! 2938 2935 !-- Prescribe top salinity flux
Note: See TracChangeset
for help on using the changeset viewer.
