<div dir="ltr"><div>HI Zilore,</div><div><br></div><div>I think the issue is that addfiles() return a *list* of filevars, where addfile returns a single filevar.  Lists are indexed with [ ] notation; perhaps something like this will work:</div><div><br></div><div>  x = a[:]->T2(0,{850},{-15},{30})</div><div><br></div><div>Rick<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sun, Apr 19, 2020 at 1:54 PM zilore mumba via ncl-talk <<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><div style="font-family:Helvetica Neue,Helvetica,Arial,sans-serif;font-size:13px"><div dir="ltr">I am trying to get time series at given locations from several WRF output files, e.g. 48H forecast temperature at 850hPa. Each file contains a day's model run.</div><div dir="ltr"> I have used ncl for years by adapting existing scripts, without really understanding. One issue I have difficulty with is getting variables from a file, especially WRF output.</div><div dir="ltr">I thought the script to get a time series at 15S, 30E would be as simple as the one below. I am able to print the varsummary, but printing x tells me that "a is not a valid file reference.</div><div dir="ltr">I would appreciate ant assistance.</div><div dir="ltr"><div>begin<br>files = systemfunc("ls -1 data/wrfout_d01_2020*") + ".nc"<br>  a = addfiles(files,"r")<br><br>  ListSetType (a, "cat")          ; concatenate (=default)<br>  tc = wrf_user_getvar(a,"tc",0)<br>  printVarSummary (tc)<br><br>  x = a->T2(0,{850},{-15},{30})<br>  asciiwrite("outfile.txt",x)<br>end<br><br></div><div dir="ltr">Zilore Mumba<br></div></div></div></div>_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a></blockquote></div>