<html><head></head><body><div class="yahoo-style-wrap" style="font-family:Helvetica Neue, Helvetica, Arial, sans-serif;font-size:13px;"><div dir="ltr" data-setdir="false">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" data-setdir="false"> 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" data-setdir="false">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" data-setdir="false">I would appreciate ant assistance.</div><div dir="ltr" data-setdir="false"><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" data-setdir="false">Zilore Mumba<br></div></div></div></body></html>