nrows_ncols=(3, 3),
                 axes_pad=0,
                 add_all=True,
                 share_all=True,
                 label_mode="L",
                 cbar_location="top",
                 cbar_mode="single",
                 cbar_size='2%',
                 aspect=True)

# homedir = '/media/raul/RauLdisk'
#homedir = '/Volumes/RauLdisk'
homedir = '/localdata'
cfsr.plot(ax=grid,
          field=['iwv_flux',range(250, 1500, 250)],
          contour=['thetaeq',None],
          dates=dates,
          homedir=homedir)

grid.axes_all[-1].axis('off')
grid.axes_all[-2].axis('off')

labticksize=11

yticks = grid.axes_all[0].yaxis.get_major_ticks()
yticks[0].label1.set_visible(False)
yticks[-1].label1.set_visible(False)
for tick in yticks:
    tick.label.set_fontsize(labticksize)

yticks = grid.axes_all[3].yaxis.get_major_ticks()
                 cbar_location="top",
                 cbar_mode="single",
                 cbar_size='2%',
                 aspect=True)

''' set tile '''
title  = 'Integrated water vapor transport $[kg m^{-1} s^{-1}]$\n '
title += 'Mean sea level pressure [hPa]\n'
title += 'Equivalent potential temperature [K] at 1000 hPa'

# homedir = '/media/raul/RauLdisk'
#homedir = '/Volumes/RauLdisk'
homedir = '/localdata'
cf = cfsr.plot(ax     = grid,
              field   = ['iwv_flux',range(250, 1500, 250)],
              contour = ['thetaeq',range(304, 340, 2)],
              dates   = dates,
              homedir = homedir,
              title   = '')

''' add title at correct position '''
ax = grid.axes_all[1]
ax.text(0.5,1.2,title,ha='center',
        transform=ax.transAxes)

''' adjust y and x labels '''
for n in [0,3]:
    yticks = grid.axes_all[n].yaxis.get_major_ticks()
    yticks[0].label1.set_visible(False)
    yticks[-1].label1.set_visible(False)

for n in [3,4,5]: