ncfile = "swan.four.nc"
lon0 = -5.1
lon1 = -4.9
lat0 = 48.35
lat1 = 48.55
# Imports
from vcmq import cdms2, data_sample, N, transect, curve2, code_file_name, os, \
add_map_lines, P, add_shadow
# Read data
f = cdms2.open(data_sample(ncfile))
hs = f('HS', squeeze=1)
f.close()
# Compute transect
hst, lons, lats = transect(hs, (lon0,lon1), (lat0,lat1), getcoords=True)
# Plot along time
s = curve2(hst, figsize=(8,3), title='Spatial transect on curved grid',
show=False, top=0.85)
# Add a small map to show the transect positions
o = add_map_lines(hs, lons, lats, map_bgcolor='w', map_bbox= [.6, .2, .3, .5], map_anchor='W')
# Save
figfile = code_file_name(ext='png')
if os.path.exists(figfile): os.remove(figfile)
s.savefig(figfile, pdf=True)
from vcmq import cdms2, data_sample, N, transect, curve2, code_file_name, os, \
add_map_lines, P, add_shadow
# Read data
f = cdms2.open(data_sample(ncfile))
hs = f('HS', squeeze=1)
f.close()
# Compute transect
hst, lons, lats = transect(hs, (lon0, lon1), (lat0, lat1), getcoords=True)
# Plot along time
s = curve2(hst,
figsize=(8, 3),
title='Spatial transect on curved grid',
show=False,
top=0.85)
# Add a small map to show the transect positions
o = add_map_lines(hs,
lons,
lats,
map_bgcolor='w',
map_bbox=[.6, .2, .3, .5],
map_anchor='W')
# Save
figfile = code_file_name(ext='png')
if os.path.exists(figfile): os.remove(figfile)
s.savefig(figfile, pdf=True)
# Read data
f = cdms2.open(data_sample(ncfile))
u = f('u')
v = f('v')
f.close()
# Transect specs
tlons, tlats = transect_specs(u.getGrid(), lon0, lat0, lon1, lat1)
ttimes = (time0, time1)
# Compute transect
tu = transect(u, tlons, tlats, times=ttimes)
tv = transect(v, tlons, tlats, times=ttimes)
# Plot along time
s = stick2(tu,
tv,
figsize=(8, 3),
title='Space-time transect of speed',
show=False,
top=0.85,
quiver_width=0.002)
# Add a small map to show the transect positions
add_map_lines(u.getGrid(), tlons, tlats, map_zoom=0.5)
# Save
figfile = code_file_name(ext='png')
if os.path.exists(figfile): os.remove(figfile)
s.savefig(figfile, pdf=True)
# Imports
from vcmq import cdms2, data_sample, N, transect, stick2, code_file_name, os, \
transect_specs, add_map_lines, create_time, IterDates, P
# Read data
f = cdms2.open(data_sample(ncfile))
u = f('u')
v = f('v')
f.close()
# Transect specs
tlons, tlats = transect_specs(u.getGrid(), lon0, lat0, lon1, lat1)
ttimes = create_time(IterDates((time0, time1), len(tlons), closed=True))
# Compute transect
tu = transect(u, tlons, tlats, times=ttimes)
tv = transect(v, tlons, tlats, times=ttimes)
# Plot along time
s = stick2(tu, tv, figsize=(8,3), title='Space-time transect of speed',
show=False, top=0.85, quiver_width=0.002)
# Add a small map to show the transect positions
add_map_lines(u.getGrid(), tlons, tlats, map_zoom=0.5)
# Save
figfile = code_file_name(ext='png')
if os.path.exists(figfile): os.remove(figfile)
s.savefig(figfile, pdf=True)
