def METHOD_NAME(self):
    """
    """
    nsp = len(self.resistivity_mod)
    # rows, columns in subplots
    # multiply rows by 2 as we are plotting res and phase
    spr = self.subplot_rows * 2
    # number of columns
    spc = 2 * int(np.ceil(float(nsp) / float(self.subplot_rows)))
    # starting subplot number
    row = -1
    for i in range(nsp):
        if (2 * i) % spc == 0:
            row += 1
        spn = 2 * i + 1 + spc * row
        for parameter, meas, mod in [['resistivity', self.resistivity_meas[i], self.resistivity_mod[i]],
                                     ['phase', self.phase_meas[i] % 180, self.phase_mod[i] % 180]]:
            print i
            for ii in range(2):
                for jj in range(2):
                    labels = ''
                    if row == int(spr / 2.) - 1:
                        if parameter == 'phase':
                            labels += 'x'
                    if (2 * i) % spc == 0:
                        if (ii != jj):
                            labels += 'y'
                    if ii == jj:
                        cpt = 0
                        plt.subplot(spr, spc, spn + 1)
                    else:
                        cpt = 1
                        plt.subplot(spr, spc, spn)
                    if ii == 0:
                        c = self.symbol_dict['colors'][0]
                    else:
                        c = self.symbol_dict['colors'][1]
                    plt.METHOD_NAME(self.period_meas[i], meas[:, ii, jj],
                             self.symbol_dict['measured'],
                             color=c)
                    plt.METHOD_NAME(self.period_mod[i], mod[:, ii, jj],
                             self.symbol_dict['modelled'],
                             color=c)
                    self.__set_axis_params(parameter,
                                           component=cpt,
                                           labels=labels)
            spn += spc
    plt.subplots_adjust(hspace=0.02, wspace=0.02)