[ncl-talk] check_for_y_lat_coord: Warning:
Noelia otero
noeli1680 at gmail.com
Fri Sep 19 03:09:47 MDT 2014
Hi Mary,
Thanks for your answer. I know that this is a common mistake, but I checked
the units and the coordinates array and they seem correct to me. I only
changed the longitude values, because they span from 347 to 34, so I
changed it to -13 34. This is why I was reattaching all of it, is it
wrong?? Otherwise, I am having problems for trying to plot between 347 and
34 lon..
The error that I get :
(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'
The script is very simple,
;***********************************************************************
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
;************************************************************************
begin
ozfil =
addfile("/work/users/nof/data/data_new/Ozone/grid1/seasons/Europe/
mda8.surfO3.EU.2000.2009.settabs_JJA.nc","r")
wks = gsn_open_wks("X11","test")
;gsn_define_colormap(wks,"BlueDarkRed18")
;ozone
oz = ozfil->MDA8_SurfO3
;change lon values -180 180
lon = ozfil ->lon
lon = where((lon.gt.180).and.(lon.lt.360),lon-360,lon)
lat = ozfil ->lat
oz&lon = lon
; oz&lat = lat ; it might be unnecessary
res = True
res at gsnAddCyclic = False
res at mpMinLonF = min(lon)
res at mpMinLatF = min(lat)
res at mpMaxLonF = max(lon)
res at mpMaxLatF = max(lat)
printVarSummary(oz)
printVarSummary(oz&lat)
printVarSummary(oz&lon)
plot = gsn_csm_contour_map(wks,oz(0,:,:),res)
end
;******************************************
The summary of the variables:
printVarSummary(oz):
Variable: oz
Type: double
Total Size: 13777920 bytes
1722240 values
Number of Dimensions: 3
Dimensions and sizes: [time | 920] x [lon | 48] x [lat | 39]
Coordinates:
time: [20000601..20090831]
lon: [ -13.. 34]
lat: [ 34.. 72]
Number Of Attributes: 2
units : parts per billion = nmol/mol
description : Maximum daily 8-hr average (MDA8) of surface ozone (ppb)
over Europe from 01/01/2000 - 12/31/2009
printVarSummary(lat)
Variable: lat (coordinate)
Type: double
Total Size: 312 bytes
39 values
Number of Dimensions: 1
Dimensions and sizes: [lat | 39]
Coordinates:
Number Of Attributes: 4
axis : Y
units : degrees_north
long_name : latitude
standard_name : latitude
printVarSummary(lon)
Variable: lon (coordinate)
Type: double
Total Size: 384 bytes
48 values
Number of Dimensions: 1
Dimensions and sizes: [lon | 48]
Coordinates:
Number Of Attributes: 4
axis : X
units : degrees_east
long_name : longitude
standard_name : longitude
Many thanks in advance for the support,
Cheers,
Noelia.
2014-09-19 1:07 GMT+02:00 Mary Haley <haley at ucar.edu>:
> Noelia,
>
> When you continue to have problems with a script, it helps to include any
> new error messages you are getting.
>
> Also, your "printVarSummary" of "oz" seems to indicate that you have
> coordinate arrays and the appropriate metadata, so I'm not sure why you are
> reattaching all of it.
>
> Please look at your data right before you plot it, and not right before
> you start adding metadata to it. Also, look at your lat/lon arrays to make
> sure they look okay.
>
> To look at your data more thoroughly, add these three lines right before
> you call gsn_csm_contour_plot:
>
> printVarSummary(oz)
> printVarSummary(oz&lat)
> printVarSummary(oz&lon)
>
> Make sure the "units" for the lat/lon arrays are correct, and that "oz"
> has all the correct coordinate information attached.
>
> --Mary
>
>
> On Thu, Sep 18, 2014 at 12:43 PM, Noelia otero <noeli1680 at gmail.com>
> wrote:
>
>> Hi,
>>
>> Thanks karin, but even with this setting I can't plot. I don't know if
>> this problem might be due to the data..actually, I don't know why this
>> warning..
>> Thanks again!
>>
>> Noelia.
>>
>> 2014-09-18 11:29 GMT+02:00 Karin Meier-Fleischer <meier-fleischer at dkrz.de
>> >:
>>
>>> Hi Noelia,
>>>
>>> maybe you have missed the setting
>>>
>>> oz&lat = lat
>>>
>>> Hope this helps,
>>> Karin
>>>
>>> Am 18.09.14 10:50, schrieb Noelia otero:
>>>
>>> Hi,
>>>
>>> I have problems to plot one field, and although I have already dealt
>>> with this error, I tried to fix it now in different ways..but I don't know
>>> what I am missing, and I have no idea what I am doing wrong. The error is
>>> always:
>>>
>>> *"check_for_y_lat_coord: Warning: Data either does not contain a valid
>>> latitude coordinate array or doesn't contain one at all......"*
>>>
>>> I would appreciate any suggestion:
>>>
>>> ;Variable
>>> oz = ozfil->MDA8_SurfO3
>>> printVarSummary(oz)
>>> ; Summary of variable: oz
>>> ;***Type: double
>>> Total Size: 13777920 bytes
>>> 1722240 values
>>> Number of Dimensions: 3
>>> Dimensions and sizes: [time | 920] x [lon | 48] x [lat | 39]
>>> Coordinates:
>>> time: [20000601..20090831]
>>> lon: [ 347.. 34]
>>> lat: [ 34.. 72]
>>>
>>> ;change lon values into -180 180
>>> lon = ozfil ->lon
>>> lon = where((lon.gt.180).and.(lon.lt.360),lon-360,lon)
>>> lat = ozfil ->lat
>>> lon at long_name = "longitude"
>>> lon at units = "degrees_east"
>>> lat at long_name = "latitude"
>>> lat at units = "degrees_north"
>>> oz&lon = lon
>>> printVarSummary(oz)
>>>
>>> Now the summary of oz is:
>>> Variable: oz
>>> Type: double
>>> Total Size: 13777920 bytes
>>> 1722240 values
>>> Number of Dimensions: 3
>>> Dimensions and sizes: [time | 920] x [lon | 48] x [lat | 39]
>>> Coordinates:
>>> time: [20000601..20090831]
>>> lon: [ -13.. 34]
>>> lat: [ 34.. 72]
>>> Number Of Attributes: 2
>>>
>>> ISo, am trying to plot:
>>> plot = gsn_csm_contour_map(wks,oz(0,:,:),res)
>>>
>>> Any idea about how can I solve this???
>>> Thank you in advance,
>>>
>>> Cheers,
>>>
>>>
>>> Noelia.
>>>
>>>
>>> _______________________________________________
>>> ncl-talk mailing list
>>> List instructions, subscriber options, unsubscribe:http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>
>>>
>>> --
>>> Dipl. Geophys. Karin Meier-Fleischer
>>> Visualization
>>> Application Support
>>>
>>> Deutsches Klimarechenzentrum GmbH (DKRZ)
>>> Bundesstrasse 45a - D20146 Hamburg - Germany
>>>
>>> Phone: +49 (0)40 460094 126
>>> Fax: +49 (0)40 460094 270
>>> E-Mail: meier-fleischer at dkrz.de
>>> URL: www.dkrz.de
>>>
>>> Geschäftsführer: Prof. Dr. Thomas Ludwig
>>> Sitz der Gesellschaft: Hamburg
>>> Amtsgericht Hamburg HRB 39784
>>>
>>>
>>> _______________________________________________
>>> 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/20140919/3f7d40b0/attachment.html
More information about the ncl-talk
mailing list