<html>
  <head>
    <meta content="text/html; charset=ISO-8859-1"
      http-equiv="Content-Type">
  </head>
  <body bgcolor="#FFFFFF" text="#000000">
    <font size="+2">hi Wendi, when I have problems like this, I break up
      the expressions and then take a look at each part, e.g.<br>
      <br>
      set t 1 2<br>
      q dim # check the dim environment<br>
      # do a 'd' after each define to see if it makes sense...<br>
      # ref is a 3-D variable?&nbsp; x,y,t?<br>
      r1=const(ref,10.0,-a)<br>
      r2=ref*0.10<br>
      p1=pow(r1,r2)<br>
      p12=p1/200.0<br>
      p22=const(ref,5.0/8.0,-a)<br>
      R=pow(p12,p22)<br>
      set t 1<br>
      precip=sum(R,t=1,t=2)<br>
      d precip<br>
      <br>
      ciao ciao Mike<br>
      <br>
      <br>
    </font>
    <div class="moz-cite-prefix">On 7/12/12 17:45 , Charles Seman wrote:<br>
    </div>
    <blockquote cite="mid:4FFF0D26.3050504@noaa.gov" type="cite">
      <pre wrap="">Wendi,
era
 From a GrADS perspective, the commands you list below look OK to me.

Chuck

On 07/12/2012 05:03 AM, wendi harjupa wrote:
</pre>
      <blockquote type="cite">
        <pre wrap="">Dear Mr Chuck

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' ;* see discussion below
'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.

Thank you very much.

Best regards,
wendi




On Thu, Jul 12, 2012 at 6:32 AM, Charles Seman &lt;<a class="moz-txt-link-abbreviated" href="mailto:charles.seman@noaa.gov">charles.seman@noaa.gov</a>
<a class="moz-txt-link-rfc2396E" href="mailto:charles.seman@noaa.gov">&lt;mailto:charles.seman@noaa.gov&gt;</a>&gt; 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
    <a class="moz-txt-link-freetext" href="http://grads.iges.org/grads/gadoc/variable.html#new">http://grads.iges.org/grads/gadoc/variable.html#new</a>); 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:
     &gt; Dear Grads users,
     &gt;
     &gt; I have problem,
     &gt; below is  my descriptor file
     &gt;
     &gt; DSET
    /misc/home/wendi/data1/radar/kototabang/XDR/data/20040410/%d2%h2.dat
     &gt; TITLE  RADAR
     &gt; OPTIONS  TEMPLATE LITTLE_ENDIAN
     &gt; UNDEF  -1000
     &gt; XDEF    321 LINEAR    99.688120996    0.004492369
     &gt; YDEF    321 LINEAR    -1.083708891    0.004522556
     &gt; ZDEF    39  LINEAR   1.0 0.5
     &gt; TDEF    15 LINEAR 21:00Z10apr2004  4mn
     &gt; VARS    2
     &gt; ref 39  99  Reflectivity(dBZ)
     &gt; vel 39  99  Doppler Velocity(m/s)
     &gt; ENDVARS
     &gt;
     &gt; As you see time is 15,
     &gt;
     &gt; I want to convert the ref (dbz) to Rain (mm/hr)
     &gt; by using Z-R correlation formula, as below :
     &gt;
     &gt; 'pow((pow(10, R/10))/200,5/8)',
     &gt;
     &gt; Then I want to accumulate the rainrate for 1 hour,
     &gt; I use this formula :
     &gt;
     &gt; 'sum(ref, t=1, t=15)'
     &gt;
     &gt; But I got strange result.
     &gt;
     &gt; Here my complete gs file,
     &gt;
     &gt; 'reinit'
     &gt; 'set grads off'
     &gt; 'set mproj latlon'
     &gt; 'set mpdset hires'
     &gt; 'set gxout shaded'
     &gt; 'open RHI.ctl'
     &gt; 'set lon 100 100.8'
     &gt; 'set lat -0.45'
     &gt; 'set Z 1 19'
     &gt; 'set parea 1 7.5 2.8 8.3'
     &gt; 'set clevs 0 50 100 150 200 250 300 350 400 450 500 550 600'
     &gt; 'set ccols 9 14 4 11 5 13 3 10 7 12 8 2 6'
     &gt; 'R=pow((pow(10, ref/10))/200,5/8)'
     &gt; 'd sum(R, t=1, t=15)'
     &gt; 'run /misc/home/wendi/grads-2.0.1/cbarn 0.7 0.7 7.9 5.2'
     &gt; 'draw string  7.7 7.8 [dBz]'
     &gt;
     &gt; Would you help me, about this,
     &gt; please let me know, which part I got wrong
     &gt;
     &gt; Thank you very much!
     &gt;
     &gt; Best regards,
     &gt; wendi
     &gt;
     &gt;
     &gt; --
     &gt; -------------
     &gt; Wendi Harjupa. ST
     &gt; ShimaneUniversityGraduateSchool ofEngineering Sciences
     &gt; Department ofElectronic ControlSystems
     &gt; RemoteSensingLaboratory
     &gt; s119447
     &gt; Cellphone : +81-080-4268-6676
     &gt; email : <a class="moz-txt-link-abbreviated" href="mailto:wendi@rslab.riko.shimane-u.ac.jp">wendi@rslab.riko.shimane-u.ac.jp</a>
    <a class="moz-txt-link-rfc2396E" href="mailto:wendi@rslab.riko.shimane-u.ac.jp">&lt;mailto:wendi@rslab.riko.shimane-u.ac.jp&gt;</a>
     &gt; &lt;<a class="moz-txt-link-freetext" href="mailto:wendi@rslab.riko.shimane-u.ac.jp">mailto:wendi@rslab.riko.shimane-u.ac.jp</a>
    <a class="moz-txt-link-rfc2396E" href="mailto:wendi@rslab.riko.shimane-u.ac.jp">&lt;mailto:wendi@rslab.riko.shimane-u.ac.jp&gt;</a>&gt;
     &gt;
     &gt;
     &gt;
     &gt; _______________________________________________
     &gt; gradsusr mailing list
     &gt; <a class="moz-txt-link-abbreviated" href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a> <a class="moz-txt-link-rfc2396E" href="mailto:gradsusr@gradsusr.org">&lt;mailto:gradsusr@gradsusr.org&gt;</a>
     &gt; <a class="moz-txt-link-freetext" href="http://gradsusr.org/mailman/listinfo/gradsusr">http://gradsusr.org/mailman/listinfo/gradsusr</a>
     &gt;

    --

    Please note that <a class="moz-txt-link-abbreviated" href="mailto:Charles.Seman@noaa.gov">Charles.Seman@noaa.gov</a>
    <a class="moz-txt-link-rfc2396E" href="mailto:Charles.Seman@noaa.gov">&lt;mailto:Charles.Seman@noaa.gov&gt;</a> should be considered my NOAA
    email address, not <a class="moz-txt-link-abbreviated" href="mailto:cjs@gfdl.noaa.gov">cjs@gfdl.noaa.gov</a> <a class="moz-txt-link-rfc2396E" href="mailto:cjs@gfdl.noaa.gov">&lt;mailto:cjs@gfdl.noaa.gov&gt;</a>.

    ********************************************************************
       Charles Seman <a class="moz-txt-link-abbreviated" href="mailto:Charles.Seman@noaa.gov">Charles.Seman@noaa.gov</a> <a class="moz-txt-link-rfc2396E" href="mailto:Charles.Seman@noaa.gov">&lt;mailto:Charles.Seman@noaa.gov&gt;</a>
       U.S. Department of Commerce / NOAA / OAR
       Geophysical Fluid Dynamics Laboratory         voice: (609)
    452-6547 <a class="moz-txt-link-rfc2396E" href="tel:%28609%29%20452-6547">&lt;tel:%28609%29%20452-6547&gt;</a>
       201 Forrestal Road                              fax: (609)
    987-5063 <a class="moz-txt-link-rfc2396E" href="tel:%28609%29%20987-5063">&lt;tel:%28609%29%20987-5063&gt;</a>
       Princeton, NJ  08540-6649 <a class="moz-txt-link-freetext" href="http://www.gfdl.noaa.gov/~cjs/">http://www.gfdl.noaa.gov/~cjs/</a>
    ********************************************************************

    "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
    <a class="moz-txt-link-abbreviated" href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a> <a class="moz-txt-link-rfc2396E" href="mailto:gradsusr@gradsusr.org">&lt;mailto:gradsusr@gradsusr.org&gt;</a>
    <a class="moz-txt-link-freetext" href="http://gradsusr.org/mailman/listinfo/gradsusr">http://gradsusr.org/mailman/listinfo/gradsusr</a>




--
-------------
Wendi Harjupa. ST
ShimaneUniversityGraduateSchool ofEngineering Sciences
Department ofElectronic ControlSystems
RemoteSensingLaboratory
s119447
Cellphone : +81-080-4268-6676
email : <a class="moz-txt-link-abbreviated" href="mailto:wendi@rslab.riko.shimane-u.ac.jp">wendi@rslab.riko.shimane-u.ac.jp</a>
<a class="moz-txt-link-rfc2396E" href="mailto:wendi@rslab.riko.shimane-u.ac.jp">&lt;mailto:wendi@rslab.riko.shimane-u.ac.jp&gt;</a>



_______________________________________________
gradsusr mailing list
<a class="moz-txt-link-abbreviated" href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>
<a class="moz-txt-link-freetext" href="http://gradsusr.org/mailman/listinfo/gradsusr">http://gradsusr.org/mailman/listinfo/gradsusr</a>

</pre>
      </blockquote>
      <pre wrap="">
</pre>
    </blockquote>
    <br>
  </body>
</html>