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