<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"><div><div>Hello once again.</div><div>I am hoping someone may have a script which they can share,or correct my code below where I am attempting to extract a time series. From the assistance I got a few days ago, I am able to write values at a given location at the surface T2) and a a fixed level (V). Though I am getting several warnings, firstly </div><div><br></div><div dir="ltr" data-setdir="false"><div><div>Dimension (west_east) and (south_north) of file (wrfout_d01_2020-02-03_00:00:00) does not have an associated coordinate variable</div><div dir="ltr" data-setdir="false">and </div><div dir="ltr" data-setdir="false"><br></div><div dir="ltr" data-setdir="false"><span>warning:Invalid stride: stride must be positive non-zero integer</span><br></div></div><div dir="ltr" data-setdir="false"><span><br></span></div><div dir="ltr" data-setdir="false"><span>But my problem is that I am not able to select the forecast time of the series, i.e e.g. daily time series of the forecast at 0600UTC. Currently I am printing 484 values from 12 files each with forecast range of 120H (5 days). I am not able to tell which value corresponds to what forecast range and why they are 484.</span></div><div dir="ltr" data-setdir="false"><span>Is there a way to select for x and y (in the code below) the specific time step for the time series?</span></div><div dir="ltr" data-setdir="false"><span><br></span></div><div dir="ltr" data-setdir="false"><span>Help will be appreciated</span></div><br></div><div dir="ltr" data-setdir="false">begin</div><div> files = systemfunc("ls -1 data/wrfout_d01_2020*") + ".nc"</div><div> a = addfiles(files,"r")</div><div><br></div><div> ListSetType (a, "cat") ; concatenate (=default)</div><div><br></div><div> varnames = getfilevarnames(a[0]) ; get var names from the first file</div><div> print(varnames)</div><div><br></div><div> x = a[:]->T2(:,{-18},{30})</div><div> y = a[:]->V(:,{2},{-18},{30})</div><div><br></div><div> asciiwrite("Lusaka.txt", x)</div><div> asciiwrite("Lusaka2.txt", y)</div><div><br></div><div>end</div><div><br></div></div><br></div><div><br></div></div></body></html>