Re: array error

Erik,

Are LAT, LON 2-dimensional by any chance? They must be 1-dimensional
and the same length in order for this kind of "new" statement to work.

--Mary

On Mon, 9 Feb 2009, Dave Allured wrote:

> Strange. Try this to see in detail what's going on. Line 6 will
> probably fail. That's okay, that is what we are trying to trace.
>
> nlats = dimsizes (LAT)
> nlons = dimsizes (LON)
>
> printVarSummary (nt)
> printVarSummary (nlats)
> printVarSummary (nlons)
>
> dims = (/nt,nlats,nlons/)
> print (dims)
>
> --Dave
>
> Erik Noble wrote:
>> Dear Debasish,
>> Thank you. ....just trying to understand, though.
>> You wrote:
>> "Instead of using
>> daily_arr = new((/29,LAT,LON/),float)
>> Please use
>> daily_arr = new((/12,dimsizes(LAT),dimsizes(LON)/),float)
>> ..."
>> I did a similar command already.
>> ncl 22> daily_arr = new((/nt,dimsizes(LAT),dimsizes(LON)/),float)
>> fatal:_NclBuildArray: each element of a literal array must have the same
>> dimension sizes, at least one item doesn't
>> fatal:Execute: Error occurred at or near line 22
>>
>> Shouldn't this work? I'm trying the other suggestions, too.
>> -Erik
>>
>>
>> On Mon, Feb 9, 2009 at 3:50 PM, Debasish Pai Mazumder
>> <debasish_at_gi.alaska.edu <mailto:debasish_at_gi.alaska.edu>> wrote:
>>
>> Dear Erik,
>>
>>
>>
>> Instead of using
>>
>> daily_arr = new((/29,LAT,LON/),float)
>>
>>
>>
>> Please use
>>
>>
>>
>> daily_arr = new((/12,dimsizes(LAT),dimsizes(LON)/),float)
>>
>>
>>
>> Then you can declare daily_arr's dimension as LAT and LON
>>
>>
>>
>> daily_arr!1="LAT"
>>
>> daily_arr!2="LON"
>>
>> daily_arr &LAT=LAT
>>
>> daily_arr &LON=LON
>>
>>
>>
>> Thanks
>>
>>
>>
>> Debasish
>>
>>
>>
>>
>>
>> ------------------------------------------------------------------------
>>
>> *From:* ncl-talk-bounces_at_ucar.edu <mailto:ncl-talk-bounces_at_ucar.edu>
>> [mailto:ncl-talk-bounces_at_ucar.edu
>> <mailto:ncl-talk-bounces_at_ucar.edu>] *On Behalf Of *Erik Noble
>> *Sent:* Monday, February 09, 2009 11:23 AM
>> *To:* ncl-talk_at_ucar.edu <mailto:ncl-talk_at_ucar.edu>
>> *Subject:* [ncl-talk] array error
>>
>>
>>
>> Hi.
>>
>> Could I have help understanding why this error is occurring when
>> creating an array in NCL?
>>
>> When I create the array with numbers, I have no problem:
>>
>> ncl 17> daily_arr = new((/29,15,21/),float)
>>
>> ncl 18>
>>
>>
>>
>> Yet, when I try to use scripting, I get an error:
>>
>> ncl 23> daily_arr = new((/nt,LAT,LAN/),float)
>>
>>
>> fatal:_NclBuildArray: each element of a literal array must have the
>> same dimension sizes, at least one item doesn't
>>
>>
>>
>> Also, where is this error described in the NCL docs?
>>
>> Thank you for your time ahead of time.
>>
>> Sincerely,
>>
>> Erik
>>
>>
>>
>> bash-3.1$ ncl
>>
>> Copyright (C) 1995-2007 - All Rights Reserved
>>
>> University Corporation for Atmospheric Research
>>
>> NCAR Command Language Version 5.0.0
>>
>> The use of this software is governed by a License Agreement.
>>
>> See http://www.ncl.ucar.edu/ for more details.
>>
>> 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> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
>>
>> ncl 3> a =
>> addfile("Daily_Region_2_WRF-GFS_SEP_1st-24-hours-forecasts-00Z.nc","r")
>>
>> ncl 4> lat2d = a->XLAT(0,0,:,:) ; WRF_grid lats
>>
>> ncl 5> lon2d = a->XLONG(0,0,:,:) ; WRF_grid lons
>>
>> ncl 6> rainnc = a->RAINNC ; Total Grid Scale Precipitation
>>
>> ncl 7> rainc = a->RAINC ; Accumulated Total Cumulus Precipitation
>>
>> ncl 8> times=a->Times
>>
>> ncl 9>
>>
>> ncl 10> nt = times(:,0,0) ; should be 29 days
>>
>> ncl 11> ;nt = 29
>>
>> ncl 12> printVarSummary(nt)
>>
>>
>>
>>
>>
>> Variable: nt
>>
>> Type: character
>>
>> Total Size: 29 bytes
>>
>> 29 values
>>
>> Number of Dimensions: 1
>>
>> Dimensions and sizes: [record | 29]
>>
>> Coordinates:
>>
>> Number Of Attributes: 0
>>
>> ncl 13> LAT = lat2d(:,0)
>>
>> ncl 14> LON = lon2d(0,:)
>>
>> ncl 15> printVarSummary(LAT)
>>
>>
>>
>>
>>
>> Variable: LAT
>>
>> Type: float
>>
>> Total Size: 60 bytes
>>
>> 15 values
>>
>> Number of Dimensions: 1
>>
>> Dimensions and sizes: [south_north | 15]
>>
>> Coordinates:
>>
>> Number Of Attributes: 5
>>
>> stagger :
>>
>> units : degree_north
>>
>> description : LATITUDE, SOUTH IS NEGATIVE
>>
>> MemoryOrder : XY
>>
>> FieldType : 104
>>
>> ncl 16> printVarSummary(LON)
>>
>>
>>
>>
>>
>> Variable: LON
>>
>> Type: float
>>
>> Total Size: 84 bytes
>>
>> 21 values
>>
>> Number of Dimensions: 1
>>
>> Dimensions and sizes: [west_east | 21]
>>
>> Coordinates:
>>
>> Number Of Attributes: 5
>>
>> stagger :
>>
>> units : degree_east
>>
>> description : LONGITUDE, WEST IS NEGATIVE
>>
>> MemoryOrder : XY
>>
>> FieldType : 104
>>
>> ncl 17> daily_arr = new((/29,15,21/),float)
>>
>> ncl 18> printVarSummary(daily_arr)
>>
>>
>>
>>
>>
>> Variable: daily_arr
>>
>> Type: float
>>
>> Total Size: 36540 bytes
>>
>> 9135 values
>>
>> Number of Dimensions: 3
>>
>> Dimensions and sizes: [29] x [15] x [21]
>>
>> Coordinates:
>>
>> Number Of Attributes: 1
>>
>> _FillValue : -999
>>
>> ncl 19> delete(daily_arr)
>>
>> ncl 20> daily_arr = new((/29,LAT,LON/),float)
>>
>> fatal:_NclBuildArray: each element of a literal array must have the
>> same dimension sizes, at least one item doesn't
>>
>> fatal:Execute: Error occurred at or near line 20
>>
>>
>>
>> ncl 21> delete(daily_arr)
>>
>> fatal:Variable (daily_arr) is undefined
>>
>> fatal:Execute: Error occurred at or near line 21
>>
>>
>>
>> ncl 22> daily_arr = new((/nt,dimsizes(LAT),dimsizes(LON)/),float)
>>
>> fatal:_NclBuildArray: each element of a literal array must have the
>> same dimension sizes, at least one item doesn't
>>
>> fatal:Execute: Error occurred at or near line 22
>>
>>
>>
>> ncl 23> daily_arr = new((/nt,LAT,LAN/),float)
>>
>>
>> fatal:_NclBuildArray: each element of a literal array must have the
>> same dimension sizes, at least one item doesn't
>>
>> fatal:Execute: Error occurred at or near line 23
>>
>>
>>
>>
>>
>> ------------------------------------------------------------------------
>>
>> _______________________________________________
>> 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 Mon Feb 09 2009 - 15:45:10 MST