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

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
>>>
>>
>
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Oct 6 09:01:42 2011