Youou could set
dataOpts@ThermoInfo = False
You will get a sounding plot with *no thermo dynamic info*
---- Actually, the underlying CAPE code requires that the input pressure levels be in descending order. Try p = histo(a:b:-1,1) tem = histo(a:b:-1,2) rh = histo(a:b:-1,3) u = histo(a:b:-1,4) v = histo(a:b:-1,5) alo, later wdir = atan2(u,v)*57.29579 + 180. On 2/15/12 11:56 AM, Mary Haley wrote: > The code appears to be bogging down in cape_thermo. I'll look into this. > > --Mary > > On Feb 14, 2012, at 8:34 AM, ahmed lasheen wrote: > >> Hello >> i am trying to produce the skew t ,but when i run the script , it freeze. >> i don't know what is wrong >> i have attached the script and the data. >> thanks in advance >> >> -- >> =============== >> Ahmed Lasheen >> Junior researcher at Cairo Numerical Weather Prediction Center (CNWPC) >> Egyptian Meteorological Authority(EMA) >> Cairo,Egypt >> =============== >> >> <wrf_nmm_skew_histogram.ncl><c7>_______________________________________________ >> ncl-talk mailing list >> List instructions, subscriber options, unsubscribe: >> http://mailman.ucar.edu/mailman/listinfo/ncl-talk > > _______________________________________________ > ncl-talk mailing list > List instructions, subscriber options, unsubscribe: > http://mailman.ucar.edu/mailman/listinfo/ncl-talk _______________________________________________ ncl-talk mailing list List instructions, subscriber options, unsubscribe: http://mailman.ucar.edu/mailman/listinfo/ncl-talkReceived on Wed Feb 15 21:34:56 2012
This archive was generated by hypermail 2.1.8 : Fri Feb 17 2012 - 08:50:10 MST