rain_czd=0.25,
nhours=2
)
try:
tta
except NameError:
tta = gapflow.get_tta_dates([2003, 2004], params)
cases = [8, 9, 11, 12, 13]
wprof_hgt = 0.5 # [km] first gate to top layer
for n, c in enumerate(cases):
gf = gapflow.run(c,
plot_theory=True,
grid=False,
wprof_hgt=wprof_hgt,
ax=axes[n],
homedir=homed,
color_surf=(0, 0, 0, 0),
color_wp=(0, 0, 0, 0)
)
subdates = get_subdates(gf, tta.index)
sub = gf.loc[subdates]
x1, y1 = sub.pdiff.values, sub.ucomp.values
x2, y2 = sub.pdiff.values, sub.wp_ucomp.values
cl_surf = (0, 0, 0.5)
cl_wprf = (0, 0.5, 0)
h1 = axes[n].scatter(x1, y1, s=90, marker='x',
color=cl_surf, lw=2)
h2 = axes[n].scatter(x2, y2, s=100, marker='+',
color=cl_wprf, lw=2)
# params = dict(wdir_surf=130, wdir_wprof=170,
# rain_czd=0.25, nhours=2)
params = dict(wdir_thres=150, wdir_layer=[0, 500], rain_czd=0.25, nhours=2)
try:
tta
except NameError:
tta = gapflow.get_tta_dates([2003, 2004], params)
cases = [8, 9, 11, 12, 13]
wprof_hgt = 0.5 # [km] first gate to top layer
for n, c in enumerate(cases):
gf = gapflow.run(c,
plot_theory=True,
grid=False,
wprof_hgt=wprof_hgt,
ax=axes[n],
homedir=homed,
color_surf=(0, 0, 0, 0),
color_wp=(0, 0, 0, 0))
subdates = get_subdates(gf, tta.index)
sub = gf.loc[subdates]
x1, y1 = sub.pdiff.values, sub.ucomp.values
x2, y2 = sub.pdiff.values, sub.wp_ucomp.values
cl_surf = (0, 0, 0.5)
cl_wprf = (0, 0.5, 0)
h1 = axes[n].scatter(x1, y1, s=90, marker='x', color=cl_surf, lw=2)
h2 = axes[n].scatter(x2, y2, s=100, marker='+', color=cl_wprf, lw=2)
axes[n].text(0.05,
0.95,
axes[n].get_gid(),
axes[3] = plt.subplot(gs0[3], gid='(d) 16-18Feb04')
params = dict(wdir_surf=130, wdir_wprof=170, rain_czd=0.25, nhours=2)
try:
tta
except NameError:
tta = np.array(xta.get_tta_dates([2003, 2004], params))
tta = tta + timedelta(minutes=55)
cases = [8, 9, 12, 11]
for n, c in enumerate(cases):
gf = gapflow.run(c,
plot=True,
grid=False,
ax=axes[n],
homedir=homed,
color_surf=(0, 0, 0, 0),
color_wp=(0, 0, 0, 0))
# tta = wp.get_tta_times(case=str(n+8), homedir=homed)
# sub = gf.loc[tta + timedelta(minutes=55)]
subdates = get_subdates(gf, tta)
sub = gf.loc[subdates]
x1, y1 = sub.pdiff.values, sub.ucomp.values
x2, y2 = sub.pdiff.values, sub.wp_ucomp.values
h1 = axes[n].scatter(x1, y1, s=90, marker='x', color=(0, 0, 0.5), lw=2)
h2 = axes[n].scatter(x2, y2, s=100, marker='+', color=(0, 0.5, 0), lw=2)
axes[n].text(0.05,
0.95,
axes[3] = plt.subplot(gs0[3],gid='(d) 16-18Feb04')
params = dict(wdir_surf=130,wdir_wprof=170,
rain_czd=0.25,nhours=2)
try:
tta
except NameError:
tta = np.array(xta.get_tta_dates([2003,2004],params))
tta = tta + timedelta(minutes=55)
cases=[8,9,12,11]
for n,c in enumerate(cases):
gf = gapflow.run(c, plot=True, grid=False,
ax=axes[n], homedir=homed,
color_surf=(0, 0, 0, 0),
color_wp=(0, 0, 0, 0))
# tta = wp.get_tta_times(case=str(n+8), homedir=homed)
# sub = gf.loc[tta + timedelta(minutes=55)]
subdates = get_subdates(gf,tta)
sub = gf.loc[subdates]
x1, y1 = sub.pdiff.values, sub.ucomp.values
x2, y2 = sub.pdiff.values, sub.wp_ucomp.values
h1=axes[n].scatter(x1, y1, s=90, marker='x', color=(0, 0, 0.5),lw=2)
h2=axes[n].scatter(x2, y2, s=100, marker='+', color=(0, 0.5, 0),lw=2)
axes[n].text(0.05,0.95,axes[n].get_gid(),size=14,va='top',
weight='bold',transform=axes[n].transAxes,
backgroundcolor='w',clip_on=True)
