Compare regridded elevation data

Lee Murray regridded a high resolution USGS elevation map to the GEOS-Chem grids (the files are in IDL save format). We compare these regridded elevation data with results from the Resample_GMTED2010_2_GCGrid notebook.

1. Imports


In [1]:
import os

import numpy as np
from scipy.io import readsav

import iris
import iris.quickplot as qiplot

2. Paths to files


In [2]:
gmted2010_awm_path = '/home/bovy/Grids/'
lm_path = '/home/bovy/Grids/'

3. Load Data and compare


In [3]:
gmted2010_2x15_cube = iris.load_cube(os.path.join(gmted2010_awm_path, 'dem_GEOS57_2x2.5_awm.nc'))
gmted2010_4x5_cube = iris.load_cube(os.path.join(gmted2010_awm_path, 'dem_GEOS57_4x5_awm.nc'))

lm_2x15 = readsav('/home/mahieu/geos.runs/elev.2x25.sav')
lm_4x5 = readsav('/home/mahieu/geos.runs/elev.4x5.sav')

lm_2x15_cube = gmted2010_2x15_cube.copy()
lm_2x15_cube.data = lm_2x15['elev2x25_m']
lm_2x15_cube.var_name = "LM-ELEV_2x25"

lm_4x5_cube = gmted2010_4x5_cube.copy()
lm_4x5_cube.data = lm_4x5['elev4x5_m']
lm_4x5_cube.var_name = "LM-ELEV_4x5"

In [4]:
diff_2x15_cube = gmted2010_2x15_cube - lm_2x15_cube
diff_2x15_cube.var_name = "Diff_GMTED2010-LM_2x25"

diff_4x5_cube = gmted2010_4x5_cube - lm_4x5_cube
diff_4x5_cube.var_name = "Diff_GMTED2010-LM_4x5"

print diff_2x15_cube
print diff_4x5_cube


Diff_GMTED2010-LM_2x25 / (m)        (latitude: 91; longitude: 144)
     Dimension coordinates:
          latitude                           x              -
          longitude                          -              x
Diff_GMTED2010-LM_4x5 / (m)         (latitude: 46; longitude: 72)
     Dimension coordinates:
          latitude                           x              -
          longitude                          -              x

4. Plot


In [5]:
%pylab inline

qiplot.pcolormesh(diff_2x15_cube, vmin=-200, vmax=200, cmap=cm.RdBu)
plt.gca().coastlines()
plt.show()


Populating the interactive namespace from numpy and matplotlib

In [6]:
qiplot.pcolormesh(diff_4x5_cube, vmin=-200, vmax=200, cmap=cm.RdBu)
plt.gca().coastlines()
plt.show()



In [6]: