Suggestions: | Users can use NCO utilities to create subsets from the "sgpallruc40" data
null
--------------------------------------------------------------
null
null
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
null
null
#!/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
null |