Re: [ coordinates : lat lon (0) check_for_y_lat_coord: Warning: Data either does not contain a valid latitude coordinate array or doesn't contain one at all...

Sorry, I had a typo. I had you do this:

> media@lat2d=f[:]->lat
> media@lon2d=f[:]->lon

when it should be:

> media@lat2d=f[0]->lat
> media@lon2d=f[0]->lon

Note, and this is important: this assumes that your lat/lon values are the same across all data files. Is this the case?

--Mary

On May 13, 2014, at 1:26 PM, Vanúcia Schumacher <vanucia-schumacher@hotmail.com> wrote:

> Hi Mary,
>
> I did what you suggested, but now new error appeared:
>
> fatal:Eq: Dimension size, for dimension number 0, of operands does not match, can't continue
>
> fatal:["Execute.c":7743]:Execute: Error occurred at or near line 887 in file $NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl
>
> fatal:["Execute.c":7743]:Execute: Error occurred at or near line 1508 in file $NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl
>
> fatal:["Execute.c":7743]:Execute: Error occurred at or near line 6845 in file $NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl
>
> fatal:["Execute.c":7743]:Execute: Error occurred at or near line 57 in file tsm.ncl
>
> Look data print:
>
> dimensions:
> time = 360 // unlimited
> j = 149
> i = 182
> bnds = 2
> vertices = 4
> variables:
> double time ( time )
> bounds : time_bnds
> units : days since 1981-01-01 00:00:00
> calendar : noleap
> axis : T
> long_name : time
> standard_name : time
>
> double time_bnds ( time, bnds )
>
> integer j ( j )
> units : 1
> long_name : cell index along second dimension
>
> integer i ( i )
> units : 1
> long_name : cell index along first dimension
>
> float lat ( j, i )
> standard_name : latitude
> long_name : latitude coordinate
> units : degrees_north
> bounds : lat_vertices
>
> float lon ( j, i )
> standard_name : longitude
> long_name : longitude coordinate
> units : degrees_east
> bounds : lon_vertices
>
> float lat_vertices ( j, i, vertices )
> units : degrees_north
>
> float lon_vertices ( j, i, vertices )
> units : degrees_east
>
> float tos ( time, j, i )
> standard_name : sea_surface_temperature
> long_name : Sea Surface Temperature
>
>
> Script:
>
> 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"
>
> begin
>
> fils = systemfunc("ls *.nc")
> f = addfiles(fils, "r")
> ListSetType(f,"join")
> t =f[:]->tos(:,0,:,:)
> tsm = t - 273.15
> copy_VarCoords(t, tsm)
> media= dim_avg_n_Wrap(tsm, 0)
>
> media@lat2d=f[:]->lat
> media@lon2d=f[:]->lon
>
> wkstype= "png"
> wkstype@wkWidth = 1000
> wkstype@wkHeight = 1000
> wks = gsn_open_wks(wkstype,"tsm")
> gsn_define_colormap(wks,"BlAqGrYeOrRe")
> res = True
> ;res@mpMinLonF = 0
> ;res@mpMaxLonF = 360
> ;res@mpCenterLonF = 180
> res@tiMainString = "TSM"
> res@gsnLeftString = "CCCma-CanCM4"
> res@cnLinesOn = False
> res@cnFillOn = True
> plot = gsn_csm_contour_map_ce(wks,media,res)
> end
>
>
>
>
> ---
> Vanúcia Schumacher
> Mestranda em Meteorologia - UFV
> Meteorologista -UFPel
> Departamento de Meteorologia Agrícola - DEA
> Cel: (31) 9978 2522
> DEA: (31) 3899 1890
>
>
> > Subject: Re: ERROR[ coordinates : lat lon (0) check_for_y_lat_coord: Warning: Data either does not contain a valid latitude coordinate array or doesn't contain one at all...
> > From: haley@ucar.edu
> > Date: Tue, 13 May 2014 09:47:28 -0600
> > CC: ncl-talk@ucar.edu
> > To: vanucia-schumacher@hotmail.com
> >
> > Vanúcia,
> >
> > It looks to me like you have a curvilinear grid and not a rectilinear grid, so you are going to be unable to use coordinate arrays, which is what you were trying to do with this code:
> >
> > > >lat = f[:]->i
> > > >lon = f[:]->j
> > > > t&i = lat
> > > >t&j = lon
> >
> > The "printVarSummary" of "t" indicates that your "t" variable has coordinates called "lat" and "lon", and hopefully these are on the file, and are two-dimensional. If so, try this (UNTESTED):
> >
> > 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"
> >
> > begin
> > fils = systemfunc("ls *.nc")
> > f = addfiles(fils, "r")
> > ListSetType(f,"join")
> > t =f[:]->tos(:,0,:,:
> >
> > tsm = t - 273.15
> > media= dim_avg_n_Wrap(tsm, 0)
> > media@lat2d = lat = f[:]->lat
> > media@lon2d = f[:]->lon
> > res = True
> > plot = gsn_csm_contour_map_ce(wks,media,res)
> > end
> >
> >
> > Note that I'm not setting any plot options via "res" here, so you are going to get a very bare-bones plot.
> >
> > For a sample of some plot options you might want to set, see the templates at:
> >
> > http://www.ncl.ucar.edu/Applications/Templates/#PlotsOverMapsTemplates
> >
> > In particular, look at "contour_map_template.ncl".
> >
> > --Mary
> >
> > On May 12, 2014, at 1:35 PM, Vanúcia Schumacher <vanucia-schumacher@hotmail.com> wrote:
> >
> > > Hi users,
> > >
> > > I'll trying to plot global map, but lat and lon have different names in variable:
> > >
> > > >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"
> > >
> > > >begin
> > > >fils = systemfunc("ls *.nc")
> > > >f = addfiles(fils, "r")
> > > >ListSetType(f,"join")
> > >
> > > >t =f[:]->tos(:,0,:,:)
> > >
> > > >lat = f[:]->i
> > > >lon = f[:]->j
> > > > t&i = lat
> > > >t&j = lon
> > > >tsm = t - 273.15
> > > >media= dim_avg_n_Wrap(tsm, 0)
> > > >plot = gsn_csm_contour_map_ce(wks,media,res)
> > > >end
> > >
> > > >printVarSummary(t)
> > >
> > >
> > > Variable: t
> > > Type: float
> > > Total Size: 108472 bytes
> > > 27118 values
> > > Number of Dimensions: 3
> > > Dimensions and sizes: [ncl_join | 1] x [j | 149] x [i | 182]
> > > Coordinates:
> > > j: [1..149]
> > > i: [1..182]
> > > Number Of Attributes: 14
> > > time : 15.5
> > > standard_name : sea_surface_temperature
> > > long_name : Sea Surface Temperature
> > > comment : "this may differ from ""surface temperature"" in regions of sea ice."
> > > units : K
> > > original_name : sosstsst
> > > original_units : degC
> > > history : 2011-11-09T14:11:22Z altered by CMOR: Converted units from 'degC' to 'K'. 2011-11-09T14:11:22Z altered by CMOR: replaced missing value flag (9.96921e+36) with standard missing value (1e+20).
> > > cell_methods : time: mean (interval: 30 minutes)
> > > cell_measures : area: areacello
> > > missing_value : 1e+20
> > > _FillValue : 1e+20
> > > associated_files : baseURL: http://cmip-pcmdi.llnl.gov/CMIP5/dataLocation gridspecFile: gridspec_ocean_fx_IPSL-CM5A-LR_decadal1980_r0i0p0.nc areacello: areacello_fx_IPSL-CM5A-LR_decadal1980_r0i0p0.nc
> > > coordinates : lat lon
> > > (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'
> > >
> > > Help me!
> > >
> > > ---
> > > Vanúcia Schumacher
> > > Mestranda em Meteorologia - UFV
> > > Meteorologista -UFPel
> > > Departamento de Meteorologia Agrícola - DEA
> > > Cel: (31) 9978 2522
> > > DEA: (31) 3899 1890
> > > _______________________________________________
> > > 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

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue May 13 13:42:23 2014