Resolution to a question regarding areas on a globe

[Karen Bice <kbice@xxxxxxxx>]

My thanks to Steve Hankin and a few others who wrote to tell me how to calculate
the area that has a particular value in a global dataset. The method for doing
this is descibed below.

Cheers.
K. Bice

Steve Hankin wrote:

> HI Karen,
>
> Did you get an answer to this?
>
> Basically, you need only define a variable
>
> LET one_1 = IF my_var EQ 1 THEN 1   ! else undefined
> LET one_2 = IF my_var EQ 2 THEN 1   ! else undefined
>
> LET area_1 = one_1[x=@din,y=@din]
> LET area_2 = one_2[x=@din,y=@din]
>
> As long as the X and Y axes are in recognized units of long/lat, then the 2D
> integration "x=@din,y=@din" will yield a result in m^2
>
> Hope this helps - steve
>
> P.S. If this leads you to a solution please consider posting a brief summary
> to ferret_users to close the question.
>
> ================
>
> Karen Bice wrote:
>
> > This must be something that is easy to do, so I apologize for needing to
> > ask:
> >
> > I have a gridded (lat, long) global data set of values (integers 1-25).
> > I want to calculate the TOTAL AREA (km^2) that has a value of 1, the
> > TOTAL AREA that has a value of 2, u.s.w.
> >
> > This sounds like something for which there must be a Ferret function or
> > script. How can I do this? I am using v. 4.9.2 on an IRIX workstation.
> >
> > Thanks. Karen
> >
> > --
> > Karen L. Bice, Assistant Scientist              Mail Stop #23
> > Dept. of Geology and Geophysics                 phone: (508) 289-3320
> > Woods Hole Oceanographic Institution            FAX:   (508) 457-2187
> > Woods Hole, MA          USA 02543-1541
>
> --
>
>                 |  NOAA/PMEL               |  ph. (206) 526-6080
> Steve Hankin    |  7600 Sand Point Way NE  |  FAX (206) 526-6744
>                 |  Seattle, WA 98115-0070  |  hankin@pmel.noaa.gov