Re: Earth Radius for CF Grid Mappings

Dear Greg

You can include earth_radius with any of the mappings. At the start of App F
it says
"The attributes which describe the ellipsoid and prime meridian may be 
included, when applicable, with any grid mapping."
and earth_radius is one of those. If it's not specified but you need it, I
think you can assume any sensible value.

Cheers

Jonathan

On Mon, Dec 15, 2008 at 06:02:10PM -0500, Greg Rappa wrote:
> Date: Mon, 15 Dec 2008 18:02:10 -0500
> From: Greg Rappa <address@hidden>
> User-Agent: Thunderbird 2.0.0.16 (X11/20080707)
> To: address@hidden
> CC: "address@hidden" <address@hidden>,
>         "'John Caron'" <address@hidden>,
>         "Newell, Oliver" <address@hidden>
> Subject: Earth Radius for CF Grid Mappings
> 
> 
>    I recently discovered that in order to resolve the location of a grid
>    element (x/y-to-lat/long conversion), using software in the [1]General
>    Cartographic Transformation Package (GCTP), a value for the
>    Earth's radius is required.  For example, the GCTP functions
>    (lamazinvint and lamazinv) for converting x & y offsets to location
>    (latitude,longitude) of a grid mapped to a Lambert Azimuthal
>    Equal Area (LAEA) projection require: Earth's radius, center
>    latitude and longitude, and offset for east and north.
>    However, the CF Conventions document Appendix F (Grid Mappings),
>    only refers to an Earth radius in the Vertical Perspective grid
>    mapping.
>    For NetCDF file readers, how are the grid element locations to be
>    resolved for mappings such as LAEA if no Earth radius is provided
>    with the file?  Is a projection-specific radius to be implied by
>    NetCDF
>    readers?  Specifically, do all grid mappings use a prescribed set of
>    Earth radius or radii (major & minor) that NetCDF readers must
>    maintain?
>    The CF Grid Mappings appendix refers readers to a Web site that
>    discusses "[2]random issues" in which it is stated that GCTP always
>    uses a
>    radius of 6370997 meters (E019 - Normal Sphere).  I suppose this
>    answers my question, regarding the LAEA projection.  However, given
>    the multitude of grid mappings and Earth radius definitions that
>    exist, I'm
>    wondering if an earth_radius attribute ought to be added to the CF
>    grid
>    mapping parameter set for all projections (perhaps earth_radius_major
>    and earth_radius_minor, if applicable).
>    Greg.
>    ---------------------------------------------------------------
>    Greg Rappa
>    MIT Lincoln Laboratory               Group 43 - Weather Sensing
>    Phone: 781-981-3660, Fax: 781-981-0632, Email: [3]address@hidden
> 
> References
> 
>    1. http://gcmd.nasa.gov/records/USGS-GCTP.html
>    2. http://www.remotesensing.org/geotiff/proj_list/random_issues.html#laea
>    3. mailto:address@hidden

From address@hidden  Tue Dec 16 16:29:15 2008
Return-Path: <address@hidden>
X-Original-To: address@hidden
Delivered-To: address@hidden
Received: by laraine.unidata.ucar.edu (Postfix, from userid 60001)
        id C6069CB94D; Tue, 16 Dec 2008 16:29:15 -0700 (MST)
X-Spam-Checker-Version: SpamAssassin 3.2.5 (2008-06-10) on
        uni3.unidata.ucar.edu
X-Spam-Level: 
X-Spam-Status: No, score=-6.6 required=5.0 tests=BAYES_00,RCVD_IN_DNSWL_MED,
        UNPARSEABLE_RELAY autolearn=ham version=3.2.5
Received: from ll.mit.edu (LLMAIL1.LL.MIT.EDU [129.55.12.41])
        by laraine.unidata.ucar.edu (Postfix) with ESMTP id CDBD4CB94B;
        Tue, 16 Dec 2008 16:28:19 -0700 (MST)
Received: (from smtp@localhost)
        by ll.mit.edu (8.12.10/8.8.8) id mBGNSBr0023085;
        Tue, 16 Dec 2008 18:28:11 -0500 (EST)
Received: from sampson.llan.ll.mit.edu(              ), claiming to be "[       
       ]"
 via SMTP by llpost, id smtpdAAAx7aOAS; Tue Dec 16 18:27:36 2008
Message-ID: <address@hidden>
Date: Tue, 16 Dec 2008 18:27:36 -0500
From: Greg Rappa <address@hidden>
User-Agent: Thunderbird 2.0.0.16 (X11/20080707)
MIME-Version: 1.0
To: Jonathan Gregory <address@hidden>
CC: "address@hidden" <address@hidden>,
        "'John Caron'" <address@hidden>,
        "Newell, Oliver" <address@hidden>,
        "address@hidden" <address@hidden>
Subject: Re: Earth Radius for CF Grid Mappings
References: <address@hidden.local> <address@hidden> <address@hidden> 
<address@hidden> <address@hidden> <address@hidden> <200812162317
In-Reply-To: <address@hidden>
Content-Type: text/plain; charset=ISO-8859-1; format=flowed
Content-Transfer-Encoding: 7bit

OK.  Thanks.  I didn't notice the definitions in Table F.1.  I'm all set 
now.

Greg.

Jonathan Gregory wrote:
> Dear Greg
>
> You can include earth_radius with any of the mappings. At the start of App F
> it says
> "The attributes which describe the ellipsoid and prime meridian may be 
> included, when applicable, with any grid mapping."
> and earth_radius is one of those. If it's not specified but you need it, I
> think you can assume any sensible value.
>
> Cheers
>
> Jonathan
>