[gradsusr] function
wendi harjupa
wendiharjupa at gmail.com
Thu Jul 12 05:31:32 EDT 2012
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>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
>> > 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
>
>
--
-------------
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/20120712/f4409909/attachment-0003.html
More information about the gradsusr
mailing list