pamtra = read_prepare(hydroset=campaign + '_all_hydro_', minhour=minhour)
pamtraTl0 = slice_data(pamtra, 'T', maxvalue=0)
lognormrule = True
#%% CFAD Sensitivity
xlim = [-60, 50]
ylim = [0, 12000]
h1, x1, y1 = hist_and_plot(pamtra,
'sensitivity cone X',
yvar='Hgt',
xvar='Z10',
xlabel='Zx [dBZ]',
ylabel='Hgt [m]',
xlim=xlim,
ylim=ylim,
lognorm=lognormrule,
savename='tripex/3f/pamtraCFAD_Zx_Hgt.png',
inverty=False,
figax=None,
bins=(100, icon150heights[::-1]),
density=False,
CFAD=False)
#p = pamtra[['Z10', 'Hgt']].dropna()
#hst, xedge, yedge = np.histogram2d(p['Z10'], p['Hgt'], bins=(100,icon150heights[::-1]))
hist_and_plot(pamtra,
'sensitivity cone Ka',
yvar='Hgt',
xvar='Z35',
xlabel='Zk [dBZ]',
plt.close('all')
addlabel = ''
addlabel = 'run'
#addlabel = 'spin'
campaign = 'tripex'
hydroset = campaign + '_all_hydro_'
data = read_prepare(hydroset, maxhour=46.0, minhour=6.0)
lognorm = True
h, x, y = hist_and_plot(data,
'CFAD X-band SNR',
'Hgt',
'N10',
'SNR X-band [dBZ]',
'Height [m]',
xlim=[-15, 70],
ylim=[0, 10000],
inverty=False,
savename='CFAD/CFAD_snrX_Hgt' + campaign + addlabel +
'.png',
lognorm=lognorm)
h, x, y = hist_and_plot(data,
'CFAD Ka-band SNR',
'Hgt',
'N35',
'SNR Ka-band [dB]',
'Height [m]',
xlim=[-15, 70],
ylim=[0, 10000],
inverty=False,
#import matplotlib.colors as colors
from statistics import hist_and_plot, make6panel_temperature
from READ import read_prepare
plt.close('all')
addlabel = ''
addlabel = 'run'
addlabel = 'spin'
hydroset='all_hydro'
data = read_prepare(hydroset, maxhour=6.0, minhour=0.0)
lognorm=True
xlim = [-10, 8]
ylim = [-5, 20]
h,x,y = hist_and_plot(data, 'DWRxk MDV X-band', 'DWRxk', 'V10',
'MDV X-band [m/s]', 'DWRxk [dB]',
xlim=xlim, ylim=ylim, inverty=False,
savename='DWRvel/DWRxk_MDV_X' + addlabel +'.png', lognorm=lognorm)
h,x,y = hist_and_plot(data, 'DWRxk MDV Ka-band', 'DWRxk', 'V35',
'MDV Ka-band [m/s]', 'DWRxk [dB]',
xlim=xlim, ylim=ylim, inverty=False,
savename='DWRvel/DWRxk_MDV_Ka' + addlabel +'.png', lognorm=lognorm)
h,x,y = hist_and_plot(data, 'DWRxk MDV W-band', 'DWRxk', 'V94',
'MDV W-band [m/s]', 'DWRxk [dB]',
xlim=xlim, ylim=ylim, inverty=False,
savename='DWRvel/DWRxk_MDV_W' + addlabel +'.png', lognorm=lognorm)
h,x,y = hist_and_plot(data, 'DWRkw MDV X-band', 'DWRkw', 'V10',
'MDV X-band [m/s]', 'DWRkw [dB]',
xlim=xlim, ylim=ylim, inverty=False,
plt.close('all')
addlabel = ''
addlabel = 'run'
#addlabel = 'spin'
campaign = 'tripex'
hydroset = campaign + '_all_hydro_'
data = read_prepare(hydroset, maxhour=6.0, minhour=0.0)
lognorm = True
h, x, y = hist_and_plot(data,
'CFAD X-band Ze',
'T',
'Z10',
'Z X-band [dBZ]',
'Temperature [deg C]',
xlim=[-70, 50],
ylim=[-60, 20],
inverty=True,
savename='CFAD/CFAD_X_temp' + campaign + addlabel +
'.png',
lognorm=lognorm)
h, x, y = hist_and_plot(data,
'CFAD Ka-band Ze',
'T',
'Z35',
'Z Ka-band [dBZ]',
'Temperature [deg C]',
xlim=[-70, 50],
ylim=[-60, 20],
inverty=True,
addlabel = ''
addlabel = 'run'
#addlabel = 'spin'
hydroset = 'only_ice'
data = read_prepare(hydroset, maxhour=40.0, minhour=6.0)
lognorm = True
xlim = [-60, 40]
ylim = [0, 1]
h, x, y = hist_and_plot(data,
hydroset + ' SW Z X-band',
'W10',
'Z10',
'Z X-band [dBZ]',
'SW X-band [m/s]',
xlim=xlim,
ylim=ylim,
inverty=False,
savename='SW_Z/SW_Z_X' + addlabel + '_' + hydroset +
'.png',
lognorm=lognorm)
h, x, y = hist_and_plot(data,
hydroset + ' SW Z Ka-band',
'W35',
'Z35',
'Z Ka-band [dBZ]',
'SW Ka-band [m/s]',
xlim=xlim,
ylim=ylim,
inverty=False,