if c in [0, 1, 2]:
matfile = source + matfs[0]
else:
matfile = source + matfs[1]
t0 = dates[c]['t0']
t1 = dates[c]['t1']
hax = of.plot(
ax[c],
matfile,
t0,
t1,
legend=leg_bool,
add=('bulk', 'upslope'),
lcolors=(cmap(0), cmap(1)),
ylims=([0, 100], [0, 35]),
second_axis='bulk',
lw=2,
legend_loc=(0.75, 0.7, 0.2, 0.2),
# ylim = [0, dates[c]['vmax']],
xtickfreq='12H',
)
if c == 5:
txt = r'$\leftarrow UTC \left[\stackrel{day}{time}\right]$'
hax[0].set_xlabel(txt, fontsize=20)
if c not in [0]:
hax[0].set_ylabel('')
if c == 0:
legend = True
else:
legend = False
matfile = source + matfs[0]
t0 = dates[c]['t0']
t1 = dates[c]['t1']
hax = of.plot(
axes[c],
matfile,
t0,
t1,
legend=legend,
add=('bulk', 'upslope'),
lcolors=(cmap(0), cmap(1)),
ylims=([0, 60], [0, 20]),
second_axis='bulk',
lw=2,
legend_loc=(0.35, 0.7, 0.2, 0.2),
xtickfreq='6H',
)
axes[c].text(0.98,
0.95,
axes[c].get_gid(),
size=14,
va='top',
ha='right',
weight='bold',
transform=axes[c].transAxes,
legend = True
else:
legend = False
if c in [0, 1, 2]:
matfile = source + matfs[0]
else:
matfile = source + matfs[1]
t0 = dates[c]['t0']
t1 = dates[c]['t1']
hax = of.plot(ax[c],
matfile,
t0,
t1,
legend=legend,
add=['bulk', 'upslope'],
legend_loc=(0.7, 0.7, 0.2, 0.2),
ylim=[0, dates[c]['vmax']],
xtickfreq='12H',
lw=2)
if c not in [0]:
hax.set_ylabel('')
if c not in [5]:
hax.set_xlabel('')
ax[c].text(0.05,
0.95,
ax[c].get_gid(),
size=14,
if c == 0:
legend=True
else:
legend=False
matfile=source+matfs[0]
t0=dates[c]['t0']
t1=dates[c]['t1']
hax = of.plot(axes[c],
matfile,t0,t1,
legend = legend,
add = ('bulk','upslope'),
lcolors = (cmap(0),cmap(1)),
ylims = ([0,60], [0,20]),
second_axis = 'bulk',
lw = 2,
legend_loc = (0.35,0.7,0.2,0.2),
xtickfreq = '6H',
)
axes[c].text(0.98,0.95,axes[c].get_gid(),size=14,
va='top',ha='right',
weight='bold',transform=axes[c].transAxes,
backgroundcolor='w',clip_on=True)
#plt.show()
fname='/home/raul/Desktop/rainfall_forcing_airborne.png'
plt.savefig(fname, dpi=300, format='png',papertype='letter',
if c == 3:
legend=True
else:
legend=False
if c in [0,1,2]:
matfile=source+matfs[0]
else:
matfile=source+matfs[1]
t0=dates[c]['t0']
t1=dates[c]['t1']
hax = of.plot(ax[c],matfile,t0,t1,
legend=legend,add=['bulk','upslope'],
legend_loc=(0.7,0.7,0.2,0.2),
ylim=[0,dates[c]['vmax']],
xtickfreq='12H',
lw=2)
if c not in [0]:
hax.set_ylabel('')
if c not in [5]:
hax.set_xlabel('')
ax[c].text(0.05,0.95,ax[c].get_gid(),size=14,va='top',
weight='bold',transform=ax[c].transAxes,
backgroundcolor='w',clip_on=True)