Custom Query (125 matches)

Filters
 
Or
 
  
 
Columns

Show under each result:


Results (67 - 69 of 125)

Ticket Resolution Summary Owner Reporter
#89 fixed standard names for vector components davidhassell markh
Description

Objective

A reinterpretation of current standard names to make the identification of vector components clear and able to meet the needs of users.

This issue is related to the proposal on #79

Proposal

To adopt the constrained standard name concept to re-interpret vector quantity standard names, without invalidating any current datasets. This would involve:

  • 'x_' type standard names being valid for all coordinate definitions:
    • '"x" indicates a vector component along the grid x-axis, positive with increasing x.'
  • 'eastward_' type standard names being valid for all 'true east' vectors:
    • '"Eastward" indicates a vector component which is positive when directed eastward (negative westward); where eastward is defined as the grid x-axis direction, this is a constrained version of the "x_" standard name';
    • this may be interpreted in two ways, as:
      1. where eastward is defined as the grid x-axis, this standard name is a constrained version of x_wind
      2. where eastward is not defined as the grid x-axis, this standard name stands independently

This enables data producers to use eastward wind in the same way they currently do, while meeting my requirements, for datasets where x may or may not be east, depending on the location and for data format interoperability with formats which do not have an explicit 'eastward_' phenomenon definition.

It enables datasets to be written where:

  • vector is x but not east
    • standard_name: x_<>
  • vector is x and may be east or eastish
    • standard_name: x_<>
  • vector is x and happens to be always east
    • standard_name: x_<>
  • vector is x constrained to always be east
    • standard_name: eastward_<>
  • vector is east but not x
    • standard_name: eastward_<>

'eastward_<>' is already interpreted in multiple ways, depending on the coordinate variable context of the dataset. 'x_<>' should also be abe to be interpreted based on coordinate variable context, to enable datasets to be encoded which currently cannot be written in a CF compliant fashion

Analogy

This approach, of constraining standard names, is analogous to qualification. For example:

  • there is a standard name of air_pressure
  • this could be defined, for a particular dataset, such that the vertical coordinate indicates that the data is at a surface
  • if the fact that the dataset is at a surface is intrinsic to the data, the qualified (constrained) standard name may be used: surface_air_pressure
#92 fixed Add oblique mercator projection davidhassell mcginnis
Description

The Oblique Mercator projection is used by at least one regional climate model, RegCM3, which is part of the NARCCAP climate modeling program. Currently we record its map projection information using the transverse_mercator projection, which I have learned is very similar but not quite the same. I propose to add this map projection so we can get it right.

Proposed text:


Oblique Mercator

grid_mapping_name = oblique_mercator

Map parameters:

  • azimuth
  • latitude_of_projection_origin
  • longitude_of_projection_origin
  • scale_factor_at_projection_origin
  • false_easting
  • false_northing

Map coordinates:

The x (abscissa) and y (ordinate) rectangular coordinates are identified by the standard_name attribute value projection_x_coordinate and projection_y_coordinate respectively.

Notes:

Notes on using the PROJ.4 software package for computing the mapping may be found at http://www.remotesensing.org/geotiff/proj_list/hotine_oblique_mercator.html . The Rotated Mercator projection is an Oblique Mercator projection with azimuth = +90.


If adding a new attribute for azimuth is problematic, this proposal could be modified to add the rotated_mercator projection instead, which is a special case of Oblique Mercator with azimuth = 90.

Note that apparently there is a subtle technical difference between an Oblique Mercator projection and a Hotine Oblique Mercator projection that depends on when the rectification from skew grid to map grid is applied. Since most mapping packages don't support a rectified grid angle parameter at all (effectively giving it a default value of 90 degrees, such that it has no effect), to avoid unnecessary proliferation of attributes I propose to omit this parameter and elide this distinction until such time as it proves necessary.

My knowledge of this topic is quite limited; I have made this proposal based on what understanding I have gleaned from the geotiff website and communications with colleagues working with our RegCM3 output in GIS. Commentary from experts would be very welcome.

#93 fixed Two new dimensionless vertical coordinate specifications for s coordinate ocean models cf-conventions@… rsignell
Description

1. Title

Two new dimensionless vertical coordinates to support ocean models

2. Moderator

Rich Signell

3. Requirement

The ocean modeling community needs two additional vertical coordinate specifications to allow modern s-coordinate model output to be CF-compliant and allow more general specification of s-coordinate model output for future development.

4. Initial Statement of Technical Proposal

The existing ocean_s_coordinate dimensionless vertical coordinate specification in CF is limited to a specific vertical stretching function and set of control parameters, while modern versions of ROMS allow for more flexible specification. Additional of these two new generalized vertical coordinate specifications would allow output from existing ROMS-derived models (and other s coordinate models) to be CF-compliant, as well as allowing more flexibility for future s_coordinate model developers and users to write and read CF-compliant model output.

Here are the proposed additions to Appendix D. Dimensionless Vertical Coordinates

Ocean s-coordinate, generic form 1

   standard_name = "ocean_s_coordinate_g1"

Definition:

   z(n,k,j,i) = S(k,j,i) + eta(n,j,i) * (1 + S(k,j,i) / depth(j,i))

    where S(k,j,i) = depth_c * s(k) + (depth(j,i) - depth_c) * C(k)

and where z(n,k,j,i) is the height, positive upwards, relative to ocean datum (e.g. mean sea level) at gridpoint (n,k,j,i); eta(n,j,i) is the height of the ocean surface, positive upwards, relative to ocean datum at gridpoint (n,j,i); s(k) is the dimensionless coordinate at vertical gridpoint (k) with a range of -1 <= s(k) <= 0, s(0) corresponds to eta(n,j,i) whereas s(-1) corresponds to depth(j,i); C(k) is the dimensionless vertical coordinate stretching function at gridpoint (k) with a range of -1 <= C(k) <= 0, C(0) corresponds to eta(n,j,i) whereas C(-1) corresponds to depth(j,i); and the constant depth_c, (positive value), is a critical depth controlling the stretching.

The format for the formula_terms attribute is

formula_terms = "s: var1 C: var2 eta: var3 depth: var4 depth_c: var5"

Ocean s-coordinate, generic form 2

   standard_name = "ocean_s_coordinate_g2"

Definition:

   z(n,k,j,i) = eta(n,j,i) + (eta(n,j,i) + depth(j,i)) * S(k,j,i)

    where S(k,j,i) = (depth_c * s(k) + depth(j,i) * C(k)) / (depth_c + depth(j,i))

and where z(n,k,j,i) is the height, positive upwards, relative to ocean datum (e.g. mean sea level) at gridpoint (n,k,j,i); eta(n,j,i) is the height of the ocean surface, positive upwards, relative to ocean datum at gridpoint (n,j,i); s(k) is the dimensionless coordinate at vertical gridpoint (k) with a range of -1 <= s(k) <= 0, s(0) corresponds to eta(n,j,i) whereas s(-1) corresponds to depth(j,i); C(k) is the dimensionless vertical coordinate stretching function at gridpoint (k) with a range of -1 <= C(k) <= 0, C(0) corresponds to eta(n,j,i) whereas C(-1) correspond to depth(j,i); and the constant depth_c, (positive value) is a critical depth controlling the stretching.

The format for the formula_terms attribute is

formula_terms = "s: var1 C: var2 eta: var3 depth: var4 depth_c: var5"

5. Benefits

The oceanographic community, especially those producing or consuming products from s coordinate models (e.g. variants and derivatives of the Regional Ocean Modeling System (ROMS)) would benefit from the addition of these two new dimensionless vertical coordinate specifications. These two new coordinates have been exercised in the Unidata Common Data Model and the Unidata NetCDF-Java Library for the last two years. It's time to add them officially to the CF Conventions.

6. Status Quo

The only option currently to store modern ROMS results as CF compliant would be to write the entire Z field as a 4D array.

Note: See TracQuery for help on using queries.