With the following code Copied from the manual for TP2
https://xhycom.readthedocs.io/en/latest/transects_transports.html
grid_xhy = xhy.open_dataset(gridfile)
bathy_xhy = xhy.open_dataset(bathyfile, grid=gridfile)
fs = xhy.Transect.named("fram_strait")
fs_on_hycom = fs.resolve(grid=grid_xhy)
print(f"T-cells along section : {fs_on_hycom.n_cells}")
print(f"C-grid faces crossed : {fs_on_hycom.n_faces}")
print(f"Section length : {fs_on_hycom.distance_km[-1]:.0f} km")
fs_on_hycom.plot(grid=grid_xhy,bathy=bathy_xhy)
plt.savefig('testsection.png')
Where the gridfile is from TP5 and the depths are from depth_TP5a0.06_05.a ( current NERSC-HYCOM-CICE) repo. These has been checked with md5sum.
I get a map that is somewhat distorted/flipped/transposed. See attachement.
With the following code Copied from the manual for TP2
https://xhycom.readthedocs.io/en/latest/transects_transports.html
grid_xhy = xhy.open_dataset(gridfile)
bathy_xhy = xhy.open_dataset(bathyfile, grid=gridfile)
fs = xhy.Transect.named("fram_strait")
fs_on_hycom = fs.resolve(grid=grid_xhy)
print(f"T-cells along section : {fs_on_hycom.n_cells}")
print(f"C-grid faces crossed : {fs_on_hycom.n_faces}")
print(f"Section length : {fs_on_hycom.distance_km[-1]:.0f} km")
fs_on_hycom.plot(grid=grid_xhy,bathy=bathy_xhy)
plt.savefig('testsection.png')
Where the gridfile is from TP5 and the depths are from depth_TP5a0.06_05.a ( current NERSC-HYCOM-CICE) repo. These has been checked with md5sum.
I get a map that is somewhat distorted/flipped/transposed. See attachement.