Bugfix/hemco grid precision (#267)

* increase precision of hemco grid vars in the hemco standalone simulation to accomodate 12km resolution

* expand base vars

src/components/prior_component/get_hemco_grid_vars.py

@@ -14,9 +14,9 @@ def get_hemco_grid_vars(sv):
     XEDGE = calc_edges(sv.lon)
 
     if option == "YMID":
-        print(" ".join([f"{num:.2f}" for num in YMID]))
+        print(" ".join([f"{num:.5f}" for num in YMID]))
     elif option == "YEDGE":
-        print(" ".join([f"{num:.3f}" for num in YEDGE]))
+        print(" ".join([f"{num:.5f}" for num in YEDGE]))
     elif option == "XMIN":
         XMIN = XEDGE[0]
         print(XMIN)

src/inversion_scripts/operators/TROPOMI_operator.py

@@ -1,3 +1,4 @@
+import os
 import numpy as np
 import xarray as xr
 import pandas as pd
@@ -102,7 +103,7 @@ def apply_average_tropomi_operator(
         jacobian_K = np.zeros([n_gridcells, n_elements], dtype=np.float32)
         jacobian_K.fill(np.nan)
 
-        pertf = (
+        pertf = os.path.expandvars(
             f'{config["OutputPath"]}/{config["RunName"]}/'
             f'archive_perturbation_sfs/pert_sf_{period_i}.npz'
         )

src/inversion_scripts/operators/operator_utilities.py

@@ -1,7 +1,8 @@
+import os
+import yaml
 import numpy as np
 import xarray as xr
 import pandas as pd
-import yaml
 from src.inversion_scripts.utils import (
     check_is_OH_element,
     check_is_BC_element
@@ -171,7 +172,7 @@ def concat_tracers(run_id, gc_date, config, sv_elems, n_elements, baserun=False)
 
     """
     prefix = config['OutputPath'] + '/' + config['RunName'] + '/jacobian_runs'
-    j_dir = f"{prefix}/{config['RunName']}_{run_id}/OutputDir"
+    j_dir = os.path.expandvars(f"{prefix}/{config['RunName']}_{run_id}/OutputDir")
     file_stub = gc_date.strftime('GEOSChem.SpeciesConc.%Y%m%d_0000z.nc4')
     dsmf = xr.open_dataset(
         '/'.join([j_dir,file_stub]),