| Description: | Users who would like access to subsetted RUC-40 data products similar to those products
currently available for RUC-20 can use NCO utilities to create subsets from the
"sgpallruc40" data products. Below is one method that can be used. This is also mentioned in the
FAQ on the ARM RUC description page at: http://www.arm.gov/xds/static/ruc.stm
--------------------------------------------------------------
Question: How can I create subsetted RUC 40 data similar to the RUC 20 data currently
available at the Archive?
Currently, the RUC 40 netCDF data available at the archive are the full area, hourly files
(sgpallruc40isobX1.c1 and sgpallruc40hybrX1.c1). Using the netCDF Operator (NCO)
utilities, it is easy to produce subsetted, daily files similar to the RUC 20 subsets available
(sgpruc20hybrX1.c1 and sgpruc20isobX1.c1).
The NCO utilities can be found at http://nco.sourceforge.net/
Using the ncks (NCO "kitchen sink") command, you can subset the ruc data as follows:
ncks -d nx,, -d ny,,
Note: xmin,xmax,ymin,ymax are index ranges.
For example to subset around the SGP site use the following index values
for latitudes 34.209 to 38.928 and longitudes -99.798 to -95.
xmin = 71
xmax = 82
ymin = 40
ymax = 53
ncks -d nx,71,82 -d ny,40,53 sgpallruc40isobX1.c1.20000314.130000.cdf
sgpruc40isobX1.c1.subset.cdf
Here is a sample C-Shell script for subsetting the ruc40 and using the NCO command ncrcat
to concatenate the hourly subsetted files into a ruc 40 subsetted daily file similar to
the sgpruc20isobX1.c1 files.
#!/bin/csh
# This script subsets the RUC 40 isobaric files
# and then concatenates hourly files into 1 daily file
# usage: subset_ruc40.csh
# where yyyymmdd is the date of the files to subset
# and concatenate. Ex. 20000314 (yyyy=2000,mm=03, dd=14)
# Example: subset_ruc40.csh 20000314
# xmin,xmax,ymin,ymax are indices for
# lat = 34.209, 38.928, lon = -99.798, -95
set xmin = 71
set xmax = 82
set ymin = 40
set ymax = 53
set date = $1
foreach file (sgpallruc40isob*$date*)
ncks -d nx,${xmin},${xmax} -d ny,${ymin},${ymax} $file $file.sub
end
ncrcat -h sgpallruc40isob*$date*sub sgpruc40isobX1.c1.$date.000000.cdf
rm *.sub
exit
For more information on the NCO commands ncks and ncrcat see:
http://nco.sourceforge.net/nco.html#ncks-netCDF-Kitchen-Sink
http://nco.sourceforge.net/nco.html#ncrcat-netCDF-Record-Concatenator
Also other commands for subsetting and analysis of netCDF data can be found at
http://www.arm.gov/data/tools.stm |
| Measurements: | sgpallruc40isobX1.c1: - Categorical Freezing Rain(catfreezerain)
- Dew Point Temperature (2 m fixed height)(tempdewpt2m)
- Snow Depth(snowdepth)
- Pressure Tendency (3-hr)(pressuretend)
- Water Equivalent Snow Depth(snowdepthwatereq)
- Altimeter Setting(altimeter)
- Earth Surface Pressure(pressuresrf)
- Storm-Relative Helicity(stormrelhelicity)
- Relative Humidity at Freezing Level(rhfrz)
- Lifted Index(indexlifted)
- v-Component of Wind (m grid)(windvm)
- Cloud Base Height(heightcldbase)
- 10 meter component u-wind(windu10m)
- v-Component of Wind (p grid)(windvp)
- Potential Temperature (2 m fixed height)(temppot2m)
- Categorical Snow(catsnow)
- Soil Moisture Availability(moistsoilavail)
- 2 meter temperature(temp2m)
- Potential Temperature at Tropopause(temppottrop)
- Pressure at Freezing Level(pressurefrz)
- u-Component of Storm Motion(stormmotionu)
- Vertical Velocity (p grid)(velocitywp)
- v-Component of Storm Motion(stormmotionv)
- Geopotential Height at Top Freezing Level(heightgpfrztop)
- Vertical Velocity (m grid)(velocitywm)
- Sub-grid Scale Precipitation (Accumulated)(precipsubgridacc)
- Dew Point Depression (2 m fixed height)(tempdewptdep2m)
- Relative Humidity (2 m fixed height)(rh2m)
- Categorical Rain(catrain)
- Precipitable water vapr, from MWR (through lssonde)(pwv)
- Altitude(altitude)
- Equilibrium Level Height(heightequillevel)
- Pressures of the p grid(pressurepgrid)
- Environmental Temperature (p grid)(tempp)
- Environmental Temperature (m grid)(tempm)
- north latitude(latitude)
- Specific Humidity (2 m fixed height)(sh2m)
- u-Component of Wind at Tropopause(windutrop)
- Relative Humidity at Top Freezing Level(rhfrztop)
- Best Lifted Index(indexliftedbest)
- Pressure at Top Freezing Level(pressurefrztop)
- Precipitation Rate at Surface(preciprate)
- Geopotential Height (p grid)(heightgpp)
- forecast period(forecast_period)
- Convective Available Potential Energy(cape)
- base time(base_time)
- Resolvable Scale Precipitation (Accumulated)(precipresolvedacc)
- Resolvable Scale Precipitation (3 hr)(precipresolved3hr)
- Equivalent Potential temperature(tempeqpot)
- Sub-surface Runoff(runoffsubsrf)
- Pressure at Tropopause(pressuretrop)
- Cloud Top Height(heightcldtop)
- Gust Wind Speed(windgustspeed)
- Convective Inhibition(cinhib)
- Sub-grid Scale Precipitation (3 hr)(precipsubgrid3hr)
- Snow Accumulation over 3-hr(snowaccum)
- Time offset of tweaks from base_time(time_offset)
- Planetary Boundary Height(heightpbl)
- u-Component of Wind (p grid)(windup)
- Categorical Ice Pellets(caticepellets)
- Surface Runoff Accumulated(runoffsrfacc)
- Pressure at Maximum Wind Level(pressurewindmax)
- east longitude(longitude)
- Geopotential Height at Freezing Level(heightgpfrz)
- Vertical visibility(visibility)
- Convective Cloud Top Height(heightcldtopconv)
- v-Component of Wind at Tropopause(windvtrop)
- u-Component of Wind (m grid)(windum)
- Relative Humidity (m grid)(rhm)
- v-Component of Wind at Maximum Wind Level(windvwindmax)
- Relative Humidity (p grid)(rhp)
- Pressure Intervals (30-0,60-30,90-60,120-90,150-120,180-150) of the m grid(pressuremgrid)
- u-Component of Wind at Maximum Wind Level(winduwindmax)
- MAPS Pressure at Mean Sea Level(pressuremsl)
- Earth Surface Height(heightsrf)
- 10 meter component v-wind(windv10m)
- Maximum Equivalent Potential Temperature Pressure(pressuretempeqpotmax)
sgpallruc40hybrX1.c1: - Precipitation Rate at Surface(preciprate)
- Water Condensation Rate near Surface(condensewaterratesrf)
- Mixing Ratio of Water Vapor (hybrid-b grid)(mxratiovapor)
- forecast period(forecast_period)
- Geopotential Height (hybrid-b grid)(heightgp)
- Pressure Vertical Velocity (hybrid-b grid)(velocityw)
- Pressure at Tropopause(pressuretrop)
- Mixing Ratio of Snow (hybrid-b grid)(mxratiosnow)
- Resolvable Scale Precipitation at Surface(precipresolved)
- north latitude(latitude)
- Mixing Ratio of Rain Water (hybrid-b grid)(mxratiorain)
- v-Component of Wind at Tropopause(windvtrop)
- Altitude(altitude)
- Sub-surface Runoff(runoffsubsrf)
- Net Longwave Irradiance at Surface(fluxlwnetsrf)
- u-Component of Wind at Tropopause(windutrop)
- Surface,Sub-surface Soil Volumetric Moisture(moistsoilvolumesrfsubsrf)
- Surface Runoff(runoffsrf)
- Rate of Water Dropping from Canopy to Ground(dropwaterratesrf)
- Evaporation of Canopy Water near Surface(evapcanopysrf)
- Latent Heat Flux at Surface(fluxlatentsrf)
- Evaporation from Bare Soil at Surface(evapsoilsrf)
- east longitude(longitude)
- Convective Available Potential Energy(cape)
- Surface/Sub-surface Heights(heightss)
- Precipitable water vapr, from MWR (through lssonde)(pwv)
- Number Concentration of Ice Particles (hybrid-b grid)(ndensice)
- Mixing Ratio of Ice Water (hybrid-b grid)(mxratioice)
- Sensible Heat Flux at Surface(fluxsensiblesrf)
- Convective Inhibition(cinhib)
- Virtual Potential Temperature (hybrid-b grid)(tempvirtpot)
- Sub-grid Scale Precipitation at Surface(precipsubgrid)
- Net Shortwave Irradiance at Surface(fluxswnetsrf)
- Plant Canopy Surface Water(watercanopy)
- Relative humidity inside the instrument enclosure(rh)
- Snow Accumulation over 3-hr(snowaccum)
- Turbulent Kinetic Energy (hybrid-b grid)(tketurb)
- v-Component of Wind (hybrid-b grid)(windv)
- u-Component of Wind (hybrid-b grid)(windu)
- Potential Temperature at Tropopause(temppottrop)
- Evapotranspiration Flux near Surface(fluxevaptranssrf)
- Snow Temperature at Snow or Top Soil Level(snowtemp)
- Mixing Ratio of Graupel (hybrid-b grid)(mxratiograupel)
- base time(base_time)
- Surface,Sub-surface Soil Moisture Content(moistsoilcontentsrfsubsrf)
- Mixing Ratio of Cloud Water (hybrid-b grid)(mxratioliquid)
- Time offset of tweaks from base_time(time_offset)
- Beam 0 Temperature(temp)
- Surface,Sub-surface Soil Temperature(tempsoilsrfsubsrf)
- MAPS Pressure at Mean Sea Level(pressuremsl)
- Snow Depth(snowdepth)
- atmospheric pressure at mean sea level and at tropopause (NGM250 predicted)(pressure)
|