snapshot...

modules/util/util.py

@@ -11,7 +11,6 @@ import h5py
 import pickle
 from netCDF4 import Dataset
 from util.setup import ancillary_path
-from util.geos_nav import get_cf_nav_parameters
 
 LatLonTuple = namedtuple('LatLonTuple', ['lat', 'lon'])
 
@@ -635,6 +634,37 @@ flt_level_ranges_str[4] = '8000_15000'
 # flt_level_ranges_str = {k: None for k in range(1)}
 # flt_level_ranges_str[0] = 'column'
 
+def get_cf_nav_parameters(satellite='GOES16', domain='FD'):
+    param_dct = None
+
+    if satellite == 'H08':  # We presently only have FD
+        param_dct = {'semi_major_axis': 6378.137,
+                     'semi_minor_axis': 6356.7523,
+                     'perspective_point_height': 35785.863,
+                     'latitude_of_projection_origin': 0.0,
+                     'longitude_of_projection_origin': 140.7,
+                     'inverse_flattening': 298.257,
+                     'sweep_angle_axis': 'y',
+                     'x_scale_factor': 5.58879902955962e-05,
+                     'x_add_offset': -0.153719917308037,
+                     'y_scale_factor': -5.58879902955962e-05,
+                     'y_add_offset': 0.153719917308037}
+    elif satellite == 'GOES16':
+        if domain == 'CONUS':
+            param_dct = {'semi_major_axis': 6378137.0,
+                         'semi_minor_axis': 6356752.31414,
+                         'perspective_point_height': 35786023.0,
+                         'latitude_of_projection_origin': 0.0,
+                         'longitude_of_projection_origin': -75,
+                         'inverse_flattening': 298.257,
+                         'sweep_angle_axis': 'x',
+                         'x_scale_factor': 5.6E-05,
+                         'x_add_offset': -0.101332,
+                         'y_scale_factor': -5.6E-05,
+                         'y_add_offset': 0.128212}
+
+    return param_dct
+
 
 def write_icing_file(clvrx_str_time, output_dir, preds_dct, probs_dct, x, y, lons, lats, elems, lines):
     outfile_name = output_dir + 'icing_prediction_'+clvrx_str_time+'.h5'