def METHOD_NAME(weatherObj, savefig=True, z1=500, z2=15000):
    '''
    Create a plot with wet refractivity and hydrostatic refractivity,
    at two different heights
    '''
    # Get the interpolator
    intFcn_w = Interpolator((weatherObj._xs, weatherObj._ys, weatherObj._zs), weatherObj._wet_refractivity.swapaxes(0, 1))
    intFcn_h = Interpolator((weatherObj._xs, weatherObj._ys, weatherObj._zs), weatherObj._hydrostatic_refractivity.swapaxes(0, 1))
    # get the points needed
    XY = np.meshgrid(weatherObj._xs, weatherObj._ys)
    x = XY[0]
    y = XY[1]
    z1a = np.zeros(x.shape) + z1
    z2a = np.zeros(x.shape) + z2
    pts1 = np.stack((x.flatten(), y.flatten(), z1a.flatten()), axis=1)
    pts2 = np.stack((x.flatten(), y.flatten(), z2a.flatten()), axis=1)
    w1 = intFcn_w(pts1)
    h1 = intFcn_h(pts1)
    w2 = intFcn_w(pts2)
    h2 = intFcn_h(pts2)
    # Now get the data to plot
    plots = [w1, h1, w2, h2]
    # titles
    titles = ('Wet refractivity {}'.format(z1),
              'Hydrostatic refractivity {}'.format(z1),
              '{}'.format(z2),
              '{}'.format(z2))
    # setup the plot
    f = plt.figure(figsize=(14, 10))
    f.suptitle(f'{weatherObj._Name} Wet and Hydrostatic refractivity at height {z1}m and {z2}m')
    # loop over each plot
    for ind, plot, title in zip(range(len(plots)), plots, titles):
        sp = f.add_subplot(2, 2, ind + 1)
        im = sp.imshow(np.reshape(plot, x.shape), cmap='viridis',
                       extent=[np.nanmin(x), np.nanmax(x), np.nanmin(y), np.nanmax(y)], origin='lower')
        divider = mal(sp)
        cax = divider.append_axes("right", size="4%", pad=0.05)
        plt.colorbar(im, cax=cax)
        sp.set_title(title)
        if ind == 0:
            sp.set_ylabel('{} m\n'.format(z1))
        if ind == 2:
            sp.set_ylabel('{} m\n'.format(z2))
    if savefig:
        wd   = os.path.dirname(os.path.dirname(weatherObj._out_name))
        f    = f'{weatherObj._Name}_refractivity_hgt{z1}_and_{z2}m.pdf'
        plt.savefig(os.path.join(wd, f))
    return f