<div dir="ltr"><div>The original code is a bit messy for debugging. Really, it is best to send a minimalist code to illustrate the problem(s).</div><div><br></div><div>--</div><div>See attachment: tst.gfsInterpolate.ncl</div><div>I used an available gfs...grib2 file.</div><div>---</div><div>The end result for 3 stations is:</div><div><br></div><div>(0)      85000   305.16   300.76   294.79<br>(1)       70000   292.31   291.38   279.63<br>(2)       60000   281.88   281.23   274.28<br>(3)       50000   269.71   269.62   264.71<br>(4)       40000   258.74   258.90   251.70<br>(5)       35000   253.37   253.66   244.32<br>(6)       30000   245.56   246.03   238.54<br>(7)       27500   240.79   241.26   236.50<br>(8)       25000   236.03   236.48   234.46<br>(9)       22500   230.07   230.59   230.70<br>(10)      20000   224.10   224.70   226.93<br>(11)      15000   210.87   210.74   216.14<br>(12)      10000   197.31   198.75   205.43</div><div><br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Dec 13, 2019 at 11:47 PM Ehsan Taghizadeh via ncl-talk <<a href="mailto:ncl-talk@ucar.edu">ncl-talk@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><div style="font-family:bookman old style,new york,times,serif;font-size:16px"><div></div>
        <div dir="ltr"><div><div>Dear Adam,</div><div>Thank you so much for your nice reply.</div><div>I changed Po to Po(::-1) in int2p_n_Wrap, like:</div><div>TMPGFS_HVInt = int2p_n_Wrap(lv0_GFS,TMPGFS_HInt,Po(::-1,:),linlog,0), also I tested followings </div><div>TMPGFS_HVInt = int2p_n_Wrap(lv0_GFS(::-1),TMPGFS_HInt,Po,linlog,0), or</div><div>TMPGFS_HVInt = int2p_n_Wrap(lv0_GFS(::-1),TMPGFS_HInt(::-1,:),Po,linlog,0).</div><div>Unfortunately, bad results continued. The temperature are wrong, relative humidity, wind speed and wind direction are sth like miss values for all points and levels. May I ask for any idea?</div><div>I hope to hear from you.</div><div>Sincerely</div><div>Ehsan</div><div><br></div></div><br></div><div><br></div>
        
        </div><div id="gmail-m_-1391252305009285122yahoo_quoted_6842090546">
            <div style="font-family:"Helvetica Neue",Helvetica,Arial,sans-serif;font-size:13px;color:rgb(38,40,42)">
                
                <div>
                    On Saturday, December 14, 2019, 01:59:21 AM GMT+3:30, Adam Phillips <<a href="mailto:asphilli@ucar.edu" target="_blank">asphilli@ucar.edu</a>> wrote:
                </div>
                <div><br></div>
                <div><br></div>
                <div><div id="gmail-m_-1391252305009285122yiv8767773030"><div><div dir="ltr">Hi Ehsan,<div>Yes, pin and pout should both either be decreasing or increasing. I see in your script that pin (lv0_GFS) is increasing and pout (Po) is decreasing. Alter one so both are behaving similarly, and hopefully that fixes your issues.</div><div>Adam</div></div><br clear="none"><div><div id="gmail-m_-1391252305009285122yiv8767773030yqt74963"><div dir="ltr">On Fri, Dec 13, 2019 at 10:34 AM Ehsan Taghizadeh via ncl-talk <<a rel="nofollow" shape="rect" href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>> wrote:<br clear="none"></div><blockquote style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><div><div></div>
        <div dir="ltr">Hi,</div><div dir="ltr">Thank you for your guidance.</div><div dir="ltr">I tried to do that (reorder, linint2_points_Wrap,int2p_n_Wrap), however it seems the result is not good at all. Relative humidity, wind speed and wind direction have written the same values. Temperature is not correct, I think.</div><div dir="ltr">May I ask does in int2p_n_Wrap(pin,xin,pout,linlog,pdim), "pin" and "pout" must both increasing or decreasing, or pin could be increasing and pout could be decreasing?</div><div dir="ltr">I've attached the script (intepgfs.ncl), input station locations (apist.dat) and the output (Interpolated.csv). I also put gfs.tar file in the ftp, inclulding two input gfs data files.</div><div dir="ltr">I'll be thankful if anyone could help me?</div><div dir="ltr"><br clear="none"></div><div dir="ltr">Sincerely</div><div dir="ltr">Ehsan</div><div><br clear="none"></div>
        
        </div><div id="gmail-m_-1391252305009285122yiv8767773030gmail-m_-6537415225624962878yahoo_quoted_6913350504">
            <div>
                
                <div>
                    On Wednesday, December 11, 2019, 09:13:15 PM GMT+3:30, Dennis Shea <<a rel="nofollow" shape="rect" href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>> wrote:
                </div>
                <div><br clear="none"></div>
                <div><br clear="none"></div>
                <div><div id="gmail-m_-1391252305009285122yiv8767773030gmail-m_-6537415225624962878yiv8130148666"><div><div dir="ltr"><div>ncl-talk does not want to see a script "with many problems." <br clear="none"></div><div>We do not have the time to sort/debug multiple issues.</div><div>---</div><div>GFS are ordered N->S</div><div>[1] reorder S->N via NCL syntax:   x(:,:,::-1,:)</div><div>[2] Use <a rel="nofollow" shape="rect" href="http://www.ncl.ucar.edu/Document/Functions/Contributed/linint2_points_Wrap.shtml" target="_blank"><b>linint2_points_Wrap</b></a> to interpolate to desired locations</div><div>[3] Use <a rel="nofollow" shape="rect" href="http://www.ncl.ucar.edu/Document/Functions/Contributed/int2p_n_Wrap.shtml" target="_blank"><b>int2p_n_Wrap</b></a> to interpolate to desired level(s)</div><div><br clear="none"></div><div><br clear="none"></div></div><br clear="none"><div><div id="gmail-m_-1391252305009285122yiv8767773030gmail-m_-6537415225624962878yiv8130148666yqt90606"><div dir="ltr">On Wed, Dec 11, 2019 at 4:42 AM Ehsan Taghizadeh via ncl-talk <<a rel="nofollow" shape="rect" href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>> wrote:<br clear="none"></div><blockquote style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><div><div dir="ltr"><div><div>Dears,</div><div>I hope you're doing great.</div><div>May I ask if anybody has a script to interpolate GFS (0.5 deg) on the desired station points and also the desired vertical level? I mean (both horizontal and vertical interpolation).</div><div dir="ltr">Should I attach my script, however it faces many problems?</div><div>I'll be thankful for any help.</div><div><br clear="none"></div><div>Sincerely</div><div>Ehsan</div></div><br clear="none"></div></div></div>_______________________________________________<br clear="none">
ncl-talk mailing list<br clear="none">
<a rel="nofollow" shape="rect" href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br clear="none">
List instructions, subscriber options, unsubscribe:<br clear="none">
<a rel="nofollow" shape="rect" href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a></blockquote></div></div></div></div></div>
            </div>
        </div></div>_______________________________________________<br clear="none">
ncl-talk mailing list<br clear="none">
<a rel="nofollow" shape="rect" href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br clear="none">
List instructions, subscriber options, unsubscribe:<br clear="none">
<a rel="nofollow" shape="rect" href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a></blockquote></div></div><br clear="all"><div><br clear="none"></div>-- <br clear="none"><div dir="ltr"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div><div><span><font style="background-color:inherit" color="#888888">Adam Phillips <br clear="none"></font></span></div><span><font style="background-color:inherit" color="#888888">Associate Scientist,  </font></span><span><font style="background-color:inherit" color="#888888">Climate and Global Dynamics Laboratory, NCAR<br clear="none"></font></span></div></div><div><span><font style="background-color:inherit" color="#888888"><a rel="nofollow" shape="rect" href="http://www.cgd.ucar.edu/staff/asphilli/" target="_blank">www.cgd.ucar.edu/staff/asphilli/</a>   </font></span><span><font style="background-color:inherit" color="#888888">303-497-1726 </font></span></div><span><font style="background-color:inherit" color="#888888"></font></span><div><div><span><font style="background-color:inherit" color="#888888"><br clear="none"></font></span><div><span><font style="background-color:inherit" color="#888888"><a rel="nofollow" shape="rect" href="http://www.cgd.ucar.edu/staff/asphilli" target="_blank"></a></font></span></div></div></div></div></div></div></div></div></div></div></div></div></div></div>
            </div>
        </div></div>_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a></blockquote></div>