[ncl-talk] unable to open multiple files
Mary Haley
haley at ucar.edu
Tue Nov 24 10:22:19 MST 2015
Geeta,
You need to do some debugging on your own to figure out these problems,
because it's getting frustrating trying to help.
The "remove primitive" type of warnings are coming from a part of the code
that you didn't show in your email.
The "fatal:ListSelect: Index out of range" error is telling you exactly
what the problem is, and that is you are trying to subscript a list
variable with a value that is out-of-range.
You said it occurred at line 40, which is:
40 wrfs = wrfa[n]
Did you check that "wrfa" indeed has at least "n-1" elements? You can use
"ListCount" to check this:
print(ListCount(wrfa))
You are looping from 0 to 4, and I know this is because you said you have 5
files, but did you make sure that trmmf and wrff actually have 5 file names
in them?
If you continue to have problems with this, then *please* try to do some
debugging on your own. You can do this with print or printVarSummary
statements.
Also, in general, it's not a good idea to "hard-code" a do loop with
something like:
do n = 0,4
It is better to do something like:
nfiles = dimsizes(wrff) ; get the number of files
do n=0,nfiles-1
But, you also need to make sure that wrff and trmmf actually have the same
number of files, since you are using the same loop counter for both of them.
If you still continue to have problems after debugging on your own, then
include the *full* script, and your data files if possible. You can use our
anonymous ftp, which I think you already have the information for.
--Mary
On Tue, Nov 24, 2015 at 8:39 AM, Geeta Geeta <geetag54 at yahoo.com> wrote:
> thanks Dennis for your reply.
> I used ls only in the systemfunc .
>
> My data is of different years. hence I wanted to see if the files are
> being in the right order in which I wanted (2011, 2013,2014 and 2015).
> Hence I wanted the date of each file too.
>
> I have modified the script slightly to include that. It is like this.
> ntrmm = dimsizes(trmm)
> ; print(ntrmm)
>
> ; do i = 0, ntr mm-1,1
> ; print(i)
> trmm0 = trmm[0]->r(0,:,:)
> trmmdate0 = trmm[0]->time
> print(trmmdate0)
>
>
> the output looks like this.
> ariable: trmm0
> Type: float
> Total Size: 2304000 bytes
> 576000 values
> Number of Dimensions: 2
> Dimensions and sizes: [latitude | 400] x [longitude | 1440]
> Coordinates:
> latitude: [-49.875..49.875]
> longitude: [0.125..359.875]
> Number Of Attributes: 8
> time : 0
> long_name : daily (0Z-21Z) rainfall total
> units : mm
> grid_name : grid-1
> grid_type : linear
> level_description : Earth surface
> _FillValue : -9999.9
> coordinates : time latitude longitude
>
>
> Variable: trmmdate0
> Type: double
> Total Size: 8 bytes
> 1 values
> Number of Dimensions: 1
> Dimensions and sizes: [time | 1]
> Coordinates:
> time: [ 0.. 0]
> Number Of Attributes: 1
> units : hours since 2011-04-22 00:00:00Z
> (0) 0
>
>
> I have 5 of TRMM files. I tried read the trmm using the do loop, But
> unsuccessful. this is what I have done now for this
>
> 18 trmmf = systemfunc("ls 3B42_daily.20*.nc") ; TRMM RF 22 May
> 2013
> 19 trmma = addfiles(trmmf,"r") ; TRMM RF 22 May 2013
> 20 trmms = trmma[0]
> 21
> 22 wrff = systemfunc("ls wrfout_d03_20*_00:00:00") ; model
> rainfall
> 23 wrfa = addfiles(wrff,"r")
> 24 wrfs = wrfa[0]
> 25
> 26 ; a = addfile("3B42_daily.2013.05.22.7.nc","r") ; TRMM RF 22
> May 2013
> 27
> 28 ; f = addfile("wrfout_d03_2013-05-21_00:00:00","r") ; model
> rainfall
> 29
> 30 ; Plotting TRMM rainfall
> 31 var = 0
> 32
> 33 plot = new(10,graphic)
> 34 wks = gsn_open_wks("ps","trmm_model")
> 35 gsn_define_colormap(wks,"prcp_2")
> 36
> 37 ; print(a)
> 38 do n = 0, 4
> 39 trmms= trmma[n]
> 40 wrfs = wrfa[n]
> 41
> 42 rf = trmms->r(0,:,:)
>
>
> rainc28 = wrfs->RAINC(28,:,:)
> 57 rainnc28 = wrfs->RAINNC(28,:,:)
> 58 rainc4 = wrfs->RAINC(4,:,:)
> 59 rainnc4 = wrfs->RAINNC(4,:,:)
> 60 ;
> 61 latitude = wrfs->XLAT(0,:,0)
> 62 longitude = wrfs->XLONG(0,0,:)
> 63 ; times_in_file = wrfs->Times
> 64 ; init = times_in_file(:,0)
> 65 ; print(init)
>
>
>
> I get an error at the LINE 40.
> (0) latitude
> (1) longitude
> warning:_NclIRemovePrimitive: First parameter is a missing value,
> returning missing values
> warning:_NclIRemovePrimitive: First parameter is a missing value,
> returning missing values
> (0) latitude
> (1) longitude
> warning:_NclIRemovePrimitive: First parameter is a missing value,
> returning missing values
> warning:_NclIRemovePrimitive: First parameter is a missing value,
> returning missing values
> (0) latitude
> (1) longitude
> warning:_NclIRemovePrimitive: First parameter is a missing value,
> returning missing values
> warning:_NclIRemovePrimitive: First parameter is a missing value,
> returning missing values
> fatal:ListSelect: Index out of range
> fatal:["Execute.c":8578]:Execute: Error occurred at or near line 40 in
> file trmm-model-
>
>
> Dont know what's going wrong.
>
>
> Geeta.
>
>
>
> On Tuesday, 24 November 2015 12:31 AM, Dennis Shea <shea at ucar.edu> wrote:
>
>
> Please do not start a new thread. Use the 'reply all' to keep all
> emails together
>
> ====
>
> [1]
> Note what I suggested
>
> trmm_files = systemfunc("ls 3B42_daily.20*.nc")
>
> I did *not* use
>
> trmm_files = systemfunc("ls -l 3B42_daily.20*.nc")
>
> the '-l' returns much more information than just the file names.
>
> ===
>
> [2]
> Had you informed ncl-talk that you wanted data for 2011 and the file names
> are
>
> 3B42_daily.2011.04.22.7.nc 3B42_daily.2014.06.05.7.nc
> 3B42_daily.2013.05.22.7.nc 3B42_daily.2015.04.23.7.nc
> 3B42_daily.2014.06.02.7.nc
>
> Then
>
> trmm_files = systemfunc("ls 3B42_daily.2011*.nc")
>
> would have been suggested. Note: no '-l'
>
> [3]
> I did *not* use
> nfiles = ListCount(trmm)
> but this should return the same as
> nfiles = dimsizes(trmm_files)
>
> [4]
>
> re: "In the printVarSummary, I dont get the time info of each file...."
>
> trmm_rf2011 = trmm[0]->r
> printVarSummary(trmm_rf2011)
> trmm_rf2013 = trmm[1]->r
> printVarSummary(trmm_rf2013)
> trmm_rf2014a = trmm[2]->r
> printVarSummary(trmm_rf2014a)
> trmm_rf2014b = trmm[3]->r
> printVarSummary(trmm_rf2014b)
> trmm_rf2015 = trmm[4]->r
> printVarSummary(trmm_rf2015)
>
>
> Since you did not include any output from 'printVarSummary', how is
> ncl-talk supposed to answer the question.
>
> On Sun, Nov 22, 2015 at 10:46 AM, Geeta Geeta <geetag54 at yahoo.com> wrote:
> > thanks Dennis.
> > I want to know if "join" will add the data???.
> > Now I have kept the default option as "cat" only.
> > But I am facing small problem.
> > when I use "ls -l" below
> >
> > trmm_files = systemfunc("ls -l 3B42_daily.20*.nc") ; with ls -l
> > printVarSummary is not shown.
> > printVarSummary(trmm_files)
> >
> > I dont get the output.
> >
> > Following are the data files.
> > aditya at agniilap:~/geeta/ncl/TS_rain$ ls 3B42*
> > 3B42_daily.2011.04.22.7.nc 3B42_daily.2014.06.05.7.nc
> > 3B42_daily.2013.05.22.7.nc 3B42_daily.2015.04.23.7.nc
> > 3B42_daily.2014.06.02.7.nc
> >
> >
> > Can I also get the time info for each file printed for me so that I know
> > that the variable trmm_rf2011 has ONLY 2011 data???
> >
> > trmm_files = systemfunc("ls 3B42_daily.20*.nc")
> > date_files = systemfunc("ncdump -h 3B42_daily.20*.nc")
> > printVarSummary(trmm_files)
> > printVarSummary(date_files)
> >
> > trmm = addfiles(trmm_files,"r")
> > printVarSummary(trmm)
> >
> > nfiles = ListCount(trmm)
> > print(nfiles)
> > ListSetType(trmm,"cat")
> >
> > trmm_rf2011 = trmm[0]->r
> > printVarSummary(trmm_rf2011)
> > trmm_rf2013 = trmm[1]->r
> > printVarSummary(trmm_rf2013)
> > trmm_rf2014a = trmm[2]->r
> > printVarSummary(trmm_rf2014a)
> > trmm_rf2014b = trmm[3]->r
> > printVarSummary(trmm_rf2014b)
> > trmm_rf2015 = trmm[4]->r
> > printVarSummary(trmm_rf2015)
> >
> >
> > In the printVarSummary, I dont get the time info of each file....
> >
> >
> > Geeta.
>
> >
> > _______________________________________________
> > ncl-talk mailing list
> > ncl-talk at ucar.edu
> > List instructions, subscriber options, unsubscribe:
> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
> >
>
>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20151124/e5b5bab6/attachment.html
More information about the ncl-talk
mailing list