<div dir="ltr">Thanks Dave for the assistance. I have made some progress. While I know that I need time to learn ncl subscripting, perhaps I could be helped to make the script work. This script is an adaptation of one of the provided wrf-ncl scripts. The script plots several variables at different levels. I have reproduced below, parts of the script which plot temperature. I have marked in<div>1. yellow, some guidance I got from Dave</div><div>2. blue some changes I made to getting the times</div><div>3. purple, getting the variables, with files_all (I assume it is equivalent to files_all[:] and gets variables from all files)</div><div>4. red I have to specify from which file to plot. This is where I have a problem, how to index files_all.</div><div>I hope I can be assisted.<br><div><br></div><div> files_all = addfiles(FILES+".nc","r")<br><br>; get time information<br> times_in_file = files_all[:]->Times<br><span style="background-color:rgb(255,255,255)"> </span><b style=""><font color="#000000" style=""> <span style="background-color:rgb(255,255,0)">dims = dimsizes (times)<br> ntimes = dimsizes (times(:,0))</span></font></b><br><br>; strip out the day and hour<br> times_to_use = new(dims(0),string)<br><br> do i=0,dims(0)-1<br> times_to_use(i) = chartostring(times_in_file(i,8:12))<br> end do<br>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<br> do it = 0,<font color="#0000ff">ntimes(0)</font>-1,2 ; TIME LOOP<br><br> print("Working on time: " +<font color="#0000ff"> times(it,0)</font> )<br> res@TimeLabel =<font color="#0000ff"> times(it,0) </font> ; Set Valid time to use on plots<br>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<br>; First get the variables we will need<br> tc = wrf_user_getvar<font color="#9900ff">(<b style="background-color:rgb(255,242,204)">files_all</b></font><b style="background-color:rgb(255,242,204)">,</b>"tc",it) ; T in C<br> u = wrf_user_getvar(files_all,"ua",it) ; u averaged to mass points<br> v = wrf_user_getvar(files_all,"va",it) ; v averaged to mass points<br> p = wrf_user_getvar(files_all, "pressure",it) ; pressure is our vertical coordinate<br> z = wrf_user_getvar(files_all, "z",it) ; grid point height<br> rh = wrf_user_getvar(files_all,"rh",it) ; relative humidity<br>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<br> do level = 0,nlevels-1 ; LOOP OVER LEVELS<br><br> pressure = pressure_levels(level)<br> tc_plane = wrf_user_intrp3d(tc,p,"h",pressure,0.,False)<br> z_plane = wrf_user_intrp3d( z,p,"h",pressure,0.,False)<br> rh_plane = wrf_user_intrp3d(rh,p,"h",pressure,0.,False)<br> u_plane = wrf_user_intrp3d( u,p,"h",pressure,0.,False)<br> v_plane = wrf_user_intrp3d( v,p,"h",pressure,0.,False)<br><br> spd = (u_plane*u_plane + v_plane*v_plane)^(0.5) ; m/sec<br> spd@description = "Wind Speed"<br> spd@units = "m/s"<br> u_plane = u_plane*1.94386 ; kts<br> v_plane = v_plane*1.94386 ; kts<br> u_plane@units = "kts"<br> v_plane@units = "kts"<br><br> wks = gsn_open_wks(type,"Plots00-" + "P" + pressure_levels(level) + "_" + times_to_use(it) + "Z")<br><br> ; Plotting options for T<br> opts = res<br> opts@cnLineColor = "Red"<br> opts@ContourParameters = (/ 5.0 /)<br> opts@cnInfoLabelOrthogonalPosF = 0.07 ; offset second label information<br> opts@gsnContourLineThicknessesScale = 2.0<br> contour_tc = wrf_contour(<b style="background-color:rgb(255,255,255)"><font color="#ff0000">files_all</font></b>,wks,tc_plane,opts)<br> delete(opts)<br><br> end do ; END OF LEVEL LOOP<br>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<br> end do ; END OF TIME LOOP<br></div></div></div><div id="DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2"><br>
<table style="border-top:1px solid #d3d4de">
<tr>
<td style="width:55px;padding-top:13px"><a href="https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail" target="_blank"><img src="https://ipmcdn.avast.com/images/icons/icon-envelope-tick-round-orange-animated-no-repeat-v1.gif" alt="" width="46" height="29" style="width: 46px; height: 29px;"></a></td>
<td style="width:470px;padding-top:12px;color:#41424e;font-size:13px;font-family:Arial,Helvetica,sans-serif;line-height:18px">Virus-free. <a href="https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail" target="_blank" style="color:#4453ea">www.avast.com</a>
</td>
</tr>
</table><a href="#DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2" width="1" height="1"></a></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Apr 2, 2022 at 5:06 PM Dave Allured - NOAA Affiliate <<a href="mailto:dave.allured@noaa.gov">dave.allured@noaa.gov</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 dir="ltr"><div>Here are two different ways to get the size of the time dimension. The second method uses dimension reduction to find the size of only one slice of the 2-D character array.</div><div><br></div><div>dims = dimsizes (times)</div><div>ntimes = dims(0)</div><div><br></div><div>ntimes = dimsizes (times(:,0))</div><div dir="ltr"><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Apr 1, 2022 at 11:04 PM Zilore Mumba <<a href="mailto:zmumba@gmail.com" target="_blank">zmumba@gmail.com</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 dir="ltr">I have 13 files, i.e. analysis plus 12 hourly files of forecasts.
print (files_all[0]->Times) and printVarSummary(times) give the following output:<div>Variable: Times (file variable)<br>Type: character<br>Total Size: 19 bytes<br> 19 values<br>Number of Dimensions: 2<br>Dimensions and sizes: [Time | 1] x [DateStrLen | 19]<br>Coordinates: <br>(0,0) 2<br>(0,1) 0<br>(0,2) 2<br>(0,3) 2<br>(0,4) -<br>(0,5) 0<br>(0,6) 3<br>(0,7) -<br>(0,8) 2<br>(0,9) 8<br>(0,10) _<br>(0,11) 0<br>(0,12) 0<br>(0,13) :<br>(0,14) 0<br>(0,15) 0<br>(0,16) :<br>(0,17) 0<br>(0,18) 0<br><br>Variable: times<br>Type: character<br>Total Size: 247 bytes<br> 247 values<br>Number of Dimensions: 2<br>Dimensions and sizes: [Time | 13] x [DateStrLen | 19]<br>Coordinates:</div><div> <br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Apr 1, 2022 at 11:53 PM Dave Allured - NOAA Affiliate <<a href="mailto:dave.allured@noaa.gov" target="_blank">dave.allured@noaa.gov</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 dir="ltr"><div>Please show print (files_all[0]->Times) and printVarSummary(times).</div><div dir="ltr"><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Apr 1, 2022 at 3:42 PM Zilore Mumba via ncl-talk <<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.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 dir="ltr">I have a number of WRF output files, one file per hour, which I open with<div>files_all = addfiles(FILES+".nc","r")<br></div><div><br></div><div>I then get the times and the total number with.</div><div> times = files_all[:]->Times ; get all times in the file<br> ntimes = dimsizes(times) ; number of times in the file<br></div><div><br></div><div>I think I am missing something here because with print(ntimes) I see two numbers.</div><div>I want to use ntimes in a loop: do it = 0,ntimes-1,2</div><div><br></div><div>This gives me the error: fatal:Loop end must be scalar, can't execute loop.</div><div>How do I get the times from files opened with addfiles. I will appreciate any assistance.</div></div></blockquote></div></div>
</blockquote></div>
</blockquote></div></div>
</blockquote></div>