<div dir="ltr">Dear Mr. Andrew<div><br></div><div>Thank you very much for your excellent supervision<br><div><br></div><div>Below is my simple script</div><div><br></div><div><div style="font-size:12.8px"> 'open ./Rain.ctl'</div><div style="font-size:12.8px"> 'set lat 33.5 35.5'</div><div style="font-size:12.8px"> 'set lon 134.3 136.6'</div><div style="font-size:12.8px"># 'd rr' ! variable</div><div style="font-size:12.8px"> 'xcbar 8.85 9.05 0.8 4 -fw 0.11 -fh 0.15 -edge triangle -fstep 1 -line on -lc 1'</div><div style="font-size:12.8px"> 'close 1'</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">##second data/resolution 0.1x0.1 (lat, lon)</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px"> 'open ./Him.ctl'</div><div style="font-size:12.8px"> 'set lat 33.5 35.5'</div><div style="font-size:12.8px"><span style="font-size:12.8px"> 'set lon 134.3 136.6'</span></div><div style="font-size:12.8px"><span style="font-size:12.8px"># 'd albedo' ! variable</span><br></div><div style="font-size:12.8px"> 'xcbar 6.3 9 0.23 0.37 -fw 0.11 -fh 0.15 -edge triangle -fs 2 -fo 1 -line on'</div><div style="font-size:12.8px"> 'close 1'</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px"> 'd lterp(rr, albedo.2, aave)'</div></div></div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">my control file as below</div><div style="font-size:12.8px">1. Rain.ctl</div><div style="font-size:12.8px"><span style="font-size:12.8px"><br></span></div><div style="font-size:12.8px"><span style="font-size:12.8px">DSET ~/data1</span></div><div style="font-size:12.8px">OPTIONS yrev</div><div style="font-size:12.8px">UNDEF -1.E35</div><div style="font-size:12.8px"> XDEF 800 LINEAR 134.375 0.003125</div><div style="font-size:12.8px"> YDEF 960 LINEAR 33.5833 0.00208333</div><div style="font-size:12.8px">ZDEF 1 LINEAR 1 1</div><div style="font-size:12.8px">TDEF 120 LINEAR 08Aug2015 1hr </div><div style="font-size:12.8px">VARS 1 </div><div style="font-size:12.8px">rr 0 98 rairate [mm/hr] </div><div style="font-size:12.8px">ENDVARS</div><div><br></div><div style="font-size:12.8px">2. Him.ctl</div><div style="font-size:12.8px"><div style="font-size:12.8px">DSET ${file_in_vor2}</div><div style="font-size:12.8px">UNDEF -1.0 </div><div style="font-size:12.8px">OPTIONS yrev</div><div style="font-size:12.8px">XDEF 280 LINEAR 121.05 0.1</div><div style="font-size:12.8px">YDEF 240 LINEAR 22.95 0.1 </div><div style="font-size:12.8px">ZDEF 1 LINEAR 1 1 </div><div style="font-size:12.8px">TDEF 1 LINEAR 13:05Z08AUG2015 5mn </div><div style="font-size:12.8px">VARS 1</div><div style="font-size:12.8px">ALBEDO 0 1 reflectivity</div><div style="font-size:12.8px">ENDVARS</div><div><br></div><div>My next question where should I locate my lterp command??</div><div>and how it would be like?</div><div>I have tried this command below, it doesnor work</div><div><br></div><div>'define try=(rr, albedo.2, aave)'</div><div>'d try'</div><div><br></div><div>Thank you very much for your help. I am very sorry for ask many questions.</div><div><br></div><div>Best regards,</div><div>wendi</div><div><br></div><div><br></div><div><br></div></div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px"><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Jan 23, 2017 at 6:48 PM, Andrew Friedman <span dir="ltr"><<a href="mailto:andfried@gmail.com" target="_blank">andfried@gmail.com</a>></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>
You can use the command as in examples 4 & 5 for the lterp documentation:<br>
'define newvar = (rr, albedo.2, aave)’<br>
What does newvar look like over specific gridboxes where you know there is data?<br>
<br>
Again, I would try first running a simplified version of your script that doesn’t call other scripts in order to simplify the problem as much as possible.<br>
<br>
Best, Andrew<br>
<div class="HOEnZb"><div class="h5"><br>
> On Jan 23, 2017, at 3:08 AM, wendi harjupa <<a href="mailto:wendiharjupa@gmail.com">wendiharjupa@gmail.com</a>> wrote:<br>
><br>
> Dear Mr. Andrew<br>
><br>
> Many thanks for your great support.<br>
><br>
> I think my problem is, how to use the lterp function. where should I write the lterp funtion<br>
> and how to write is. I just found some sample in using lterp command, but It did work.<br>
><br>
> Based on the sample script that I have sent in the previous email, where should I locate my lterp<br>
> command?<br>
><br>
> Thank you very much.<br>
><br>
> Best regards,<br>
> wendi<br>
><br>
> On Sun, Jan 22, 2017 at 11:43 PM, Andrew Friedman <<a href="mailto:andfried@gmail.com">andfried@gmail.com</a>> wrote:<br>
> 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>
><br>
> Andrew<br>
><br>
> > On Jan 21, 2017, at 4:11 AM, wendi harjupa <<a href="mailto:wendiharjupa@gmail.com">wendiharjupa@gmail.com</a>> wrote:<br>
> ><br>
> > Dear Mr. Andrew,<br>
> ><br>
> > Many thanks for your great support,<br>
> ><br>
> > Below is q dims and q file of my data,<br>
> ><br>
> > ## Data 1 (large scale : 0.1x0.1 deg)<br>
> ><br>
> > q dims 1<br>
> > X is varying Lon = 121.05 to 148.95 X = 1 to 280<br>
> > Y is varying Lat = 22.95 to 46.85 Y = 1 to 240<br>
> > Z is fixed Lev = 1 Z = 1<br>
> > T is fixed Time = 13:05Z08AUG2015 T = 1<br>
> > E is fixed Ens = 1 E = 1<br>
> ><br>
> > q file 1<br>
> ><br>
> > Descriptor: Him.ctl<br>
> > Binary: /TAKAYASU/wendi/RDCA/rdca_<wbr>201508080520.dat<br>
> > Type = Gridded<br>
> > Xsize = 280 Ysize = 240 Zsize = 1 Tsize = 1 Esize = 1<br>
> > Number of Variables = 1<br>
> > albedo 0 1 reflectivity<br>
> ><br>
> ><br>
> > ## data 2 (small scale : 0.002x0.003 deg)<br>
> ><br>
> > q dims 2<br>
> ><br>
> > X is varying Lon = 134.375 to 136.872 X = 1 to 800<br>
> > Y is varying Lat = 33.5833 to 35.5812 Y = 1 to 960<br>
> > Z is fixed Lev = 1 Z = 1<br>
> > T is fixed Time = 00Z08AUG2015 T = 1<br>
> > E is fixed Ens = 1 E = 1<br>
> ><br>
> > q file 2<br>
> ><br>
> > Binary: /TAKAYASU/wendi/xband/G000.<wbr>grads/KINKI/20150808_0900-<wbr>1700/data/KINKI00001- 20150808-1420-G000-EL000000.<wbr>bin<br>
> > Type = Gridded<br>
> > Xsize = 800 Ysize = 960 Zsize = 1 Tsize = 120 Esize = 1<br>
> > Number of Variables = 1<br>
> > rr 0 98 rairate [mm/hr]<br>
> ><br>
> > is there any mistakes on that information?<br>
> ><br>
> > Best regards,<br>
> > wendi<br>
> ><br>
> ><br>
> ><br>
> ><br>
> > On Fri, Jan 20, 2017 at 6:25 PM, Andrew Friedman <<a href="mailto:andfried@gmail.com">andfried@gmail.com</a>> wrote:<br>
> > Hi Wendi,<br>
> > What is your dimension environment?<br>
> > Can you post the output of ‘q dims’, ‘q file 1’ and ‘q file 2'<br>
> > Andrew<br>
> ><br>
> > > On Jan 20, 2017, at 1:28 AM, wendi harjupa <<a href="mailto:wendiharjupa@gmail.com">wendiharjupa@gmail.com</a>> wrote:<br>
> > ><br>
> > > Dear Mr. Andrew<br>
> > ><br>
> > > Thank you very much for your explanation,<br>
> > > I am still struggling to solve it.<br>
> > ><br>
> > > Best regards,<br>
> > > wendi<br>
> > ><br>
> > > On Wed, Jan 18, 2017 at 10:16 PM, Andrew Friedman <<a href="mailto:andfried@gmail.com">andfried@gmail.com</a>> wrote:<br>
> > > Hi Wendi,<br>
> > ><br>
> > > You can check ‘q dims’ to see what your time dimensions are.<br>
> > ><br>
> > > 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>
> > > As they show, you often need to specify time or dfile, or some combination. Another possibility could be wrapped longitude coordinates.<br>
> > ><br>
> > > If none of these steps work, can you post a simplified version of your problem, without the anciallary scripts in your code?<br>
> > ><br>
> > > Andrew<br>
> > ><br>
> > > > On Jan 18, 2017, at 2:56 AM, wendi harjupa <<a href="mailto:wendiharjupa@gmail.com">wendiharjupa@gmail.com</a>> wrote:<br>
> > > ><br>
> > > > Dear Mr Andrew,<br>
> > > ><br>
> > > > Many thanks for your response<br>
> > > ><br>
> > > > There is no specific error, but the is no change of result when I use<br>
> > > > 'd lterp command'.<br>
> > > ><br>
> > > > I set my dimension environment as below :<br>
> > > ><br>
> > > > 'set lat 33.000 36.000'<br>
> > > > 'set lon 134.000 137.000'<br>
> > > ><br>
> > > > I am little confused how to use the 'd lterp' function in my script program.<br>
> > > ><br>
> > > > Would you mind to give me more clue?<br>
> > > ><br>
> > > > Best regards,<br>
> > > > wendi<br>
> > > ><br>
> > > ><br>
> > > ><br>
> > > ><br>
> > > ><br>
> > > > On Wed, Jan 18, 2017 at 5:55 AM, Andrew Friedman <<a href="mailto:andfried@gmail.com">andfried@gmail.com</a>> wrote:<br>
> > > > Hi Wendi,<br>
> > > > What specific error message are you getting?<br>
> > > > 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>
> > > > Andrew<br>
> > > ><br>
> > > > > On Jan 17, 2017, at 1:53 PM, wendi harjupa <<a href="mailto:wendiharjupa@gmail.com">wendiharjupa@gmail.com</a>> wrote:<br>
> > > > ><br>
> > > > > Dear Mr. Andrew and colleagues<br>
> > > > ><br>
> > > > > Please tell me what wrong with this script. It does not work to interpolate the coarse grid to fine grid<br>
> > > > ><br>
> > > > > ##first data/resolution 0.002x0.003 (lat, lon)<br>
> > > > ><br>
> > > > > 'open ./Rain.ctl'<br>
> > > > > 'setXrain.gs'<br>
> > > > > 'setCOMMON.gs' ! set parameter (lon, lat, parea, grid off, etc)<br>
> > > > > # 'd rr' ! variable<br>
> > > > > 'xcbar 8.85 9.05 0.8 4 -fw 0.11 -fh 0.15 -edge triangle -fstep 1 -line on -lc 1'<br>
> > > > > 'close 1'<br>
> > > > ><br>
> > > > > ##second data/resolution 0.1x0.1 (lat, lon)<br>
> > > > ><br>
> > > > > 'open ./Him.ctl'<br>
> > > > > 'setCOMMON.gs'<br>
> > > > > 'setHimawari.gs'<br>
> > > > > # 'd albedo' ! variable<br>
> > > > > 'xcbar 6.3 9 0.23 0.37 -fw 0.11 -fh 0.15 -edge triangle -fs 2 -fo 1 -line on'<br>
> > > > > 'close 1'<br>
> > > > ><br>
> > > > > 'd lterp(rr, albedo.2, aave)'<br>
> > > > ><br>
> > > > > Any help will be very appreciated.<br>
> > > > > Thank you very much.<br>
> > > > ><br>
> > > > > Best regards,<br>
> > > > > wendi<br>
> > > > ><br>
> > > > ><br>
> > > > ><br>
> > > > ><br>
> > > > ><br>
> > > > ><br>
> > > > ><br>
> > > > ><br>
> > > > ><br>
> > > > ><br>
> > > > ><br>
> > > > ><br>
> > > > > On Tue, Jan 17, 2017 at 7:12 PM, wendi harjupa <<a href="mailto:wendiharjupa@gmail.com">wendiharjupa@gmail.com</a>> wrote:<br>
> > > > > Dear Mr. Andrew..<br>
> > > > ><br>
> > > > > Many thanks for your valuable information.<br>
> > > > ><br>
> > > > > Best regards,<br>
> > > > > wendi<br>
> > > > ><br>
> > > > > On Tue, Jan 17, 2017 at 6:45 PM, Andrew Friedman <<a href="mailto:andfried@gmail.com">andfried@gmail.com</a>> wrote:<br>
> > > > > Hi Wendi,<br>
> > > > > 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>
> > > > > -Andrew<br>
> > > > ><br>
> > > > > > On Jan 17, 2017, at 10:29 AM, wendi harjupa <<a href="mailto:wendiharjupa@gmail.com">wendiharjupa@gmail.com</a>> wrote:<br>
> > > > > ><br>
> > > > > > Dear All<br>
> > > > > ><br>
> > > > > > I got problem in overlaying two type of data.<br>
> > > > > > As I have attached the plot of data in this email.<br>
> > > > > ><br>
> > > > > > First data is rain data with resolution is about 0.002x0.003 deg (lon ,lat).<br>
> > > > > > The rain data is indicated by rainbow color.<br>
> > > > > > and second data is cloud data with resolution is 0.1x0.1 deg (lon, lat)<br>
> > > > > > The cloud data is indicated by gray color<br>
> > > > > ><br>
> > > > > > The plot looks not good. Is there any tools in grads to smooth the cloud data?<br>
> > > > > > Then, we can interpolate cloud data to become 0.005 or 0.001, and plot data will be more nice.<br>
> > > > > ><br>
> > > > > > Thank you very much for your attention and help.<br>
> > > > > ><br>
> > > > > > Best regards,<br>
> > > > > > wendi<br>
> > > > > ><br>
> > > > > > <large skale.png>____________________<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>
> > > > ><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>
> > > > ><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>
> > > ><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>
> > > ><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>
> > ><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>
> > ><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>
> ><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>
> ><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>
><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>
><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>
<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>