added humidity or avg rhmin and rhmax for new onhm

nhm-usgs · Jul 21, 2020 · d4f4891 · d4f4891
1 parent adf8ff7
commit d4f4891
Showing 1 changed file with 9 additions and 2 deletions.
diff --git a/gridmetetl/etl.py b/gridmetetl/etl.py
@@ -521,12 +521,19 @@ def getxy(pt):
 
             elif var == 'srad':
                 srad = ncfile.createVariable('srad', dtype(float32).char, ('time', 'hruid'))
-                srad.long_name = 'surface_downwelling_shortwave_flux_in_air '
+                srad.long_name = 'surface_downwelling_shortwave_flux_in_air'
                 srad.units = 'W m-2'
-                srad.standard_name = 'srad '
+                srad.standard_name = 'srad'
                 srad.fill_value = default_fillvals['f8']
                 srad[:, :] = self.np_srad[:, :]
 
+        hum = ncfile.createVariable('humidity', dtype(float32).char, ('time', 'hruid'))
+        hum.long_name = 'Daily mean relative humidity'
+        hum.units = 'percent'
+        hum.standard_name = 'rhavg'
+        hum.fill_value = default_fillvals['f8']
+        hum[:, :] = (self.np_rhmax[:, :] + self.np_rhmin[:, :])/2.0
+
         ncfile.close()
         print("dataset is closed", flush=True)