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()