Dear Mr. Chuck<div><br></div><div>Thank you very much for checking my script,</div><div><br></div><div>Here I attached 2 figures in this email. Figure by using log command (log.gif)</div><div>and figure without log command (without_log.gif)</div>
<div>Would you tell me about those figures, what should I do to make good figure</div><div>and easy to understant.</div><div><br></div><div>Thank you very much.</div><div><br></div><div>Best regards,</div><div>wendi</div>
<div><br><br><div class="gmail_quote">On Fri, Jul 13, 2012 at 2:53 AM, Charles Seman <span dir="ltr"><<a href="mailto:charles.seman@noaa.gov" target="_blank">charles.seman@noaa.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Wendi,<br>
<br>
I'm sorry but I don't know what you mean by "strange result"?<br>
<br>
Chuck<br>
<div class="im"><br>
On 07/12/2012 05:31 AM, wendi harjupa wrote:<br>
> Dear Mr Chuck<br>
><br>
> I am sorry I resend my email:<br>
><br>
> Thank you very much for your explanation,<br>
> I am sorry, I want to ask you more...<br>
><br>
> If I have time series data is 4485, so can I make the script like this:<br>
><br>
> 'set t 1 4485'<br>
> 'R=pow((pow(10, ref/10))/200,5/8)'<br>
> 'set t 1'<br>
> 'd sum(R/15, t=1, t=4485)'<br>
><br>
> R=pow((pow(10, ref/10))/200,5/8) is formula to convert from Z (dBz)<br>
> to R (mm/hr)<br>
> to get value mm/hr, I have to use R/15.<br>
><br>
><br>
> if I just want to make calculation and no need to convert to R,<br>
> can I use this script :<br>
><br>
> 'set t 1 4485'<br>
> 'set t 1'<br>
> 'R1=sum(ref, t=1, t=4485)'<br>
> 'd log(R1)'<br>
><br>
> I got strange result if I donot use log!<br>
><br>
><br>
> Thank you very much.<br>
><br>
> Best regards,<br>
> wendi<br>
><br>
><br>
><br>
><br>
> On Thu, Jul 12, 2012 at 6:32 AM, Charles Seman<br>
</div><div><div class="h5">> <<a href="mailto:charles.seman@noaa.gov">charles.seman@noaa.gov</a> <mailto:<a href="mailto:charles.seman@noaa.gov">charles.seman@noaa.gov</a>>> wrote:<br>
><br>
> Wendi,<br>
><br>
> An idea; set the time dimension to span t=1,15 for the<br>
> definition of "R":<br>
><br>
> 'set t 1 15'<br>
> 'R=pow((pow(10, ref/10))/200,5/8)'<br>
> 'set t 1' ;* see discussion below<br>
> 'd sum(R, t=1, t=15)'<br>
><br>
> 'set t 1': for displaying sum; this is OK to go from time range to a<br>
> single time level to do the sum (see "wild cards" discussion at<br>
> <a href="http://grads.iges.org/grads/gadoc/variable.html#new" target="_blank">http://grads.iges.org/grads/gadoc/variable.html#new</a>); if left as<br>
> 'set t<br>
> 1 15', I think you would get 15 time levels with same "sum"<br>
> result for<br>
> the display command... you could try it both ways to see what<br>
> you get...<br>
><br>
> Hope this helps,<br>
> Chuck<br>
><br>
> On 07/06/2012 05:53 AM, wendi harjupa wrote:<br>
> > Dear Grads users,<br>
> ><br>
> > I have problem,<br>
> > below is my descriptor file<br>
> ><br>
> > DSET<br>
> /misc/home/wendi/data1/radar/kototabang/XDR/data/20040410/%d2%h2.dat<br>
> > TITLE RADAR<br>
> > OPTIONS TEMPLATE LITTLE_ENDIAN<br>
> > UNDEF -1000<br>
> > XDEF 321 LINEAR 99.688120996 0.004492369<br>
> > YDEF 321 LINEAR -1.083708891 0.004522556<br>
> > ZDEF 39 LINEAR 1.0 0.5<br>
> > TDEF 15 LINEAR 21:00Z10apr2004 4mn<br>
> > VARS 2<br>
> > ref 39 99 Reflectivity(dBZ)<br>
> > vel 39 99 Doppler Velocity(m/s)<br>
> > ENDVARS<br>
> ><br>
> > As you see time is 15,<br>
> ><br>
> > I want to convert the ref (dbz) to Rain (mm/hr)<br>
> > by using Z-R correlation formula, as below :<br>
> ><br>
> > 'pow((pow(10, R/10))/200,5/8)',<br>
> ><br>
> > Then I want to accumulate the rainrate for 1 hour,<br>
> > I use this formula :<br>
> ><br>
> > 'sum(ref, t=1, t=15)'<br>
> ><br>
> > But I got strange result.<br>
> ><br>
> > Here my complete gs file,<br>
> ><br>
> > 'reinit'<br>
> > 'set grads off'<br>
> > 'set mproj latlon'<br>
> > 'set mpdset hires'<br>
> > 'set gxout shaded'<br>
> > 'open RHI.ctl'<br>
> > 'set lon 100 100.8'<br>
> > 'set lat -0.45'<br>
> > 'set Z 1 19'<br>
> > 'set parea 1 7.5 2.8 8.3'<br>
> > 'set clevs 0 50 100 150 200 250 300 350 400 450 500 550 600'<br>
> > 'set ccols 9 14 4 11 5 13 3 10 7 12 8 2 6'<br>
> > 'R=pow((pow(10, ref/10))/200,5/8)'<br>
> > 'd sum(R, t=1, t=15)'<br>
> > 'run /misc/home/wendi/grads-2.0.1/cbarn 0.7 0.7 7.9 5.2'<br>
> > 'draw string 7.7 7.8 [dBz]'<br>
> ><br>
> > Would you help me, about this,<br>
> > please let me know, which part I got wrong<br>
> ><br>
> > Thank you very much!<br>
> ><br>
> > Best regards,<br>
> > wendi<br>
> ><br>
> ><br>
> > --<br>
> > -------------<br>
> > Wendi Harjupa. ST<br>
> > ShimaneUniversityGraduateSchool ofEngineering Sciences<br>
> > Department ofElectronic ControlSystems<br>
> > RemoteSensingLaboratory<br>
> > s119447<br>
</div></div>> > Cellphone : <a href="tel:%2B81-080-4268-6676" value="+818042686676">+81-080-4268-6676</a> <tel:%2B81-080-4268-6676><br>
> > email : <a href="mailto:wendi@rslab.riko.shimane-u.ac.jp">wendi@rslab.riko.shimane-u.ac.jp</a><br>
> <mailto:<a href="mailto:wendi@rslab.riko.shimane-u.ac.jp">wendi@rslab.riko.shimane-u.ac.jp</a>><br>
> > <mailto:<a href="mailto:wendi@rslab.riko.shimane-u.ac.jp">wendi@rslab.riko.shimane-u.ac.jp</a><br>
<div class="im">> <mailto:<a href="mailto:wendi@rslab.riko.shimane-u.ac.jp">wendi@rslab.riko.shimane-u.ac.jp</a>>><br>
> ><br>
> ><br>
> ><br>
> > _______________________________________________<br>
> > gradsusr mailing list<br>
</div>> > <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a> <mailto:<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>><br>
<div class="im">> > <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
> ><br>
><br>
> --<br>
><br>
> Please note that <a href="mailto:Charles.Seman@noaa.gov">Charles.Seman@noaa.gov</a><br>
</div><div class="im">> <mailto:<a href="mailto:Charles.Seman@noaa.gov">Charles.Seman@noaa.gov</a>> should be considered my NOAA<br>
> email address, not <a href="mailto:cjs@gfdl.noaa.gov">cjs@gfdl.noaa.gov</a> <mailto:<a href="mailto:cjs@gfdl.noaa.gov">cjs@gfdl.noaa.gov</a>>.<br>
><br>
> ********************************************************************<br>
> Charles Seman <a href="mailto:Charles.Seman@noaa.gov">Charles.Seman@noaa.gov</a><br>
</div>> <mailto:<a href="mailto:Charles.Seman@noaa.gov">Charles.Seman@noaa.gov</a>><br>
<div class="im">> U.S. Department of Commerce / NOAA / OAR<br>
> Geophysical Fluid Dynamics Laboratory voice: (609)<br>
</div>> 452-6547 <tel:%28609%29%20452-6547><br>
<div class="im">> 201 Forrestal Road fax: (609)<br>
</div>> 987-5063 <tel:%28609%29%20987-5063><br>
<div class="im">> Princeton, NJ 08540-6649 <a href="http://www.gfdl.noaa.gov/~cjs/" target="_blank">http://www.gfdl.noaa.gov/~cjs/</a><br>
> ********************************************************************<br>
><br>
> "The contents of this message are mine personally and do not<br>
> reflect any<br>
> official or unofficial position of the United States Federal<br>
> Government,<br>
> the United States Department of Commerce, or NOAA."<br>
><br>
><br>
> _______________________________________________<br>
> gradsusr mailing list<br>
</div>> <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a> <mailto:<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>><br>
<div class="im">> <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
><br>
><br>
><br>
><br>
> --<br>
> -------------<br>
> Wendi Harjupa. ST<br>
</div><div class="im">> ShimaneUniversityGraduateSchool ofEngineering Sciences<br>
> Department ofElectronic ControlSystems<br>
> RemoteSensingLaboratory<br>
> s119447<br>
</div>> Cellphone : <a href="tel:%2B81-080-4268-6676" value="+818042686676">+81-080-4268-6676</a> <tel:%2B81-080-4268-6676><br>
> email : <a href="mailto:wendi@rslab.riko.shimane-u.ac.jp">wendi@rslab.riko.shimane-u.ac.jp</a><br>
> <mailto:<a href="mailto:wendi@rslab.riko.shimane-u.ac.jp">wendi@rslab.riko.shimane-u.ac.jp</a>><br>
<div class="im HOEnZb">><br>
><br>
><br>
><br>
> --<br>
> -------------<br>
> Wendi Harjupa. ST<br>
</div><div class="HOEnZb"><div class="h5">> ShimaneUniversityGraduateSchool ofEngineering Sciences<br>
> Department ofElectronic ControlSystems<br>
> RemoteSensingLaboratory<br>
> s119447<br>
> Cellphone : <a href="tel:%2B81-080-4268-6676" value="+818042686676">+81-080-4268-6676</a><br>
> email : <a href="mailto:wendi@rslab.riko.shimane-u.ac.jp">wendi@rslab.riko.shimane-u.ac.jp</a><br>
> <mailto:<a href="mailto:wendi@rslab.riko.shimane-u.ac.jp">wendi@rslab.riko.shimane-u.ac.jp</a>><br>
><br>
><br>
><br>
> _______________________________________________<br>
> gradsusr mailing list<br>
> <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
> <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
><br>
<br>
--<br>
<br>
Please note that <a href="mailto:Charles.Seman@noaa.gov">Charles.Seman@noaa.gov</a> should be considered my NOAA<br>
email address, not <a href="mailto:cjs@gfdl.noaa.gov">cjs@gfdl.noaa.gov</a>.<br>
<br>
********************************************************************<br>
Charles Seman <a href="mailto:Charles.Seman@noaa.gov">Charles.Seman@noaa.gov</a><br>
U.S. Department of Commerce / NOAA / OAR<br>
Geophysical Fluid Dynamics Laboratory voice: <a href="tel:%28609%29%20452-6547" value="+16094526547">(609) 452-6547</a><br>
201 Forrestal Road fax: <a href="tel:%28609%29%20987-5063" value="+16099875063">(609) 987-5063</a><br>
Princeton, NJ 08540-6649 <a href="http://www.gfdl.noaa.gov/~cjs/" target="_blank">http://www.gfdl.noaa.gov/~cjs/</a><br>
********************************************************************<br>
<br>
"The contents of this message are mine personally and do not reflect any<br>
official or unofficial position of the United States Federal Government,<br>
the United States Department of Commerce, or NOAA."<br>
<br>
<br>
_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br><span style="font-family:arial,helvetica,sans-serif;font-size:16px;background-color:rgb(255,255,255)">-------------</span><br style="font-family:arial,helvetica,sans-serif;font-size:16px;background-color:rgb(255,255,255)">
<span style="font-family:arial,helvetica,sans-serif;font-size:16px;background-color:rgb(255,255,255)">Wendi Harjupa. ST </span><br style="font-family:arial,helvetica,sans-serif;font-size:16px;background-color:rgb(255,255,255)">
<span style="color:rgb(51,51,51);font-family:arial,sans-serif;font-size:16px;background-color:rgb(245,245,245)">Shimane</span><span style="color:rgb(51,51,51);font-family:arial,sans-serif;font-size:16px;background-color:rgb(245,245,245)"> </span><span style="color:rgb(51,51,51);font-family:arial,sans-serif;font-size:16px;background-color:rgb(245,245,245)">University</span><span style="color:rgb(51,51,51);font-family:arial,sans-serif;font-size:16px;background-color:rgb(245,245,245)"> </span><span style="color:rgb(51,51,51);font-family:arial,sans-serif;font-size:16px;background-color:rgb(245,245,245)">Graduate</span><span style="color:rgb(51,51,51);font-family:arial,sans-serif;font-size:16px;background-color:rgb(245,245,245)"> </span><span style="color:rgb(51,51,51);font-family:arial,sans-serif;font-size:16px;background-color:rgb(245,245,245)">School of</span><span style="color:rgb(51,51,51);font-family:arial,sans-serif;font-size:16px;background-color:rgb(245,245,245)"> </span><span style="color:rgb(51,51,51);font-family:arial,sans-serif;font-size:16px;background-color:rgb(245,245,245)">Engineering Sciences</span><br style="color:rgb(51,51,51);font-family:arial,sans-serif;font-size:16px;background-color:rgb(245,245,245)">
<span style="color:rgb(51,51,51);font-family:arial,sans-serif;font-size:16px;background-color:rgb(245,245,245)">Department of</span><span style="color:rgb(51,51,51);font-family:arial,sans-serif;font-size:16px;background-color:rgb(245,245,245)"> </span><span style="color:rgb(51,51,51);font-family:arial,sans-serif;font-size:16px;background-color:rgb(245,245,245)">Electronic Control</span><span style="color:rgb(51,51,51);font-family:arial,sans-serif;font-size:16px;background-color:rgb(245,245,245)"> </span><span style="color:rgb(51,51,51);font-family:arial,sans-serif;font-size:16px;background-color:rgb(245,245,245)">Systems</span><br style="color:rgb(51,51,51);font-family:arial,sans-serif;font-size:16px;background-color:rgb(245,245,245)">
<span style="color:rgb(51,51,51);font-family:arial,sans-serif;font-size:16px;background-color:rgb(245,245,245)">Remote</span><span style="color:rgb(51,51,51);font-family:arial,sans-serif;font-size:16px;background-color:rgb(245,245,245)"> </span><span style="color:rgb(51,51,51);font-family:arial,sans-serif;font-size:16px;background-color:rgb(245,245,245)">Sensing</span><span style="color:rgb(51,51,51);font-family:arial,sans-serif;font-size:16px;background-color:rgb(245,245,245)"> </span><span style="color:rgb(51,51,51);font-family:arial,sans-serif;font-size:16px;background-color:rgb(245,245,245)">Laboratory</span><div>
<span style="color:rgb(51,51,51);font-family:arial,sans-serif;font-size:16px;background-color:rgb(245,245,245)">s119447</span></div><div><font color="#333333" face="arial, sans-serif" size="3">Cellphone : +81-080-4268-6676</font></div>
<div><font color="#333333" face="arial, sans-serif" size="3">email : <a href="mailto:wendi@rslab.riko.shimane-u.ac.jp" target="_blank">wendi@rslab.riko.shimane-u.ac.jp</a></font></div><br>
</div>