[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: Fwd: Polar stereo projection parameter name problem

--1971325068-1010076504-1192724011=:29009
Content-Type: TEXT/PLAIN; charset=X-UNKNOWN; format=flowed
Content-Transfer-Encoding: QUOTED-PRINTABLE

On Thu, 18 Oct 2007, Ben Domenico wrote:

> Hi John and Robb,
>
> This is a follow up on an error in a variable name in a file served by th=
e
> 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=20
wasn't 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=20
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 be=
=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=20
done since Grib uses numbers internally so the parameters will be=20
correct.



Thanks,
Robb...



The THREDDS catalog for such datasets is
> at:
>
>
> http://nomads.ncdc.noaa.gov:8085/thredds/catalogServices?catalog=3Dhttp:/=
/nomads.ncdc.noaa.gov:8085/thredds/catalog/radar/StIV/2006/catalog.xml
>
> The NCDC folks are talking about using ncML to correct the problem, but i=
t
> could be that the grib iosp is putting the wrong name in.  I thought this
> was worth checking because, if it is the iosp, it will pop up on any TDS
> that serves grib files in polar stereographic.
>
> -- Ben
>
> ---------- 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 <address@hidden.e=
du>,
> Nomads NCDC <address@hidden>, Brian Nelson <address@hidden>,
> Rich Baldwin <address@hidden>
>
> Hello all,
>
> The STIV (Stage-4 Multisensor Precip from NCEP) are in GRIB files from
> NCEP.
>
> John,
> Is this a parameter that should be different in the GRIB file?  Can we us=
e
> NcML to describe an 'alias' for this parameter?  If so, can you send a
> sample.
>
> Thanks,
> Steve
>
>
>
> Glenn.Rutledge said the following on 10/18/2007 8:24 AM:
>
> Hi Ben,
> This server data is populated and maintained (mostly) by Steve Ansari and
> Brian Nelson.  They are cc;ed to this and I'm certain they'll look into i=
t.
>
>
> It turns out that this data, and some others (GOES),k and projected grids
> have some " incompatibility" issues for their datum description while bei=
ng
> 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.
>
> Regards, Glenn
>
> Ben Domenico wrote the following on 10/17/2007 3:36 PM:
>
> Hi Glenn,
>
> You may not be the right person for this note, but I'll start with you wi=
th
> the idea you can forward it to the appropriate person.  As the note below
> from John Caron indicates, there is a problem with the name being used to
> 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, th=
is
> is not an uncommon error, but it would be good to fix it.
>
> Thanks.
>
> -- Ben
>
> ---------- 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.e=
du
>
> Hi Noman:
>
> Unfortunately its common for these files to not follow the conventions
> exactly.
>
> "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 for
> both names.
>
> Ben, can we notify the writer? Also, can I get a sample of the file for m=
y
> tests?
>
> John
>
> Ben Domenico wrote:
>> Noman,
>>
>> I'm forwarding this message to John Caron who can help you with the
>> details of the projection in that netCDF file.
>>
>> -- Ben
>>
>> 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 email
>>     below) and the parameters listed on the CF Convention page at _
>>     http://cf-pcmdi.llnl.gov/documents/cf-conventions/1.0/cf-conventions=
=2Ehtml#appendix-grid-mappings_
>
>>     Example F.5.
>>
>>     *Example F.5. Polar stereographic*
>>
>>     grid_mapping_name =3D polar_stereographic
>>
>>
>>
>>     */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 from
>>
>>     the IDV.  In each case, you can see at the bottom of the screen, the
>>
>>     lat/lon values associated with the cursor position.  I'll also appen=
d
>>
>>     the ncDump output which shows the projection information -- in netCD=
F
>>
>>     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) subcent=
er
>>     =3D 4";
>>
>>      :Generating_Process_or_Model =3D "River Forecast Center Quantitativ=
e
>>
>>     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>> wrot=
e:
>>
>>    >
>>
>>    >  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 whether
>>
>>    >  this comes through OK.
>>
>>    >
>>
>>    >  -- Ben
>>
>>    >
>>
>>    >
>>
>>
>
>
> --=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
>
> NOMADS: http://nomads.ncdc.noaa.gov/
>
>
> --=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
>

=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=
=3D=3D=3D=3D=3D=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-1010076504-1192724011=:29009--