[GRASS-user] Import CORDEX data in GRASS

[Ken Mankoff]

Hi Johannes,

On 2019-10-23 at 10:49 -07, Johannes Radinger <johannesradinger at gmail.com> wrote...
> So maybe another way (if rotated pole is not possible with GRASS/gdal)

Rotated pole works fine in GRASS. Here are my notes from when I did this a few years ago.

 -k.

https://lists.osgeo.org/pipermail/grass-user/2011-October/062180.html

#+BEGIN_SRC bash :results verbatim

# Normal and rotated GRASS locations
trash Gnorm Grot

grass74 -c EPSG:4326 Gnorm

# Import something into the normal location
v.import input=~/data/Zwally_2012/sectors output=Z

cat << EOF > ./Gnorm/PERMANENT/PROJ_INFO
name: General Oblique Transformation
datum: wgs84
towgs84: 0.000,0.000,0.000
proj: ob_tran
o_proj: latlon
ellps: wgs84
a: 6378137.0000000000
es: 0.0066943800
f: 298.2572235630
lat_0: 0.0000000000
lon_0: 180.0000000000
o_lat_p: 18.0
o_lon_p: -200.0
EOF

# rotated_pole:grid_north_pole_latitude = 18. ;
# rotated_pole:grid_north_pole_longitude = -200.

cat << EOF > ./Gnorm/PERMANENT/PROJ_UNITS
unit: degree
units: degrees
meters: .0174532925
EOF

grass74 -e -c EPSG:4326 Grot
grass74 ./Grot/PERMANENT

v.proj location=Gnorm input=Z
g.region vector=Z

d.mon start=wx0
d.erase
d.vect Z
d.grid 1:0:0
#+END_SRC