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

Hi Gerry,
When you have 2D latitudes/longitudes, you need to assign the 2D lat/lon
arrays as lat2d/lon2d attributes of your data array.

I tried this with the data file you sent in:

ncl 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> a = addfile("fcast_201009.nc","r")
ncl 3> t = a->t(:,0,:,:)
ncl 4> t@lat2d = a->alat
ncl 5> t@lon2d = a->alon
ncl 6> wks = gsn_open_wks("x11","test")
ncl 7> plot = gsn_csm_contour_map(wks,t(0,:,:),False)

and it worked. Let us know if you continue to have a problem plotting
the data.. Adam

On 10/06/2011 09:10 AM, Gerald Creager wrote:
> I just tried but passive mode was refused (I've had this problem
> recently here, too!). Donna or I will get a subset of one of the model
> files uploaded, as well as the script shortly.
>
> Thanks, Gerry
>
>
> On 10/06/2011 10:01 AM, Dennis Shea wrote:
>> GC .. Send 2 or 3 sample files and a minimalistic script that isolates
>> the issues you are having.
>>
>> ftp ftp.cgd.ucar.edu
>> anonymous
>> Use your email address for the password
>> cd incoming
>> put ...
>> put ...
>> :
>> quit
>>
>> Send email to ncl-talk with the following information:
>>
>> The names of the files after you have successfully completed the transfer.
>> We cannot look at the ftp/incoming directory for security reasons.
>> The type of system you are on (uname -a).
>>
>> THX
>> On 10/6/11 8:33 AM, Gerald Creager wrote:
>>> 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
>>>>>
>>>>
>>>
>

-- 
______________________________________________________________
Adam Phillips                                asphilli@ucar.edu
NCAR/Climate and Global Dynamics Division       (303) 497-1726
P.O. Box 3000				
Boulder, CO 80307-3000    http://www.cgd.ucar.edu/cas/asphilli
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk