snapshot...

modules/util/util.py

@@ -802,17 +802,13 @@ def write_icing_file_nc4(clvrx_str_time, output_dir, preds_dct, probs_dct, x, y,
     icing_pred_ds.setncattr('missing', -1)
     icing_pred_ds[:,] = max_lvl
 
-    # lon_ds = h5f_out.create_dataset('longitude', data=lons, dtype='f4')
-    # lon_ds.attrs.create('units', data='degrees_east')
-    # lon_ds.attrs.create('long_name', data='icing prediction longitude')
-    # lon_ds.dims[0].label = dim_0_name
-    # lon_ds.dims[1].label = dim_1_name
-    #
-    # lat_ds = h5f_out.create_dataset('latitude', data=lats, dtype='f4')
-    # lat_ds.attrs.create('units', data='degrees_north')
-    # lat_ds.attrs.create('long_name', data='icing prediction latitude')
-    # lat_ds.dims[0].label = dim_0_name
-    # lat_ds.dims[1].label = dim_1_name
+    lon_ds = rootgrp.createVariable('longitude', 'f4', [dim_1_name, dim_0_name])
+    lon_ds.units = 'degrees_east'
+    lon_ds[:,] = lons
+
+    lat_ds = rootgrp.createVariable('latitude', 'f4', [dim_1_name, dim_0_name])
+    lat_ds.units = 'degrees_north'
+    lat_ds[:,] = lats
 
     # proj_ds = h5f_out.create_dataset('Projection', data=0, dtype='b')
     # proj_ds.attrs.create('long_name', data='Himawari Imagery Projection')