From 66f1689b69d1d76064951726455c7328fd222084 Mon Sep 17 00:00:00 2001 From: David Jean du Preez Date: Fri, 13 May 2022 16:01:28 +0200 Subject: [PATCH] organics --- configurations/mbees_les_test_10.conf | 4 +- configurations/mbees_rans_test_100.conf | 61 ++++++++++++++++++ .../palm_dynamic_aerosol.cpython-38.pyc | Bin 6842 -> 6842 bytes .../__pycache__/palm_wrf_utils.cpython-38.pyc | Bin 17439 -> 17439 bytes dynamic/palm_dynamic.py | 2 + 5 files changed, 65 insertions(+), 2 deletions(-) create mode 100644 configurations/mbees_rans_test_100.conf diff --git a/configurations/mbees_les_test_10.conf b/configurations/mbees_les_test_10.conf index e1387f4..d23566e 100644 --- a/configurations/mbees_les_test_10.conf +++ b/configurations/mbees_les_test_10.conf @@ -8,7 +8,7 @@ scenario = 'les_test' nested_domain = False # dynamic driver output -dynamic_driver_file = "/cfs/home/d/u/dupreeda/MBEES/PALM/palm_model_system-v21.10/JOBS/augs3_1_3/INPUT/augs3_1_3_dynamic" +dynamic_driver_file = "/cfs/home/d/u/dupreeda/MBEES/PALM/palm_model_system-v21.10/JOBS/augs3_1_3/INPUT/augs3_1_dynamic_12hr" # import grid parameters for dynamic driver from static driver grid_from_static = True @@ -29,7 +29,7 @@ dz_max = 10.0 # time origin and extent of the simulation (format YYYY-MM-DD hh:mm:ss) origin_time = "2022-03-16 06:00:00" -simulation_hours = 3 +simulation_hours = 12 # wrf_chem input dir and mask wrf_dir_name = '/cfs/home/d/u/dupreeda/MBEES/PALM/wrf_chem_data' diff --git a/configurations/mbees_rans_test_100.conf b/configurations/mbees_rans_test_100.conf new file mode 100644 index 0000000..5f77829 --- /dev/null +++ b/configurations/mbees_rans_test_100.conf @@ -0,0 +1,61 @@ +# Test case for developing WRF-CHEM for PALM +# see GMD: https://doi.org/10.5194/gmd-2020-175 + +# PALM case, domain, and configuration parameters +domain = 'mbees' +resolution = '100' +scenario = 'rans_test' +nested_domain = False + +# dynamic driver output +dynamic_driver_file = "/cfs/home/d/u/dupreeda/MBEES/PALM/palm_model_system-v21.10/JOBS/mbees/mbees_dynamic12hr" + +# import grid parameters for dynamic driver from static driver +grid_from_static = True +# static driver input +static_driver_file = "/cfs/home/d/u/dupreeda/MBEES/PALM/palm_model_system-v21.10/JOBS/mbees/mbees_static" + +# reference coordinate system of PALM simulation +proj_palm = "EPSG:32633" +# projection lon-lat +proj_wgs84 = 'EPSG:4326' + +# vertical grid +dz = 10.0 +nz = 120 # z in grids +dz_stretch_level = 3000.0 # in meters +dz_stretch_factor = 1.00 +dz_max = 10.0 + +# time origin and extent of the simulation (format YYYY-MM-DD hh:mm:ss) +origin_time = "2022-03-16 06:00:00" +simulation_hours = 12 + +# wrf_chem input dir and mask +wrf_dir_name = '/cfs/home/d/u/dupreeda/MBEES/PALM/wrf_chem_data' + +wrf_file_mask = "wrfout_d01_2022*" +wrf_hybrid_levs = False +vinterp_terrain_smoothing = None +# interpolated files +interp_dir_name = '/cfs/home/d/u/dupreeda/MBEES/PALM/wrf_chem_data/interp' + +# WRF-chem gaseous species +# If salsa is used, h2so4, hno3, nh3, ocnv, ocsv must be included +#wrfchem_spec = ['no', 'no2'] +wrfchem_spec = ['h2so4','hno3','nh3','ocnv','ocsv'] + +# Salsa parameters (only aerosol species) +aerosol_wrfchem = True +wrfchem_bin_limits = [3.9e-8, 1.56e-7, 6.25e-7, 2.5e-6, 1.0e-5] +nbin = [1,7] +reglim = [3.9e-8, 5.0e-8, 2.5e-6] +nf2a = 1.0 +# Either SO4 or OC must be included in listspec +listspec = ['DU', 'SO4'] +#listspec = ['DU','SO4', 'OC','BC','SS','NH','NO'] + +# radiation +radiation_from_wrf = False +wrf_rad_file_mask = "auxhist6_*" +radiation_smoothing_distance = 10000.0 diff --git a/dynamic/__pycache__/palm_dynamic_aerosol.cpython-38.pyc b/dynamic/__pycache__/palm_dynamic_aerosol.cpython-38.pyc index af1c1b1dbf1f47229aa634c3b081e036c7ff24cd..5ba851d32a195cd2f8b301e60b414610d857795a 100644 GIT binary patch delta 19 ZcmdmGy33R+l$V!_0SLY>+{m>-3IH{81$_Vj delta 19 ZcmdmGy33R+l$V!_0SN8~Z{*q_1pqT?1sMPU diff --git a/dynamic/__pycache__/palm_wrf_utils.cpython-38.pyc b/dynamic/__pycache__/palm_wrf_utils.cpython-38.pyc index a9540272c3abb20c2cb20b3a6722bb675094b682..29708878012c30a5e0f8765bbe65275d635f4fa8 100644 GIT binary patch delta 21 bcmbQ=!8pHzkt>v!mx}=izAoIzCF%kIKRX3} delta 21 bcmbQ=!8pHzkt>v!mx}=i>=HL}iMjv)Ik5zX diff --git a/dynamic/palm_dynamic.py b/dynamic/palm_dynamic.py index c3707a4..d8f4dca 100644 --- a/dynamic/palm_dynamic.py +++ b/dynamic/palm_dynamic.py @@ -415,6 +415,7 @@ for wrf_file in wrf_files_proc: v_wrf_valout = N_avr * (1/inv_den) * v_wrf_val/mol_w v_out = f_out.createVariable(varname, 'f4', v_wrf.dimensions) v_out[:] = regridder.regrid(v_wrf_valout[...,regridder.ys,regridder.xs]) + # smpa, smpbb, biog1_c, biog1_o # NVOCs elif varname == 'ocnv': nvoc_list = ['gly', 'n2o5'] @@ -431,6 +432,7 @@ for wrf_file in wrf_files_proc: # SVOCs elif varname == 'ocsv': svoc_list = ['hcho', 'ch3oh'] + mol_w = [30.03, 32.04] vsize = f_wrf.variables[svoc_list[0]].shape v_wrf = np.zeros(vsize) cnt = 0 -- GitLab