Changeset 3632 for palm/trunk/TESTS

Timestamp:

Dec 17, 2018 2:52:50 PM (6 years ago)

Author:

kanani

Message:

Fix for masked output, change default chemistry mechanism

Location:

palm/trunk/TESTS/cases

Files:

• urban_environment/INPUT/urban_environment_p3d (modified) (4 diffs)
• urban_environment/OUTPUT (added)
• urban_environment/OUTPUT/urban_environment_pr.nc (added)
• urban_environment/OUTPUT/urban_environment_ts.nc (added)
• urban_environment/case_config.yml (modified) (1 diff)
• urban_environment_restart/INPUT/urban_environment_restart_p3dr (modified) (3 diffs)
• urban_environment_restart/OUTPUT (added)
• urban_environment_restart/OUTPUT/urban_environment_restart_pr.001.nc (added)
• urban_environment_restart/OUTPUT/urban_environment_restart_ts.001.nc (added)
• urban_environment_restart/case_config.yml (modified) (1 diff)

palm/trunk/TESTS/cases/urban_environment/INPUT/urban_environment_p3d

@@ -30 +30 @@
 !-- Initialization
 !------------------------------------------------------------------------------
-    initializing_actions = 'inifor',   !'set_constant_profiles', or 'inifor',
+    initializing_actions = 'inifor set_constant_profiles',   !'set_constant_profiles', or 'inifor', 
+                                                             ! or both (met. from inifor, chem constant profiles)
 
     pt_surface           = 277.15,
@@ -90 +91 @@
     data_output_pr        = '#u',  'u*2',  'wu',  'w*u*',  'w"u"',
                             '#v',  'v*2',  'wv',  'w*v*',  'w"v"',
-                            'w',   'w*2',
+!                            'w',   'w*2',
                             '#theta', 'theta*2', 'wtheta', 'w*theta*', 'w"theta"',
                             '#q',  'q*2',
                             'e', 'e*', '#km', '#l',
                             '#t_soil', '#m_soil', 'rad_lw_in', 'rad_lw_out',
+                            'kc_NO', 'kc_NO2', 'kc_O3', 'kc_PM10', 
 !
 !-- 2D/3D output settings
@@ -127 +129 @@
                   'tsurf*_xy', 'tsurf*_xy_av',
                   'pcm_lad', 'pcm_heatrate',
-                  'kc_PM10', 'kc_PM25', 'kc_PM10_av', 'kc_PM25_av',
-
-
+                  'kc_NO', 'kc_NO2', 'kc_O3', 'kc_PM10', 
+                  'kc_NO_av', 'kc_NO2_av', 'kc_O3_av', 'kc_PM10_av', 
 
  / ! end of d3par namelist
@@ -206 +207 @@
  &chemistry_parameters
 
-    chem_gasphase_on           = .F.,
+    chem_gasphase_on           = .T.,
+    call_chem_at_all_substeps  = .FALSE.,
     
-    do_emis                    = .T.,
-    mode_emis                  = 'PARAMETERIZED',
-    surface_csflux_name        = 'PM10','PM25',  !"cs" stands for chemical species
-    surface_csflux             = 0.75, 0.375,    !surface emission of PM10, PM25 (arbitrary units)
-    cs_name                    = 'PM10','PM25',
-    cs_surface                 = 2.0, 1.0,       !surface concentration of PM10, PM25 (arbitrary units)
+    photolysis_scheme          = 'simple',
+
+    cs_name                    = 'NO', 'NO2', 'O3', 'PM10',
+    cs_surface                 = 0.02,  0.03, 0.02, 25.0e-9,   ! surface concentration of PM10, PM25 in kg m-3, 
+                                                               ! units for gases in ppm   
+!-- initial profiles (also for lateral boundaries if 'set_constant_profiles' is chosen)
+!-- units for PM: kg m-3, units for gases: ppm
+    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,
+    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,
+    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,
+    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, 
+
+    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,
+    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,
+    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,
+    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,
+
     bc_cs_b                    = 'neumann',
     bc_cs_t                    = 'neumann',
     
-!-- for emission depending on street type in static driver
-    emiss_factor_main = 1.667, 1.667,   !emiss_factor_main*surface_csflux for emission on main streets
-    emiss_factor_side = 0.334, 0.334,   !emiss_factor_side*surface_csflux for emission on side streets
-    main_street_id = 12,          !id is for identifying main/side streets (see chem_emissions in chemistry_model_mod.f90)
+    do_emis                    = .T.,
+    mode_emis                  = 'PARAMETERIZED',  ! emission depending on street type from static driver
+    surface_csflux_name        = 'NO',   'NO2',  'PM10',   !"cs" stands for chemical species
+    surface_csflux             = 4745.0, 1326.0, 2700.e-09,! surface emission of PM10, PM25 in kg m-2 d-1, 
+                                                           ! for gases in umol m-2 d-1 for PARAMETERIZED 
+    emiss_factor_main          = 1.667,  1.667,  1.667,    ! emiss_factor_main*surface_csflux for emission on main streets
+    emiss_factor_side          = 0.334,  0.334,  0.334,    !emiss_factor_side*surface_csflux for emission on side streets
+    main_street_id = 12,       ! id is for identifying main/side streets (see chem_emissions in chemistry_model_mod.f90)
     side_street_id =  8,
-    max_street_id = 19,
+    max_street_id  = 19,
+
+! -- KPP solver steering. For details  see comments in chem_gasphase_mod.f90, SUBROUTINE rosenbrock 
+!                         (or other solver subroutine, if another solver is chosen)
+    icntrl(3)                  = 1,   ! solver ros2
+    icntrl(4)                  = 500, ! max. number of chem-substeps
+    rcntrl(3)                  = 0.1, ! Hstart,starting value for the integration step size
 
  / ! end of chemistry_par namelist

palm/trunk/TESTS/cases/urban_environment/case_config.yml

@@ -16 +16 @@
 
 significant_digits_for_netcdf_checks:
-  timeseries: 6
-  profiles: 6
-  other: 8
+  timeseries: 5
+  profiles: 5
+  other: 2

palm/trunk/TESTS/cases/urban_environment_restart/INPUT/urban_environment_restart_p3dr

@@ -90 +90 @@
     data_output_pr        = '#u',  'u*2',  'wu',  'w*u*',  'w"u"',
                             '#v',  'v*2',  'wv',  'w*v*',  'w"v"',
-                            'w',   'w*2',
+!                            'w',   'w*2',
                             '#theta', 'theta*2', 'wtheta', 'w*theta*', 'w"theta"',
                             '#q',  'q*2',
                             'e', 'e*', '#km', '#l',
                             '#t_soil', '#m_soil', 'rad_lw_in', 'rad_lw_out',
+                            'kc_NO', 'kc_NO2', 'kc_O3', 'kc_PM10', 
 !
 !-- 2D/3D output settings
@@ -127 +128 @@
                   'tsurf*_xy', 'tsurf*_xy_av',
                   'pcm_lad', 'pcm_heatrate',
-                  'kc_PM10', 'kc_PM25', 'kc_PM10_av', 'kc_PM25_av',
-
-
+                  'kc_NO', 'kc_NO2', 'kc_O3', 'kc_PM10', 
+                  'kc_NO_av', 'kc_NO2_av', 'kc_O3_av', 'kc_PM10_av', 
 
  / ! end of d3par namelist
@@ -206 +206 @@
  &chemistry_parameters
 
-    chem_gasphase_on           = .F.,
+    chem_gasphase_on           = .T.,
+    call_chem_at_all_substeps  = .FALSE.,
     
-    do_emis                    = .T.,
-    mode_emis                  = 'PARAMETERIZED',
-    surface_csflux_name        = 'PM10','PM25',  !"cs" stands for chemical species
-    surface_csflux             = 0.75, 0.375,    !surface emission of PM10, PM25 (arbitrary units)
-    cs_name                    = 'PM10','PM25',
-    cs_surface                 = 2.0, 1.0,       !surface concentration of PM10, PM25 (arbitrary units)
+    photolysis_scheme          = 'simple',
+
+    cs_name                    = 'NO', 'NO2', 'O3', 'PM10',
+    cs_surface                 = 0.02,  0.03, 0.02, 25.0e-9,   ! surface concentration of PM10, PM25 in kg m-3, 
+                                                               ! units for gases in ppm   
+!-- initial profiles (also for lateral boundaries if 'set_constant_profiles' is chosen)
+!-- units for PM: kg m-3, units for gases: ppm
+    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,
+    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,
+    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,
+    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, 
+
+    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,
+    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,
+    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,
+    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,
+
     bc_cs_b                    = 'neumann',
     bc_cs_t                    = 'neumann',
     
-!-- for emission depending on street type in static driver
-    emiss_factor_main = 1.667, 1.667,   !emiss_factor_main*surface_csflux for emission on main streets
-    emiss_factor_side = 0.334, 0.334,   !emiss_factor_side*surface_csflux for emission on side streets
-    main_street_id = 12,          !id is for identifying main/side streets (see chem_emissions in chemistry_model_mod.f90)
+    do_emis                    = .T.,
+    mode_emis                  = 'PARAMETERIZED',  ! emission depending on street type from static driver
+    surface_csflux_name        = 'NO',   'NO2',  'PM10',   !"cs" stands for chemical species
+    surface_csflux             = 4745.0, 1326.0, 2700.e-09,! surface emission of PM10, PM25 in kg m-2 d-1, 
+                                                           ! for gases in umol m-2 d-1 for PARAMETERIZED 
+    emiss_factor_main          = 1.667,  1.667,  1.667,    ! emiss_factor_main*surface_csflux for emission on main streets
+    emiss_factor_side          = 0.334,  0.334,  0.334,    !emiss_factor_side*surface_csflux for emission on side streets
+    main_street_id = 12,       ! id is for identifying main/side streets (see chem_emissions in chemistry_model_mod.f90)
     side_street_id =  8,
-    max_street_id = 19,
+    max_street_id  = 19,
+
+! -- KPP solver steering. For details  see comments in chem_gasphase_mod.f90, SUBROUTINE rosenbrock 
+!                         (or other solver subroutine, if another solver is chosen)
+    icntrl(3)                  = 1,   ! solver ros2
+    icntrl(4)                  = 500, ! max. number of chem-substeps
+    rcntrl(3)                  = 0.1, ! Hstart,starting value for the integration step size
 
  / ! end of chemistry_par namelist

palm/trunk/TESTS/cases/urban_environment_restart/case_config.yml

@@ -18 +18 @@
 
 significant_digits_for_netcdf_checks:
-  timeseries: 6
-  profiles: 6
+  timeseries: 4
+  profiles: 4
   other: 8