[ncl-talk] check_for_y_lat_coord: Warning:
Noelia otero
noeli1680 at gmail.com
Sat Sep 20 02:28:14 MDT 2014
Hi Mary,
It works now.
Thanks a lot for your help and your time,
cheers,
Noelia.
2014-09-19 20:06 GMT+02:00 Mary Haley <haley at ucar.edu>:
> 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/20140920/d61796be/attachment.html
More information about the ncl-talk
mailing list