if c == 14:
cbar = ax
else:
cbar = False
wspdMerid=-wspd*cosd(wdir);
ax, hcbar = wp.plot_time_height(ax=ax,
wspd=wspdMerid,
time=time, height=hgt,
spd_range=[0, 30], spd_delta=4,
cmap='nipy_spectral',
cbar=cbar
)
wp.add_windstaff(wspd, wdir, time, hgt, color='k',ax=ax,
vdensity=2, hdensity=2)
''' add arrow annotations '''
scale = 4.1 # use for adjust png output
alpha = 0.6
if None not in times[c]['xpol']:
st = times[c]['xpol'][0]
en = times[c]['xpol'][1]
# h = times[c]['xpol'][2]*scale
h = np.abs(st-en)*(-66.6667*scale)+(15.3333*scale)
connectst = 'bar,armA={},armB={}'.format(h,h),
ax.annotate('',
xy=(st, -0.1 ),
xycoords='axes fraction',
ax=ax,
wspd=wspd_target,
time=time,
height=hgt,
spd_range=[-4, 32],
spd_delta=2,
cmap='nipy_spectral',
cbar=(ax, cbarinvi),
kind='pcolormesh',
)
wp.add_windstaff(wspd,
wdir,
time,
hgt,
color=(0.4, 0.4, 0.4),
ax=ax,
vdensity=2,
hdensity=1,
head_size=0.08,
tail_length=5)
''' determine xticks '''
xtl = ax.get_xticklabels()
xt = ax.get_xticks()
nmod = 0
''' add xtick labels '''
newxtl = []
for i, lb in enumerate(xtl):
if np.mod(i, 6) in [nmod]:
newxtl.append(lb.get_text())
else:
wspd=wspdMerid,
time=time,
height=hgt,
spd_range=[0, 30],
spd_delta=2,
cmap='jet',
cbar=(ax, cbar_inv),
cbar_label='meridional wind $['
'm\ s^{-1}]$')
wp.add_windstaff(
wspd,
wdir,
time,
hgt,
# color=(0.6,0.6,0.6),
color='k',
ax=ax,
vdensity=1,
hdensity=0,
head_size=0.08,
tail_length=5)
''' add arrow annotations '''
vpos1 = -3.11
vpos2 = -2.3
arrstyle = '|-|,widthA = 0.5,widthB = 0.5'
ttacolor = (0, 0, 0)
xplcolor = (0.7, 0.7, 0.7)
for t, p in zip(times, params):
if t['tta'][0] is None:
vpos2 += 2.3
ax, hcbar, img = wp.plot_time_height(ax=ax,
wspd=wspdMerid,
time=time,
height=hgt,
spd_range=[-5, 30],
spd_delta=4,
cmap='nipy_spectral',
cbar=ax,
cbarinvi=cbar_inv
)
wp.add_windstaff(wspd, wdir, time, hgt,
color=(0.0, 0.0, 0.0),
ax=ax,
vdensity=2,
hdensity=2,
head_size=0.07,
tail_length=4.5,
)
if c == 8:
ax.text(0.82, -0.04, '*', color='r',
fontsize=35, transform=ax.transAxes,
zorder=100000)
''' add arrow annotations '''
vpos1 = -3.11
vpos2 = -2.0
arrstyle = '|-|,widthA = 0.5,widthB = 0.5'
ttacolor = (0, 0, 0)
xplcolor = (0.7, 0.7, 0.7)
ax, hcbar = wp.plot_time_height(ax=ax,
wspd=wspdMerid,
time=time,
height=hgt,
spd_range=[-5, 30],
spd_delta=4,
cmap='nipy_spectral',
cbar=(ax, cbar_inv)
)
wp.add_windstaff(wspd, wdir, time, hgt,
# color = (0.5,0.5,0.5),
color=(0.0, 0.0, 0.0),
ax=ax,
vdensity=2,
hdensity=2,
head_size=0.07,
tail_length=4.5,
)
''' add arrow annotations '''
vpos1 = -3.11
vpos2 = -2.0
arrstyle = '|-|,widthA = 0.5,widthB = 0.5'
ttacolor = (0, 0, 0)
xplcolor = (0.7, 0.7, 0.7)
if None not in times[c]['xpol']:
st = times[c]['xpol'][0] - 0.3
en = times[c]['xpol'][1]
ax, hcbar, im = wp.plot_time_height(ax=ax,
wspd=wspdMerid,
time=time,
height=hgt,
spd_range=[0, 30],
spd_delta=2,
cmap='jet',
cbar=(ax, cbar_inv)
)
wp.add_windstaff(wspd, wdir, time, hgt,
# color=(0.6,0.6,0.6),
color='k',
ax=ax,
vdensity=1,
hdensity=0,
head_size=0.08,
tail_length=5
)
" add arrow annotations "
vpos1 = -3.11
vpos2 = -2.3
arrstyle = '|-|,widthA = 0.5,widthB = 0.5'
ttacolor = (0, 0, 0)
xplcolor = (0.7, 0.7, 0.7)
for t,p in zip(times, params):
if t['tta'][0] is None:
vpos2 += 2.3
ax, hcbar = wp.plot_time_height(ax = ax,
wspd = wspd_target,
time = time,
height = hgt,
spd_range = [-4, 32],
spd_delta = 2,
cmap = 'nipy_spectral',
cbar = (ax,cbarinvi),
kind = 'pcolormesh',
)
wp.add_windstaff(wspd, wdir, time, hgt,
color = (0.4,0.4,0.4),
ax = ax,
vdensity = 2,
hdensity = 1,
head_size = 0.08,
tail_length = 5)
''' determine xticks '''
xtl = ax.get_xticklabels()
xt = ax.get_xticks()
nmod = 0
''' add xtick labels '''
newxtl = []
for i, lb in enumerate(xtl):
if np.mod(i, 6) in [nmod]:
wspdMerid = -wspd * cosd(wdir)
ax, hcbar = wp.plot_time_height(ax=ax,
wspd=wspdMerid,
time=time,
height=hgt,
spd_range=[0, 30],
spd_delta=4,
cmap='nipy_spectral',
cbar=cbar)
wp.add_windstaff(wspd,
wdir,
time,
hgt,
color='k',
ax=ax,
vdensity=2,
hdensity=2)
''' add arrow annotations '''
scale = 4.1 # use for adjust png output
alpha = 0.6
if None not in times[c]['xpol']:
st = times[c]['xpol'][0]
en = times[c]['xpol'][1]
# h = times[c]['xpol'][2]*scale
h = np.abs(st - en) * (-66.6667 * scale) + (15.3333 * scale)
connectst = 'bar,armA={},armB={}'.format(h, h),
ax.annotate('',
xy=(st, -0.1),
xycoords='axes fraction',