diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index 60adf4ff1d3e2c5cc04e254f23b03f4652113910..249d1d6a341abb30264ac2b388a5128bb363926a 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -99,7 +99,7 @@ def get_grid_values(h5f, grid_name, j_c, i_c, half_width, scale_factor_name='sca
     return grd_vals
 
 
-def create_file(filename, data_dct, ds_list, lon_c, lat_c, time_s, fl_alt_s):
+def create_file(filename, data_dct, ds_list, lon_c, lat_c, time_s, fl_alt_s, icing_intensity=None):
     h5f_expl = h5py.File(a_clvr_file, 'r')
     h5f = h5py.File(filename, 'w')
 
@@ -123,6 +123,10 @@ def create_file(filename, data_dct, ds_list, lon_c, lat_c, time_s, fl_alt_s):
     ice_alt_ds.dims[0].label = 'time'
     ice_alt_ds.attrs.create('units', data='m')
 
+    if icing_intensity is not None:
+        icing_int_ds = h5f.create_dataset('icing_intensity', data=icing_intensity)
+        icing_int_ds.attrs.create('long_name', data='0:No intensity report, 1:Trace, 2:Light, 3:Light Moderate, 4:Moderate, 5:Moderate Severe, 6:Severe')
+
     # copy relevant attributes
     for ds_name in ds_list:
         h5f_ds = h5f[ds_name]