Project

General

Profile

Rotated pole gird

Added by Felix Emlinger over 3 years ago

Hello,

I am using climate projection data out of the Climex_project to calculate changes in regional precipitation extremes. In the data set description I cant find any information about the used coordinate system of the data. I just need the precipitation data out of a single grid sell. I used the lon lat coordinates of the center point of my watershed to extract the nearest gridcell (using: remapnn,lon=8.730004_lat=51.491288).

Applying the cdo griddes command gives me following information:

scanningMode = 64
gridtype = projection
gridsize = 78400
xsize = 280
ysize = 280
xname = rlon
xlongname = "longitude in rotated pole grid"
xunits = "degrees"
yname = rlat
ylongname = "latitude in rotated pole grid"
yunits = "degrees"
xfirst = -21.6649932861328
xinc = 0.11000000875056
yfirst = -16.7749977111816
yinc = 0.10999999507781
scanningMode = 64
grid_mapping = rotated_pole
grid_mapping_name = rotated_latitude_longitude
grid_north_pole_latitude = 39.25f
grid_north_pole_longitude = 198.f
north_pole_grid_longitude = 0.f

I don't understand weather I have to reproject my whole files first (in total 360 with 78000 gridpoints) before I can extract the required coordinates or if the rotated_pole is just important when I do use the rlon and rlat values. I already did cut the files and merged them together, right now I am just not 100% sure if I extracted the correct location and want to be sure before I continue with the statistical calculations.

Here the information printed out when using the cdo sinfo command:

File format : NetCDF4 classic zip
-1 : Institut Source T Steptype Levels Num Points Num Dtype : Parameter ID
1 : unknown unknown v instant 1 1 78400 1 F32z : -1
Grid coordinates :
1 : curvilinear : points=78400 (280x280)
lon : -24.98154 to 37.91363 degrees_east
lat : 30.42796 to 64.665 degrees_north
mapping : rotated_latitude_longitude
rlon : -21.66499 to 9.025009 by 0.11 degrees
rlat : -16.775 to 13.915 by 0.11 degrees

Thank you in advance for your help! I just want to be sure that I don't find out that I did all my calculations for a wrong grid cell because I didn't reproject it.

Greetings:
Felix


Replies (2)

RE: Rotated pole gird - Added by Karin Meier-Fleischer over 3 years ago

Hello Felix,

I would guess that the grid cell is the one you are looking for.

Using the following for the first timestep

cdo -remapnn,lon=8.730004_lat=51.491288 -seltimestep,1 pr_EUR-11_CCCma-CanESM2_rcp85_r1-r1i1p1_OURANOS-CRCM5_kba_1h_207205.nc outfile.nc

the latitude value (of the original regional non-rotated grid) is set to lat=51.491288 which is above the maximum value of rlat (of the rotated grid).
Looking at the outfile.nc shows that there is a reasonable value for it

cdo infon outfile.nc 
    -1 :       Date     Time   Level Gridsize    Miss :     Minimum        Mean     Maximum : Parameter name
     1 : 2072-05-01 00:30:00       0        1       0 :              0.00026834             : pr            
cdo    infon: Processed 1 variable over 1 timestep [0.01s 11MB].

-Karin

RE: Rotated pole gird - Added by Felix Emlinger over 3 years ago

Perfect, thank you alot. :) So I will just continuing working with my lon lat values and hope that the precipitation output will be resonable.

Greetings: Felix

    (1-2/2)