netcdf ngm-q { // NCEP Early Domestic Product Set from NGM dimensions: record = UNLIMITED ; // (reference time, forecast time) level = 19 ; // isobaric levels for most parameters level_c = 4 ; // isobaric levels for absolute vorticity level_o = 9 ; // isobaric levels for omega liso = 1 ; // layer between two isobaric levels ls_all = 1 ; // whole atmosphere layer x = 93 ; y = 65 ; datetime_len = 21 ; // string length for datetime strings accum = 2 ; // time range for accumulations nav = 1 ; // For navigation. Variables that use // this dimension define a mapping between // (x, y) indices and (lat, lon) coords. nav_len = 100 ; // max length for navigation strings variables: double reftime(record) ; // reference time of the model reftime:long_name = "reference time" ; reftime:units = "hours since 1992-1-1" ; double valtime(record) ; // forecast time ("valid" time) valtime:long_name = "valid time" ; valtime:units = "hours since 1992-1-1" ; :record = "reftime, valtime" ; // "dimension attribute" -- means // (reftime, valtime) uniquely // determine record char datetime(record, datetime_len) ; // derived from reftime datetime:long_name = "reference date and time" ; // units YYYY-MM-DD hh:mm:ssZ (ISO 8601) float valtime_offset(record) ; // derived as valtime-reftime valtime_offset:long_name = "hours from reference time" ; valtime_offset:units = "hours" ; float level(level) ; level:long_name = "level" ; level:units = "hectopascals" ; float level_c(level_c) ; level_c:long_name = "level subset" ; level_c:units = "hectopascals" ; float level_o(level_o) ; level_o:long_name = "level subset" ; level_o:units = "hectopascals" ; :liso = "liso_bot, liso_top" ; // (liso_bot, liso_top) uniquely // determines liso float liso_bot(liso) ; liso_bot:long_name = "bottom level of layer between 2 isobaric levels" ; liso_bot:units = "hectopascals" ; float liso_top(liso) ; liso_top:long_name = "top level of layer between 2 isobaric levels" ; liso_top:units = "hectopascals" ; :ls_all = "ls_all_bot, ls_all_top" ; float ls_all_bot(ls_all) ; ls_all_bot:long_name = "bottom level of atmosphere between 2 sigma levels" ; ls_all_bot:units = "" ; float ls_all_top(ls_all) ; ls_all_top:long_name = "top level of atmosphere between 2 sigma levels" ; ls_all_top:units = "" ; long model_id ; model_id:long_name = "generating process ID number" ; // navigation variables all use nav dimension char nav_model(nav, nav_len) ; // navigation parameterization nav_model:long_name = "navigation model name" ; int grid_type_code(nav) ; grid_type_code:long_name = "GRIB-1 GDS data representation type" ; char grid_type(nav, nav_len) ; grid_type:long_name = "GRIB-1 grid type" ; char grid_name(nav, nav_len) ; grid_name:long_name = "grid name" ; int grid_center(nav) ; grid_center:long_name = "GRIB-1 originating center ID" ; int grid_number(nav) ; grid_number:long_name = "GRIB-1 catalogued grid number" ; char earth_shape(nav, nav_len) ; earth_shape:long_name = "assumed earth shape" ; char x_dim(nav, nav_len) ; x_dim:long_name = "x dimension name" ; char y_dim(nav, nav_len) ; y_dim:long_name = "y dimension name" ; int Nx(nav) ; Nx:long_name = "number of points along x-axis" ; int Ny(nav) ; Ny:long_name = "number of points along y-axis" ; float La1(nav) ; La1:long_name = "latitude of first grid point" ; La1:units = "degrees_north" ; float Lo1(nav) ; Lo1:long_name = "longitude of first grid point" ; Lo1:units = "degrees_east" ; byte ResCompFlag(nav) ; ResCompFlag:long_name = "resolution and component flags" ; float Lov(nav) ; Lov:long_name = "orientation of the grid" ; Lov:units = "degrees_east" ; float Dx(nav) ; Dx:long_name = "x-direction grid length" ; Dx:units = "m" ; float Dy(nav) ; Dy:long_name = "y-direction grid length" ; Dy:units = "m" ; byte ProjFlag(nav) ; ProjFlag:long_name = "projection center flag" ; float Latin1(nav) ; Latin1:long_name = "first intersecting latitude" ; Latin1:units = "degrees_north" ; float Latin2(nav) ; Latin2:long_name = "second intersecting latitude" ; Latin2:units = "degrees_north" ; float SpLat(nav) ; SpLat:long_name = "Latitude of the southern pole" ; SpLat:units = "degrees_north" ; float SpLon(nav) ; SpLon:long_name = "Longitude of the southern pole" ; SpLon:units = "degrees_east" ; // end of navigation variables float LI(record, liso, y, x) ; // Lifted index between two levels LI:long_name = "lifted index" ; LI:units = "degK" ; LI:_FillValue = -9999.f ; // To fill grid corners LI:navigation = "nav" ; float P_trop(record, y, x) ; P_trop:long_name = "pressure at tropopause" ; P_trop:units = "Pa" ; P_trop:_FillValue = -9999.f ; P_trop:navigation = "nav" ; float P_sfc(record, y, x) ; P_sfc:long_name = "pressure at surface" ; P_sfc:units = "Pa" ; P_sfc:_FillValue = -9999.f ; P_sfc:navigation = "nav" ; float PRECIP(record, y, x) ; PRECIP:long_name = "total precipitation over accumulation interval" ; PRECIP:units = "kg/m2" ; PRECIP:begin_or_end = "end" ; PRECIP:_FillValue = -9999.f ; PRECIP:navigation = "nav" ; float PRECIP_accum_times(record, accum) ; PRECIP_accum_times:long_name = "precipitation accumulation interval" ; PRECIP_accum_times:units = "hours" ; PRECIP_accum_times:_FillValue = -9999.f ; float precip_cn(record, y, x) ; precip_cn:long_name = "convective precipitation" ; precip_cn:units = "kg/m2" ; precip_cn:_FillValue = -9999.f ; precip_cn:navigation = "nav" ; float precip_cn_accum_times(record, accum) ; precip_cn_accum_times:long_name = "convective precipitation accumulation interval" ; precip_cn_accum_times:units = "hours" ; precip_cn_accum_times:_FillValue = -9999.f ; float P_msl(record, y, x) ; P_msl:long_name = "pressure reduced to MSL" ; P_msl:units = "Pa" ; P_msl:_FillValue = -9999.f ; P_msl:navigation = "nav" ; float RH(record, level, y, x) ; RH:long_name = "relative humidity" ; RH:units = "percent" ; RH:_FillValue = -9999.f ; RH:navigation = "nav" ; float T(record, level, y, x) ; T:long_name = "temperature" ; T:units = "degK" ; T:_FillValue = -9999.f ; T:navigation = "nav" ; float T_trop(record, y, x) ; T_trop:long_name = "temperature at tropopause" ; T_trop:units = "degK" ; T_trop:_FillValue = -9999.f ; T_trop:navigation = "nav" ; float Z(record, level, y, x) ; Z:long_name = "geopotential height" ; Z:units = "gp m" ; Z:_FillValue = -9999.f ; Z:navigation = "nav" ; // georeference info float Z_sfc(y, x) ; Z_sfc:long_name = "terrain" ; Z_sfc:units = "gp m" ; Z_sfc:_FillValue = -9999.f ; Z_sfc:navigation = "nav" ; float omega(record, level_o, y, x) ; omega:long_name = "pressure vertical velocity" ; omega:units = "Pa/s" ; omega:_FillValue = -9999.f ; omega:navigation = "nav" ; float u(record, level, y, x) ; u:long_name = "u-component of wind" ; u:units = "meters/second" ; u:_FillValue = -9999.f ; u:navigation = "nav" ; float u_trop(record, y, x) ; // Eastward wind at tropopause u_trop:long_name = "u-component of wind at tropopause" ; u_trop:units = "meters/second" ; u_trop:_FillValue = -9999.f ; // To fill grid corners u_trop:navigation = "nav" ; float v(record, level, y, x) ; v:long_name = "v-component of wind" ; v:units = "meters/second" ; v:_FillValue = -9999.f ; v:navigation = "nav" ; float v_trop(record, y, x) ; // Northward wind at tropopause v_trop:long_name = "v-component of wind at tropopause" ; v_trop:units = "meters/second" ; v_trop:_FillValue = -9999.f ; // To fill grid corners v_trop:navigation = "nav" ; float pr_water_ls(record, ls_all, y, x) ; pr_water_ls:long_name = "precipitable water" ; pr_water_ls:units = "kg/m2" ; pr_water_ls:_FillValue = -9999.f ; pr_water_ls:navigation = "nav" ; float absvor(record, level_c, y, x) ; absvor:long_name = "absolute vorticity" ; absvor:units = "1/s" ; absvor:_FillValue = -9999.f ; absvor:navigation = "nav" ; // global attributes: :history = "created by gribtonc from HDS broadcast" ; :title = "NCEP Early Domestic Product Set" ; :Conventions = "NUWG" ; :version = 0.0 ; // still just a draft data: level = 1000, 950, 900, 850, 800, 750, 700, 650, 600, 550, 500, 450, 400, 350, 300, 250, 200, 150, 100 ; level_c = 850, 700, 500, 250 ; level_o = 850, 700, 500, 400, 300, 250, 200, 150, 100 ; liso_bot = 0 ; liso_top = 500 ; ls_all_top = 0.0 ; ls_all_bot = 1.0 ; model_id = 39 ; // Navigation nav_model = "GRIB1" ; grid_type_code = 3 ; grid_type = "Lambert conformal" ; grid_name = "AWIPS grid 211: Regional CONUS" ; grid_center = 7 ; grid_number = 211 ; earth_shape = "oblate spheroid (IAU 1965)" ; x_dim = "x" ; y_dim = "y" ; Nx = 93 ; Ny = 65 ; La1 = 12.190 ; Lo1 = -133.459 ; ResCompFlag = 8 ; Lov = -95.0 ; Dx = 81271 ; Dy = 81271 ; ProjFlag = 0 ; Latin1 = 25.0 ; Latin2 = 25.0 ; SpLat = 0.0 ; SpLon = 0.0 ; }