Ignore:
Timestamp:
Dec 17, 2018 2:52:50 PM (6 years ago)
Author:
kanani
Message:

Fix for masked output, change default chemistry mechanism

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/TESTS/cases/urban_environment/INPUT/urban_environment_p3d

    r3421 r3632  
    3030!-- Initialization
    3131!------------------------------------------------------------------------------
    32     initializing_actions = 'inifor',   !'set_constant_profiles', or 'inifor',
     32    initializing_actions = 'inifor set_constant_profiles',   !'set_constant_profiles', or 'inifor',
     33                                                             ! or both (met. from inifor, chem constant profiles)
    3334
    3435    pt_surface           = 277.15,
     
    9091    data_output_pr        = '#u',  'u*2',  'wu',  'w*u*',  'w"u"',
    9192                            '#v',  'v*2',  'wv',  'w*v*',  'w"v"',
    92                             'w',   'w*2',
     93!                            'w',   'w*2',
    9394                            '#theta', 'theta*2', 'wtheta', 'w*theta*', 'w"theta"',
    9495                            '#q',  'q*2',
    9596                            'e', 'e*', '#km', '#l',
    9697                            '#t_soil', '#m_soil', 'rad_lw_in', 'rad_lw_out',
     98                            'kc_NO', 'kc_NO2', 'kc_O3', 'kc_PM10',
    9799!
    98100!-- 2D/3D output settings
     
    127129                  'tsurf*_xy', 'tsurf*_xy_av',
    128130                  'pcm_lad', 'pcm_heatrate',
    129                   'kc_PM10', 'kc_PM25', 'kc_PM10_av', 'kc_PM25_av',
    130 
    131 
     131                  'kc_NO', 'kc_NO2', 'kc_O3', 'kc_PM10',
     132                  'kc_NO_av', 'kc_NO2_av', 'kc_O3_av', 'kc_PM10_av',
    132133
    133134 / ! end of d3par namelist
     
    206207 &chemistry_parameters
    207208
    208     chem_gasphase_on           = .F.,
     209    chem_gasphase_on           = .T.,
     210    call_chem_at_all_substeps  = .FALSE.,
    209211   
    210     do_emis                    = .T.,
    211     mode_emis                  = 'PARAMETERIZED',
    212     surface_csflux_name        = 'PM10','PM25',  !"cs" stands for chemical species
    213     surface_csflux             = 0.75, 0.375,    !surface emission of PM10, PM25 (arbitrary units)
    214     cs_name                    = 'PM10','PM25',
    215     cs_surface                 = 2.0, 1.0,       !surface concentration of PM10, PM25 (arbitrary units)
     212    photolysis_scheme          = 'simple',
     213
     214    cs_name                    = 'NO', 'NO2', 'O3', 'PM10',
     215    cs_surface                 = 0.02,  0.03, 0.02, 25.0e-9,   ! surface concentration of PM10, PM25 in kg m-3,
     216                                                               ! units for gases in ppm   
     217!-- initial profiles (also for lateral boundaries if 'set_constant_profiles' is chosen)
     218!-- units for PM: kg m-3, units for gases: ppm
     219    cs_profile(1,:)            =  0.020, 0.020, 0.015, 0.015, 0.015, 0.010, 0.010, 0.010, 0.010, 0.010, 0.007,
     220    cs_profile(2,:)            =  0.030, 0.030, 0.020, 0.020, 0.020, 0.015, 0.015, 0.015, 0.015, 0.015, 0.010,
     221    cs_profile(3,:)            =  0.020, 0.020, 0.050, 0.050, 0.050, 0.050, 0.050, 0.050, 0.050, 0.050, 0.050,
     222    cs_profile(1,:)            =  25.e-09, 25.e-09, 23.0e-09, 23.e-09, 23.e-09, 23.e-09, 23.e-09, 23.e-09, 22.e-09, 22.e-09, 22.0e-09,
     223
     224    cs_heights(1,:)            =  0.0,   5.0,   15.0,   25.0,  35.0,  45.0,  55.0,  65.0,  75.0,  85.0, 95.0,
     225    cs_heights(2,:)            =  0.0,   5.0,   15.0,   25.0,  35.0,  45.0,  55.0,  65.0,  75.0,  85.0, 95.0,
     226    cs_heights(3,:)            =  0.0,   5.0,   15.0,   25.0,  35.0,  45.0,  55.0,  65.0,  75.0,  85.0, 95.0,
     227    cs_heights(4,:)            =  0.0,   5.0,   15.0,   25.0,  35.0,  45.0,  55.0,  65.0,  75.0,  85.0, 95.0,
     228
    216229    bc_cs_b                    = 'neumann',
    217230    bc_cs_t                    = 'neumann',
    218231   
    219 !-- for emission depending on street type in static driver
    220     emiss_factor_main = 1.667, 1.667,   !emiss_factor_main*surface_csflux for emission on main streets
    221     emiss_factor_side = 0.334, 0.334,   !emiss_factor_side*surface_csflux for emission on side streets
    222     main_street_id = 12,          !id is for identifying main/side streets (see chem_emissions in chemistry_model_mod.f90)
     232    do_emis                    = .T.,
     233    mode_emis                  = 'PARAMETERIZED',  ! emission depending on street type from static driver
     234    surface_csflux_name        = 'NO',   'NO2',  'PM10',   !"cs" stands for chemical species
     235    surface_csflux             = 4745.0, 1326.0, 2700.e-09,! surface emission of PM10, PM25 in kg m-2 d-1,
     236                                                           ! for gases in umol m-2 d-1 for PARAMETERIZED
     237    emiss_factor_main          = 1.667,  1.667,  1.667,    ! emiss_factor_main*surface_csflux for emission on main streets
     238    emiss_factor_side          = 0.334,  0.334,  0.334,    !emiss_factor_side*surface_csflux for emission on side streets
     239    main_street_id = 12,       ! id is for identifying main/side streets (see chem_emissions in chemistry_model_mod.f90)
    223240    side_street_id =  8,
    224     max_street_id = 19,
     241    max_street_id  = 19,
     242
     243! -- KPP solver steering. For details  see comments in chem_gasphase_mod.f90, SUBROUTINE rosenbrock
     244!                         (or other solver subroutine, if another solver is chosen)
     245    icntrl(3)                  = 1,   ! solver ros2
     246    icntrl(4)                  = 500, ! max. number of chem-substeps
     247    rcntrl(3)                  = 0.1, ! Hstart,starting value for the integration step size
    225248
    226249 / ! end of chemistry_par namelist
Note: See TracChangeset for help on using the changeset viewer.