Gauss grid resolutionΒΆ
[1]:
%matplotlib inline
# for figures in notebook
# import & initialize epygram
import epygram
epygram.init_env()
import os
INPUTS_DIR = os.path.join('..', 'inputs')
[2]:
# open resource in read mode
g = epygram.formats.resource(os.path.join(INPUTS_DIR, 'analysis.full-arpege.tl149-c24.fa'), 'r')
f = g.readfield('SURFTEMPERATURE')
[3]:
print(f.geometry.resolution_ll(2,45))
55350.67187324146
[4]:
r = f.geometry.resolution_field_from_stretching()
fig, ax = r.cartoplot(colorsnumber=10)
[ ]: