Re: [error]fatal:(lon) is not a named dimension in variable (var).

Hi Dennis,
I added these lines to change the name of i and j for lat and lon, and th=
e following error appears:
var = f[:]->sitvar!0="time"var!1="lat"var!2="lon"i=var!1j=var!2
Variable: varType: floatTotal Size: 39049920 bytes 9762480 value=
sNumber of Dimensions: 3Dimensions and sizes: [time | 360] x [lat | 149] x =
[lon | 182]Coordinates: time: [15.5..10934.5] lat: [=
1..149] lon: [1..182]Number Of Attributes: 12
fatal:linint2: The rightmost dimensions of fi must be nyi x nxi, where ny=
i and nxi are the lengths of yi and xi respectivelyfatal:["Execute.c":8567]=
:Execute: Error occurred at or near line 40 in file resolucao_2g.ncl
But if I not change values of i and j, and change the line to new_var =
= linint2_Wrap(var&i,var&j,var,True,new_lat,new_lon,0) it's =
right (run).
However, does not plot any data when using the new data generated (with n=
ew resolution).

---Vanúcia SchumacherMestranda em Meteorologia - UFVMeteorologista -UFPel
Departamento de Meteorologia Agrícola - DEACel: (31) 9978 2522 DEA: (31) =
3899 1890

Date: Fri, 13 Jun 2014 09:03:07 -0600
Subject: Re: [error]fatal:(lon) is not a named dimension in vari=
able (var).
From: shea@ucar.edu
To: vanucia-schumacher@hotmail.com
CC: ncl-talk@ucar.edu

The most important rule in data processing is *look at your data*. From pri=
ntVarSummary(var)

Variable: varDimensions and sizes: [time | 360] x [j | 149] x [i | 182] =
Coordinates: time: [15.5..10934.5] j: [1..149] =
      i: [1..182]

There is no 'lat' dimension name. In fact, there are no latitude/longitud=
e coordinates on the file at all.

There is an attribute that tells you where you can get the file that contai=
ns the needed information.

associated_files : baseURL: http://cmip-pcmdi.llnl.gov/CMIP5/dataLocation g=
ridspecFile: gridspec_seaIce_fx_IPSL-CM5A-LR_decadal1980_r0i0p0.nc areacell=
o: areacello_fx_IPSL-CM5A-LR_decadal1980_r0i0p0.nc

This is your responsibility to get this file. Unfortunately, this site ap=
pears to be obsolete.





On Thu, Jun 12, 2014 at 1:27 PM, Vanúcia Schumacher <vanucia-schuma=
cher@hotmail.com> wrote:

Variable: varType: floatTotal Size: 39049920 bytes 9762480 value=
sNumber of Dimensions: 3 =
Dimensions and sizes: [time | 360] x [j | 149] x [i | 182]Coordinates: =
        time: [15.5..10934.5] j: [1..149] i: [1..182]=

Number Of Attributes: 12 standard_name : sea_ice_thickness long_name : Se=
a Ice Thickness comment : the mean thickness of sea ice in the ocean porti=
on of the grid cell (averaging over the entire ocean portion, including t=
he ice-free fraction). Reported as 0.0 in regions free of sea ice. =
  units : m original_name : iicethic cell_methods : time: mean (interval:=
 30 minutes) area: mean where sea =
  cell_measures : area: areacello history : 2011-11-11T12:43:09Z altered b=
y CMOR: replaced missing value flag (9.96921e+36) with standard missing val=
ue (1e+20). =
  missing_value : 1e+20 _FillValue : 1e+20 associated_files : baseURL: ht=
tp://cmip-pcmdi.llnl.gov/CMIP5/dataLocation gridspecFile: gridspec_seaIce_f=
x_IPSL-CM5A-LR_decadal1980_r0i0p0.nc areacello: areacello_fx_IPSL-CM5A-LR_d=
ecadal1980_r0i0p0.nc =
  coordinates : lat lon


From: asphilli@ucar.edu
Date: Thu, 12 Jun 2014 09:27:59 -0600

Subject: Re: [error]fatal:(lon) is not a named dimension in vari=
able (var).
To: shea@ucar.edu
CC: vanucia-schumacher@hotmail.com; ncl-talk@ucar.edu

Hi Vanúcia,
To add on to what Dennis is saying: From the printVarSummary you should be =
able to see if lon is a named dimension of var. I am guessing that it is no=
t, and that NCL is failing at this line because of it:

new_var = linint2_Wrap(var&lon,var&lat,var,True,new_lon,new_lat=
,0)

Adam

On Thu, Jun 12, 2014 at 8:16 AM, Dennis Shea <shea@ucar.edu> wrote:

It would be better if you showed the output from printVarSummary

var = f[:]->sit

printVarSummary(var) ; <==============
=

newfile = "sit1.nc"
system("rm -f "+newfile) =

file_create = addfile(newfile, "c") ;criando arquivo.

new_lon = fspan(0, 358, 180) ; usually no cyclic point
new_lon@units = "degrees_east"

new_lon!0 = "lon"

new_lat = fspan(-87.86, 87.86, 90)
new_lat@units = "degrees_north"

new_lat!0 = "lat"

new_var = linint2_Wrap(var&lon,var&lat,var,True,new_lon,new_lat=
,0)printVarSummary(new_var)

file_create->media = new_var

On Thu, Jun 12, 2014 at 7:19 AM, Vanúcia Schumacher <vanucia-schuma=
cher@hotmail.com> wrote:

Hi users,
I need to interpolate the data for resolution of 2 degrees, and I'm havi=
ng this error: =

[error]fatal:(lon) is not a named dimension in variable (var).

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"load "$NCARG_ROOT/l=
ib/ncarg/nclscripts/csm/gsn_csm.ncl" =

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"load "$NCARG_ROO=
T/lib/ncarg/nclscripts/csm/shea_util.ncl"

beginfils = systemfunc("ls *.nc") f = addfiles(fils, "r") =
                 =

ListSetType(f,"join") t =f[:]->sitvar = f[:]->sittime=f[:]->t=
ime =

system("rm -f teste2.nc") file_create = addfile("sit1.nc", "c") ;c=
riando arquivo. =

new_lon = fspan(0, 360, 180) new_lat = fspan(-87.86, 87.86=
, 90) new_var = linint2_Wrap(var&lon,var&lat,var,True,ne=
w_lon,new_lat,0) =

new_var!0="time"new_var!1 = "lon"new_var!2 ="lat"new_var&time=time=

new_var&lat = new_latnew_var&lon = new_lonnew_var&lat@units = "degree=
s_north"new_var&lon@units = "degrees_east" =

printVarSummary(new_var)file_create->media = new_varend

print data: =

 dimensions: =

      time = 120 // unlimited lat = 64 =

      lon = 128 =

      bnds = 2 variables: =

      double time ( time ) bounds : time_bnds =

         units : days since 1850-1-1 =

         calendar : 365_day =

         axis : T =

         long_name : time =

         standard_name : time =

      double time_bnds ( time, bnds ) =

      double lat ( lat ) =

         bounds : lat_bnds =

         units : degrees_north =

         axis : Y =

         long_name : latitude =

         standard_name : latitude =

      double lat_bnds ( lat, bnds ) =

      double lon ( lon ) bounds : lon_bnds =

         units : degrees_east =

         axis : X =

         long_name : longitude =

         standard_name : longitude =

      double lon_bnds ( lon, bnds ) =

      float sit ( time, lat, lon ) =

         standard_name : sea_ice_thickness =

         long_name : Sea Ice Thickness

---Vanúcia Schumacher =

Mestranda em Meteorologia - UFVMeteorologista -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

--
Adam Phillips
Associate Scientist, Climate and Global Dynamics Division, NCAR

www.cgd.ucar.edu/staff/asphilli/ 303-497-1726 =

                                                =

_______________________________________________

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 Fri Jun 13 03:55:24 2014