fix units for vector hh output

swensosc · Apr 16, 2024 · cad8d42 · cad8d42
1 parent a65df14
commit cad8d42
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/src/main/histFileMod.F90 b/src/main/histFileMod.F90
@@ -3147,16 +3147,16 @@ subroutine htape_timeconst(t, mode)
                   units='m', ncid=nfid(t))             
              call ncd_defvar(varname='hillslope_area', xtype=ncd_double, &
                   dim1name=namec, long_name='hillslope column area', &
-                  units='m', ncid=nfid(t))             
+                  units='m2', ncid=nfid(t))
              call ncd_defvar(varname='hillslope_elev', xtype=ncd_double, &
                   dim1name=namec, long_name='hillslope column elevation', &
                   units='m', ncid=nfid(t))             
              call ncd_defvar(varname='hillslope_slope', xtype=ncd_double, &
                   dim1name=namec, long_name='hillslope column slope', &
-                  units='m', ncid=nfid(t))             
+                  units='m/m', ncid=nfid(t))
              call ncd_defvar(varname='hillslope_aspect', xtype=ncd_double, &
                   dim1name=namec, long_name='hillslope column aspect', &
-                  units='m', ncid=nfid(t))             
+                  units='radians', ncid=nfid(t))
              call ncd_defvar(varname='hillslope_index', xtype=ncd_int, &
                   dim1name=namec, long_name='hillslope index', &
                   ncid=nfid(t))