<div dir="ltr">Dear Mr. Andrew<div><br></div><div>Thank you very much for your instructions.</div><div><br></div><div>I am very sorry for asking you many questions.</div><div><br></div><div>But, I got error message as below :</div><div><br></div><div><div>Data Request Warning: Request is completely outside file limits</div><div> Entire grid contents are set to missing data </div><div> Grid limits of file: X = 1 280 Y = 1 240 Z = 1 1 T = 1 1 E = 1 1 </div><div> Grid limits of request: X = 134 160 Y = 107 128 Z = 1 1 T = -156 -156 E =1 1</div></div><div><br></div><div>Would you mind to explain me the error meaning?</div><div><br></div><div>Thank you very much.</div><div><br></div><div>Best regards,</div><div>wendi</div><div><br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Jan 24, 2017 at 12:33 AM, 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 Wendy,<br>
<br>
I would place the command at the end of your script.<br>
<br>
'define try=(rr, albedo.2, aave)’<br>
<br>
Do you get an error? Based on your script, I’m wondering if you might be closing your files too soon.<br>
Assuming no error message, you can go to a specific lat/lon point where you know there is data in the original file and type<br>
'd try’<br>
<br>
You can also try ‘set gxout grid’ to look at the values in specific gridboxes.<br>
<br>
Hope this helps.<br>
<span class="HOEnZb"><font color="#888888">Andrew<br>
</font></span><div class="HOEnZb"><div class="h5"><br>
> On Jan 23, 2017, at 1:20 PM, 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 excellent supervision<br>
><br>
> Below is my simple script<br>
><br>
> 'open ./Rain.ctl'<br>
> 'set lat 33.5 35.5'<br>
> 'set lon 134.3 136.6'<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>
> 'set lat 33.5 35.5'<br>
> 'set lon 134.3 136.6'<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>
> my control file as below<br>
> 1. Rain.ctl<br>
><br>
> DSET ~/data1<br>
> OPTIONS yrev<br>
> UNDEF -1.E35<br>
> XDEF 800 LINEAR 134.375 0.003125<br>
> YDEF 960 LINEAR 33.5833 0.00208333<br>
> ZDEF 1 LINEAR 1 1<br>
> TDEF 120 LINEAR 08Aug2015 1hr<br>
> VARS 1<br>
> rr 0 98 rairate [mm/hr]<br>
> ENDVARS<br>
><br>
> 2. Him.ctl<br>
> DSET ${file_in_vor2}<br>
> UNDEF -1.0<br>
> OPTIONS yrev<br>
> XDEF 280 LINEAR 121.05 0.1<br>
> YDEF 240 LINEAR 22.95 0.1<br>
> ZDEF 1 LINEAR 1 1<br>
> TDEF 1 LINEAR 13:05Z08AUG2015 5mn<br>
> VARS 1<br>
> ALBEDO 0 1 reflectivity<br>
> ENDVARS<br>
><br>
> My next question where should I locate my lterp command??<br>
> and how it would be like?<br>
> I have tried this command below, it doesnor work<br>
><br>
> 'define try=(rr, albedo.2, aave)'<br>
> 'd try'<br>
><br>
> Thank you very much for your help. I am very sorry for ask many questions.<br>
><br>
> Best regards,<br>
> wendi<br>
><br>
><br>
><br>
><br>
><br>
><br>
> On Mon, Jan 23, 2017 at 6:48 PM, Andrew Friedman <<a href="mailto:andfried@gmail.com">andfried@gmail.com</a>> wrote:<br>
> 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>
><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>
><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>