search

kmaterna / Strain_2D

Strain rate modeling from GNSS velocity fields
MIT License
43 stars 17 forks source link

compare strain methods writes gridline-node-registered files #17

Closed kmaterna closed 2 years ago

kmaterna commented 2 years ago

I had to manually use grdedit -T to switch them to pixel-node-registered. it happens at compare_strain_grids.py line 50: 

mean_ds.to_netcdf(os.path.join(MyParams.outdir, "means_"+plot_type+".nc"))
stds_ds.to_netcdf(os.path.join(MyParams.outdir, "stds_" + plot_type + ".nc"))
kmaterna commented 2 years ago

Fixed this by saving means and standard deviations into two layers of the same netcdf, so user will have to extract them and impose grid registration separately. https://github.com/kmaterna/Strain_2D/commit/604d2967ca572ff14b23d039743428488b2af57b