[ncl-talk] check_for_y_lat_coord: Warning:
Noelia otero
noeli1680 at gmail.com
Fri Sep 19 08:26:42 MDT 2014
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/136cdf08/attachment.html
More information about the ncl-talk
mailing list