[gradsusr] gradsusr Digest, Vol 44, Issue 25

Jeff Duda jeffduda319 at gmail.com
Thu Oct 10 01:14:10 EDT 2013


WSM3 is single moment microphysics.  It doesn't predict number
concentration for graupel.  I'm not even sure WDM3 does since there are
only 3 classes.  Try WDM5, WDM6, Morrison, or M-Y instead if you really
want to get the graupel number concentration out.

Jeff


On Thu, Oct 10, 2013 at 12:07 AM, Geeta Geeta <geeta124 at hotmail.com> wrote:

> Thanks Jeff.
>
> here is my namelist file. I am using mp_phy=3 This is WRF single Moment 3
> class scheme that has ICe, snow and graupel processes!!!!.
>
> Here is my namelist.input file...........
> &time_control
>  run_minutes            = 3240 ,
>  start_year             = 2013, 2013, 2010,
>  start_month            = 05, 05, 03,
>  start_day              = 21, 21, 04,
>  start_hour             = 00, 00, 00,
>  end_year               = 2013, 2013, 2010,
>  end_month              = 05, 05, 03,
>  end_day                = 23, 23, 05,
>  end_hour               = 06, 06, 12,
>  interval_seconds       = 21600,
>  input_from_file        = .true.,.true.,.true.,
>  history_interval       = 60, 60, 60,
>  frames_per_outfile     = 1000, 1000, 1,
>  restart                = .false.,
>  restart_interval       = 5000,
>  io_form_history        = 2
>  io_form_restart        = 2
>  io_form_input          = 2
>  io_form_boundary       = 2
>  debug_level            = 1
> /
>
>
> &domains
>  time_step              = 45 ,
>  time_step_fract_num    = 0,
>  time_step_fract_den    = 1,
>  max_dom                = 1,
>  s_we                   = 1, 1,
>  e_we                   = 148, 259, ,
>  s_sn                   = 1, 1, 1,
>  e_sn                   = 160, 223, ,
>  s_vert                 = 1,     1,     1,
>  e_vert                 = 38,    38,    ,
>  num_metgrid_levels     = 27,
>  dx                     = 09000, 3000, ,
>  dy                     = 09000, 3000, ,
>  grid_id                = 1,     2,     3,
>  parent_id              = 1,     1,     2,
>  i_parent_start         = 1,     37,    ,
>  j_parent_start         = 1,     37,    ,
>  parent_grid_ratio      = 1,     3,     3,
>  parent_time_step_ratio = 1,     3,     3,
>  p_top_requested        = 5000,
>  eta_levels             = 1.000, 0.994, 0.983, 0.968, 0.950, 0.930, 0.908,
> 0.882, 0.853, 0.821, 0.788, 0.752, 0.715, 0.677, 0.637, 0.597, 0.557,
> 0.517, 0.477, 0.438, 0.401, 0.365, 0.332, 0.302, 0.274, 0.248, 0.224,
> 0.201, 0.179, 0.158, 0.138, 0.118, 0.098, 0.078, 0.058, 0.038, 0.018, 0.000,
>  feedback               = 1,
>  smooth_option          = 1
> /
>
>
> &physics
>  mp_physics             = 3,     3,     3,
>  ra_lw_physics          = 1,     1,     1,
>  ra_sw_physics          = 1,     1,     1,
>  radt                   = 30,    30,    30,
>  sf_sfclay_physics      = 1,     1,     1,
>  sf_surface_physics     = 2,     2,     2,
>  bl_pbl_physics         = 1,     1,     1,
>  bldt                   = 0,     0,     0,
>  cu_physics             = 0,     3,     0,
>  cudt                   = 5,     5,     5,
>  isfflx                 = 1,
>  ifsnow                 = 0,
>  icloud                 = 1,
>  surface_input_source   = 1,
>  num_soil_layers        = 4,
>  ucmcall                = 0,
>  mp_zero_out            = 0,
>
>
> can you Pls suggest.
>
>
> geeta
>
> > From: gradsusr-request at gradsusr.org
> > Subject: gradsusr Digest, Vol 44, Issue 25
> > To: gradsusr at gradsusr.org
> > Date: Wed, 9 Oct 2013 06:36:32 -0400
> >
> > Send gradsusr mailing list submissions to
> > gradsusr at gradsusr.org
> >
> > To subscribe or unsubscribe via the World Wide Web, visit
> > http://gradsusr.org/mailman/listinfo/gradsusr
> > or, via email, send a message with subject or body 'help' to
> > gradsusr-request at gradsusr.org
> >
> > You can reach the person managing the list at
> > gradsusr-owner at gradsusr.org
> >
> > When replying, please edit your Subject line so it is more specific
> > than "Re: Contents of gradsusr digest..."
> >
> >
> > Today's Topics:
> >
> > 1. Re: unable to get GRAUPELNC from WRF (Jeff Duda)
> > 2. Re: Cloudiness cut (Jeff Duda)
> > 3. Select by Attribute (Ryan Danielson)
> > 4. How to sum variables?? (uzair manan)
> > 5. request (Rajesh Maurya)
> > 6. plot precipitation from a 3hourly data to hourly
> > (Lyndon Mark Olaguera)
> > 7. Re: How to calculate "quartile" like value usinggradsscript?
> > (Brandon Brown)
> >
> >
> > ----------------------------------------------------------------------
> >
> > Message: 1
> > Date: Tue, 8 Oct 2013 14:53:40 -0500
> > From: Jeff Duda <jeffduda319 at gmail.com>
> > Subject: Re: [gradsusr] unable to get GRAUPELNC from WRF
> > To: GrADS Users Forum <gradsusr at gradsusr.org>
> > Message-ID:
> > <CAAig09A18wmd0SA8EjunFvCZW79AFBZzuxXGOoj04tdFT=qzAg at mail.gmail.com>
> > Content-Type: text/plain; charset="iso-8859-1"
> >
> > It's generated by the microphysics scheme. You have to be using a
> > microphysics scheme that outputs number concentration for graupel. As far
> > as I know, only the M-Y scheme, Morrison scheme, and possibly the NSSL
> > 2-moment scheme actually do that (as of version 3.4.1).
> >
> > Jeff Duda
> >
> >
> > On Tue, Oct 8, 2013 at 5:04 AM, Geeta Geeta <geeta124 at hotmail.com>
> wrote:
> >
> > > Hi,
> > > I am getting GRAUPELNC as constant field 0 from WRF output.
> > > How is the GRAUPELNC generated by the model???.
> > >
> > > ga-> c
> > > ga-> reinit
> > > No hardcopy metafile open
> > > All files closed; all defined objects released
> > > All GrADS attributes have been reinitialized
> > > ga-> open wrf_d01_2013-04-01_00.ctl
> > > Scanning description file: wrf_d01_2013-04-01_00.ctl
> > > Data file wrf_d01_2013-04-01_00.dat is open as file 1
> > > LON set to 72.0905 83.9095
> > > LAT set to 7.21236 19.6264
> > > LEV set to 1000 1000
> > > Time values set: 2013:4:1:0 2013:4:1:0
> > > ga-> set t 6
> > > Time values set: 2013:4:1:5 2013:4:1:5
> > > ga-> d GRAUPEL
> > > Syntax Error: Invalid Operand
> > > 'graupel' not a variable or function name
> > > Error ocurred at column 1
> > > DISPLAY error: Invalid expression
> > > Expression = graupel
> > > ga-> d GRAUPELNC
> > > Constant field. Value = 0
> > > ga->
> > >
> > >
> > > Pls suggest.
> > > Geeta
> > >
> > > geeta
> > >
> > > _______________________________________________
> > > gradsusr mailing list
> > > gradsusr at gradsusr.org
> > > http://gradsusr.org/mailman/listinfo/gradsusr
> > >
> > >
> >
> >
> > --
> > Jeff Duda
> > Graduate research assistant
> > University of Oklahoma School of Meteorology
> > Center for Analysis and Prediction of Storms
> > -------------- next part --------------
> > An HTML attachment was scrubbed...
> > URL:
> http://gradsusr.org/pipermail/gradsusr/attachments/20131008/b59d6e24/attachment-0001.html
> >
> > ------------------------------
> >
> > Message: 2
> > Date: Tue, 8 Oct 2013 14:54:40 -0500
> > From: Jeff Duda <jeffduda319 at gmail.com>
> > Subject: Re: [gradsusr] Cloudiness cut
> > To: GrADS Users Forum <gradsusr at gradsusr.org>
> > Message-ID:
> > <CAAig09A7E9f35Xp1gCRzDXkF-wALwFsNqsvohTdXF9j68LHMzg at mail.gmail.com>
> > Content-Type: text/plain; charset="iso-8859-1"
> >
> > I would suggest obtaining the Unified Post Processor, available on the
> WRF
> > site. It can diagnose cloud cover percentage from WRF output. If you are
> > using a different data set, then I don't know.
> >
> > Jeff Duda
> >
> >
> > On Tue, Oct 8, 2013 at 5:24 AM, avi kojokro <kojokroavi at gmail.com>
> wrote:
> >
> > > Can anyone help?
> > >
> > >
> > > 2013/10/6 avi kojokro <kojokroavi at gmail.com>
> > >
> > >> Hello to all users of Grads.
> > >>
> > >> Does anyone know a script that displays a percentage cloudiness data
> section
> > >> of a single point of geographical latitude and longitude, or graph
> paints
> > >> the clouds at the height of the layer in which the clouds, as in the
> > >> attached example?
> > >>
> > >> thank you,
> > >>
> > >> Ragards,
> > >>
> > >> Avi.
> > >>
> > >>
> > > _______________________________________________
> > > gradsusr mailing list
> > > gradsusr at gradsusr.org
> > > http://gradsusr.org/mailman/listinfo/gradsusr
> > >
> > >
> >
> >
> > --
> > Jeff Duda
> > Graduate research assistant
> > University of Oklahoma School of Meteorology
> > Center for Analysis and Prediction of Storms
> > -------------- next part --------------
> > An HTML attachment was scrubbed...
> > URL:
> http://gradsusr.org/pipermail/gradsusr/attachments/20131008/ad5e50f2/attachment-0001.html
> >
> > ------------------------------
> >
> > Message: 3
> > Date: Tue, 8 Oct 2013 15:15:32 -0700
> > From: Ryan Danielson <rydan418 at gmail.com>
> > Subject: [gradsusr] Select by Attribute
> > To: GrADS Users Forum <gradsusr at gradsusr.org>
> > Message-ID:
> > <CAATaBjtAD+akXUt6K3Nzbcn8T_VQAHxTS0fuqosgdNEZ=+eKfQ at mail.gmail.com>
> > Content-Type: text/plain; charset="iso-8859-1"
> >
> > Hi,
> >
> > I am running GrADS v1.9.0-rc1. I am running an IsoGSM model and am
> > currently trying to select data by attributes. I have isotope points
> > graphed on a scatter plot and most of it is trending linearly, similar to
> > the line y=8x+10. I would like to set a range of values above and below
> > this line, and any points that fall outside this range I would like to
> > select and determine their geographic location. Is there an available
> > command or set of commands for this? Thanks.
> >
> > Ryan
> > -------------- next part --------------
> > An HTML attachment was scrubbed...
> > URL:
> http://gradsusr.org/pipermail/gradsusr/attachments/20131008/3b95380a/attachment-0001.html
> >
> > ------------------------------
> >
> > Message: 4
> > Date: Wed, 9 Oct 2013 11:10:43 +0500
> > From: uzair manan <uzair_july at hotmail.com>
> > Subject: [gradsusr] How to sum variables??
> > To: grads <gradsusr at gradsusr.org>
> > Message-ID: <COL125-W341506AA5C774648510AFDEF1D0 at phx.gbl>
> > Content-Type: text/plain; charset="iso-8859-1"
> >
> > Hi everyone
> >
> > I've a netcdf output file from a regional model.. It has a precipiation
> (pr) variable at four different intervals during a day i.e
> >
> > Time 00, 06, 12, 18 for Day 1
> >
> > I've to set the time in grads as:-
> >
> > ga-> set t 1 4
> >
> > in order to set the environment from Time 00 to 18....
> >
> > How can we add the variable precipitation (pr) in above four intervals
> by a single command?
> >
> >
> >
> >
> >
> >
> >
> > Best regards
> >
> > Uzair Mannan
> > Cell: +923323489950
> > NED University of Engineering & Technology
> > Karachi, Pakistan
> >
> >
> > -------------- next part --------------
> > An HTML attachment was scrubbed...
> > URL:
> http://gradsusr.org/pipermail/gradsusr/attachments/20131009/e9f1278d/attachment-0001.html
> >
> > ------------------------------
> >
> > Message: 5
> > Date: Wed, 9 Oct 2013 13:38:07 +0530
> > From: Rajesh Maurya <rksmauryagpbhu at gmail.com>
> > Subject: [gradsusr] request
> > To: gradsusr at gradsusr.org
> > Message-ID:
> > <CAGn9A+JzY4Y=0_GCEVca8AQWEjgH6dsu1DP02p2ktQ-K+aScVg at mail.gmail.com>
> > Content-Type: text/plain; charset="iso-8859-1"
> >
> > thanks
> >
> > --
> > Rajesh Kumar Maurya
> > M.Sc(Tech)- Geophysics ; Junior Research Fellow ;Research Scholar
> > Department of Geophysics
> > Banaras Hindu University
> > Varanasi 221005
> > -------------- next part --------------
> > An HTML attachment was scrubbed...
> > URL:
> http://gradsusr.org/pipermail/gradsusr/attachments/20131009/dbf8f0c0/attachment-0001.html
> >
> > ------------------------------
> >
> > Message: 6
> > Date: Wed, 9 Oct 2013 16:33:08 +0800
> > From: Lyndon Mark Olaguera <olagueralyndonmark429 at gmail.com>
> > Subject: [gradsusr] plot precipitation from a 3hourly data to hourly
> > To: Amulya Chevuturi <amulya.chevuturi at gmail.com>, GrADS Users Forum
> > <gradsusr at gradsusr.org>
> > Message-ID:
> > <CA+V+8dWk5veoOZDd9d_Uw8be0WRxbfz3wusXisXgc7WvZGtsig at mail.gmail.com>
> > Content-Type: text/plain; charset="iso-8859-1"
> >
> > hi all,
> >
> > Im having problem in plotting the precipitation from TRMM. its a three
> > hourly data but I'd like to output an hourly precipitation. how do i do
> > this?
> >
> > thanks for any help.
> >
> > *Lyndon Mark P. Olaguera*
> > -------------- next part --------------
> > An HTML attachment was scrubbed...
> > URL:
> http://gradsusr.org/pipermail/gradsusr/attachments/20131009/11cfc81d/attachment-0001.html
> >
> > ------------------------------
> >
> > Message: 7
> > Date: Wed, 9 Oct 2013 18:36:32 +0800
> > From: Brandon Brown <whicily1988 at gmail.com>
> > Subject: Re: [gradsusr] How to calculate "quartile" like value
> > usinggradsscript?
> > To: GrADS Users Forum <gradsusr at gradsusr.org>
> > Message-ID:
> > <CAOtqrS3XTQxqO0wzM80nL7GMt1+_TAsSmpAGoMRrjd-nyym8BA at mail.gmail.com>
> > Content-Type: text/plain; charset="iso-8859-1"
> >
> > Hi, Dave.
> >
> > Thanks very much for showing me your code, It helps me learning GrADS
> > Script.
> >
> > I had used python interface to grads to implement this, with a few lines
> of
> > python code, very easy and straightforward.
> >
> > Brandon.
> >
> >
> > On Tue, Oct 1, 2013 at 3:28 AM, Ryglicki, David CIV FNMOC <
> > david.ryglicki at navy.mil> wrote:
> >
> > > And then I went ahead and attached my script, anyway.
> > >
> > > Anyway, Brandon, the key trick I learned in the past for plume plots
> is a
> > > clunky work-around by defining the percentile variables in t and y, but
> > > setting x-size as t-size, which you'll see in my define statements.
> Later
> > > on, after a bubble sort, you can set defval, chunking through the time
> > > variable, and that'll eventually get you enough information to plot the
> > > percentiles as a faux x-variable.
> > >
> > > It's extremely clunky and convoluted, but I got it to work to make
> plume
> > > plot percentiles.
> > >
> > > David Ryglicki, Ph.D.
> > > N3 Models
> > > Fleet Numerical Meteorology and Oceanography Center
> > > Room 702, Room 146
> > > 831|656|4519
> > >
> > >
> > >
> > > -----Original Message-----
> > > From: gradsusr-bounces at gradsusr.org on behalf of Ryglicki, David CIV
> FNMOC
> > > Sent: Mon 9/30/2013 11:35 AM
> > > To: GrADS Users Forum
> > > Subject: Re: [gradsusr] How to calculate "quartile" like value
> > > usinggradsscript?
> > >
> > > Hi, Brandon.
> > >
> > > Are you trying to make a plume plot or 2-D spatial plots?
> > >
> > > ~Dave
> > >
> > > -----Original Message-----
> > > From: gradsusr-bounces at gradsusr.org on behalf of Brandon Brown
> > > Sent: Mon 9/30/2013 4:09 AM
> > > To: GrADS Users Forum
> > > Subject: Re: [gradsusr] How to calculate "quartile" like value using
> > > gradsscript?
> > >
> > > I have known how to this.
> > >
> > > Thanks.
> > >
> > > On 9/30/13, Brandon Brown <whicily1988 at gmail.com> wrote:
> > > > Hi Jennifer
> > > >
> > > > Thanks for pointing out this, it saves my time.
> > > >
> > > > I am planning to use pygrads instead, but pygrads does not support
> > > > ensemble, I can't import result back to grads.
> > > >
> > > > After some reseach , it seems that I can write binary data with vars
> > > > and dims, but all samples use fortran witch I am not familiar with,
> ;(
> > > > , is there some specs to do this, thus I can python.
> > > >
> > > > Great thanks.
> > > >
> > > >
> > > > On 9/30/13, Jennifer Adams <jma at cola.iges.org> wrote:
> > > >> Unfortunately, GrADS does not have sorting built in (yet), so there
> is
> > > no
> > > >> straightforward way to do the quartile-like calculations you are
> asking
> > > >> for.
> > > >> Perhaps other members of the forum can offer their workarounds.
> > > >> --Jennifer
> > > >>
> > > >> On Sep 26, 2013, at 3:00 AM, Brandon Brown wrote:
> > > >>
> > > >>> Hi,
> > > >>>
> > > >>> I am newbie, now I want to draw a boxplot, but I don't now how to
> > > >>> calculate the values which ares quartile like.
> > > >>>
> > > >>> As you know, quartile calculate q1(25%), q2(50%), q3(75%), and now
> I
> > > >>> want
> > > >>> to calculate q1(15%), q2(25%), q3(50%), q4(75%), q5(85%).
> > > >>>
> > > >>> Could any one help me, thanks.
> > > >>>
> > > >>> Brown.
> > > >>> _______________________________________________
> > > >>> gradsusr mailing list
> > > >>> gradsusr at gradsusr.org
> > > >>> http://gradsusr.org/mailman/listinfo/gradsusr
> > > >>
> > > >> --
> > > >> Jennifer M. Adams
> > > >> Center for Ocean-Land-Atmosphere Studies (COLA)
> > > >> 111 Research Hall, Mail Stop 2B3
> > > >> George Mason University
> > > >> 4400 University Drive
> > > >> Fairfax, VA 22030
> > > >>
> > > >>
> > > >>
> > > >>
> > > >>
> > > >>
> > > >
> > > _______________________________________________
> > > gradsusr mailing list
> > > gradsusr at gradsusr.org
> > > http://gradsusr.org/mailman/listinfo/gradsusr
> > >
> > >
> > >
> > >
> > >
> > > _______________________________________________
> > > gradsusr mailing list
> > > gradsusr at gradsusr.org
> > > http://gradsusr.org/mailman/listinfo/gradsusr
> > >
> > -------------- next part --------------
> > An HTML attachment was scrubbed...
> > URL:
> http://gradsusr.org/pipermail/gradsusr/attachments/20131009/116ef699/attachment.html
> >
> > ------------------------------
> >
> > _______________________________________________
> > gradsusr mailing list
> > gradsusr at gradsusr.org
> > http://gradsusr.org/mailman/listinfo/gradsusr
> >
> >
> > End of gradsusr Digest, Vol 44, Issue 25
> > ****************************************
>



-- 
Jeff Duda
Graduate research assistant
University of Oklahoma School of Meteorology
Center for Analysis and Prediction of Storms
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20131010/6b347669/attachment-0001.html 


More information about the gradsusr mailing list