I’m 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 can’t 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]]