[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: Fwd: Polar stereo projection parameter name problem
- Subject: Re: Fwd: Polar stereo projection parameter name problem
- Date: Thu, 18 Oct 2007 10:22:34 -0600 (MDT)
- Boundary="1971325068-2114347592-1192724554=: 29009"
--1971325068-2114347592-1192724554=:29009
Content-Type: TEXT/PLAIN; charset=X-UNKNOWN; format=flowed
Content-Transfer-Encoding: QUOTED-PRINTABLE
Please disregard the last message, I miss interpreted the question because=
=20
the message was long and I assumed it was a parameter question.
Sorry,
Robb...
On Thu, 18 Oct 2007, Robb Kambic wrote:
> On Thu, 18 Oct 2007, Ben Domenico wrote:
>
>> Hi John and Robb,
>>=20
>> This is a follow up on an error in a variable name in a file served by t=
he
>> NCDC NOMADS TDS. It turns out the source file is from ncep in grib.
>> Perhaps the GRIB iosp for these files is putting the wrong variable name=
in
>> for the polar stereographic case.
>
>
> All,
>
> This is a GRIB table problem and it occurs because the desired table wasn=
't=20
> available, so the code uses the closes table available.
>
> dumping a radar file gives table info:
> center =3D 7
> sub_center =3D 4
> table_version =3D 2
>
> table used was:
>
> center =3D 7
> sub_center =3D -1 default sub_center table_version =3D 2
>
> The fix is to get table 7 4 2 and to include in the resources and make a=
=20
> release or if you don't want to change your release, then the table can b=
e=20
> added and referenced. This is a more time consuming process.
>
> Does anyone have table 7 4 2 in NWS format ?
>
> Once the correct table is added to the code, no modifications need to be =
done=20
> since Grib uses numbers internally so the parameters will be correct.
>
>
>
> Thanks,
> Robb...
>
>
>
> The THREDDS catalog for such datasets is
>> at:
>>=20
>>=20
>> http://nomads.ncdc.noaa.gov:8085/thredds/catalogServices?catalog=3Dhttp:=
//nomads.ncdc.noaa.gov:8085/thredds/catalog/radar/StIV/2006/catalog.xml
>>=20
>> The NCDC folks are talking about using ncML to correct the problem, but =
it
>> could be that the grib iosp is putting the wrong name in. I thought thi=
s
>> was worth checking because, if it is the iosp, it will pop up on any TDS
>> that serves grib files in polar stereographic.
>>=20
>> -- Ben
>>=20
>> ---------- Forwarded message ----------
>> From: Steve.Ansari <address@hidden>
>> Date: Oct 18, 2007 7:09 AM
>> Subject: Re: Polar stereo projection parameter name problem
>> To: "Glenn.Rutledge" <address@hidden>
>> Cc: Ben Domenico <address@hidden>, John Caron=20
>> <address@hidden>,
>> Nomads NCDC <address@hidden>, Brian Nelson <address@hidden>=
,
>> Rich Baldwin <address@hidden>
>>=20
>> Hello all,
>>=20
>> The STIV (Stage-4 Multisensor Precip from NCEP) are in GRIB files from
>> NCEP.
>>=20
>> John,
>> Is this a parameter that should be different in the GRIB file? Can we u=
se
>> NcML to describe an 'alias' for this parameter? If so, can you send a
>> sample.
>>=20
>> Thanks,
>> Steve
>>=20
>>=20
>>=20
>> Glenn.Rutledge said the following on 10/18/2007 8:24 AM:
>>=20
>> Hi Ben,
>> This server data is populated and maintained (mostly) by Steve Ansari an=
d
>> Brian Nelson. They are cc;ed to this and I'm certain they'll look into =
it.
>>=20
>>=20
>> It turns out that this data, and some others (GOES),k and projected grid=
s
>> have some " incompatibility" issues for their datum description while be=
ing
>> access with the GeoPortal Tooklit by ESRI (aka our ACCESS conversions of
>> late)- as reported by Rich Baldwin here and he's also cc:ed just in case
>> this might be germane.
>>=20
>> Regards, Glenn
>>=20
>> Ben Domenico wrote the following on 10/17/2007 3:36 PM:
>>=20
>> Hi Glenn,
>>=20
>> You may not be the right person for this note, but I'll start with you w=
ith
>> the idea you can forward it to the appropriate person. As the note belo=
w
>> from John Caron indicates, there is a problem with the name being used t=
o
>> specify the latitude coordinate for the polar stereographic coordinate
>> system being used to serve STIV radar data on the nomads server. The
>> prescribed CF name is "straight_vertical_longitude_from_pole" but the
>> netCDF file uses "longitude_of_projection_origin." According to John, t=
his
>> is not an uncommon error, but it would be good to fix it.
>>=20
>> Thanks.
>>=20
>> -- Ben
>>=20
>> ---------- Forwarded message ----------
>> From: John Caron <address@hidden >
>> Date: Oct 17, 2007 1:08 PM
>> Subject: Re: attached file
>> To: Ben Domenico <address@hidden>
>> Cc: Nawajish Noman <address@hidden >, address@hidden.=
edu
>>=20
>> Hi Noman:
>>=20
>> Unfortunately its common for these files to not follow the conventions
>> exactly.
>>=20
>> "longitude_of_projection_origin" should be
>> "straight_vertical_longitude_from_pole". Its a common error, and we tend=
to
>> try to be lenient in these cases, so our software would probably look fo=
r
>> both names.
>>=20
>> Ben, can we notify the writer? Also, can I get a sample of the file for =
my
>> tests?
>>=20
>> John
>>=20
>> Ben Domenico wrote:
>>> Noman,
>>>=20
>>> I'm forwarding this message to John Caron who can help you with the
>>> details of the projection in that netCDF file.
>>>=20
>>> -- Ben
>>>=20
>>> On 10/16/07, * Nawajish Noman* <address@hidden
>>> <mailto: address@hidden>> wrote:
>>>
>>> Hi Ben,
>>>
>>> We are working on the problem with Polar stereographic projection. =
I
>>> am comparing the grid_mapping parameters in the file (see your emai=
l
>>> below) and the parameters listed on the CF Convention page at _
>>> http://cf-pcmdi.llnl.gov/documents/cf-conventions/1.0/cf-convention=
s.html#appendix-grid-mappings_
>>
>>> Example F.5.
>>>
>>> *Example F.5. Polar stereographic*
>>>
>>> grid_mapping_name =3D polar_stereographic
>>>=20
>>>=20
>>>
>>> */Map parameters:/*
>>>
>>> straight_vertical_longitude_from_pole
>>>
>>> latitude_of_projection_origin - Either +90. or -90.
>>>
>>> Either standard_parallel or scale_factor_at_projection_origin
>>>
>>> false_easting
>>>
>>> false_northing
>>>
>>> *Parameters in file are:***
>>>
>>> char Polar_Stereographic_projection_Grid;
>>>
>>> :grid_mapping_name =3D "polar_stereographic";
>>>
>>> :longitude_of_projection_origin =3D - 105.0; // double
>>>
>>> :scale_factor_at_projection_origin =3D 0.933; // double
>>>
>>> Here are my questions:
>>>
>>> Q1: Some of the parameters are missing in file. What do we do about
>>> missing parameters?
>>>
>>> Q2: You have " longitude_of_projection_origin " in file but the CF
>>> convention says " latitude_of_projection_origin" - What do with
>>> "longitude_of_projection_origin "?
>>>
>>> We look for parameters listed in CF Conventions and map them to our
>>> Spatial reference parameters. If parameters are missing in file we
>>> assume default value. If incorrect parameters are specified - we
>>> ignore them.
>>>
>>> I want to make sure we are reading the parameters from your file
>>> correctly. Please comment=85
>>>
>>> Thanks,
>>>
>>> Noman
>>>
>>> -----Original Message-----
>>> From: address@hidden <mailto:address@hidden> [
>>> mailto: address@hidden] On Behalf Of Ben Domenico
>>> Sent: Tuesday, June 05, 2007 1:56 PM
>>> To: Steve Kopp
>>> Cc: Nawajish Noman
>>> Subject: Re: attached file
>>>
>>> Steve,
>>>
>>> I'll attach jpegs of the dataset as displayed in arcMap and one fro=
m
>>>
>>> the IDV. In each case, you can see at the bottom of the screen, th=
e
>>>
>>> lat/lon values associated with the cursor position. I'll also appe=
nd
>>>
>>> the ncDump output which shows the projection information -- in netC=
DF
>>>
>>> terms (i.e., Polar_Stereographic_projection_Grid) anyway.
>>>
>>> -- Ben
>>>
>>> ncDump output
>>>
>>> =3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=
=3D=3D=3D=3D
>>>
>>> netcdf C:/Documents and
>>> Settings/ben/Desktop/ST4.2006010100.06hWithLatLon.nc {
>>>
>>> dimensions:
>>>
>>> time =3D 1; // (has coord.var)
>>>
>>> y =3D 815; // (has coord.var)
>>>
>>> x =3D 1121; // (has coord.var)
>>>
>>> variables:
>>>
>>> float Total_precipitation(time=3D1, y=3D815, x=3D1121);
>>>
>>> :units =3D "kg/m^2";
>>>
>>> :long_name =3D "Total precipitation @ surface";
>>>
>>> :missing_value =3D -9999.0f; // float
>>>
>>> :grid_mapping =3D "Polar_Stereographic_projection_Grid";
>>>
>>> :GRIB_param_name =3D "Total precipitation";
>>>
>>> :GRIB_center_id =3D 7; // int
>>>
>>> :GRIB_table_id =3D 2; // int
>>>
>>> :GRIB_param_number =3D 61; // int
>>>
>>> :GRIB_param_id =3D 1, 7, 2, 61; // int
>>>
>>> :GRIB_product_definition_type =3D "Accumulation";
>>>
>>> :GRIB_level_type =3D 1; // int
>>>
>>> int time(time=3D1);
>>>
>>> :long_name =3D "forecast time";
>>>
>>> :units =3D "hour since 2005-12-31T12:00:00Z";
>>>
>>> :GRIB_orgReferenceTime =3D "2005-12-31T12:00:00Z";
>>>
>>> :GRIB2_significanceOfRTName =3D "Start of forecast";
>>>
>>> :_CoordinateAxisType =3D "Time";
>>>
>>> double y(y=3D815);
>>>
>>> :units =3D "km";
>>>
>>> :long_name =3D "y coordinate of projection";
>>>
>>> :standard_name =3D "projection_y_coordinate";
>>>
>>> :grid_spacing =3D "4.763 km";
>>>
>>> :_CoordinateAxisType =3D "GeoY";
>>>
>>> double x(x=3D1121);
>>>
>>> :units =3D "km";
>>>
>>> :long_name =3D "x coordinate of projection";
>>>
>>> :standard_name =3D "projection_x_coordinate";
>>>
>>> :grid_spacing =3D "4.763 km";
>>>
>>> :_CoordinateAxisType =3D "GeoX";
>>>
>>> char Polar_Stereographic_projection_Grid;
>>>
>>> :grid_mapping_name =3D "polar_stereographic";
>>>
>>> :longitude_of_projection_origin =3D -105.0; // double
>>>
>>> :scale_factor_at_projection_origin =3D 0.933 ; // double
>>>
>>> :earth_shape =3D "spherical";
>>>
>>> :GRIB_param_Dx =3D 4763.0; // double
>>>
>>> :GRIB_param_Dy =3D 4763.0; // double
>>>
>>> :GRIB_param_GDSkey =3D 1294122487; // int
>>>
>>> :GRIB_param_La1 =3D 23.117; // double
>>>
>>> :GRIB_param_Lo1 =3D -119.023; // double
>>>
>>> :GRIB_param_LoV =3D -105.0; // double
>>>
>>> :GRIB_param_Nx =3D 1121; // int
>>>
>>> :GRIB_param_Ny =3D 881; // int
>>>
>>> :GRIB_param_ProjFlag =3D 0; // int
>>>
>>> :GRIB_param_ResCompFlag =3D 8; // int
>>>
>>> :GRIB_param_ScanningMode =3D 64; // int
>>>
>>> :GRIB_param_grid_name =3D "Polar Stereographic projection Grid=
";
>>>
>>> :GRIB_param_grid_radius_spherical_earth =3D 6367.47; // double
>>>
>>> :GRIB_param_grid_shape =3D "spherical";
>>>
>>> :GRIB_param_grid_shape_code =3D 0; // int
>>>
>>> :GRIB_param_grid_type =3D 5; // int
>>>
>>> :_CoordinateTransformType =3D "Projection";
>>>
>>> :_CoordinateAxisTypes =3D "GeoX GeoY";
>>>
>>> double lat(y=3D815, x=3D1121);
>>>
>>> :units =3D "degrees_north";
>>>
>>> :long_name =3D "latitude coordinate";
>>>
>>> :standard_name =3D "latitude";
>>>
>>> :_CoordinateAxisType =3D "Lat";
>>>
>>> double lon(y=3D815, x=3D1121);
>>>
>>> :units =3D "degrees_east";
>>>
>>> :long_name =3D "longitude coordinate";
>>>
>>> :standard_name =3D "longitude";
>>>
>>> :_CoordinateAxisType =3D "Lon";
>>>
>>> :History =3D "GridDatatype extracted from dataset
>>>
>>> /eclipse1a/ftp/pub/radar/StIV/2006/01/01/hr/ST4.2006010100.06h";
>>>
>>> :Conventions =3D "CF-1.0";
>>>
>>> :Originating_center =3D "US National Weather Service (NCEP) subcen=
ter
>>> =3D 4";
>>>
>>> :Generating_Process_or_Model =3D "River Forecast Center Quantitati=
ve
>>>
>>> Precipitation estimate mosaic generated by NCEP";
>>>
>>> :Product_Type =3D "Accumulation";
>>>
>>> :cdm_data_type =3D "Grid";
>>>
>>> :creator_name =3D "US National Weather Service (NCEP) subcenter =
=3D 4";
>>>
>>> :file_format =3D "GRIB-1";
>>>
>>> :location =3D
>>> "/eclipse1a/ftp/pub/radar/StIV/2006/01/01/hr/ST4.2006010100.06h";
>>>
>>> :history =3D "Direct read of GRIB into NetCDF-Java 2.2 API";
>>>
>>> :_CoordinateModelRunDate =3D "2005-12-31T12:00:00Z";
>>>
>>> =3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=
=3D=3D=3D=3D=3D=3D
>>>
>>> On 6/5/07, Steve Kopp < address@hidden <mailto:address@hidden>> wro=
te:
>>>
>>> >
>>>
>>> > Ben, the dataset contains no projection information. Without this=
,
>> we
>>>
>>> > won't know where else to draw anything.
>>>
>>> >
>>>
>>> > Could you send us a screen shot of what you are seeing?
>>>
>>> >
>>>
>>> > Thanks,
>>>
>>> >
>>>
>>> > Steve
>>>
>>> >
>>>
>>> >
>>>
>>> > -----Original Message-----
>>>
>>> > From: address@hidden <mailto:address@hidden> [
>>> mailto:address@hidden] On Behalf Of
>>>
>>> > Ben Domenico
>>>
>>> > Sent: Tuesday, June 05, 2007 11:32 AM
>>>
>>> > To: Steve Kopp
>>>
>>> > Subject: attached file
>>>
>>> >
>>>
>>> > Hi Steve,
>>>
>>> >
>>>
>>> > Don't know if this will work, but I'm attaching one of the
>> CF-netCDF
>>>
>>> > files I got from that NCDC radar data server. Let me know whethe=
r
>>>
>>> > this comes through OK.
>>>
>>> >
>>>
>>> > -- Ben
>>>
>>> >
>>>
>>> >
>>>=20
>>>=20
>>=20
>>=20
>> --=20
>> Glenn K. Rutledge
>> Services Team Leader
>> Remote Sensing and Applications Division
>> NOMADS Project Manager
>> National Oceanic and Atmospheric Administration
>> National Climatic Data Center
>> Asheville NC 28801
>> Phone: (828) 271-4097
>> Fax: (828) 271-4328
>>=20
>> NOMADS: http://nomads.ncdc.noaa.gov/
>>=20
>>=20
>> --=20
>> Steve Ansari
>> Physical Scientist
>> NOAA's National Climatic Data Center
>> Veach-Baley Federal Building
>> 151 Patton Avenue
>> Asheville, NC 28801
>> Ph: 828-271-4611
>> Fax: 828-271-4022
>>=20
>
> =3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=
=3D=3D=3D=3D=3D
> Robb Kambic =09=09=09 Unidata Program Center
> Software Engineer III=09=09=09 Univ. Corp for Atmospheric=20
> Research
> address@hidden=09=09 WWW: http://www.unidata.ucar.edu/
> =3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=
=3D=3D=3D=3D=3D
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=
=3D=3D=3D=3D
Robb Kambic =09=09=09 Unidata Program Center
Software Engineer III=09=09=09 Univ. Corp for Atmospheric Research
address@hidden=09=09 WWW: http://www.unidata.ucar.edu/
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=
=3D=3D=3D=3D
--1971325068-2114347592-1192724554=:29009--