Re: Lat/Lon Question

From: Winston Hao <wxhao_at_nyahnyahspammersnyahnyah>
Date: Thu Jul 10 2014 - 08:03:58 MDT

Hi Greg,
 
         On the left side of HRRR page ( http://ruc.noaa.gov/hrrr/ ), you =
can follow a link called "HRRR static fields inc lat/lon" to download =
geo_em.d01.nc Perhaps the information in that file would be useful to =
you?
 
                                   Winston

>>> Greg Garrett - NOAA Federal <greg.garrett@noaa.gov> 7/9/2014 4:11 PM =
>>>
Rick,

What I am really trying to do is extract a subset of data from a conus =
native lambert conformal projected domain. The bounding box of my sub =
domain is in lat/lon. Since lat/lon is not a grid in my grib2 file, I have =
no way that I can figure out how to tell ncl to only pull information out =
for my area of interest.

The grib2 file is the 3km HRRR model for the conus which is extremely =
large. I am using the opendap version of ncl in an effort to only extract =
parameter data for just my immediate area after opening the url of the =
grib2 file.

I have a ncl_filedump of one of the grib files below for reference. Any =
ideas?

Greg

Copyright (C) 1995-2014 - All Rights Reserved
University Corporation for Atmospheric Research
NCAR Command Language Version 6.2.0
The use of this software is governed by a License Agreement.
See http://www.ncl.ucar.edu/ for more details.


Variable: f
Type: file
filename: wrftwo_subh_RUN_2014-07-09T16:00:00Z
path: http://esrl.noaa.gov/gsd/thredds/dodsC/hrrr/conus/wrftwo_subh/runs/wr=
ftwo_subh_RUN_2014-07-09T16:00:00Z
file global attributes:
Originating_or_generating_Center : The NOAA Forecast Systems Laboratory, =
Boulder, CO, USA
Originating_or_generating_Subcenter : 0
GRIB_table_version : 2,1
Type_of_generating_process : Forecast
Analysis_or_forecast_generating_process_identifier_defined_by_originating_c=
entre : High-Resolution Rapid Refresh
Conventions : CF-1.4, _Coordinates
history : Read using CDM IOSP Grib2Collection ;
FMRC Run 2014-07-09T16:00:00Z Dataset
featureType : GRID
file_format : GRIB-2
_CoordSysBuilder : ucar.nc2.dataset.conv.CF1Convention
cdm_data_type : GRID
location : Proto fmrc:wrftwo_subh
dimensions:
ncl_scalar = 1
height_above_ground = 1
height_above_ground1 = 1
height_above_ground2 = 1
height_above_ground3 = 1
height_above_ground_layer = 1
height_above_ground_layer1 = 1
height_above_ground_layer1_bounds_1 = 2
height_above_ground_layer_bounds_1 = 2
time = 61
x = 1799
y = 1059
variables:
integer LambertConformal_Projection ( ncl_scalar )
long_name :
grid_mapping_name : lambert_conformal_conic
latitude_of_projection_origin : 38.5
longitude_of_central_meridian : 262.5
standard_parallel : 38.5
earth_radius : 6367470
_CoordinateTransformType : Projection
_CoordinateAxisTypes : GeoX GeoY


float height_above_ground_layer_bounds ( height_above_ground_layer, =
height_above_ground_layer_bounds_1 )
units : m
long_name : bounds for height_above_ground_layer


float height_above_ground_layer1_bounds ( height_above_ground_layer1, =
height_above_ground_layer1_bounds_1 )
units : m
long_name : bounds for height_above_ground_layer1


double time_offset ( time )
long_name : offset hour from start of run for coordinate = time
standard_name : forecast_period
units : hours since 2014-07-08T18:00:00Z
missing_value : nan
_FillValue : nan


float x ( x )
units : km
long_name : projection_x_coordinate
standard_name : projection_x_coordinate
_CoordinateAxisType : GeoX


float y ( y )
units : km
long_name : projection_y_coordinate
standard_name : projection_y_coordinate
_CoordinateAxisType : GeoY


float height_above_ground ( height_above_ground )
units : m
long_name : Specified height level above ground
positive : up
Grib2_level_type : 103
datum : ground
_CoordinateAxisType : Height
_CoordinateZisPositive : up


float height_above_ground1 ( height_above_ground1 )
units : m
long_name : Specified height level above ground
positive : up
Grib2_level_type : 103
datum : ground
_CoordinateAxisType : Height
_CoordinateZisPositive : up


float height_above_ground_layer ( height_above_ground_layer )
units : m
long_name : Specified height level above ground
positive : up
Grib2_level_type : 103
datum : ground
bounds : height_above_ground_layer_bounds
_CoordinateAxisType : Height
_CoordinateZisPositive : up


float height_above_ground_layer1 ( height_above_ground_layer1 )
units : m
long_name : Specified height level above ground
positive : up
Grib2_level_type : 103
datum : ground
bounds : height_above_ground_layer1_bounds
_CoordinateAxisType : Height
_CoordinateZisPositive : up


float height_above_ground2 ( height_above_ground2 )
units : m
long_name : Specified height level above ground
positive : up
Grib2_level_type : 103
datum : ground
_CoordinateAxisType : Height
_CoordinateZisPositive : up


float height_above_ground3 ( height_above_ground3 )
units : m
long_name : Specified height level above ground
positive : up
Grib2_level_type : 103
datum : ground
_CoordinateAxisType : Height
_CoordinateZisPositive : up


double time ( time )
long_name : Forecast time for ForecastModelRunCollection
standard_name : time
units : hours since 2014-07-08T18:00:00Z
missing_value : nan
_CoordinateAxisType : Time
_FillValue : nan


float Temperature_height_above_ground ( time, height_above_ground, y, x )
units : K
long_name : Temperature @ Specified height level above ground
missing_value : nan
description : Mean temperature in lowest 30 mb
grid_mapping : LambertConformal_Projection
Grib_Variable_Id : VAR_0-0-0_L103
Grib2_Parameter : ( 0, 0, 0 )
Grib2_Parameter_Discipline : Meteorological products
Grib2_Parameter_Category : Temperature
Grib2_Parameter_Name : Temperature
Grib2_Level_Type : 103
Grib2_Generating_Process_Type : Forecast
coordinates : time_run time height_above_ground y x
_FillValue : nan


float Dew-point_temperature_height_above_ground ( time, height_above_ground=
, y, x )
units : K
long_name : Dew-point temperature @ Specified height level above ground
missing_value : nan
description : Mean dewpoint temperature in lowest 30 mb
grid_mapping : LambertConformal_Projection
Grib_Variable_Id : VAR_0-0-6_L103
Grib2_Parameter : ( 0, 0, 6 )
Grib2_Parameter_Discipline : Meteorological products
Grib2_Parameter_Category : Temperature
Grib2_Parameter_Name : Dew-point temperature
Grib2_Level_Type : 103
Grib2_Generating_Process_Type : Forecast
coordinates : time_run time height_above_ground y x
_FillValue : nan


float Convective_Initiation1_Entire_atmosphere ( time, y, x )
units : minutes
long_name : Convective Initiation1 @ null
missing_value : nan
description : Time since last convective initation (> 1 hr) based upon =
lightning threat 3 and no adjacent activity
grid_mapping : LambertConformal_Projection
Grib_Variable_Id : VAR_0-0-195_L200
Grib2_Parameter : ( 0, 0, 195 )
Grib2_Parameter_Discipline : Meteorological products
Grib2_Parameter_Category : Temperature
Grib2_Parameter_Name : Convective Initiation1
Grib2_Level_Type : 200
Grib2_Generating_Process_Type : Forecast
coordinates : time_run time y x
_FillValue : nan


float Convective_Acitivity1_Entire_atmosphere ( time, y, x )
units : minutes
long_name : Convective Acitivity1 @ null
missing_value : nan
description : Time since last convective activity based upon lightning =
threat 3
grid_mapping : LambertConformal_Projection
Grib_Variable_Id : VAR_0-0-196_L200
Grib2_Parameter : ( 0, 0, 196 )
Grib2_Parameter_Discipline : Meteorological products
Grib2_Parameter_Category : Temperature
Grib2_Parameter_Name : Convective Acitivity1
Grib2_Level_Type : 200
Grib2_Generating_Process_Type : Forecast
coordinates : time_run time y x
_FillValue : nan


float Convective_Activitiy2_Entire_atmosphere ( time, y, x )
units : minutes
long_name : Convective Activitiy2 @ null
missing_value : nan
description : Time since last convective activity based upon vertical flux =
of rain and graupel
grid_mapping : LambertConformal_Projection
Grib_Variable_Id : VAR_0-0-201_L200
Grib2_Parameter : ( 0, 0, 201 )
Grib2_Parameter_Discipline : Meteorological products
Grib2_Parameter_Category : Temperature
Grib2_Parameter_Name : Convective Activitiy2
Grib2_Level_Type : 200
Grib2_Generating_Process_Type : Forecast
coordinates : time_run time y x
_FillValue : nan


float Convective_Activity3_Entire_atmosphere ( time, y, x )
units : minutes
long_name : Convective Activity3 @ null
missing_value : nan
description : Time since last convective activitiy based upon 35 dBZ =
reflectivity at -10C
grid_mapping : LambertConformal_Projection
Grib_Variable_Id : VAR_0-0-202_L200
Grib2_Parameter : ( 0, 0, 202 )
Grib2_Parameter_Discipline : Meteorological products
Grib2_Parameter_Category : Temperature
Grib2_Parameter_Name : Convective Activity3
Grib2_Level_Type : 200
Grib2_Generating_Process_Type : Forecast
coordinates : time_run time y x
_FillValue : nan


float Specific_humidity_height_above_ground ( time, height_above_ground, =
y, x )
units : kg/kg
long_name : Specific humidity @ Specified height level above ground
missing_value : nan
description : Specific humidity at 2 m above ground level
grid_mapping : LambertConformal_Projection
Grib_Variable_Id : VAR_0-1-0_L103
Grib2_Parameter : ( 0, 1, 0 )
Grib2_Parameter_Discipline : Meteorological products
Grib2_Parameter_Category : Moisture
Grib2_Parameter_Name : Specific humidity
Grib2_Level_Type : 103
Grib2_Generating_Process_Type : Forecast
coordinates : time_run time height_above_ground y x
_FillValue : nan


float Precipitation_rate_surface ( time, y, x )
units : kg.m-2.s-1
long_name : Precipitation rate @ Ground or water surface
missing_value : nan
description : Precipitation rate
grid_mapping : LambertConformal_Projection
Grib_Variable_Id : VAR_0-1-7_L1
Grib2_Parameter : ( 0, 1, 7 )
Grib2_Parameter_Discipline : Meteorological products
Grib2_Parameter_Category : Moisture
Grib2_Parameter_Name : Precipitation rate
Grib2_Level_Type : 1
Grib2_Generating_Process_Type : Forecast
coordinates : time_run time y x
_FillValue : nan


float Precipitation_Type4_surface ( time, y, x )
units : flag
long_name : Precipitation Type4 @ Ground or water surface
missing_value : nan
description : Categorical rain (yes=1 no=0)
grid_mapping : LambertConformal_Projection
Grib_Variable_Id : VAR_0-1-192_L1
Grib2_Parameter : ( 0, 1, 192 )
Grib2_Parameter_Discipline : Meteorological products
Grib2_Parameter_Category : Moisture
Grib2_Parameter_Name : Precipitation Type4
Grib2_Level_Type : 1
Grib2_Generating_Process_Type : Forecast
coordinates : time_run time y x
_FillValue : nan


float Precipitation_Type3_surface ( time, y, x )
units : flag
long_name : Precipitation Type3 @ Ground or water surface
missing_value : nan
description : Categorical freezing rain (yes=1 no=0)
grid_mapping : LambertConformal_Projection
Grib_Variable_Id : VAR_0-1-193_L1
Grib2_Parameter : ( 0, 1, 193 )
Grib2_Parameter_Discipline : Meteorological products
Grib2_Parameter_Category : Moisture
Grib2_Parameter_Name : Precipitation Type3
Grib2_Level_Type : 1
Grib2_Generating_Process_Type : Forecast
coordinates : time_run time y x
_FillValue : nan


float Precipitation_Type2_surface ( time, y, x )
units : flag
long_name : Precipitation Type2 @ Ground or water surface
missing_value : nan
description : Categorical ice pellets (yes=1 no=0)
grid_mapping : LambertConformal_Projection
Grib_Variable_Id : VAR_0-1-194_L1
Grib2_Parameter : ( 0, 1, 194 )
Grib2_Parameter_Discipline : Meteorological products
Grib2_Parameter_Category : Moisture
Grib2_Parameter_Name : Precipitation Type2
Grib2_Level_Type : 1
Grib2_Generating_Process_Type : Forecast
coordinates : time_run time y x
_FillValue : nan


float Precipitation_Type1_surface ( time, y, x )
units : flag
long_name : Precipitation Type1 @ Ground or water surface
missing_value : nan
description : Categorical snow (yes=1 no=0)
grid_mapping : LambertConformal_Projection
Grib_Variable_Id : VAR_0-1-195_L1
Grib2_Parameter : ( 0, 1, 195 )
Grib2_Parameter_Discipline : Meteorological products
Grib2_Parameter_Category : Moisture
Grib2_Parameter_Name : Precipitation Type1
Grib2_Level_Type : 1
Grib2_Generating_Process_Type : Forecast
coordinates : time_run time y x
_FillValue : nan


float Convective_Initiation2_Entire_atmosphere ( time, y, x )
units : minutes
long_name : Convective Initiation2 @ null
missing_value : nan
description : Time since last convective initation (> 1 hr) based upon =
flux of rain/graupel and no adjacent activity
grid_mapping : LambertConformal_Projection
Grib_Variable_Id : VAR_0-1-213_L200
Grib2_Parameter : ( 0, 1, 213 )
Grib2_Parameter_Discipline : Meteorological products
Grib2_Parameter_Category : Moisture
Grib2_Parameter_Name : Convective Initiation2
Grib2_Level_Type : 200
Grib2_Generating_Process_Type : Forecast
coordinates : time_run time y x
_FillValue : nan


float Convective_Initiation3_Entire_atmosphere ( time, y, x )
units : minutes
long_name : Convective Initiation3 @ null
missing_value : nan
description : Time since last convective initation (> 1 hr) based upon 35 =
dBZ reflectivity at -10C and no adjacent activity
grid_mapping : LambertConformal_Projection
Grib_Variable_Id : VAR_0-1-214_L200
Grib2_Parameter : ( 0, 1, 214 )
Grib2_Parameter_Discipline : Meteorological products
Grib2_Parameter_Category : Moisture
Grib2_Parameter_Name : Convective Initiation3
Grib2_Level_Type : 200
Grib2_Generating_Process_Type : Forecast
coordinates : time_run time y x
_FillValue : nan


float u-component_of_wind_height_above_ground ( time, height_above_ground3,=
 y, x )
units : m/s
long_name : u-component of wind @ Specified height level above ground
missing_value : nan
description : U-component of wind at 10 m above ground level
grid_mapping : LambertConformal_Projection
Grib_Variable_Id : VAR_0-2-2_L103
Grib2_Parameter : ( 0, 2, 2 )
Grib2_Parameter_Discipline : Meteorological products
Grib2_Parameter_Category : Momentum
Grib2_Parameter_Name : u-component of wind
Grib2_Level_Type : 103
Grib2_Generating_Process_Type : Forecast
coordinates : time_run time height_above_ground3 y x
_FillValue : nan


float v-component_of_wind_height_above_ground ( time, height_above_ground3,=
 y, x )
units : m/s
long_name : v-component of wind @ Specified height level above ground
missing_value : nan
description : V-component of wind at 10 m above ground level
grid_mapping : LambertConformal_Projection
Grib_Variable_Id : VAR_0-2-3_L103
Grib2_Parameter : ( 0, 2, 3 )
Grib2_Parameter_Discipline : Meteorological products
Grib2_Parameter_Category : Momentum
Grib2_Parameter_Name : v-component of wind
Grib2_Level_Type : 103
Grib2_Generating_Process_Type : Forecast
coordinates : time_run time height_above_ground3 y x
_FillValue : nan


float Wind_mixing_energy_height_above_ground ( time, height_above_ground2, =
y, x )
units : J
long_name : Wind mixing energy @ Specified height level above ground
missing_value : nan
description : Wind energy generation potential using 0.5*density*speed^3 =
at 80 m above ground level
grid_mapping : LambertConformal_Projection
Grib_Variable_Id : VAR_0-2-19_L103
Grib2_Parameter : ( 0, 2, 19 )
Grib2_Parameter_Discipline : Meteorological products
Grib2_Parameter_Category : Momentum
Grib2_Parameter_Name : Wind mixing energy
Grib2_Level_Type : 103
Grib2_Generating_Process_Type : Forecast
coordinates : time_run time height_above_ground2 y x
_FillValue : nan


float Wind_speed_gust_surface ( time, y, x )
units : m/s
long_name : Wind speed (gust) @ Ground or water surface
missing_value : nan
description : Wind gust speed at 10 m above ground level
grid_mapping : LambertConformal_Projection
Grib_Variable_Id : VAR_0-2-22_L1
Grib2_Parameter : ( 0, 2, 22 )
Grib2_Parameter_Discipline : Meteorological products
Grib2_Parameter_Category : Momentum
Grib2_Parameter_Name : Wind speed (gust)
Grib2_Level_Type : 1
Grib2_Generating_Process_Type : Forecast
coordinates : time_run time y x
_FillValue : nan


float Pressure_surface ( time, y, x )
units : Pa
long_name : Pressure @ Ground or water surface
missing_value : nan
description : Surface pressure
grid_mapping : LambertConformal_Projection
Grib_Variable_Id : VAR_0-3-0_L1
Grib2_Parameter : ( 0, 3, 0 )
Grib2_Parameter_Discipline : Meteorological products
Grib2_Parameter_Category : Mass
Grib2_Parameter_Name : Pressure
Grib2_Level_Type : 1
Grib2_Generating_Process_Type : Forecast
coordinates : time_run time y x
_FillValue : nan


float Geopotential_height_cloud_base ( time, y, x )
units : gpm
long_name : Geopotential height @ Cloud base level
missing_value : nan
description : Adiabatic condensation level above ground (LCL)
grid_mapping : LambertConformal_Projection
Grib_Variable_Id : VAR_0-3-5_L2
Grib2_Parameter : ( 0, 3, 5 )
Grib2_Parameter_Discipline : Meteorological products
Grib2_Parameter_Category : Mass
Grib2_Parameter_Name : Geopotential height
Grib2_Level_Type : 2
Grib2_Generating_Process_Type : Forecast
coordinates : time_run time y x
_FillValue : nan


float Geopotential_height_cloud_tops ( time, y, x )
units : gpm
long_name : Geopotential height @ Level of cloud tops
missing_value : nan
description : Adiabatic condensation level above ground (LCL)
grid_mapping : LambertConformal_Projection
Grib_Variable_Id : VAR_0-3-5_L3
Grib2_Parameter : ( 0, 3, 5 )
Grib2_Parameter_Discipline : Meteorological products
Grib2_Parameter_Category : Mass
Grib2_Parameter_Name : Geopotential height
Grib2_Level_Type : 3
Grib2_Generating_Process_Type : Forecast
coordinates : time_run time y x
_FillValue : nan


float Echo_Top_cloud_tops ( time, y, x )
units : m
long_name : Echo Top @ Level of cloud tops
missing_value : nan
description : Echo top height of 18dBZ surface
grid_mapping : LambertConformal_Projection
Grib_Variable_Id : VAR_0-3-193_L3
Grib2_Parameter : ( 0, 3, 193 )
Grib2_Parameter_Discipline : Meteorological products
Grib2_Parameter_Category : Mass
Grib2_Parameter_Name : Echo Top
Grib2_Level_Type : 3
Grib2_Generating_Process_Type : Forecast
coordinates : time_run time y x
_FillValue : nan


float VIL_Entire_atmosphere ( time, y, x )
units : kg/m2
long_name : VIL @ null
missing_value : nan
description : Vertically integrated liquid
grid_mapping : LambertConformal_Projection
Grib_Variable_Id : VAR_0-3-194_L200
Grib2_Parameter : ( 0, 3, 194 )
Grib2_Parameter_Discipline : Meteorological products
Grib2_Parameter_Category : Mass
Grib2_Parameter_Name : VIL
Grib2_Level_Type : 200
Grib2_Generating_Process_Type : Forecast
coordinates : time_run time y x
_FillValue : nan


float Radar_VIL_Entire_atmosphere ( time, y, x )
units : kg/m2
long_name : Radar VIL @ null
missing_value : nan
description : Radar-derived vertically integrated liquid
grid_mapping : LambertConformal_Projection
Grib_Variable_Id : VAR_0-3-195_L200
Grib2_Parameter : ( 0, 3, 195 )
Grib2_Parameter_Discipline : Meteorological products
Grib2_Parameter_Category : Mass
Grib2_Parameter_Name : Radar VIL
Grib2_Level_Type : 200
Grib2_Generating_Process_Type : Forecast
coordinates : time_run time y x
_FillValue : nan


float Radiation_surface ( time, y, x )
units : W/m2
long_name : Radiation @ Ground or water surface
missing_value : nan
description : Incoming shortwave radiation at the surface
grid_mapping : LambertConformal_Projection
Grib_Variable_Id : VAR_0-4-192_L1
Grib2_Parameter : ( 0, 4, 192 )
Grib2_Parameter_Discipline : Meteorological products
Grib2_Parameter_Category : Short-wave radiation
Grib2_Parameter_Name : Radiation
Grib2_Level_Type : 1
Grib2_Generating_Process_Type : Forecast
coordinates : time_run time y x
_FillValue : nan


float Radiation_atmosphere_top ( time, y, x )
units : W/m2
long_name : Radiation @ Nominal top of the atmosphere
missing_value : nan
description : Outgoing longwave radiation at top of atmosphere
grid_mapping : LambertConformal_Projection
Grib_Variable_Id : VAR_0-5-193_L8
Grib2_Parameter : ( 0, 5, 193 )
Grib2_Parameter_Discipline : Meteorological products
Grib2_Parameter_Category : Long-wave radiation
Grib2_Parameter_Name : Radiation
Grib2_Level_Type : 8
Grib2_Generating_Process_Type : Forecast
coordinates : time_run time y x
_FillValue : nan


float Updraft_Helicity_height_above_ground_layer ( time, height_above_groun=
d_layer, y, x )
units : m2/s2
long_name : Updraft Helicity @ Specified height level above ground layer
missing_value : nan
description : Updraft helicity in layer 2-5 km above ground level
grid_mapping : LambertConformal_Projection
Grib_Variable_Id : VAR_0-7-197_L103_layer
Grib2_Parameter : ( 0, 7, 197 )
Grib2_Parameter_Discipline : Meteorological products
Grib2_Parameter_Category : Thermodynamic stability indices
Grib2_Parameter_Name : Updraft Helicity
Grib2_Level_Type : 103
Grib2_Generating_Process_Type : Forecast
coordinates : time_run time height_above_ground_layer y x
_FillValue : nan


float Reflectivity_height_above_ground ( time, height_above_ground1, y, x =
)
units : dBZ
long_name : Reflectivity @ Specified height level above ground
missing_value : nan
description : Reflectivity at 1 km above ground level
grid_mapping : LambertConformal_Projection
Grib_Variable_Id : VAR_0-16-195_L103
Grib2_Parameter : ( 0, 16, 195 )
Grib2_Parameter_Discipline : Meteorological products
Grib2_Parameter_Category : Forecast radar imagery
Grib2_Parameter_Name : Reflectivity
Grib2_Level_Type : 103
Grib2_Generating_Process_Type : Forecast
coordinates : time_run time height_above_ground1 y x
_FillValue : nan


float Composite_Reflectivity_Entire_atmosphere ( time, y, x )
units : dBZ
long_name : Composite Reflectivity @ null
missing_value : nan
description : Maximum reflectivity in model column
grid_mapping : LambertConformal_Projection
Grib_Variable_Id : VAR_0-16-196_L200
Grib2_Parameter : ( 0, 16, 196 )
Grib2_Parameter_Discipline : Meteorological products
Grib2_Parameter_Category : Forecast radar imagery
Grib2_Parameter_Name : Composite Reflectivity
Grib2_Level_Type : 200
Grib2_Generating_Process_Type : Forecast
coordinates : time_run time y x
_FillValue : nan


float Visibility_surface ( time, y, x )
units : m
long_name : Visibility @ Ground or water surface
missing_value : nan
description : Surface visibility
grid_mapping : LambertConformal_Projection
Grib_Variable_Id : VAR_0-19-0_L1
Grib2_Parameter : ( 0, 19, 0 )
Grib2_Parameter_Discipline : Meteorological products
Grib2_Parameter_Category : Physical atmospheric properties
Grib2_Parameter_Name : Visibility
Grib2_Level_Type : 1
Grib2_Generating_Process_Type : Forecast
coordinates : time_run time y x
_FillValue : nan


float SBT123_height_above_ground_layer ( time, height_above_ground_layer1, =
y, x )
units : m2/s2
long_name : SBT123 @ Specified height level above ground layer
missing_value : nan
description : Updraft helicity in layer 1-6 km above ground level
grid_mapping : LambertConformal_Projection
Grib_Variable_Id : VAR_3-192-1_L103_layer
Grib2_Parameter : ( 3, 192, 1 )
Grib2_Parameter_Discipline : Space products
Grib2_Parameter_Name : SBT123
Grib2_Level_Type : 103
Grib2_Generating_Process_Type : Forecast
coordinates : time_run time height_above_ground_layer1 y x
_FillValue : nan


float u-component_of_current_height_above_ground ( time, height_above_groun=
d2, y, x )
units : m/s
long_name : u-component of current @ Specified height level above ground
missing_value : nan
description : U-component of wind at 80 m above ground level
grid_mapping : LambertConformal_Projection
Grib_Variable_Id : VAR_10-1-2_L103
Grib2_Parameter : ( 10, 1, 2 )
Grib2_Parameter_Discipline : Oceanographic products
Grib2_Parameter_Category : Currents
Grib2_Parameter_Name : u-component of current
Grib2_Level_Type : 103
Grib2_Generating_Process_Type : Forecast
coordinates : time_run time height_above_ground2 y x
_FillValue : nan


float v-component_of_current_height_above_ground ( time, height_above_groun=
d2, y, x )
units : m/s
long_name : v-component of current @ Specified height level above ground
missing_value : nan
description : V-component of wind at 80 m above ground level
grid_mapping : LambertConformal_Projection
Grib_Variable_Id : VAR_10-1-3_L103
Grib2_Parameter : ( 10, 1, 3 )
Grib2_Parameter_Discipline : Oceanographic products
Grib2_Parameter_Category : Currents
Grib2_Parameter_Name : v-component of current
Grib2_Level_Type : 103
Grib2_Generating_Process_Type : Forecast
coordinates : time_run time height_above_ground2 y x
_FillValue : nan







On Wed, Jul 9, 2014 at 12:03 PM, Rick Brownrigg <brownrig@ucar.edu> wrote:


Hi Greg,

I'm not aware of a means to do this from within NCL. If I had to do it, =
I'd use the Proj4 cartographic projection library, which would entail:

- writing the X/Ys to a temp file,
- issue a system() call with a proper command-line containing the L.C. =
projection details to the "proj" command, which would then read the temp =
file and write a new file containing lat/lons
- read the lat/lon

Its a bit clunky, but not too hard *if* you already have Proj4 and are =
familiar with it (http://trac.osgeo.org/proj/)

Hope that helps...
Rick



On Wed, Jul 9, 2014 at 9:51 AM, Greg Garrett - NOAA Federal <greg.garrett@n=
oaa.gov> wrote:


I am working with a Lambert Conformal projected dataset with GeoX and GeoY =
x and y coordinates (excerpt below). There are no latitude or longitude =
coordinates.
Can anyone tell me if it is possible to compute lat/lon from the x,y, and =
projection information and if so, how?

Thanks,

Greg


x = 1799
y = 1059
variables:
integer LambertConformal_Projection ( ncl_scalar )
long_name :
grid_mapping_name : lambert_conformal_conic
latitude_of_projection_origin : 38.5
longitude_of_central_meridian : 262.5
standard_parallel : 38.5
earth_radius : 6367470
_CoordinateTransformType : Projection
_CoordinateAxisTypes : GeoX GeoY

float x ( x )
units : km
long_name : projection_x_coordinate
standard_name : projection_x_coordinate
_CoordinateAxisType : GeoX


float y ( y )
units : km
long_name : projection_y_coordinate
standard_name : projection_y_coordinate
_CoordinateAxisType : GeoY

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk







_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk


Received on Thu Jul 10 02:04:11 2014

This archive was generated by hypermail 2.1.8 : Wed Jul 23 2014 - 15:33:46 MDT