[ncl-talk] GFS interpolation (vertical and horizontal)

[Ehsan Taghizadeh]

 Dear Adam,Thank you so much for your nice reply.I changed Po to Po(::-1) in int2p_n_Wrap, like:TMPGFS_HVInt = int2p_n_Wrap(lv0_GFS,TMPGFS_HInt,Po(::-1,:),linlog,0), also I tested followings TMPGFS_HVInt = int2p_n_Wrap(lv0_GFS(::-1),TMPGFS_HInt,Po,linlog,0), orTMPGFS_HVInt = int2p_n_Wrap(lv0_GFS(::-1),TMPGFS_HInt(::-1,:),Po,linlog,0).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?I hope to hear from you.SincerelyEhsan


    On Saturday, December 14, 2019, 01:59:21 AM GMT+3:30, Adam Phillips <asphilli at ucar.edu> wrote:  
 
 Hi Ehsan,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.Adam
On Fri, Dec 13, 2019 at 10:34 AM Ehsan Taghizadeh via ncl-talk <ncl-talk at ucar.edu> wrote:

 Hi,Thank you for your guidance.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.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?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.I'll be thankful if anyone could help me?
SincerelyEhsan
    On Wednesday, December 11, 2019, 09:13:15 PM GMT+3:30, Dennis Shea <shea at ucar.edu> wrote:  
 
 ncl-talk does not want to see a script "with many problems." 
We do not have the time to sort/debug multiple issues.---GFS are ordered N->S[1] reorder S->N via NCL syntax:   x(:,:,::-1,:)[2] Use linint2_points_Wrap to interpolate to desired locations[3] Use int2p_n_Wrap to interpolate to desired level(s)


On Wed, Dec 11, 2019 at 4:42 AM Ehsan Taghizadeh via ncl-talk <ncl-talk at ucar.edu> wrote:

Dears,I hope you're doing great.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).Should I attach my script, however it faces many problems?I'll be thankful for any help.
SincerelyEhsan
_______________________________________________
ncl-talk mailing list
ncl-talk at ucar.edu
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
  _______________________________________________
ncl-talk mailing list
ncl-talk at ucar.edu
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk


-- 
Adam Phillips 
Associate Scientist,  Climate and Global Dynamics Laboratory, NCAR
www.cgd.ucar.edu/staff/asphilli/   303-497-1726 
  
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191214/0423b9e3/attachment.html>