Skip to content
Snippets Groups Projects
Commit 4d166471 authored by tomrink's avatar tomrink
Browse files

snapshot...

parent e7341110
No related branches found
No related tags found
No related merge requests found
......@@ -1067,7 +1067,7 @@ def write_icing_file_nc4_viirs(clvrx_str_time, output_dir, preds_dct, probs_dct,
dim_0_name = 'x'
dim_1_name = 'y'
time_dim_name = 'time'
geo_coords = 'time y x'
geo_coords = 'longitude latitude'
dim_1_len, dim_0_len = lons.shape
......@@ -1092,7 +1092,7 @@ def write_icing_file_nc4_viirs(clvrx_str_time, output_dir, preds_dct, probs_dct,
icing_pred_ds = rootgrp.createVariable('icing_prediction_level_'+flt_level_ranges_str[flvl], 'i2', var_dim_list)
icing_pred_ds.setncattr('coordinates', geo_coords)
#icing_pred_ds.setncattr('grid_mapping', 'Projection')
icing_pred_ds.setncattr('grid_mapping', 'Projection')
icing_pred_ds.setncattr('missing', -1)
if has_time:
preds = preds.reshape((1, dim_1_len, dim_0_len))
......@@ -1104,7 +1104,7 @@ def write_icing_file_nc4_viirs(clvrx_str_time, output_dir, preds_dct, probs_dct,
icing_prob_ds = rootgrp.createVariable('icing_probability_level_'+flt_level_ranges_str[flvl], 'f4', var_dim_list)
icing_prob_ds.setncattr('coordinates', geo_coords)
#icing_prob_ds.setncattr('grid_mapping', 'Projection')
icing_prob_ds.setncattr('grid_mapping', 'Projection')
if not use_nan:
icing_prob_ds.setncattr('missing', -1.0)
else:
......@@ -1124,7 +1124,7 @@ def write_icing_file_nc4_viirs(clvrx_str_time, output_dir, preds_dct, probs_dct,
icing_prob_ds = rootgrp.createVariable('max_icing_probability_column', 'f4', var_dim_list)
icing_prob_ds.setncattr('coordinates', geo_coords)
#icing_prob_ds.setncattr('grid_mapping', 'Projection')
icing_prob_ds.setncattr('grid_mapping', 'Projection')
if not use_nan:
icing_prob_ds.setncattr('missing', -1.0)
else:
......@@ -1140,14 +1140,14 @@ def write_icing_file_nc4_viirs(clvrx_str_time, output_dir, preds_dct, probs_dct,
icing_pred_ds = rootgrp.createVariable('max_icing_probability_level', 'i2', var_dim_list)
icing_pred_ds.setncattr('coordinates', geo_coords)
#icing_pred_ds.setncattr('grid_mapping', 'Projection')
icing_pred_ds.setncattr('grid_mapping', 'Projection')
icing_pred_ds.setncattr('missing', -1)
icing_pred_ds[:,] = max_lvl
if bt_10_4 is not None:
bt_ds = rootgrp.createVariable('bt_10_4', 'f4', var_dim_list)
bt_ds.setncattr('coordinates', geo_coords)
#bt_ds.setncattr('grid_mapping', 'Projection')
bt_ds.setncattr('grid_mapping', 'Projection')
bt_ds[:,] = bt_10_4
lon_ds = rootgrp.createVariable('longitude', 'f4', [dim_1_name, dim_0_name])
......@@ -1158,15 +1158,7 @@ def write_icing_file_nc4_viirs(clvrx_str_time, output_dir, preds_dct, probs_dct,
lat_ds.units = 'degrees_north'
lat_ds[:,] = lats
# proj_ds = rootgrp.createVariable('Projection', 'b')
# proj_ds.setncattr('long_name', long_name)
# proj_ds.setncattr('grid_mapping_name', 'geostationary')
# proj_ds.setncattr('sweep_angle_axis', cf_nav_dct['sweep_angle_axis'])
# proj_ds.setncattr('semi_major_axis', cf_nav_dct['semi_major_axis'])
# proj_ds.setncattr('semi_minor_axis', cf_nav_dct['semi_minor_axis'])
# proj_ds.setncattr('inverse_flattening', cf_nav_dct['inverse_flattening'])
# proj_ds.setncattr('perspective_point_height', cf_nav_dct['perspective_point_height'])
# proj_ds.setncattr('latitude_of_projection_origin', cf_nav_dct['latitude_of_projection_origin'])
# proj_ds.setncattr('longitude_of_projection_origin', cf_nav_dct['longitude_of_projection_origin'])
proj_ds = rootgrp.createVariable('Projection', 'b')
proj_ds.setncattr('grid_mapping_name', 'latitude_longitude')
rootgrp.close()
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment