I am fairly new at using R and I am currently trying to use it to extract data for a specific region of a netCDF file and convert it to a csv. Ultimately I would like to have a timeseries for each variable over my area of interest.
I've downloaded the packages ncdf
and clim.pact
for R, and as my netCDF file is quite large (daily data for 25 years) I think I should be able to extract a the data for my region of interest using cdfextract, and then I will convert that smaller file to a csv. Howver, I keep getting the same error when I try to run it. This is the script I've been using and the error that comes up:
cdfextract("file.nc","tas", x.rng=c(-93, -85), y.rng=c(16, 23),t.rng=c(1, 9125))
[1] "time.origin= 01 Jan 0850" Error in system("rm -f cdfcont.txt", intern = TRUE) : 'rm' not found
Where my file is file.nc
and the variable of interest is tas
. I think it may have something to do with the time format, but when I check the time units of the file it says it is in "days since 01 Jan 0850" so I think the units are right.
Anyways any comments on this issue would be very welcome! Thank you!
If you don't have to use R, you can use
ncks
from the NCO package: http://nco.sourceforge.net/nco.html#ncks-netCDF-Kitchen-Sinkncks
is a wonderful command-line tool for extracting data from netCDF files. In general, you can extract a variable over certain lat/lon bands and then output that data into a .csv file using something like:I couldn't tell if you wanted an area-average of
tas
first and then a time series of that area-averaged quantity in the .csv file. If so, you'll need to usencwa
to perform the area averaging first: http://nco.sourceforge.net/nco.html#Averaging