<div dir="ltr">Hi:<div><br></div><div>Here is some sample code.</div><div><br></div><div>external GET_IJ "./get_ij.so"   ; convert get_ij.f with WRAPIT get_ij.f <br></div><div>begin</div><div>  xlat2d=wrf_user_getvar(a,"XLAT",0)<br>  xlon2d=wrf_user_getvar(a,"XLONG",0)<br>  dims2d=dimsizes(xlat2d)<br>  sta_i  = new(n_sta,integer)   ; Create array to hold text objects<br>  sta_j  = new(n_sta,integer)   ; Create array to hold text objects<br>  do i_sta=0,n_sta-1<br>  i_loc = 1<br>  j_loc = 1<br>  GET_IJ::get_ij(xlat2d,xlon2d,Latitude(i_sta),Longitude(i_sta),i_loc,j_loc,dims2d(0),dims2d(1))<br>  sta_i(i_sta) = i_loc-1<br>  sta_j(i_sta) = j_loc-1<br>  end do<br></div><div><br></div><div>   t2=wrf_user_getvar(a,"T2",it)<br></div><div>   q2=wrf_user_getvar(a,"Q2",it)<br>   rh   = relhum(t2, q2, psfc)<br>   u10 = wrf_user_getvar(a,"U10",it)<br>   v10 = wrf_user_getvar(a,"V10",it)<br>   rainnc = wrf_user_getvar(a,"RAINNC",it)<br>   u_plane := u10*3.6<br>   v_plane := v10*3.6<br>   wspd = wind_speed(u_plane,v_plane)<br>   wdir = wind_direction(u_plane,v_plane,0)<br></div><div><br></div><div>   do i_sta = 0,n_sta-1<br>     sta_t2(i_sta,time)=t2(sta_j(i_sta),sta_i(i_sta))-273.15<br>     sta_rh(i_sta,time)=rh(sta_j(i_sta),sta_i(i_sta))<br>     sta_speed(i_sta,time)=wspd(sta_j(i_sta),sta_i(i_sta))<br>     sta_wdir(i_sta,time)=wdir(sta_j(i_sta),sta_i(i_sta))<br>     sta_prec(i_sta,time)=rainnc(sta_j(i_sta),sta_i(i_sta))<br>   end do<br>end</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Oct 23, 2021 at 8:38 PM Barry Lynn <<a href="mailto:barry.h.lynn@gmail.com">barry.h.lynn@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">Hi:<div><br></div><div><a href="https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_getvar.shtml" target="_blank">https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_getvar.shtml</a><br></div><div><br></div><div>Did you get this far?</div><div><br></div><div>If  so, I can send you a program to find the location of points with specific latitude and longitude.</div><div><br></div><div>Barry</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Oct 23, 2021 at 7:50 PM Setareh Rahimi 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>Dear NCL users,</div><div>I have a WRF output file in nc format. I need to have point values of wind, temperature,and humidity  at specific lat/lon and specific times. I do not need a map plot or time series. I just want point values for each variable at specific times ( 00,03,06,09,12,15,18,21 UTC) in ascii format .For example, for wind variable at lat:y and lon:x</div><div>Time:             00   ,  03   , 06   , 09   ,12</div><div>wind value :   x1   ,  x2   , x3   , x4   , x5, <br></div><div>Would  you please kindly help me in this regard?</div><div>Best wishes,<br></div><div><div><div><div><br>-- <br><div dir="ltr">S.Rahimi<br><br></div></div></div></div></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></blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr"><div dir="ltr"><div><div dir="ltr"><br></div><div dir="ltr">Barry H. Lynn, Ph.D<div><div>Senior <font color="#000000">Scientist, Lecturer,</font></div><div><div><font color="#000000">The Institute of Earth Sciences, <br>The Hebrew University of Jerusalem, <br>Givat Ram, Jerusalem 91904, Israel <br></font></div><font color="#000000">Tel: 972 547 231 170<br>Fax: (972)-25662581</font></div></div><div><span style="color:rgb(136,136,136)"><br></span></div><div>Weather It Is, LTD<br>Weather and Climate Focus<br><a href="http://weather-it-is.com" target="_blank">https://weather-it-is.com</a><br>Jerusalem, Israel<br>Local: 02 930 9525<br>Cell: 054 7 231 170<br>Int-IS: x972 2 930 9525<br><br></div></div></div></div></div>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><br></div><div dir="ltr">Barry H. Lynn, Ph.D<div><div>Senior <font color="#000000">Scientist, Lecturer,</font></div><div><div><font color="#000000">The Institute of Earth Sciences, <br>The Hebrew University of Jerusalem, <br>Givat Ram, Jerusalem 91904, Israel <br></font></div><font color="#000000">Tel: 972 547 231 170<br>Fax: (972)-25662581</font></div></div><div><span style="color:rgb(136,136,136)"><br></span></div><div>Weather It Is, LTD<br>Weather and Climate Focus<br><a href="http://weather-it-is.com" target="_blank">https://weather-it-is.com</a><br>Jerusalem, Israel<br>Local: 02 930 9525<br>Cell: 054 7 231 170<br>Int-IS: x972 2 930 9525<br><br></div></div></div></div></div>