Re: vertical cross section problem

Dear Shea, Thanks very much for your helpful advise, the problem was solved thanks very much.   Mohamed Ibrahim El-Sayed OMAR Egyptian Meteorological Authority Scientific Research Department Climate Studies Centre Kobry El-Quobba, P.O.Box. 11784 Cairo, EGYPT Mobil phone  : 002 0160802805  Home Phone : 002 0482370261 ________________________________ From: Dennis Shea <shea_at_ucar.edu> To: Mohamad Ibrahim OMAR <egy_mohamad_at_yahoo.com> Cc: "ncl-talk_at_ucar.edu" <ncl-talk_at_ucar.edu> Sent: Sunday, March 11, 2012 7:15 PM Subject: Re: vertical cross section problem The error message is telling you the nature of the problem. > float level(level) ; > level:long_name = "Height_Index" ; > level:units = "level" ; > level:actual_range = 1050.f, 0.f ; the units of the variable "level" on the file are not one of the recognized units (see message). varave=tmpAve(0,{1000:100},:,{30}) printVarSummary(varave) print("varave&level_at_units="+varave&level_at_units) The plot function only sees the level information associated with only variable 'varave' Before the plot set varave&level_at_units = "hPa" On 3/10/12 11:43 AM, Mohamad Ibrahim OMAR wrote: > Hi all, > I have a problem in plotting the vertical cross section for U wind. > > ******** the error is ; > Variable: level > Type: float > Total Size: 44 bytes > 11 values > Number of Dimensions: 1 > Dimensions and sizes: [level | 11] > Coordinates: > level: [1000..100] > Number Of Attributes: 3 > long_name : Height_Index > units : level > actual_range : ( 1050, 0 ) > (0) 1000 > (1) 925 > (2) 850 > (3) 700 > (4) 500 > (5) 400 > (6) 300 > (7) 250 > (8) 200 > (9) 150 > (10) 100 > (0) gsn_csm_pres_hgt: Fatal: The coordinate array for the first > dimension of the input data must be in Pascals, Hecto-pascals, or millibars > (0) and it must contain the attribute 'units' set to one of the > following strings (depending on your units): > (0) 'mb' 'Mb' 'MB' 'millibar' 'millibars' 'MILLIBARS' > 'hybrid_sigma_pressure' 'Pa' 'pa' 'PA' 'Pascals' 'pascals' 'PASCALS' > 'hpa' 'hPa' 'HPA' 'hecto-pascals' 'HECTO-PASCALS' > (0) Cannot create plot. > fatal:Illegal right-hand side type for assignment > fatal:Execute: Error occurred at or near line 104 in file > U_cross_section_test_10yr.ncl > > ******and the ncdump for the file is ; > > netcdf ATMc-monthlyAVG { > dimensions: > lon = 158 ; > lat = 138 ; > level = 11 ; > time = UNLIMITED ; // (119 currently) > variables: > float lon(lon) ; > lon:long_name = "Longitude" ; > lon:units = "degrees_east" ; > lon:actual_range = -18.04229f, 58.50375f ; > float lat(lat) ; > lat:long_name = "Latitude" ; > lat:units = "degrees_north" ; > lat:actual_range = -7.783508f, 52.90308f ; > float level(level) ; > level:long_name = "Height_Index" ; > level:units = "level" ; > level:actual_range = 1050.f, 0.f ; > double time(time) ; > time:long_name = "Time" ; > time:units = "hours since 1900-1-1 00:00:0.0" ; > time:actual_range = 841848., 928056. ; > float U(time, level, lat, lon) ; > U:long_name = "Zonal Wind" ; > U:units = "m/s" ; > U:missing_value = -1.e+30f ; > float V(time, level, lat, lon) ; > V:long_name = "Meridional Wind" ; > V:units = "m/s" ; > V:missing_value = -1.e+30f ; > float W(time, level, lat, lon) ; > W:long_name = "Vertical Velocity" ; > W:units = "m/s" ; > W:missing_value = -1.e+30f ; > float TK(time, level, lat, lon) ; > TK:long_name = "Temperature" ; > TK:units = "K" ; > TK:missing_value = -1.e+30f ; > float QD(time, level, lat, lon) ; > QD:long_name = "Mixing Ratio" ; > QD:units = "kg/kg" ; > QD:missing_value = -1.e+30f ; > float QC(time, level, lat, lon) ; > QC:long_name = "Cloud Mixing Ratio" ; > QC:units = "kg/kg" ; > QC:missing_value = -1.e+30f ; > float RH(time, level, lat, lon) ; > RH:long_name = "Relative Humidity" ; > RH:units = "fraction" ; > RH:missing_value = -1.e+30f ; > float HGT(time, level, lat, lon) ; > HGT:long_name = "Geopotential Height" ; > HGT:units = "m" ; > HGT:missing_value = -1.e+30f ; > float TH(time, level, lat, lon) ; > TH:long_name = "Potential Temperature" ; > TH:units = "K" ; > TH:missing_value = -1.e+30f ; > float TD(time, level, lat, lon) ; > TD:long_name = "Dew Point Temperature" ; > TD:units = "K" ; > TD:missing_value = -1.e+30f ; > float VOR(time, level, lat, lon) ; > VOR:long_name = "Horz Vorticity" ; > VOR:units = "1/s" ; > VOR:missing_value = -1.e+30f ; > float DIV(time, level, lat, lon) ; > DIV:long_name = "Horz Divergence" ; > DIV:units = "1/s" ; > DIV:missing_value = -1.e+30f ; > float MSE(time, level, lat, lon) ; > MSE:long_name = "Moist Static Energy" ; > MSE:units = "m2/s2" ; > MSE:missing_value = -1.e+30f ; > float PS(time, lat, lon) ; > PS:long_name = "Surface Pressure" ; > PS:units = "hPa" ; > PS:missing_value = -1.e+30f ; > float RT(time, lat, lon) ; > RT:long_name = "Total Precip" ; > RT:units = "mm/day" ; > RT:missing_value = -1.e+30f ; > float TGRND(time, lat, lon) ; > TGRND:long_name = "Ground Temperature" ; > TGRND:units = "K" ; > TGRND:missing_value = -1.e+30f ; > float SMT(time, lat, lon) ; > SMT:long_name = "Total Soil Water" ; > SMT:units = "mm" ; > SMT:missing_value = -1.e+30f ; > float BS(time, lat, lon) ; > BS:long_name = "Base Flow" ; > BS:units = "mm/day" ; > BS:missing_value = -1.e+30f ; > float SLP1(time, lat, lon) ; > SLP1:long_name = "Sea Level Pressure" ; > SLP1:units = "hPa" ; > SLP1:missing_value = -1.e+30f ; > float SLP2(time, lat, lon) ; > SLP2:long_name = "Sea Level Pressure" ; > SLP2:units = "hPa" ; > SLP2:missing_value = -1.e+30f ; > > // global attributes: > :domxmin = -18.04229f ; > :domxmax = 58.50375f ; > :domymin = -7.783508f ; > :domymax = 52.90308f ; > :domzmin = 1050.f ; > :domzmax = 0.f ; > } > > so is there any thing cal help? > THANKS > > > Mohamed Ibrahim El-Sayed OMAR > Egyptian Meteorological Authority > Scientific Research Department > Climate Studies Centre > Kobry El-Quobba, P.O.Box. 11784 > Cairo, EGYPT > Mobil phone : 002 0160802805 > Home Phone : 002 0482370261 > > > _______________________________________________ > 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
Received on Sun Mar 11 12:01:32 2012