def main():
''' raf files are created using matlab function
W:/MATLAB/p3_matlab/convert_to_raf_ascii.m '''
f = '/home/rvalenzuela/P3_stdtape/ascii/010123I.std.ascii.raf'
ini, end = get_sounding_times()
m = 0
for i, e in zip(ini, end):
if m >= 0:
sound = mf.parse_acft_sounding(f, i, e, return_interp=True)
# print sound
temp = sound.AIR_TEMP.values #[C]
dewp = sound.DEW_POINT.values #[C]
wspd = sound.WIND_SPD.values
wdir = sound.WIND_DIR.values
u, v = get_components(wspd, wdir)
pres = sound.AIR_PRESS.values
lats = sound.LAT.values
lons = sound.LON.values
hgt = sound.index.values
theta = sound.theta.values
thetaeq = sound.thetaeq.values
bvf_dry = sound.bvf_dry.values
bvf_moist = sound.bvf_moist.values
location = [np.average(lats), np.average(lons)]
# plot_skew1(temp=temp, dewp=dewp, u=u, v=v, press=pres, date=[i, e], loc=location)
plot_skew2(temp=temp,
dewp=dewp,
u=u,
v=v,
press=pres,
date=[i, e],
loc=location,
hgt=hgt)
# plot_thermo(temp=temp, dewp=dewp, u=u, v=v, hgt=hgt,theta=theta, thetaeq=thetaeq,
# bvf_dry=bvf_dry, bvf_moist=bvf_moist,press=pres, date=[i, e], loc=location,top=5000)
m += 1
# compare_potential_temp(df,date)
# break
plt.show()
def main(): ''' raf files are created using matlab function W:/MATLAB/p3_matlab/convert_to_raf_ascii.m ''' f='/home/rvalenzuela/P3_stdtape/ascii/010123I.std.ascii.raf' ini,end = get_sounding_times() m=0 for i,e in zip(ini,end): if m>=0: sound=mf.parse_acft_sounding(f, i, e, return_interp=True) # print sound temp=sound.AIR_TEMP.values #[C] dewp=sound.DEW_POINT.values #[C] wspd=sound.WIND_SPD.values wdir=sound.WIND_DIR.values u,v = get_components(wspd,wdir) pres=sound.AIR_PRESS.values lats=sound.LAT.values lons=sound.LON.values hgt=sound.index.values theta=sound.theta.values thetaeq=sound.thetaeq.values bvf_dry=sound.bvf_dry.values bvf_moist=sound.bvf_moist.values location=[np.average(lats),np.average(lons)] # plot_skew1(temp=temp, dewp=dewp, u=u, v=v, press=pres, date=[i, e], loc=location) plot_skew2(temp=temp, dewp=dewp, u=u, v=v, press=pres, date=[i, e], loc=location, hgt=hgt) # plot_thermo(temp=temp, dewp=dewp, u=u, v=v, hgt=hgt,theta=theta, thetaeq=thetaeq, # bvf_dry=bvf_dry, bvf_moist=bvf_moist,press=pres, date=[i, e], loc=location,top=5000) m+=1 # compare_potential_temp(df,date) # break plt.show()
