[ncl-talk] check_for_y_lat_coord: Warning:
Noelia otero
noeli1680 at gmail.com
Fri Sep 19 07:59:33 MDT 2014
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/0bbd2463/attachment.html
More information about the ncl-talk
mailing list