Re: Variable input error: what am I doing wrong?

I checked last night. We'*ARE* running v6.0.0:
[root@scoopdata gomexstep3]# ncl
  Copyright (C) 1995-2011 - All Rights Reserved
  University Corporation for Atmospheric Research
  NCAR Command Language Version 6.0.0
  The use of this software is governed by a License Agreement.
  See http://www.ncl.ucar.edu/ for more details.

but the problem persisted. Not a real problem. I'm looking for SST and
SSH at this point, so I snagged the top-level temperature and I'm good
on this.

Now, however, another problem with the same netcdf input. The snippet of
code where I'm assigning variables for lat/lon:

model2Lat=model2[:]->alat ; Degrees NORTH
model2Lon=model2[:]->alon ; Degrees WEST (?)

And, a print statement to see what the vars look like:
print(model2Lat)

Variable: model2Lat
Type: float
Total Size: 3767796 bytes
             941949 values
Number of Dimensions: 2
Dimensions and sizes: [y | 2349] x [x | 401]
Coordinates:
Number Of Attributes: 3
   long_name : latitude
   units : deg-N
   coordinates : alon alat
(0,0) 18
(0,1) 18
(0,2) 18
(0,3) 18
(0,4) 18
...

print(model2Lon)

Variable: model2Lon
Type: float
Total Size: 3767796 bytes
             941949 values
Number of Dimensions: 2
Dimensions and sizes: [y | 2349] x [x | 401]
Coordinates:
Number Of Attributes: 3
   long_name : longitude
   units : deg-W
   coordinates : alon alat
(0,0) -98
(0,1) -97.95
(0,2) -97.9
(0,3) -97.85
(0,4) -97.8
(0,5) -97.75
(0,6) -97.7
...

OK, units are screwy:
model2Lat@units = "degrees_north"
model2Lon@units = "degrees_east"
print(model2Lat)

Variable: model2Lat
Type: float
Total Size: 3767796 bytes
             941949 values
Number of Dimensions: 2
Dimensions and sizes: [y | 2349] x [x | 401]
Coordinates:
Number Of Attributes: 3
   long_name : latitude
   units : degrees_north
   coordinates : alon alat
(0,0) 18
(0,1) 18
(0,2) 18
(0,3) 18
...
print(model2Lon)

Variable: model2Lon
Type: float
Total Size: 3767796 bytes
             941949 values
Number of Dimensions: 2
Dimensions and sizes: [y | 2349] x [x | 401]
Coordinates:
Number Of Attributes: 3
   long_name : longitude
   units : degrees_east
   coordinates : alon alat
(0,0) -98
(0,1) -97.95
(0,2) -97.9
(0,3) -97.85
...

When I plot, I get a plot that isn't at the resolution I'm expecting,
and I get:
(0) check_for_y_lat_coord: Warning: Data either does not contain a valid
latitude coordinate array or doesn't co
ntain 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' 'd
egreeN' 'degreesN' 'deg north'
(0) check_for_lon_coord: Warning: Data either does not contain a valid
longitude coordinate array or doesn't con
tain 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' 'degree
E' 'degreesE' 'deg east'

As near as I can tell, I've reassigned units in an appropriate manner.
I'm guessing this is is associate with gsn_csm_contour_map, but I've not
found where I've gone wrong.

Thanks,
gerry

On 10/05/2011 06:18 PM, Gerry Creager wrote:
> I'll give it a try in a little while (parental duties right now!)
>
> gc
>
> Adam Phillips wrote:
>> Hi Gerry,
>> The newest version of ncl (v6.0.0) allows array sizes larger than 2GB.
>> Previous versions do not. I am guessing you are using a pre-6.0.0
>> version and your aggregated T array is (potentially) larger than 2GB..
>>
>> If I'm off on my guess let us know.. Adam
>>
>> On 10/05/2011 04:50 PM, Gerald Creager wrote:
>>> Thanks, Gentlemen. I didn't realize that (but I should have)...
>>>
>>> OK, so I tried this:
>>> model2Sst=model2[:]->t ; potential temperature with
>>> coordinates : time, zz, y, x
>>> fatal:Aggregating variable t from file list variable model2 as specified
>>> would exceed maximum NCL variable size
>>>
>>> gerry
>>>
>>> On 10/05/2011 04:45 PM, David Brown wrote:
>>>> Hi Gerry,
>>>> When you use 'addfiles' rather than 'addfile', you get a file list variable rather than a file variable. It requires a different syntax.
>>>> Try
>>>> model2Sst=model2[:]->t
>>>> and
>>>> Time = model2[:]->time
>>>>
>>>> You can say 'model2[0]' in order to "reference a file".
>>>> -dave
>>>>
>>>> On Oct 5, 2011, at 3:09 PM, Gerald Creager wrote:
>>>>
>>>>> I'm trying to assign values to a variable from a netcdf file (ncl_filedump output for the file is attached). A sample snippet of code is attached. Note the error... And, note that line 6 does, in fact load the data into 'model2'.
>>>>>
>>>>> ncl 30> model2Sst=model2->t
>>>>> fatal:(model2) does not reference a filencl 0> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
>>>>> ncl 1> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
>>>>> ncl 2> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
>>>>> ncl 3> load "$NCARG_ROOT/lib/ncarg/nclscripts/contrib/calendar_decode2.ncl"
>>>>> ncl 4>
>>>>> ncl 5> fnameModel2=systemfunc ("ls fcast_20100*.nc") ; files location
>>>>> ncl 6> model2=addfiles(fnameModel2, "r")
>>>>> ncl 7> type="x11"
>>>>> ncl 8> Time=model2->time
>>>>> fatal:(model2) does not reference a file
>>>>> fatal:Execute: Error occurred at or near line 8
>>>>>
>>>>> fatal:Execute: Error occurred at or near line 30
>>>>>
>>>>> Thanks, Gerry
>>>>> -
>>>>> Gerry Creager -- gerry.creager@tamu.edu
>>>>> Texas Mesonet -- AATLT, Texas A&M University
>>>>> Cell: 979.229.5301 Office: 979.458.4020 FAX: 979.862.3983
>>>>> Office: 1700 Research Parkway Ste 160, TAMU, College Station, TX 77843
>>>>> <nclfiledump-model2.txt>_______________________________________________
>>>>> ncl-talk mailing list
>>>>> List instructions, subscriber options, unsubscribe:
>>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>

-- 
Gerry Creager -- gerry.creager@tamu.edu
Texas Mesonet -- AATLT, Texas A&M University
Cell: 979.229.5301 Office: 979.458.4020 FAX: 979.862.3983
Office: 1700 Research Parkway Ste 160, TAMU, College Station, TX 77843
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk