Im trying to read some data from a NetCDF file with a rotated coordinate
system using the function retrieve.nc() from the clim.pact package. I know,
that the rotated pole is supposed to be at 198.0 E and 39.25 N, the top left
corner is at 331.79 E and 21.67 N. The grid resolution is 0.44°x0.44° and the
model operates over an equatorial domain with a quasi- uniform resolution of
approximately 50km.
As you can see in the example below, I would like to extract the cells
(non-rotated) from 7.75 to 8.3 E and 47.75 to 48.3 N (Southern Germany).
I simply cant figure out, how to proceed.
>library(RNetCDF)
>library(ncdf)
>library(date)
>library(ncdf)
>library(clim.pact)
>library(survival)
>var = 'tas'
>file_in <-
>"tas_EUR-44_MPI-M-MPI-ESM-LR_rcp45_r1i1p1_SMHI-RCA4_v1_day_20210101-20251231.nc"
>tas_data <- retrieve.nc(filename=file_in, v.nam= "tas", x.nam= "lon", y.nam=
>"lat", x.rng= c(7.75, 8.30), y.rng= c(47.75, 48.30))
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
