<html><head></head><body><div class="ydpef9c9addyahoo-style-wrap" style="font-family:Helvetica Neue, Helvetica, Arial, sans-serif;font-size:13px;"><div></div>
<div dir="ltr" data-setdir="false">Sorry Rick I put in a surface variable T2 which only 3 coordinates, so I put <span>"x = a[:]->T2(0,{-15},{30})"</span>. Now I am getting error of coordinates west-east and south-north <span>not having associated coordinate variables. I have seen this error before. I will find out how to get out of it </span><br></div><div><br></div>
</div><div id="yahoo_quoted_7471934217" class="yahoo_quoted">
<div style="font-family:'Helvetica Neue', Helvetica, Arial, sans-serif;font-size:13px;color:#26282a;">
<div>
On Sunday, April 19, 2020, 11:53:33 PM GMT+2, Rick Brownrigg <brownrig@ucar.edu> wrote:
</div>
<div><br></div>
<div><br></div>
<div><div id="yiv0600931807"><div><div dir="ltr"><div>Perhaps do a printVarSummary(a[0]->T2) ; it sounds like its not the shape you think it might/should be; i.e. 3D vs 4D</div><div><br clear="none"></div><div><br clear="none"></div></div><br clear="none"><div class="yiv0600931807yqt0411640310" id="yiv0600931807yqt33899"><div class="yiv0600931807gmail_quote"><div class="yiv0600931807gmail_attr" dir="ltr">On Sun, Apr 19, 2020 at 3:29 PM zilore mumba <<a rel="nofollow" shape="rect" ymailto="mailto:zmumba@yahoo.com" target="_blank" href="mailto:zmumba@yahoo.com">zmumba@yahoo.com</a>> wrote:<br clear="none"></div><blockquote class="yiv0600931807gmail_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></div>
<div><br clear="none"></div><div><br clear="none"></div>
</div><div id="yiv0600931807gmail-m_4060975555631802753yahoo_quoted_7983023808">
<div style="font-family:Helvetica, Arial, sans-serif;font-size:13px;color:rgb(38,40,42);">
<div>
On Sunday, April 19, 2020, 11:27:47 PM GMT+2, zilore mumba <<a rel="nofollow" shape="rect" ymailto="mailto:zmumba@yahoo.com" target="_blank" href="mailto:zmumba@yahoo.com">zmumba@yahoo.com</a>> wrote:
</div>
<div><br clear="none"></div>
<div><br clear="none"></div>
<div><div id="yiv0600931807gmail-m_4060975555631802753yiv7781110550"><div><div style="font-family:Helvetica Neue, Helvetica, Arial, sans-serif;font-size:13px;"><div></div>
<div dir="ltr">Thanks Rick, that is close. The error am getting now is:</div><div dir="ltr"> <span>fatal:Number of subscripts on rhs do not match number of dimensions of aggregated cat type variable, (4) Subscripts used, (3) Subscripts expected</span><br clear="none"></div><div><br clear="none"></div>
</div><div id="yiv0600931807gmail-m_4060975555631802753yiv7781110550yqt05792"><div id="yiv0600931807gmail-m_4060975555631802753yiv7781110550yahoo_quoted_8222480418">
<div style="font-family:Helvetica, Arial, sans-serif;font-size:13px;color:rgb(38,40,42);">
<div>
On Sunday, April 19, 2020, 10:17:12 PM GMT+2, Rick Brownrigg <<a rel="nofollow" shape="rect" ymailto="mailto:brownrig@ucar.edu" target="_blank" href="mailto:brownrig@ucar.edu">brownrig@ucar.edu</a>> wrote:
</div>
<div><br clear="none"></div>
<div><br clear="none"></div>
<div><div id="yiv0600931807gmail-m_4060975555631802753yiv7781110550"><div><div dir="ltr"><div>HI Zilore,</div><div><br clear="none"></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 clear="none"></div><div> x = a[:]->T2(0,{850},{-15},{30})</div><div><br clear="none"></div><div>Rick<br clear="none"></div></div><br clear="none"><div><div id="yiv0600931807gmail-m_4060975555631802753yiv7781110550yqt67703"><div dir="ltr">On Sun, Apr 19, 2020 at 1:54 PM zilore mumba via ncl-talk <<a rel="nofollow" shape="rect" ymailto="mailto:ncl-talk@ucar.edu" target="_blank" href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>> wrote:<br clear="none"></div><blockquote 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 clear="none">files = systemfunc("ls -1 data/wrfout_d01_2020*") + ".nc"<br clear="none"> a = addfiles(files,"r")<br clear="none"><br clear="none"> ListSetType (a, "cat") ; concatenate (=default)<br clear="none"> tc = wrf_user_getvar(a,"tc",0)<br clear="none"> printVarSummary (tc)<br clear="none"><br clear="none"> x = a->T2(0,{850},{-15},{30})<br clear="none"> asciiwrite("outfile.txt",x)<br clear="none">end<br clear="none"><br clear="none"></div><div dir="ltr">Zilore Mumba<br clear="none"></div></div></div></div>_______________________________________________<br clear="none">
ncl-talk mailing list<br clear="none">
<a rel="nofollow" shape="rect" ymailto="mailto:ncl-talk@ucar.edu" target="_blank" href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br clear="none">
List instructions, subscriber options, unsubscribe:<br clear="none">
<a rel="nofollow" shape="rect" target="_blank" href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a></blockquote></div></div></div></div></div>
</div>
</div></div></div></div></div>
</div>
</div></div></blockquote></div></div></div></div></div>
</div>
</div></body></html>