Skip to content

Wrong builtin UWCW-projection information #31

Description

@mpvginde

The current built in grid/projection information for the UWCW (DINI) domain:

uwcw=dict(
projection="lcc",
kws_globe=dict(
semimajor_axis=6371229.0,
semiminor_axis=6371229.0,
),
kws_projection=dict(
central_longitude=-1.96590281,
central_latitude=55.5164337,
standard_parallels=[55.499996, 55.499996],
),
kws_grid=dict(
lon_ll=-25.4470005,
lat_ll=39.6389999,
lon_ur=40.1508102,
lat_ur=62.6713715,
dx=2000.0,
dy=2000.0,
nx=1909,
ny=1609,
),
),

is incorrect.

Looking directly at a UWCW-reforecast grib file:

  # Lambert conformal (Can be secant or tangent, conical or bipolar)  (grib2/tables/21/3.1.table)  
  gridDefinitionTemplateNumber = 30;
  # Earth assumed spherical with radius of 6 371 229.0 m (grib2/tables/21/3.2.table)  
  shapeOfTheEarth = 6;
  Nx = 1909;
  Ny = 1609;
  latitudeOfFirstGridPointInDegrees = 39.639;
  longitudeOfFirstGridPointInDegrees = 334.553;
  LaDInDegrees = 55.5;
  LoVInDegrees = 352;
  DxInMetres = 2000;
  DyInMetres = 2000;
  # (1=0) North Pole is on the projection plane;(2=0) Only one projection centre is used:grib2/tables/21/3.5.table 
  # flags: 00000000
  projectionCentreFlag = 0;
  iScansNegatively = 0;
  jScansPositively = 1;
  jPointsAreConsecutive = 0;
  alternativeRowScanning = 0;
  Latin1InDegrees = 55.5;
  Latin2InDegrees = 55.5;
  latitudeOfSouthernPoleInDegrees = -90;
  longitudeOfSouthernPoleInDegrees = 0;
  gridType = lambert;
  NV = 182;

the correct description should rather be:

    uwcw=dict(
        projection="lcc",
        kws_globe=dict(
            semimajor_axis=6371229.0,
            semiminor_axis=6371229.0,
        ),
        kws_projection=dict(
            central_longitude=-8,
            central_latitude=55.5,
            standard_parallels=[55.5, 55.5],
        ),
        kws_grid=dict(
            lon_ll=-25.447,
            lat_ll=39.639,
            lon_ur=40.1508102,
            lat_ur=62.6713715,
            dx=2000.0,
            dy=2000.0,
            nx=1909,
            ny=1609,
        ),
    ),

Tagging @bastien-francois, @Sooph-knmi for visibility.

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type
    No fields configured for issues without a type.

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions