Ignore:
Timestamp:
Dec 17, 2018 2:52:50 PM (5 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_restart/INPUT/urban_environment_restart_p3dr

    r3421 r3632  
    9090    data_output_pr        = '#u',  'u*2',  'wu',  'w*u*',  'w"u"',
    9191                            '#v',  'v*2',  'wv',  'w*v*',  'w"v"',
    92                             'w',   'w*2',
     92!                            'w',   'w*2',
    9393                            '#theta', 'theta*2', 'wtheta', 'w*theta*', 'w"theta"',
    9494                            '#q',  'q*2',
    9595                            'e', 'e*', '#km', '#l',
    9696                            '#t_soil', '#m_soil', 'rad_lw_in', 'rad_lw_out',
     97                            'kc_NO', 'kc_NO2', 'kc_O3', 'kc_PM10',
    9798!
    9899!-- 2D/3D output settings
     
    127128                  'tsurf*_xy', 'tsurf*_xy_av',
    128129                  'pcm_lad', 'pcm_heatrate',
    129                   'kc_PM10', 'kc_PM25', 'kc_PM10_av', 'kc_PM25_av',
    130 
    131 
     130                  'kc_NO', 'kc_NO2', 'kc_O3', 'kc_PM10',
     131                  'kc_NO_av', 'kc_NO2_av', 'kc_O3_av', 'kc_PM10_av',
    132132
    133133 / ! end of d3par namelist
     
    206206 &chemistry_parameters
    207207
    208     chem_gasphase_on           = .F.,
     208    chem_gasphase_on           = .T.,
     209    call_chem_at_all_substeps  = .FALSE.,
    209210   
    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)
     211    photolysis_scheme          = 'simple',
     212
     213    cs_name                    = 'NO', 'NO2', 'O3', 'PM10',
     214    cs_surface                 = 0.02,  0.03, 0.02, 25.0e-9,   ! surface concentration of PM10, PM25 in kg m-3,
     215                                                               ! units for gases in ppm   
     216!-- initial profiles (also for lateral boundaries if 'set_constant_profiles' is chosen)
     217!-- units for PM: kg m-3, units for gases: ppm
     218    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,
     219    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,
     220    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,
     221    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,
     222
     223    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,
     224    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,
     225    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,
     226    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,
     227
    216228    bc_cs_b                    = 'neumann',
    217229    bc_cs_t                    = 'neumann',
    218230   
    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)
     231    do_emis                    = .T.,
     232    mode_emis                  = 'PARAMETERIZED',  ! emission depending on street type from static driver
     233    surface_csflux_name        = 'NO',   'NO2',  'PM10',   !"cs" stands for chemical species
     234    surface_csflux             = 4745.0, 1326.0, 2700.e-09,! surface emission of PM10, PM25 in kg m-2 d-1,
     235                                                           ! for gases in umol m-2 d-1 for PARAMETERIZED
     236    emiss_factor_main          = 1.667,  1.667,  1.667,    ! emiss_factor_main*surface_csflux for emission on main streets
     237    emiss_factor_side          = 0.334,  0.334,  0.334,    !emiss_factor_side*surface_csflux for emission on side streets
     238    main_street_id = 12,       ! id is for identifying main/side streets (see chem_emissions in chemistry_model_mod.f90)
    223239    side_street_id =  8,
    224     max_street_id = 19,
     240    max_street_id  = 19,
     241
     242! -- KPP solver steering. For details  see comments in chem_gasphase_mod.f90, SUBROUTINE rosenbrock
     243!                         (or other solver subroutine, if another solver is chosen)
     244    icntrl(3)                  = 1,   ! solver ros2
     245    icntrl(4)                  = 500, ! max. number of chem-substeps
     246    rcntrl(3)                  = 0.1, ! Hstart,starting value for the integration step size
    225247
    226248 / ! end of chemistry_par namelist
Note: See TracChangeset for help on using the changeset viewer.