[ncl-talk] check_for_y_lat_coord: Warning:
Mary Haley
haley at ucar.edu
Fri Sep 19 08:04:26 MDT 2014
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/9d6ce745/attachment.html
More information about the ncl-talk
mailing list