[gradsusr] function
Charles Seman
charles.seman at noaa.gov
Thu Jul 12 13:53:46 EDT 2012
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>> 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>
> > 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 <tel:%2B81-080-4268-6676>
> email : 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
> 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."
More information about the gradsusr
mailing list