Compute vor/div from u/v

# for figures in notebook
%matplotlib inline

# import & initialize epygram
import epygram
epygram.init_env()

# [2025/01/16-17:42:22][epygram.formats][<module>:0072][INFO]: Format: HDF5SAF is deactivated at runtime (Error: No module named 'h5py'). Please deactivate from config.implemented_formats or fix error.

import os
INPUTS_DIR = os.path.join('..', 'inputs')
# open resource in read mode
r = epygram.open(os.path.join(INPUTS_DIR, 'ICMSHAROM+0022'), 'r')

u = r.readfield('S080WIND.U.PHYS')
v = r.readfield('S080WIND.V.PHYS')
# re-compute u/v from psi/khi
uv = epygram.fields.make_vector_field(u, v)
spgeom = u.spectral_geometry
uv.sp2gp()

fig, ax = uv.cartoplot(subsampling=5, vector_plot_method='quiver',
                       components_are_projected_on='grid')

# fields need to be spectral to compute vor/div : spectral derivatives are used
uv.gp2sp(spgeom)
vor, div = uv.compute_vordiv()

# map factor due to stretching:
fig, ax = vor.cartoplot()
fig, ax = div.cartoplot()