Ignore:
Timestamp:
Nov 2, 2020 10:31:45 AM (2 years ago)
Author:
Giersch
Message:

Tutorials updated for r4761

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/TUTORIALS/cases/chem_phstatp/INPUT/chem_phstatp_p3d

    r4409 r4765  
    2020!-- pressure solver
    2121!--------------------------------------
    22     fft_method =  'fftw',
    23 
    24 !     psolver = 'multigrid',
    25 !     cycle_mg  = 'w',
    26 !     mg_cycles = 2,
    27 !     mg_switch_to_pe0_level = -1,
     22!   fft_method =  'fftw',
     23
     24      psolver = 'multigrid',
     25      cycle_mg  = 'w',
     26      mg_cycles = 2,
     27      mg_switch_to_pe0_level = -1,
    2828
    2929!
     
    3232    initializing_actions = 'set_constant_profiles',  ! options: 'set_constant_profiles', or 'inifor',
    3333                                                     ! or both, e.g. 'inifor set_constant_profiles'
    34     constant_flux_layer = .TRUE.,
     34    constant_flux_layer = .T.,
    3535!
    3636!-- wind speed
     
    4848!-- humidity
    4949!--------------------------------------
    50     humidity   = .TRUE.,
     50    humidity   = .T.,
    5151    q_surface  = 0.001,
    5252    bc_q_b     = 'dirichlet',
     
    5454!-- time
    5555!--------------------------------------
    56     origin_date_time     = '2017-07-31 03:00:00 +00'
     56    origin_date_time     = '2017-07-31 06:00:00 +00'
    5757!
    5858!-- wall/soil spinup
     
    6262     spinup_pt_amplitude = 20.0,
    6363     dt_spinup           = 120.0,           ! time step during spinup
    64      data_output_during_spinup = .FALSE.,       ! set .TRUE. if you want output
     64     data_output_during_spinup = .F.,       ! set .T. if you want output
    6565!
    6666!-- building topology
     
    8787!-- 2d/3d output
    8888!--------------------------------------
    89     do3d_at_begin         = .TRUE.,
    90     do2d_at_begin         = .TRUE.,
     89    do3d_at_begin         = .T.,
     90    do2d_at_begin         = .T.,
    9191    dt_data_output        = 300.0,
    9292    dt_data_output_av     = 600.0,
     
    9696
    9797     
    98     data_output = 'theta', 'q', 'u', 'v', 'w','e', 't_soil', 'm_soil', 'tsurf*_xy', 'shf*_xy', 'qsws*_xy', 'ghf*_xy', 'rad_net*_xy', 'us*_xy', 't*_xy', 'r_a*_xy', 'm_liq*_xy', 'kc_NO2','kc_O3','kc_NO', 'kc_PM10', 'theta_av', 'q_av', 'u_av', 'v_av', 'w_av','e_av', 'tsurf*_xy_av', 'shf*_xy_av', 'qsws*_xy_av', 'ghf*_xy_av', 'kc_NO2_av','kc_O3_av','kc_NO_av', 'kc_PM10_av',
     98    data_output = 'theta', 'q', 'u', 'v', 'w','e',
     99                   't_soil', 'm_soil',
     100                   'tsurf*_xy', 'shf*_xy', 'qsws*_xy', 'ghf*_xy',
     101                   'rad_net*_xy', 'us*_xy', 't*_xy', 'r_a*_xy', 'm_liq*_xy',
     102                   'kc_NO2','kc_O3','kc_NO', 'kc_PM10',
     103                   'theta_av', 'q_av', 'u_av', 'v_av', 'w_av','e_av',
     104                   'tsurf*_xy_av', 'shf*_xy_av', 'qsws*_xy_av', 'ghf*_xy_av',
     105                   'kc_NO2_av','kc_O3_av','kc_NO_av', 'kc_PM10_av',
    99106!
    100107! Masked output - here for first layer above ground
     
    115122    dt_averaging_input_pr =  600.0,
    116123
    117     data_output_pr = 'e', 'e*', '#u', 'w"u"', 'w*u*', 'wu', 'u*2', '#v', 'w"v"', 'w*v*', 'wv', 'v*2', 'w',  'w*2', '#theta', 'w"theta"', 'w*theta*','wtheta', 'theta*2', 'kc_NO2','kc_O3','kc_NO',
     124    data_output_pr = 'e', 'e*',
     125                     '#u', 'w"u"', 'w*u*', 'wu', 'u*2',
     126                     '#v', 'w"v"', 'w*v*', 'wv', 'v*2',
     127                     'w',  'w*2',
     128                     '#theta', 'w"theta"', 'w*theta*','wtheta', 'theta*2',
     129                     'kc_NO2','kc_O3','kc_NO',
    118130 /
    119131
     
    127139
    128140    albedo_type      = 5,
    129     constant_albedo  = .FALSE.,
    130 
    131     surface_reflections = .TRUE.,
     141    constant_albedo  = .F.,
     142
     143    surface_reflections = .T.,
    132144    nrefsteps = 5,
    133 
    134     rad_angular_discretization = .FALSE.,
    135     raytrace_mpi_rma = .FALSE.,
    136145
    137146 / ! end of radiation_parameters namelist
     
    143152 &urban_surface_parameters
    144153
    145     usm_material_model = .TRUE.,
    146154 / ! end of urban_surface_parameters namelist
    147155
     
    160168    deep_soil_temperature = 285.,
    161169
    162     aero_resist_kray       = .TRUE.,
    163     constant_roughness     = .TRUE.,
    164     conserve_water_content = .TRUE.,
     170    aero_resist_kray       = .T.,
     171    constant_roughness     = .T.,
     172    conserve_water_content = .T.,
    165173 / ! end of land_surface_parameters namelist
    166174
     
    182190
    183191    chem_gasphase_on           = .TRUE.,
    184     emissions_anthropogenic    = .TRUE.,
    185     emiss_lod                  = 0,          ! 0 "PARAMETERIZED", 1 "DEFAULT", 2 "PRE-PROCESSED",
    186192    chem_mechanism             = "phstatp",
    187193    deposition_dry             = .TRUE.,
    188194    photolysis_scheme          = 'simple',
    189195
    190     surface_csflux_name        = 'NO',   'NO2', 'PM10',     !"cs" stands for chemical species, NOT to be used when running
    191                                                                            ! with mode_emis = "PRE-PROCESSED" or mode_emis = "DEFAULT"
    192     surface_csflux             = 4745.0, 1326.0, 2700e-09,  ! surface emission of NO, NO2 in umol m-2 d-1 (PMs in kg m-2 d-1),
    193                                                                            ! NOT to be used when running with mode_emis = "PRE-PROCESSED" or mode_emis = "DEFAULT"
    194     emiss_factor_main          = 1.667, 1.667, 1.667,      ! emiss_factor_main*surface_csflux for emission on main streets, one factor per surface_csflux component needed!
    195     emiss_factor_side          = 0.334, 0.334, 0.334,      ! emiss_factor_side*surface_csflux for emission on side streets, one factor per surface_csflux component needed!
    196     main_street_id             = 11,                                       ! id is for identifying main/side streets (see chem_emissions in chemistry_model_mod.f90)
     196! initial profiles
     197    cs_name                    = 'O3',  'NO',  'NO2',   'PM10',
     198    cs_surface                 = 0.025, 0.000, 0.010,    2.0e-09,       !units- chem spcs in ppm, and PM10 in ug/m**2
     199
     200    cs_profile(1,:)            = 0.025, 0.040, 0.050,
     201    cs_profile(2,:)            = 0.000, 0.000,
     202    cs_profile(3,:)            = 0.010, 0.002,
     203    cs_profile(4,:)            = 2.0e-09, 2.0e-09, 2.0e-09, 2.0e-09, 2.0e-09, 2.0e-09, 2.0e-09, 2.0e-09, 2.0e-09, 2.0e-09, 1.0e-09,
     204
     205    cs_heights(1,:)            = 0.0,  305.0,  495.0,
     206    cs_heights(2,:)            = 0.0,   295.0,
     207    cs_heights(3,:)            = 0.0,   295.0,
     208    cs_heights(4,:)            = 0.0,   15.0,   15.0,   45.0,  65.0,  85.0, 105.0, 125.0, 155.0,  185.0, 295.0,
     209
     210! Emissions
     211    emissions_anthropogenic    = .TRUE.,
     212    emiss_lod                  = 0,                 ! 0 "PARAMETERIZED", 1 "DEFAULT", 2 "PRE-PROCESSED",
     213! Parameters for parameterized traffic emissions
     214    surface_csflux_name        = 'NO',   'NO2', 'PM10',     !"cs" stands for chemical species,
     215                                                            ! only for parameterized traffic emissions
     216                                                            ! NOT to be used when running with emiss_lod>0
     217    surface_csflux             = 4745.0, 1326.0, 2700e-09,  ! surface emission of NO, NO2 in umol m-2 d-1 (PMs in kg m-2 d-1),only for parameterized traffic emissions
     218    emiss_factor_main          = 1.667, 1.667, 1.667,       ! emiss_factor_main*surface_csflux for emission on main streets, one factor per surface_csflux component needed!
     219    emiss_factor_side          = 0.334, 0.334, 0.334,       ! emiss_factor_side*surface_csflux for emission on side streets, one factor per surface_csflux component needed!
     220    main_street_id             = 11,                        ! id is for identifying main/side streets (see chem_emissions in chemistry_model_mod.f90)
    197221    side_street_id             =  8,
    198222    max_street_id              = 19,
    199223
    200     cs_name                    = 'O3',  'NO',  'NO2',   'PM10',
    201     cs_surface                 = 0.020, 0.000, 0.010,    2.0e-09,       !units- chem spcs in ppm, and PM10 in ug/m**2
    202 
    203     cs_profile(1,:)            =  0.015, 0.018, 0.020, 0.021, 0.022, 0.023, 0.024, 0.025, 0.027, 0.030, 0.050,    !units for gasphas ppm
    204     cs_profile(2,:)            =  0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,    !units for gasphas ppm
    205     cs_profile(3,:)            =  0.010, 0.008, 0.006, 0.005, 0.004, 0.003, 0.002, 0.002, 0.002, 0.002, 0.002,
    206     cs_profile(4,:)            =  2.0e-09, 2.0e-09, 2.0e-09, 2.0e-09, 2.0e-09, 2.0e-09, 2.0e-09, 2.0e-09, 2.0e-09, 2.0e-09, 1.0e-09,
    207 
    208 
    209     cs_heights(1,:)            =  0.0,   15.0,   25.0,   45.0,  65.0,  85.0, 105.0, 125.0, 155.0,  185.0, 295.0,
    210     cs_heights(2,:)            =  0.0,   15.0,   25.0,   45.0,  65.0,  85.0, 105.0, 125.0, 155.0,  185.0, 295.0,
    211     cs_heights(3,:)            =  0.0,   15.0,   25.0,   45.0,  65.0,  85.0, 105.0, 125.0, 155.0,  185.0, 295.0,
    212     cs_heights(4,:)            =  0.0,   15.0,   15.0,   45.0,  65.0,  85.0, 105.0, 125.0, 155.0,  185.0, 295.0,
    213 
    214     bc_cs_b                    = 'neumann',                      !'dirichlet'
     224! Bottom and top boundary conditions (lateral BCs are cyclic, which is the default)
     225    bc_cs_b                    = 'neumann',
    215226    bc_cs_t                    = 'neumann',
    216227
    217 
    218     icntrl(3)                  = 1,    !solver ros2
    219     icntrl(4)                  = 500,  !max number of chem-substeps
     228! Parameters chemistry solver (see KPP description)
     229    icntrl(3)                  = 1,    ! solver ros2
     230    icntrl(4)                  = 500,  ! max number of chem-substeps
    220231    rcntrl(3)                  = 0.1,  ! Hstart, starting value for the integration step size
    221232
Note: See TracChangeset for help on using the changeset viewer.