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]]

