diff --git a/modules/util/util.py b/modules/util/util.py index 94a2a8dd8a1d2cbdad7870cb38002121237c6165..c4475ff4ed7a39c277757ec716d01d37c08eb72c 100644 --- a/modules/util/util.py +++ b/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'