[gradsusr] area

[Jeff Duda]

Fernanda,
The best I can come up with is to define a new field that satisfies the
condition, but use the const function to set all values as 1.  Then take
the sum of that field across a larger area that contains your new field.
For your example, you could define a new field as the binary of temperature
within the 10 C contour (you'd have to clarify whether above or below 10 C)
as

'define above = maskout(temperature,10)'
'define aboveone = const(above,1)'

or nest the above formulas into one step: 'define aboveone =
const(maskout(temperature,10),1)'

Then: 'd asumg(aboveone,x=x1,x=x2,y=y1,y=y2)', where you'd have to insert
specific values for x1, x2, y1, y2.

If your temperature field contained values both above and below a closed
10C contour, you'd have to define a second field for the other direction,
then sum the two final values together.

Jeff Duda

On Tue, Aug 28, 2012 at 9:28 AM, Fernanda Casagrande <fenc10 at hotmail.com>wrote:

>  Hello,
> Does anyone know how to calculate an area (Kmˆ2) inside a contour
> interval? example, the area inside the  10oC contour?
> best,
> Fernanda
>
>
> > Date: Tue, 28 Aug 2012 15:06:35 +0100
> > From: bernd.becker at metoffice.gov.uk
> > To: gradsusr at gradsusr.org
> > Subject: Re: [gradsusr] Dimension Environment
> >
> >
> > That is your nearest grid point.
> > Do you think you create extra information by interpolation?
> >
> > I doubt it.
> >
> > --
> > Bernd Becker Climate Impacts Product Development
> > Met Office Hadley Centre FitzRoy Road Exeter Devon EX1 3PB United
> > Kingdom
> > Tel.: +44 (0) 1392 884511 Fax: +44 (0)870 900 5050
> > E-mail:bernd.becker at metoffice.gov.uk -
> > http://www.metoffice.gov.uk/research/people/bernd-becker
> >
> > -----Original Message-----
> > From: gradsusr-bounces at gradsusr.org
> > [mailto:gradsusr-bounces at gradsusr.org] On Behalf Of one weather
> > Sent: 28 August 2012 14:49
> > To: gradsusr at gradsusr.org
> > Subject: [gradsusr] Dimension Environment
> >
> > Hello, everyone.
> >
> > I'm working with forecast from GFS 0.5x0.5 Degree (I'm working with
> > GrADS 2.0.a9), but I've got some problems with environment dimensions.
> > When I write
> >
> > ga-> set lat 45.7
> >
> > GrADS replay:
> >
> > LAT set to 46 46
> >
> > Someone knows what should I do to get the forecast at lat 45.7? There
> > are any posibility to interpolate data?
> >
> > Thank you so much!!
> > _______________________________________________
> > 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
>
> _______________________________________________
> 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/20120828/cd46ff26/attachment-0003.html