[ncl-talk] check_for_y_lat_coord: Warning:
Mary Haley
haley at ucar.edu
Fri Sep 19 12:06:06 MDT 2014
Hi,
The problem is that the "oz" array is ordered time x lon x lat, and it
needs to be time x lat x lon in order for NCL to be able to plot it.
Also, your data has NaNs in it, so I used replace_ieeenan to fix these.
Please see attached.
--Mary
On Fri, Sep 19, 2014 at 8:26 AM, Noelia otero <noeli1680 at gmail.com> wrote:
> Hi Mary,
>
> Sure, I put the data in your ftp, in incoming, it is a file with the name *mda8.surfO3.EU.2000.2009.settabs_JJA.nc
> <http://mda8.surfO3.EU.2000.2009.settabs_JJA.nc>*
>
> I am using the 6.2.0 ncl version and Debian 3.2.60-1+deb7u1 x86_64
> GNU/Linux for the system.
>
> Again, thanks a lot for the help.
>
> Noelia
>
> 2014-09-19 16:04 GMT+02:00 Mary Haley <haley at ucar.edu>:
>
>> You're right about lonFlip, sorry about that.
>>
>> Can you provide your data on our ftp so I can take a look? You can give
>> me the information offline if you don't want to share it on ncl-talk:
>>
>> http://www.ncl.ucar.edu/report_bug.shtml#HowToFTP
>>
>> --Mary
>>
>>
>> On Fri, Sep 19, 2014 at 7:59 AM, Noelia otero <noeli1680 at gmail.com>
>> wrote:
>>
>>> Hi Mary,
>>>
>>> Thanks for the suggestion (I had the idea that lonFlip only works for
>>> global longitudes).
>>> Unfortunately, when I try to use lonFlip, the ncl process stops and I
>>> receive the message:
>>>
>>> *(0) lonflip: longitude dimension size must be even: mlon=39*
>>>
>>> I don't understand much the message.
>>>
>>> Many thanks again,
>>>
>>> Noelia.
>>>
>>> 2014-09-19 15:02 GMT+02:00 Mary Haley <haley at ucar.edu>:
>>>
>>>> Hi Noelia,
>>>>
>>>> Thanks for the additional information. I admit I'm stumped, because it
>>>> looks like you are doing everything right.
>>>>
>>>> However, there's an easier way to fix the longitudes, using the
>>>> "lonFlip" function. Try this code instead:
>>>>
>>>> ;***********************************************************************
>>>> 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
>>>> oz = lonFlip(oz)
>>>>
>>>> res = True
>>>> res at gsnAddCyclic = False
>>>> res at mpMinLonF = min(oz&lon)
>>>> res at mpMinLatF = min(oz&lat)
>>>> res at mpMaxLonF = max(oz&lon)
>>>> res at mpMaxLatF = max(oz&lat)
>>>>
>>>>
>>>> printVarSummary(oz)
>>>> printVarSummary(oz&lat)
>>>> printVarSummary(oz&lon)
>>>> plot = gsn_csm_contour_map(wks,oz(0,:,:),res)
>>>>
>>>> end
>>>>
>>>> --Mary
>>>>
>>>>
>>>> On Fri, Sep 19, 2014 at 3:09 AM, Noelia otero <noeli1680 at gmail.com>
>>>> wrote:
>>>>
>>>>> 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/58a1135b/attachment.html
-------------- next part --------------
A non-text attachment was scrubbed...
Name: mda.png
Type: image/png
Size: 175458 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20140919/58a1135b/attachment.png
-------------- next part --------------
A non-text attachment was scrubbed...
Name: mda.ncl
Type: application/octet-stream
Size: 1837 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20140919/58a1135b/attachment.obj
More information about the ncl-talk
mailing list