[ncl-talk] error in plot from data set of different number of grid points
tms_l
littithomas at gmail.com
Thu Dec 18 16:22:49 MST 2014
hi,
variable with dimension of lat=145 .the number of lat points is 145 in one
variable..its not a typo.
however, the other two variables (ud & vd) having number of points , lat
=144..
but the number of longitude points are same..I think this will also make an
error..can i get any suggestion to deal with the data sets of different
number of latitude points???ultimately , i would like to plot these three
data sets together..
thank you so much
On Fri, Dec 19, 2014 at 1:19 AM, Dennis Shea <shea at ucar.edu> wrote:
>
> There are no coordinate variables associated with your variables.
> How would NCL know where to plot values??
>
> For example
>
> printVarSummary(u)
>
> Variable: u
> Type: float
> Number of Dimensions: 2
> Dimensions and sizes: [lat | 96] x [lon | 144]
> Coordinates:
> lat: [ -90..90] <===== These tell ncl where to ploat
> lon: [ 0..357.5] <=====
> Number Of Attributes: 3
> units : m/s
> long_name : Zonal wind
>
> =================
> Are there lat/lon values available? If so, one approach would be
>
> u!0 = "lat"
> u!1 = "lon"
> u&lat = lat
> u&lon = lon
> printVarSummary(u)
>
> --
> Also, the variable with dimension of lon=145 ... Does that have the cyclic
> point
> contained or was that a typo?
>
>
> On Wed, Dec 17, 2014 at 5:17 PM, tms_l <littithomas at gmail.com> wrote:
>>
>> Hi,
>>
>> I am plotting ud,vd & psl_clim using the following plot function .
>>
>> *plot=gsn_csm_vector_scalar_map(wks,ud,vd,psl_clim,res)*
>>
>> Variable: *ud*
>> Type: float
>> Total Size: 110592 bytes
>> 27648 values
>> Number of Dimensions: 2
>> Dimensions and sizes: [144] x [192]
>> Coordinates:
>> Number Of Attributes: 1
>> _FillValue : -1.073742e+09
>>
>>
>> Variable: *vd*
>> Type: float
>> Total Size: 110592 bytes
>> 27648 values
>> Number of Dimensions: 2
>> Dimensions and sizes: [144] x [192]
>> Coordinates:
>> Number Of Attributes: 1
>> _FillValue : -1.073742e+09
>>
>>
>> Variable: *psl_clim*
>> Type: float
>> Total Size: 111360 bytes
>> 27840 values
>> Number of Dimensions: 2
>> Dimensions and sizes: [145] x [192]
>> Coordinates:
>> Number Of Attributes: 3
>> long_name : MSLP(hPa)& wind-850hPa (a)
>> _FillValue : 9.96921e+36
>>
>>
>>
>> Thanks in advance,
>>
>> Thomas
>>
>>
>>
>> On Thu, Dec 18, 2014 at 6:26 AM, Mary Haley <haley at ucar.edu> wrote:
>>>
>>> Hi Thomas,
>>>
>>> In order to plot data over a map, either as vectors or contours, you
>>> have to provide latitude and longitude locations for your data.
>>>
>>> If you have a rectilinear grid, where the lat / lon are attached
>>> directly to the data as 1-dimensional coordinate variables, then the
>>> plotting routines will use these automatically.
>>>
>>> If you have 2-dimensional lat / lon coordinate points, then you either
>>> need to attach these to your data as special "lat2d" / "lon2d" attributes,
>>> or set the sfYArray / sfXArray resources.
>>>
>>> From your description, it sounds like you have a rectilinear grid, but
>>> your data variable may not be constructed properly. The error message is
>>> saying that it can't determine if you have valid lat / lon coordinate
>>> variables, which can happen if your "units" attribute is not one of the
>>> recognized values, like "degrees_north" or "degrees_east", and/or the
>>> coordinate variables somehow got stripped off your data variable. This can
>>> happen if you do a calculation on your data and put it in a new array,
>>> causing the metadata to be stripped off.
>>>
>>> Please do a "printVarSummary" on the data variables you are trying to
>>> plot, and send it back to ncl-talk so we can see what it looks like.
>>>
>>> Thanks,
>>>
>>> --Mary
>>>
>>>
>>> On Tue, Dec 16, 2014 at 8:23 PM, tms_l <littithomas at gmail.com> wrote:
>>>
>>>> Dear ncl users
>>>>
>>>> I would like to plot u & v components of wind along with mean sea level
>>>> pressure data.The wind data (time,level,lat_1,lon_1) is 4 dimensional and
>>>> mean sea level pressure data (time,lat,lon) 3 dimensional.
>>>>
>>>> Also the number of grids points are arranged differently in each data
>>>> sets.
>>>>
>>>> For wind data latitude , lat_1 = 144 points
>>>> longitude , lon_1 =192 points
>>>>
>>>> For mean sea level pressure latitude , lat =145
>>>> longitude, lon =192
>>>>
>>>> I am getting the error message as follows and the mean sea level
>>>> pressure data seems to plotted reversely..Is there any function to deal
>>>> with data sets of different grid points or how can solve this issue? Any
>>>> help will be appreciated.
>>>>
>>>>
>>>> error message :::::::
>>>> (0) check_for_y_lat_coord: Warning: Data either does not contain a
>>>> valid latitude coordinate array or doesn't contain one at all.
>>>> (0) A valid latitude coordinate array should have a 'units'
>>>> attribute equal to one of the following values:
>>>> (0) 'degrees_north' 'degrees-north' 'degree_north' 'degrees
>>>> north' 'degrees_N' 'Degrees_north' 'degree_N' 'degreeN' 'degreesN' 'deg
>>>> north'
>>>> (0) check_for_lon_coord: Warning: Data either does not contain a
>>>> valid longitude coordinate array or doesn't contain one at all.
>>>> (0) A valid longitude coordinate array should have a 'units'
>>>> attribute equal to one of the following values:
>>>> (0) 'degrees_east' 'degrees-east' 'degree_east' 'degrees east'
>>>> 'degrees_E' 'Degrees_east' 'degree_E' 'degreeE' 'degreesE' 'deg east'
>>>> (0) check_for_y_lat_coord: Warning: Data either does not contain a
>>>> valid latitude coordinate array or doesn't contain one at all.
>>>> (0) A valid latitude coordinate array should have a 'units'
>>>> attribute equal to one of the following values:
>>>> (0) 'degrees_north' 'degrees-north' 'degree_north' 'degrees
>>>> north' 'degrees_N' 'Degrees_north' 'degree_N' 'degreeN' 'degreesN' 'deg
>>>> north'
>>>> (0) check_for_lon_coord: Warning: Data either does not contain a
>>>> valid longitude coordinate array or doesn't contain one at all.
>>>> (0) A valid longitude coordinate array should have a 'units'
>>>> attribute equal to one of the following values:
>>>> (0) 'degrees_east' 'degrees-east' 'degree_east' 'degrees east'
>>>> 'degrees_E' 'Degrees_east' 'degree_E' 'degreeE' 'degreesE' 'deg east'
>>>>
>>>>
>>>>
>>>> Regards
>>>>
>>>> Thomas
>>>>
>>>> _______________________________________________
>>>> 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
>>
>>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20141219/41861a1c/attachment.html
More information about the ncl-talk
mailing list