Expedition Expedition Navigation Software

You are not logged in. Would you like to login or register?



2/16/2022 3:22 pm  #1


Unsupported grib2 file

I am trying to open a historical RTOFS data set in expedition, downloaded from Noaa in grib2 format, expedition is not able to open it, would like to understand if that is something that we could see in the future.

Noaa homepage
Output fields from the NOAA ... on NOAA OneStop

Download page: Index of /data/oceans/ncep/rtofs (noaa.gov) (tried multiple)

Details of the dataset when opened in Panoply:netcdf file:/C:/Users/admin/OneDrive/Desktop/ofs_atl.t00z.f001.20170321.grb.grib2 {
  dimensions:
    lon = 1200;
    lat = 1684;
    reftime = 1;
    time = 1;
    time1 = 1;
    hybrid_layer = 1;
  variables:
    double reftime;
      :units = "Hour since 2017-03-21T00:00:00Z";
      :standard_name = "forecast_reference_time";
      :long_name = "GRIB reference time";
      :calendar = "proleptic_gregorian";

    double time(time=1);
      :units = "Hour since 2017-03-21T00:00:00Z";
      :standard_name = "time";
      :long_name = "GRIB forecast or observation time";
      :calendar = "proleptic_gregorian";

    double time1(time1=1);
      :units = "Hour since 2017-03-21T00:00:00Z";
      :standard_name = "time";
      :long_name = "GRIB forecast or observation time";
      :calendar = "proleptic_gregorian";

    float hybrid_layer(hybrid_layer=1);
      :units = "sigma";
      :long_name = "Hybrid level";
      ositive = "down";
      :Grib_level_type = 105; // int
      :bounds = "hybrid_layer_bounds";

    float hybrid_layer_bounds(hybrid_layer=1, 2);
      :units = "sigma";
      :long_name = "bounds for hybrid_layer";

    float Mixed_layer_depth_surface(time1=1, lat=1684, lon=1200);
      :long_name = "Mixed layer depth @ Ground or water surface";
      :units = "m";
      :abbreviation = "MIXHT";
      :missing_value = NaNf; // float
      :coordinates = "reftime time1 lat lon Latitude_of_Presure_Point_surface  Longitude_of_Presure_Point_surface ";
      :Grib_Variable_Id = "VAR_0-19-3_L1";
      :Grib2_Parameter = 0, 19, 3; // int
      :Grib2_Parameter_Discipline = "Meteorological products";
      :Grib2_Parameter_Category = "Physical atmospheric Properties";
      :Grib2_Parameter_Name = "Mixed layer depth";
      :Grib2_Level_Type = 1; // int
      :Grib2_Level_Desc = "Ground or water surface";
      :Grib2_Generating_Process_Type = "Forecast";
      :Grib2_Statistical_Process_Type = "UnknownStatType--1";

    float Montgomery_stream_function_surface(time1=1, lat=1684, lon=1200);
      :long_name = "Montgomery stream function @ Ground or water surface";
      :units = "m2.s-2";
      :abbreviation = "MNTSF";
      :missing_value = NaNf; // float
      :coordinates = "reftime time1 lat lon Latitude_of_Presure_Point_surface  Longitude_of_Presure_Point_surface ";
      :Grib_Variable_Id = "VAR_0-2-6_L1";
      :Grib2_Parameter = 0, 2, 6; // int
      :Grib2_Parameter_Discipline = "Meteorological products";
      :Grib2_Parameter_Category = "Momentum";
      :Grib2_Parameter_Name = "Montgomery stream function";
      :Grib2_Level_Type = 1; // int
      :Grib2_Level_Desc = "Ground or water surface";
      :Grib2_Generating_Process_Type = "Forecast";
      :Grib2_Statistical_Process_Type = "UnknownStatType--1";

    float Latitude_of_U_Wind_Component_of_Velocity_surface(lat=1684, lon=1200);
      :long_name = "Latitude of U Wind Component of Velocity @ Ground or water surface";
      :units = "degrees_north";
      :abbreviation = "LAUV";
      :missing_value = NaNf; // float
      :coordinates = "reftime time lat lon ";
      :Grib_Variable_Id = "VAR_0-2-198_L1";
      :Grib2_Parameter = 0, 2, 198; // int
      :Grib2_Parameter_Discipline = "Meteorological products";
      :Grib2_Parameter_Category = "Momentum";
      :Grib2_Parameter_Name = "Latitude of U Wind Component of Velocity";
      :Grib2_Level_Type = 1; // int
      :Grib2_Level_Desc = "Ground or water surface";
      :Grib2_Generating_Process_Type = "Forecast";
      :Grib2_Statistical_Process_Type = "UnknownStatType--1";

    float Longitude_of_U_Wind_Component_of_Velocity_surface(lat=1684, lon=1200);
      :long_name = "Longitude of U Wind Component of Velocity @ Ground or water surface";
      :units = "degrees_east";
      :abbreviation = "LOUV";
      :missing_value = NaNf; // float
      :coordinates = "reftime time lat lon ";
      :Grib_Variable_Id = "VAR_0-2-199_L1";
      :Grib2_Parameter = 0, 2, 199; // int
      :Grib2_Parameter_Discipline = "Meteorological products";
      :Grib2_Parameter_Category = "Momentum";
      :Grib2_Parameter_Name = "Longitude of U Wind Component of Velocity";
      :Grib2_Level_Type = 1; // int
      :Grib2_Level_Desc = "Ground or water surface";
      :Grib2_Generating_Process_Type = "Forecast";
      :Grib2_Statistical_Process_Type = "UnknownStatType--1";

    float Latitude_of_V_Wind_Component_of_Velocity_surface(lat=1684, lon=1200);
      :long_name = "Latitude of V Wind Component of Velocity @ Ground or water surface";
      :units = "degrees_north";
      :abbreviation = "LAVV";
      :missing_value = NaNf; // float
      :coordinates = "reftime time lat lon ";
      :Grib_Variable_Id = "VAR_0-2-200_L1";
      :Grib2_Parameter = 0, 2, 200; // int
      :Grib2_Parameter_Discipline = "Meteorological products";
      :Grib2_Parameter_Category = "Momentum";
      :Grib2_Parameter_Name = "Latitude of V Wind Component of Velocity";
      :Grib2_Level_Type = 1; // int
      :Grib2_Level_Desc = "Ground or water surface";
      :Grib2_Generating_Process_Type = "Forecast";
      :Grib2_Statistical_Process_Type = "UnknownStatType--1";

    float Longitude_of_V_Wind_Component_of_Velocity_surface(lat=1684, lon=1200);
      :long_name = "Longitude of V Wind Component of Velocity @ Ground or water surface";
      :units = "degrees_east";
      :abbreviation = "LOVV";
      :missing_value = NaNf; // float
      :coordinates = "reftime time lat lon ";
      :Grib_Variable_Id = "VAR_0-2-201_L1";
      :Grib2_Parameter = 0, 2, 201; // int
      :Grib2_Parameter_Discipline = "Meteorological products";
      :Grib2_Parameter_Category = "Momentum";
      :Grib2_Parameter_Name = "Longitude of V Wind Component of Velocity";
      :Grib2_Level_Type = 1; // int
      :Grib2_Level_Desc = "Ground or water surface";
      :Grib2_Generating_Process_Type = "Forecast";
      :Grib2_Statistical_Process_Type = "UnknownStatType--1";

    float Latitude_of_Presure_Point_surface(lat=1684, lon=1200);
      :long_name = "Latitude of Presure Point @ Ground or water surface";
      :units = "degrees_north";
      :abbreviation = "LAPP";
      :missing_value = NaNf; // float
      :coordinates = "reftime time lat lon ";
      :Grib_Variable_Id = "VAR_0-2-202_L1";
      :Grib2_Parameter = 0, 2, 202; // int
      :Grib2_Parameter_Discipline = "Meteorological products";
      :Grib2_Parameter_Category = "Momentum";
      :Grib2_Parameter_Name = "Latitude of Presure Point";
      :Grib2_Level_Type = 1; // int
      :Grib2_Level_Desc = "Ground or water surface";
      :Grib2_Generating_Process_Type = "Forecast";
      :Grib2_Statistical_Process_Type = "UnknownStatType--1";

    float Longitude_of_Presure_Point_surface(lat=1684, lon=1200);
      :long_name = "Longitude of Presure Point @ Ground or water surface";
      :units = "degrees_east";
      :abbreviation = "LOPP";
      :missing_value = NaNf; // float
      :coordinates = "reftime time lat lon ";
      :Grib_Variable_Id = "VAR_0-2-203_L1";
      :Grib2_Parameter = 0, 2, 203; // int
      :Grib2_Parameter_Discipline = "Meteorological products";
      :Grib2_Parameter_Category = "Momentum";
      :Grib2_Parameter_Name = "Longitude of Presure Point";
      :Grib2_Level_Type = 1; // int
      :Grib2_Level_Desc = "Ground or water surface";
      :Grib2_Generating_Process_Type = "Forecast";
      :Grib2_Statistical_Process_Type = "UnknownStatType--1";

    float Barotropic_U_velocity_entire_ocean_single_layer(time1=1, lat=1684, lon=1200);
      :long_name = "Barotropic U velocity @ Entire ocean layer";
      :units = "m.s-1";
      :abbreviation = "UBARO";
      :missing_value = NaNf; // float
      :coordinates = "reftime time1 lat lon Latitude_of_Presure_Point_surface  Longitude_of_Presure_Point_surface ";
      :Grib_Variable_Id = "VAR_10-1-194_L201";
      :Grib2_Parameter = 10, 1, 194; // int
      :Grib2_Parameter_Discipline = "Oceanographic products";
      :Grib2_Parameter_Category = "Currents";
      :Grib2_Parameter_Name = "Barotropic U velocity";
      :Grib2_Level_Type = 201; // int
      :Grib2_Level_Desc = "Entire ocean layer";
      :Grib2_Generating_Process_Type = "Forecast";
      :Grib2_Statistical_Process_Type = "UnknownStatType--1";

    float Barotropic_V_velocity_entire_ocean_single_layer(time1=1, lat=1684, lon=1200);
      :long_name = "Barotropic V velocity @ Entire ocean layer";
      :units = "m.s-1";
      :abbreviation = "VBARO";
      :missing_value = NaNf; // float
      :coordinates = "reftime time1 lat lon Latitude_of_Presure_Point_surface  Longitude_of_Presure_Point_surface ";
      :Grib_Variable_Id = "VAR_10-1-195_L201";
      :Grib2_Parameter = 10, 1, 195; // int
      :Grib2_Parameter_Discipline = "Oceanographic products";
      :Grib2_Parameter_Category = "Currents";
      :Grib2_Parameter_Name = "Barotropic V velocity";
      :Grib2_Level_Type = 201; // int
      :Grib2_Level_Desc = "Entire ocean layer";
      :Grib2_Generating_Process_Type = "Forecast";
      :Grib2_Statistical_Process_Type = "UnknownStatType--1";

    float Sea_Surface_Height_Relative_to_Geoid_surface(time1=1, lat=1684, lon=1200);
      :long_name = "Sea Surface Height Relative to Geoid @ Ground or water surface";
      :units = "m";
      :abbreviation = "SSHG";
      :missing_value = NaNf; // float
      :coordinates = "reftime time1 lat lon Latitude_of_Presure_Point_surface  Longitude_of_Presure_Point_surface ";
      :Grib_Variable_Id = "VAR_10-3-195_L1";
      :Grib2_Parameter = 10, 3, 195; // int
      :Grib2_Parameter_Discipline = "Oceanographic products";
      :Grib2_Parameter_Category = "Surface Properties";
      :Grib2_Parameter_Name = "Sea Surface Height Relative to Geoid";
      :Grib2_Level_Type = 1; // int
      :Grib2_Level_Desc = "Ground or water surface";
      :Grib2_Generating_Process_Type = "Forecast";
      :Grib2_Statistical_Process_Type = "UnknownStatType--1";

    float 3-D_Temperature_hybrid_layer(time1=1, hybrid_layer=1, lat=1684, lon=1200);
      :long_name = "3-D Temperature @ Hybrid level layer";
      :units = "c";
      :abbreviation = "WTMPC";
      :missing_value = NaNf; // float
      :coordinates = "reftime time1 hybrid_layer lat lon Latitude_of_Presure_Point_surface  Longitude_of_Presure_Point_surface ";
      :Grib_Variable_Id = "VAR_10-4-192_L105_layer";
      :Grib2_Parameter = 10, 4, 192; // int
      :Grib2_Parameter_Discipline = "Oceanographic products";
      :Grib2_Parameter_Category = "Sub-surface Properties";
      :Grib2_Parameter_Name = "3-D Temperature";
      :Grib2_Level_Type = 105; // int
      :Grib2_Level_Desc = "Hybrid level";
      :Grib2_Generating_Process_Type = "Forecast";
      :Grib2_Statistical_Process_Type = "UnknownStatType--1";

    float 3-D_Salinity_hybrid_layer(time1=1, hybrid_layer=1, lat=1684, lon=1200);
      :long_name = "3-D Salinity @ Hybrid level layer";
      :units = "psu";
      :abbreviation = "SALIN";
      :missing_value = NaNf; // float
      :coordinates = "reftime time1 hybrid_layer lat lon Latitude_of_Presure_Point_surface  Longitude_of_Presure_Point_surface ";
      :Grib_Variable_Id = "VAR_10-4-193_L105_layer";
      :Grib2_Parameter = 10, 4, 193; // int
      :Grib2_Parameter_Discipline = "Oceanographic products";
      :Grib2_Parameter_Category = "Sub-surface Properties";
      :Grib2_Parameter_Name = "3-D Salinity";
      :Grib2_Level_Type = 105; // int
      :Grib2_Level_Desc = "Hybrid level";
      :Grib2_Generating_Process_Type = "Forecast";
      :Grib2_Statistical_Process_Type = "UnknownStatType--1";

    float u-component_of_current_hybrid_layer(time1=1, hybrid_layer=1, lat=1684, lon=1200);
      :long_name = "u-component of current @ Hybrid level layer";
      :units = "m/s";
      :abbreviation = "UOGRD";
      :missing_value = NaNf; // float
      :coordinates = "reftime time1 hybrid_layer lat lon Latitude_of_U_Wind_Component_of_Velocity_surface Longitude_of_U_Wind_Component_of_Velocity_surface ";
      :Grib_Variable_Id = "VAR_10-1-2_L105_layer";
      :Grib2_Parameter = 10, 1, 2; // int
      :Grib2_Parameter_Discipline = "Oceanographic products";
      :Grib2_Parameter_Category = "Currents";
      :Grib2_Parameter_Name = "u-component of current";
      :Grib2_Level_Type = 105; // int
      :Grib2_Level_Desc = "Hybrid level";
      :Grib2_Generating_Process_Type = "Forecast";
      :Grib2_Statistical_Process_Type = "UnknownStatType--1";

    float v-component_of_current_hybrid_layer(time1=1, hybrid_layer=1, lat=1684, lon=1200);
      :long_name = "v-component of current @ Hybrid level layer";
      :units = "m/s";
      :abbreviation = "VOGRD";
      :missing_value = NaNf; // float
      :coordinates = "reftime time1 hybrid_layer lat lon Latitude_of_V_Wind_Component_of_Velocity_surface Longitude_of_V_Wind_Component_of_Velocity_surface ";
      :Grib_Variable_Id = "VAR_10-1-3_L105_layer";
      :Grib2_Parameter = 10, 1, 3; // int
      :Grib2_Parameter_Discipline = "Oceanographic products";
      :Grib2_Parameter_Category = "Currents";
      :Grib2_Parameter_Name = "v-component of current";
      :Grib2_Level_Type = 105; // int
      :Grib2_Level_Desc = "Hybrid level";
      :Grib2_Generating_Process_Type = "Forecast";
      :Grib2_Statistical_Process_Type = "UnknownStatType--1";

  // global attributes:
  :Originating_or_generating_Center = "US National Weather Service, National Centres for Environmental Prediction (NCEP)";
  :Originating_or_generating_Subcenter = "0";
  :GRIB_table_version = "2,1";
  :Type_of_generating_process = "Forecast";
  :Analysis_or_forecast_generating_process_identifier_defined_by_originating_centre = "HYCOM - North Pacific basin";
  :file_format = "GRIB-2";
  :Conventions = "CF-1.6";
  :history = "Read using CDM IOSP GribCollection v3";
  :featureType = "GRID";
}

 

2/16/2022 3:59 pm  #2


Re: Unsupported grib2 file

It isn't on a projection that Exp suports. I think it is a curvilinear orthogonal grid, but with the points defined in the lat & lon fields in the file.

You might be able to use wgrib2 (-lola etc) to convert it to a projection Exp supports (equirectangular, Mercator, Lambert, polar stereo, Gaussian etc).

 

Board footera

 

Powered by Boardhost. Create a Free Forum


Interested in advertising here? Over a thousand active navigators and Expedition users visit this forum. Click here to contact the administrator.