<div dir="ltr"><div>Oops, it was pointed out that I did not attach the script.</div><div>Attached<br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Dec 31, 2020 at 8:33 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>The attached prints the values of grid points surrounding the 940.504 values.</div><div><br></div><div>(0) ==================<br>(0) n=2 m=2 LAT2D(n,m)=31.9429 LON2D(n,m)=51.921 EMFLX(n,m)=940.504<br>(0) ==================<br>(0) 116.882 267.648 93.284 13.129 0.019<br>(0) 406.997 579.607 570.014 403.070 34.221<br>(0) 0.000 724.817 <b>940.504</b> 1244.555 415.407<br>(0) 0.000 0.000 485.039 1189.545 1862.362<br>(0) 0.000 0.000 551.603 889.335 2226.263</div><div><br></div><div><br></div><div>as you can readily see, the flux is not 'smoothly varying' due to 0's. The 0s are treated a real values. <br></div><div>What should the plotted value be? The raster and cell use multiple values for the plotted value.</div><div>If you want, try the following<br></div><div><br></div><div><br></div><div>The variable(s) have no _FillValue attribute, <br></div><div><br></div><div> emflx = f->emflx</div><div> emflx = where(emflx.eq.0.0, -999.0, emflx) ; change 0.0 to -999.0<br></div><div> emflx@_FillValue = -999.0 ; assign _FillValue attribute<br></div><div> printVarSummary(emflx)</div><div> printMinMax(emflx,0)</div><div><br></div><div>Try this and plot the emflx <br></div><div><br></div><div><br></div><div> <br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Dec 30, 2020 at 11:39 AM Setareh Rahimi <<a href="mailto:setareh.rahimi@gmail.com" target="_blank">setareh.rahimi@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">That attached was 'dry deposition flux', please ignore that, and have a look at the attached file.<div>Thanks,</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Dec 30, 2020 at 10:01 PM Setareh Rahimi <<a href="mailto:setareh.rahimi@gmail.com" target="_blank">setareh.rahimi@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">Dear Dennis,<div>This is something weird to me, too.I attached a plot of output model for the same date. Comparing the plot with the value for that lat/lon, results do not match together. What is your idea?</div><div>Many thanks again,</div><div>Best wishes,</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Dec 30, 2020 at 7:09 PM Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">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>???? <span style="background-color:rgb(255,255,255)">does not return t</span><span style="background-color:rgb(255,255,255)">he value of the emflx variable for the desired point (lat/lon)</span></div><div><span style="background-color:rgb(255,255,255)"><br></span></div><div><span style="background-color:rgb(255,255,255)">Well, what is the "correct value" ? <br></span></div><div><span style="background-color:rgb(255,255,255)">======<br></span></div><div><br></div><div> lat = 32.0 <br> lon = 52.0 <br> </div><div>(0) ==================</div>(0) n=21 m=8 lat2d(n,m)=31.9429 lon2d(n,m)=51.921 <b><====This is the closest grid point ... pretty close!!!</b><br><div>(0) ==================</div><div><br></div><div>(5) 605628 ymdh=2019010212<br>(6) 605634 ymdh=2019010218<br><b>(7)</b> 605640 ymdh=<b>2019010300 <=== NT=7</b><br>(8) 605646 ymdh=2019010306<br>(9) 605652 ymdh=2019010312<br></div><div><br></div><div> <b>YMDH = 2019010300 </b> ; user specified<br> print("YMDH="+YMDH)<br> NT = ind(ymdh .eq. YMDH)<br> print("NT="+NT+" ymdh(NT)="+ymdh(NT) )</div><div><br></div><div><b>(0) YMDH=2019010300<br>(0) NT=7 ymdh(NT)=2019010300</b><br>(0) ==================<br></div><div><br></div><div> EMFLX = f->emflx<br> emflx = EMFLX(NT,n,m) ; scalar<br> print("==================")<br> print( emflx ) ; ====</div><div><br></div><div>(0) 940.5043</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Dec 30, 2020 at 8:15 AM Setareh Rahimi <<a href="mailto:setareh.rahimi@gmail.com" target="_blank">setareh.rahimi@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">Dear Dennis,<div>Thanks again for your guide. However, after running the code I noticed tha<span style="background-color:rgb(243,243,243)">t </span><span style="background-color:rgb(0,255,0)"> Y = f->emflx(NT,n,m) </span><span style="background-color:rgb(255,255,255)"> does not return the value of the emflx variable for the desired point (lat/lon). </span></div><div><span style="background-color:rgb(255,255,255)">May I ask you please to give me more advice in this regard?</span></div><div><span style="background-color:rgb(255,255,255)">Best wishes,</span></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Dec 29, 2020 at 3:22 AM Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">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>I think we have been through this before. There are many <a href="https://www.ncl.ucar.edu/Document/Functions/date.shtml" target="_blank"><b>time/date functions.</b></a><br></div><div><br></div><div> time = f->time ; hours since 1949-12-01 00:00:00 UTC" </div><div><div>; [1]<br></div><div> TIME = <a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/cd_calendar.shtml" target="_blank"><b>cd_calendar</b></a>(time, 0) ; PLEASE read the documentation <br> year = toint( TIME(:,0) ) ; toint strips meta data<br> month = toint( TIME(:,1) )<br> day = toint( TIME(:,2) ) ; day of month <br> hour = toint( TIME(:,3) ) ; hour of day</div><div> minute = toint( TIME(:,4) ) <br></div><div> second = toint( TIME(:,5) ) <br></div><div> <br></div><div><br></div><div>;[2] <br></div><div> ymdh = <a href="https://www.ncl.ucar.edu/Document/Functions/Contributed/yyyymmddhh_time.shtml" target="_blank"><b>yyyymmddhh_time</b></a>(time, -3) ; PLEASE read the documentation</div><div><br></div><div> YMDH = 2019010300</div><div> nt =<a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/ind.shtml" target="_blank"><b> ind</b></a>(ymdh .eq. YMDH)</div>================================================================</div><div>A simple script is attached.<br></div><div>================================================================</div><div><br></div><div>%> ncl setareh.ncl</div><div><br></div><div>Please Look-at/READ the documentation of the functions used<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Dec 28, 2020 at 3:04 PM Setareh Rahimi <<a href="mailto:setareh.rahimi@gmail.com" target="_blank">setareh.rahimi@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><br></div><div dir="auto">Dear Dennis,</div><div dir="auto">I am not looking for interpotation, temperature was just an example. Emflx is one of the varibles in the file. I need to know for say Jan-03-00:00UTC , what is the value of emflx for a specific point .</div><div dir="auto">Best wishes,<br><div class="gmail_quote" dir="auto"><div dir="ltr" class="gmail_attr">On Tue, Dec 29, 2020 at 01:23 Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">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"><pre>As noted by Eshan, you want interpolation. <br>Please read the documentation for <a href="https://www.ncl.ucar.edu/Document/Functions/Contributed/rcm2points_Wrap.shtml" target="_blank"><b>rcm2points_Wrap</b></a><br><br>Note: Temperature is NOT in the file.<br></pre><pre><br> f = <a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/addfile.shtml" target="_blank"><strong>addfile</strong></a>("<a href="http://Yazd-2019_DUST01.2019010100.nc" target="_blank">Yazd-2019_DUST01.2019010100.nc</a>" , "r")
lat2d = f->xlat
lon2d = f->xlon
<a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/printVarSummary.shtml" target="_blank"><strong>printVarSummary</strong></a>(lat2d)
<a href="https://www.ncl.ucar.edu/Document/Functions/Contributed/printMinMax.shtml" target="_blank"><strong>printMinMax</strong></a>(lat2d, 0)
<a href="https://www.ncl.ucar.edu/Document/Functions/Contributed/printMinMax.shtml" target="_blank"><strong>printMinMax</strong></a>(lon2d, 0)<br><br> stalat = 32.0 ;(/ 32.0 , 17.31, 24.05 /) ; user specified coordinate pairs
stalon = 52.0 ; (/ 52.0 ,-101.00,-92.46 /)<br><br> X = f->mixrat ; (time, kz, iy, jx)<br> Y = d->ddflx ; (time, iy, jx)<br><br> x = <a href="https://www.ncl.ucar.edu/Document/Functions/Contributed/rcm2points_Wrap.shtml" target="_blank"><b>rcm2points_Wrap</b></a>(lat2d, lon2d, X, stalat, stalon, 0)<br> y = <a href="https://www.ncl.ucar.edu/Document/Functions/Contributed/rcm2points_Wrap.shtml" target="_blank"><b>rcm2points_Wrap</b></a>(lat2d, lon2d, Y, stalat, stalon, 0)<br><br></pre><pre> print(x)<br></pre><pre> print(y)<br></pre></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Dec 28, 2020 at 12:10 PM Setareh Rahimi <<a href="mailto:setareh.rahimi@gmail.com" target="_blank">setareh.rahimi@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">Dear Dennis,<div>Many thanks for your response. What I exactly mean is how to extract for example "Temperature variable" value for a specific lat/lon from the output result. For example: what is temperature value for a point at (lat=32, lon 52)?</div><div>Thanks again</div><div>Best wishes </div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Dec 28, 2020 at 10:17 PM Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">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/Contributed/getind_latlon2d.shtml" target="_blank"><b>getind_latlon2d </b></a></div><div><br></div><div>A modification of Example 1</div><div><br></div><div><pre> f = <a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/addfile.shtml" target="_blank"><strong>addfile</strong></a>("<a href="http://Yazd-2019_DUST01.2019010100.nc" target="_blank">Yazd-2019_DUST01.2019010100.nc</a>" , "r")
lat2d = f->xlat
lon2d = f->xlon
<a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/printVarSummary.shtml" target="_blank"><strong>printVarSummary</strong></a>(lat2d)
<a href="https://www.ncl.ucar.edu/Document/Functions/Contributed/printMinMax.shtml" target="_blank"><strong>printMinMax</strong></a>(lat2d, 0)
<a href="https://www.ncl.ucar.edu/Document/Functions/Contributed/printMinMax.shtml" target="_blank"><strong>printMinMax</strong></a>(lon2d, 0)
lat = (/ 31.0 , 17.31, 24.05 /) ; user specified coordinate pairs
lon = (/ -86.45,-101.00,-92.46 /)
; return 2d subscripts
nm = <strong>getind_latlon2d</strong> (lat2d,lon2d, lat, lon)
<a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/print.shtml" target="_blank"><strong>print</strong></a>(nm)<br><br></pre><pre> X = f->mixrat ; (time, kz, iy, jx)<br> Y = d->ddflx ; (time, iy, jx)<br><br></pre><pre> nt = .... the index some specified time<br></pre><pre> kl = ... index<br></pre><pre> do k=0,<a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/dimsizes.shtml" target="_blank"><strong>dimsizes</strong></a>(lat)-1
n = nm(k,0)
m = nm(k,1)<br></pre><pre> x = X(nt,kl,n,m)<br> y = Y(nt,n,m)<br></pre><pre> <a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/print.shtml" target="_blank"><strong>print</strong></a>(lat2d(n,m)+" "+lon2d(n,m)+" "+x+" "+y)
end do<br><br></pre><pre> <br></pre></div><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Dec 28, 2020 at 10:16 AM 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">Dear all,<div><br></div><div>I need to extract data (it could be any variable like temperature, wind ...) for a specific point (say synoptic stations) at a specific day and time </div><div class="gmail_chip gmail_drive_chip" style="width:396px;height:18px;max-height:18px;background-color:rgb(245,245,245);padding:5px;font-family:arial;font-weight:bold;font-size:13px;border:1px solid rgb(221,221,221);line-height:1"><a href="https://drive.google.com/file/d/105vU4Hx2sRjvWlsmRBKOxa97kXp8Cx1_/view?usp=drive_web" style="display:inline-block;max-width:366px;overflow:hidden;text-overflow:ellipsis;white-space:nowrap;text-decoration-line:none;padding:1px 0px;border:medium none" target="_blank"><img style="vertical-align: bottom; border: medium none;" src="https://ssl.gstatic.com/docs/doclist/images/icon_10_generic_list.png"> <span dir="ltr" style="vertical-align:bottom;text-decoration:none">Yazd-2019_DUST01.2019010100.nc</span></a><img style="opacity: 0.55; float: right; display: none;"></div><div>( Jan-03- 00:00 UTC), from a model output (file attached). I wonder how it could be achieved using NCL (any example...) .</div><div><br clear="all"><div>I appreciate any suggestion,</div><div>Best wishes,</div>-- <br><div dir="ltr">S.Rahimi<br><br></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>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr">S.Rahimi<br><br></div>
</blockquote></div>
</blockquote></div></div>-- <br><div dir="ltr">S.Rahimi<br><br></div>
</blockquote></div>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr">S.Rahimi<br><br></div>
</blockquote></div>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr">S.Rahimi<br><br></div>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr">S.Rahimi<br><br></div>
</blockquote></div>
</blockquote></div>