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@units="+varave&level@units)

The plot function only sees the level information
associated with only variable 'varave'

Before the plot set

varave&level@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 07:15:39 2012