<div dir="ltr"><div>Dear Dennis,</div><div><br></div><div>Thanks for your advice. <br></div><div>I prepare a script as follow:</div><div>----------------------------------------------------------------------------------------------------------------------------------------------------------------</div><div>lon_pts  = 51<br> lat_pts  = 31<br><br>   f     = addfile ("wrfout_d03_2021-01-14_00:00:00", "r")<br>   lat2d  = f->XLAT(0,:,:)                  ; size = (nlat,nlon)<br>   lon2d = f->XLONG(0,:,:)              ; size = (nlat,nlon)<br>   time  = wrf_user_getvar(f,"times",-1)<br><br>;=========================================<br><br>   u10 = f->U10<br>   v10 = f->V10<br>   ;u_plane := u10*3.6<br>   ;v_plane := v10*3.6<br>   wspd = wind_speed(u10,v10)<br>   wdir = wind_direction(u10,v10,0)<br><br>;===========================================<br>  printVarSummary(time)<br>  printVarSummary(rh)<br>   x1     = f->T2<br>   x=x1-273.15<br><br>  nm = getind_latlon2d (lat2d,lon2d, lat_pts, lon_pts)<br><br>  print(nm)<br><br>  do k=0,dimsizes(lat_pts)-1    ; loop over locations<br>     n = nm(k,0)<br>     m = nm(k,1)<br><br>     x_nm =  x(:,n,m)     ;  x(n,m)<br><br>     wspd_nm = wspd (:,n,m)<br>     print("---")<br>     print(lat2d(n,m)+"   "+lon2d(n,m))<br>     print(time)<br>     print(x_nm)<br>     print(wspd_nm)<br><br>     asciiwrite ("temperature", x_nm)<br>     asciiwrite ("wind", wspd_nm)<br>    asciiwrite ("time", time)<br><br><br>  end do<br></div><div>----------------------------------------------------------------------------------------------------------------------------------------------------------------</div><div>Now my questions are that:</div><div>1. Is my script correct?</div><div>2. How can I have time and one of variable (say temperature) net together in an ascii file please?</div><div><br></div><div>like:</div><div><br></div><div>time     temp<br></div><div>00         2<br></div><div>06         2.3<br></div><div>12         4</div><div>Look forward to hearing from you.</div><div>Best wishes,<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Jul 15, 2023 at 6:22 AM Dennis Shea <<a href="mailto:shea@ucar.edu">shea@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"><div><a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/rcm2points.shtml" target="_blank">rcm2points</a></div><div>or to include meta data<br></div><div><a href="https://www.ncl.ucar.edu/Document/Functions/Contributed/rcm2points_Wrap.shtml" target="_blank"><b>rcm2points_Wrap</b></a><br></div><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Jul 14, 2023 at 8:31 PM zoe jacobs 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"><div>Hi dear NCL users,</div><div>I need to extract values for some specific lat/lon from my WRF output files ( for example temperature value fat lat=23 and lon= 36 at 00:00  UTC) . But I am completely new in NCL, and do not know how I can d o this. May I ask you please to help me to achieve this?</div><div>Looking forward to hearing from you.</div><div>Kind regards.<br></div></div>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="https://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">https://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
</blockquote></div>
</blockquote></div>