<div dir="ltr">Dear Mr. Andrew<div><br></div><div>Many thanks for your great support.</div><div><br></div><div>I think my problem is, how to use the lterp function. where should I write the lterp funtion</div><div>and how to write is. I just found some sample in using lterp command, but It did work.</div><div><br></div><div>Based on the sample script that I have sent in the previous email, where should I locate my lterp </div><div>command?</div><div><br></div><div>Thank you very much.</div><div><br></div><div>Best regards,</div><div>wendi</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sun, Jan 22, 2017 at 11:43 PM, Andrew Friedman <span dir="ltr">&lt;<a href="mailto:andfried@gmail.com" target="_blank">andfried@gmail.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Wendi,<br>
<br>
I don’t see an error offhand with your dimensions.<br>
What does the numerical output of lterp look like in the regions where you have data? Do you get any values?<br>
<span class="HOEnZb"><font color="#888888"><br>
Andrew<br>
</font></span><div class="HOEnZb"><div class="h5"><br>
&gt; On Jan 21, 2017, at 4:11 AM, wendi harjupa &lt;<a href="mailto:wendiharjupa@gmail.com">wendiharjupa@gmail.com</a>&gt; wrote:<br>
&gt;<br>
&gt; Dear Mr. Andrew,<br>
&gt;<br>
&gt; Many thanks for your great support,<br>
&gt;<br>
&gt; Below is q dims and q file of my data,<br>
&gt;<br>
&gt; ## Data 1 (large scale : 0.1x0.1 deg)<br>
&gt;<br>
&gt; q dims 1<br>
&gt; X is varying   Lon = 121.05 to 148.95   X = 1 to 280<br>
&gt; Y is varying   Lat = 22.95 to 46.85   Y = 1 to 240<br>
&gt; Z is fixed     Lev = 1  Z = 1<br>
&gt; T is fixed     Time = 13:05Z08AUG2015  T = 1<br>
&gt; E is fixed     Ens = 1  E = 1<br>
&gt;<br>
&gt; q file 1<br>
&gt;<br>
&gt;   Descriptor: Him.ctl<br>
&gt;   Binary: /TAKAYASU/wendi/RDCA/rdca_<wbr>201508080520.dat<br>
&gt;   Type = Gridded<br>
&gt;   Xsize = 280  Ysize = 240  Zsize = 1  Tsize = 1  Esize = 1<br>
&gt;   Number of Variables = 1<br>
&gt;      albedo  0  1  reflectivity<br>
&gt;<br>
&gt;<br>
&gt; ## data 2 (small scale : 0.002x0.003 deg)<br>
&gt;<br>
&gt; q dims 2<br>
&gt;<br>
&gt; X is varying   Lon = 134.375 to 136.872   X = 1 to 800<br>
&gt; Y is varying   Lat = 33.5833 to 35.5812   Y = 1 to 960<br>
&gt; Z is fixed     Lev = 1  Z = 1<br>
&gt; T is fixed     Time = 00Z08AUG2015  T = 1<br>
&gt; E is fixed     Ens = 1  E = 1<br>
&gt;<br>
&gt; q file 2<br>
&gt;<br>
&gt;   Binary: /TAKAYASU/wendi/xband/G000.<wbr>grads/KINKI/20150808_0900-<wbr>1700/data/KINKI00001-    20150808-1420-G000-EL000000.<wbr>bin<br>
&gt;   Type = Gridded<br>
&gt;   Xsize = 800  Ysize = 960  Zsize = 1  Tsize = 120  Esize = 1<br>
&gt;   Number of Variables = 1<br>
&gt;      rr  0  98  rairate [mm/hr]<br>
&gt;<br>
&gt; is there any mistakes on that information?<br>
&gt;<br>
&gt; Best regards,<br>
&gt; wendi<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt; On Fri, Jan 20, 2017 at 6:25 PM, Andrew Friedman &lt;<a href="mailto:andfried@gmail.com">andfried@gmail.com</a>&gt; wrote:<br>
&gt; Hi Wendi,<br>
&gt; What is your dimension environment?<br>
&gt; Can you post the output of ‘q dims’, ‘q file 1’ and ‘q file 2&#39;<br>
&gt; Andrew<br>
&gt;<br>
&gt; &gt; On Jan 20, 2017, at 1:28 AM, wendi harjupa &lt;<a href="mailto:wendiharjupa@gmail.com">wendiharjupa@gmail.com</a>&gt; wrote:<br>
&gt; &gt;<br>
&gt; &gt; Dear Mr. Andrew<br>
&gt; &gt;<br>
&gt; &gt; Thank you very much for your explanation,<br>
&gt; &gt; I am still struggling to solve it.<br>
&gt; &gt;<br>
&gt; &gt; Best regards,<br>
&gt; &gt; wendi<br>
&gt; &gt;<br>
&gt; &gt; On Wed, Jan 18, 2017 at 10:16 PM, Andrew Friedman &lt;<a href="mailto:andfried@gmail.com">andfried@gmail.com</a>&gt; wrote:<br>
&gt; &gt; Hi Wendi,<br>
&gt; &gt;<br>
&gt; &gt; You can check ‘q dims’ to see what your time dimensions are.<br>
&gt; &gt;<br>
&gt; &gt; I recommend going through the website examples for <a href="http://lterp.gs" rel="noreferrer" target="_blank">lterp.gs</a>, particularly #3-5 - these show situations that often come up about mismatched dimensions.<br>
&gt; &gt; As they show, you often need to specify time or dfile, or some combination. Another possibility could be wrapped longitude coordinates.<br>
&gt; &gt;<br>
&gt; &gt; If none of these steps work, can you post a simplified version of your problem, without the anciallary scripts in your code?<br>
&gt; &gt;<br>
&gt; &gt; Andrew<br>
&gt; &gt;<br>
&gt; &gt; &gt; On Jan 18, 2017, at 2:56 AM, wendi harjupa &lt;<a href="mailto:wendiharjupa@gmail.com">wendiharjupa@gmail.com</a>&gt; wrote:<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; Dear Mr Andrew,<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; Many thanks for your response<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; There is no specific error, but the is no change of result when I use<br>
&gt; &gt; &gt; &#39;d lterp command&#39;.<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; I set my dimension environment as below :<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;     &#39;set lat 33.000  36.000&#39;<br>
&gt; &gt; &gt;     &#39;set lon 134.000  137.000&#39;<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; I am little confused how to use the &#39;d lterp&#39;  function in my script program.<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; Would you mind to give me more clue?<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; Best regards,<br>
&gt; &gt; &gt; wendi<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; On Wed, Jan 18, 2017 at 5:55 AM, Andrew Friedman &lt;<a href="mailto:andfried@gmail.com">andfried@gmail.com</a>&gt; wrote:<br>
&gt; &gt; &gt; Hi Wendi,<br>
&gt; &gt; &gt; What specific error message are you getting?<br>
&gt; &gt; &gt; A typical problem is that the time axes of the two files have to be the same - you may want to check your dimension environment.<br>
&gt; &gt; &gt; Andrew<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; &gt; On Jan 17, 2017, at 1:53 PM, wendi harjupa &lt;<a href="mailto:wendiharjupa@gmail.com">wendiharjupa@gmail.com</a>&gt; wrote:<br>
&gt; &gt; &gt; &gt;<br>
&gt; &gt; &gt; &gt; Dear Mr. Andrew and colleagues<br>
&gt; &gt; &gt; &gt;<br>
&gt; &gt; &gt; &gt; Please tell me what wrong with this script. It does not work to interpolate the coarse grid to fine grid<br>
&gt; &gt; &gt; &gt;<br>
&gt; &gt; &gt; &gt; ##first data/resolution 0.002x0.003 (lat, lon)<br>
&gt; &gt; &gt; &gt;<br>
&gt; &gt; &gt; &gt;    &#39;open ./Rain.ctl&#39;<br>
&gt; &gt; &gt; &gt;    &#39;setXrain.gs&#39;<br>
&gt; &gt; &gt; &gt;    &#39;setCOMMON.gs&#39;    ! set parameter (lon, lat, parea, grid off, etc)<br>
&gt; &gt; &gt; &gt; #   &#39;d rr&#39;      ! variable<br>
&gt; &gt; &gt; &gt;    &#39;xcbar  8.85 9.05 0.8 4 -fw 0.11 -fh 0.15 -edge triangle -fstep 1 -line on -lc 1&#39;<br>
&gt; &gt; &gt; &gt;    &#39;close 1&#39;<br>
&gt; &gt; &gt; &gt;<br>
&gt; &gt; &gt; &gt; ##second data/resolution 0.1x0.1 (lat, lon)<br>
&gt; &gt; &gt; &gt;<br>
&gt; &gt; &gt; &gt;      &#39;open ./Him.ctl&#39;<br>
&gt; &gt; &gt; &gt;     &#39;setCOMMON.gs&#39;<br>
&gt; &gt; &gt; &gt;     &#39;setHimawari.gs&#39;<br>
&gt; &gt; &gt; &gt; #    &#39;d albedo&#39;    ! variable<br>
&gt; &gt; &gt; &gt;     &#39;xcbar 6.3 9 0.23 0.37 -fw 0.11 -fh 0.15 -edge triangle -fs 2 -fo 1 -line on&#39;<br>
&gt; &gt; &gt; &gt;      &#39;close 1&#39;<br>
&gt; &gt; &gt; &gt;<br>
&gt; &gt; &gt; &gt;      &#39;d lterp(rr, albedo.2, aave)&#39;<br>
&gt; &gt; &gt; &gt;<br>
&gt; &gt; &gt; &gt; Any help will be very appreciated.<br>
&gt; &gt; &gt; &gt; Thank you very much.<br>
&gt; &gt; &gt; &gt;<br>
&gt; &gt; &gt; &gt; Best regards,<br>
&gt; &gt; &gt; &gt; wendi<br>
&gt; &gt; &gt; &gt;<br>
&gt; &gt; &gt; &gt;<br>
&gt; &gt; &gt; &gt;<br>
&gt; &gt; &gt; &gt;<br>
&gt; &gt; &gt; &gt;<br>
&gt; &gt; &gt; &gt;<br>
&gt; &gt; &gt; &gt;<br>
&gt; &gt; &gt; &gt;<br>
&gt; &gt; &gt; &gt;<br>
&gt; &gt; &gt; &gt;<br>
&gt; &gt; &gt; &gt;<br>
&gt; &gt; &gt; &gt;<br>
&gt; &gt; &gt; &gt; On Tue, Jan 17, 2017 at 7:12 PM, wendi harjupa &lt;<a href="mailto:wendiharjupa@gmail.com">wendiharjupa@gmail.com</a>&gt; wrote:<br>
&gt; &gt; &gt; &gt; Dear Mr. Andrew..<br>
&gt; &gt; &gt; &gt;<br>
&gt; &gt; &gt; &gt; Many thanks for your valuable information.<br>
&gt; &gt; &gt; &gt;<br>
&gt; &gt; &gt; &gt; Best regards,<br>
&gt; &gt; &gt; &gt; wendi<br>
&gt; &gt; &gt; &gt;<br>
&gt; &gt; &gt; &gt; On Tue, Jan 17, 2017 at 6:45 PM, Andrew Friedman &lt;<a href="mailto:andfried@gmail.com">andfried@gmail.com</a>&gt; wrote:<br>
&gt; &gt; &gt; &gt; Hi Wendi,<br>
&gt; &gt; &gt; &gt; You can use the function lterp for spatial interpolation: <a href="http://cola.gmu.edu/grads/gadoc/gradfunclterp.html" rel="noreferrer" target="_blank">http://cola.gmu.edu/grads/<wbr>gadoc/gradfunclterp.html</a><br>
&gt; &gt; &gt; &gt; -Andrew<br>
&gt; &gt; &gt; &gt;<br>
&gt; &gt; &gt; &gt; &gt; On Jan 17, 2017, at 10:29 AM, wendi harjupa &lt;<a href="mailto:wendiharjupa@gmail.com">wendiharjupa@gmail.com</a>&gt; wrote:<br>
&gt; &gt; &gt; &gt; &gt;<br>
&gt; &gt; &gt; &gt; &gt; Dear All<br>
&gt; &gt; &gt; &gt; &gt;<br>
&gt; &gt; &gt; &gt; &gt; I got problem in overlaying two type of data.<br>
&gt; &gt; &gt; &gt; &gt; As I have attached the plot of data in this email.<br>
&gt; &gt; &gt; &gt; &gt;<br>
&gt; &gt; &gt; &gt; &gt; First data is rain data with resolution is about 0.002x0.003 deg (lon ,lat).<br>
&gt; &gt; &gt; &gt; &gt; The rain data is indicated by rainbow color.<br>
&gt; &gt; &gt; &gt; &gt; and second data is cloud data with resolution is 0.1x0.1 deg (lon, lat)<br>
&gt; &gt; &gt; &gt; &gt; The cloud data is indicated by gray color<br>
&gt; &gt; &gt; &gt; &gt;<br>
&gt; &gt; &gt; &gt; &gt; The plot looks not good. Is there any tools in grads to smooth the cloud data?<br>
&gt; &gt; &gt; &gt; &gt; Then, we can interpolate cloud data to become 0.005 or 0.001, and plot data will be more nice.<br>
&gt; &gt; &gt; &gt; &gt;<br>
&gt; &gt; &gt; &gt; &gt; Thank you very much for your attention and help.<br>
&gt; &gt; &gt; &gt; &gt;<br>
&gt; &gt; &gt; &gt; &gt; Best regards,<br>
&gt; &gt; &gt; &gt; &gt; wendi<br>
&gt; &gt; &gt; &gt; &gt;<br>
&gt; &gt; &gt; &gt; &gt; &lt;large skale.png&gt;____________________<wbr>___________________________<br>
&gt; &gt; &gt; &gt; &gt; gradsusr mailing list<br>
&gt; &gt; &gt; &gt; &gt; <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
&gt; &gt; &gt; &gt; &gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/<wbr>listinfo/gradsusr</a><br>
&gt; &gt; &gt; &gt;<br>
&gt; &gt; &gt; &gt;<br>
&gt; &gt; &gt; &gt; ______________________________<wbr>_________________<br>
&gt; &gt; &gt; &gt; gradsusr mailing list<br>
&gt; &gt; &gt; &gt; <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
&gt; &gt; &gt; &gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/<wbr>listinfo/gradsusr</a><br>
&gt; &gt; &gt; &gt;<br>
&gt; &gt; &gt; &gt;<br>
&gt; &gt; &gt; &gt; ______________________________<wbr>_________________<br>
&gt; &gt; &gt; &gt; gradsusr mailing list<br>
&gt; &gt; &gt; &gt; <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
&gt; &gt; &gt; &gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/<wbr>listinfo/gradsusr</a><br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; ______________________________<wbr>_________________<br>
&gt; &gt; &gt; gradsusr mailing list<br>
&gt; &gt; &gt; <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
&gt; &gt; &gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/<wbr>listinfo/gradsusr</a><br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; ______________________________<wbr>_________________<br>
&gt; &gt; &gt; gradsusr mailing list<br>
&gt; &gt; &gt; <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
&gt; &gt; &gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/<wbr>listinfo/gradsusr</a><br>
&gt; &gt;<br>
&gt; &gt;<br>
&gt; &gt; ______________________________<wbr>_________________<br>
&gt; &gt; gradsusr mailing list<br>
&gt; &gt; <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
&gt; &gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/<wbr>listinfo/gradsusr</a><br>
&gt; &gt;<br>
&gt; &gt; ______________________________<wbr>_________________<br>
&gt; &gt; gradsusr mailing list<br>
&gt; &gt; <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
&gt; &gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/<wbr>listinfo/gradsusr</a><br>
&gt;<br>
&gt;<br>
&gt; ______________________________<wbr>_________________<br>
&gt; gradsusr mailing list<br>
&gt; <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
&gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/<wbr>listinfo/gradsusr</a><br>
&gt;<br>
&gt; ______________________________<wbr>_________________<br>
&gt; gradsusr mailing list<br>
&gt; <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
&gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/<wbr>listinfo/gradsusr</a><br>
<br>
<br>
______________________________<wbr>_________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/<wbr>listinfo/gradsusr</a><br>
</div></div></blockquote></div><br></div>