Copy dataset units in geocode.py

[EJFielding]

Description of the desired feature

The new addition of UNIT in the velocity.h5 output files is very helpful when we have run timeseries2velocity.py with a step function fit. It would be great if these UNIT were copied to the geocoded output when we run geocode.py velocity.h5.

Is your feature request related to a problem? Please describe

The velocity.h5 file has units.

HDF5 dataset "/step20210908        ": shape (1155, 2053)        , dtype <float32>
  UNIT    m
HDF5 dataset "/step20210908Std     ": shape (1155, 2053)        , dtype <float32>
  UNIT    m
HDF5 dataset "/velocity            ": shape (1155, 2053)        , dtype <float32>
  UNIT    m/year
HDF5 dataset "/velocityStd         ": shape (1155, 2053)        , dtype <float32>
  UNIT    m/year

After running geocode.py on it, the units are lost in the geo_velocity.h5 file:

HDF5 dataset "/step20210908        ": shape (1155, 2053)        , dtype <float32>
HDF5 dataset "/step20210908Std     ": shape (1155, 2053)        , dtype <float32>
HDF5 dataset "/velocity            ": shape (1155, 2053)        , dtype <float32>
HDF5 dataset "/velocityStd         ": shape (1155, 2053)        , dtype <float32>

Describe the solution you'd like

It would be great if the UNIT values were copied over by geocode.py.

Describe alternatives you have considered

Additional context

Are you willing to help implement and maintain this feature?

• [ ] Yes
• [ x] No