Re: Error spanning multiple files, addfiles_GetVar

Erik,

Please start reading your error messages more closely before posting to
this list. Your error message is:

     fatal:Subscript out of range, error in subscript #0

If you look at the line where the error message is coming from:

      finarr(ifil*4+2,:,:,:) = addfiles_GetVar(b,ALLFILES(ifil*4+2),"VV")

Subscript #0 refers to the leftmost subscript. You are subscripting
finarr with "ifil*4+2", where "ifil" goes from 0 to (numFILES-1)/4. The
leftmost dimension of finarr is dimensioned (numFILES-1)/4. If you do
the arithmetic, "ifil*4+2" is eventually going to be greater than
(numFILES-1)/4, and hence you will get a subscript out of range error
message. I think maybe you meant to do something like:

      finarr(ifil,:,:,:) = addfiles_GetVar(b,ALLFILES(ifil*4+2),"VV")

However, your do loop will probably need to be fixed, since remember,
you are starting at index 0, and not index 1. You may need to do this:

  do ifil = 0, ((numFILES-1)/4)-1

--Mary

On Feb 5, 2008, at 6:05 PM, Erik Noble wrote:

> Hi.
> May I have some help?
>
> I am trying to span multiple files, return a variable, and put it into
> an
> array.
> I have 245 4xdaily data files that I am trying to read "VV" from ( yet
> only
> at 12 UTC). So this means (245-1)/4 = 61 timesteps tat I am trying to
> span.
>
> I tried making the code below using the concatenate "ListSetType (b,
> "cat"),
> addfiles, and addfiles_GetVar" functions similar to the example at
> http://www.ncl.ucar.edu/Applications/addfiles.shtml
>
> I am not sure why I keep getting the "out of range" error below, even
> though
> I allocated enough space.
>
> Could I have a suggestion here?
> Humbly yours,
> Erik
>
> fatal:Subscript out of range, error in subscript #0
> fatal:Execute: Error occurred at or near line 22 in file
> test_script3.ncl
>
>
> ; Script to Read FNL Reanalysis 4 x Daily data (245 files)
> ; Read in V data for every 12 UTC timestep
> ; Create a new array
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
>
> begin
> ;========================
> ; get list of all met_em files and open as "one big file" for array
> ;========================
> DATADir = "/Volumes/Data_and_Models/NAMA-SOP3/" ; path to your
> data -
> must have trailing "/"
> ALLFILES = systemfunc ("csh -c ' ls " + DATADir +
> "met_em.d01.2006*nc'")
> ; root of data file you want to plot
> numFILES = dimsizes(ALLFILES) ; 245 files present, 4 hourly data
> finarr = new((/(numFILES-1)/4,27,110,139/),"float") ;244/4=61,27
> levels,
> 110x139 grid
>
> do ifil = 0, (numFILES-1)/4
> b = addfiles(ALLFILES(ifil*4+2),"r")
> ListSetType (b, "cat")
> finarr(ifil*4+2,:,:,:) = addfiles_GetVar(b,ALLFILES(ifil*4+2),"VV")
> end do
> printVarSummary(finarr)
> end
>
>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk_at_ucar.edu
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Feb 05 2008 - 21:45:54 MST