[gradsusr] function

wendi harjupa wendiharjupa at gmail.com
Fri Jul 13 06:07:41 EDT 2012


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>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>> 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."
>
>
> _______________________________________________
> 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/20120713/52260f93/attachment-0003.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: log.gif
Type: image/gif
Size: 23514 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20120713/52260f93/attachment-0006.gif 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: without_log.gif
Type: image/gif
Size: 15030 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20120713/52260f93/attachment-0007.gif 


More information about the gradsusr mailing list