def plot(self, filename_base=None, file_extension=None, **kwargs):
        import flopy.plot.plotutil as pu

        if file_extension is not None:
            fext = file_extension
        else:
            fext = 'png'

        if 'kper' in kwargs:
            kk = kwargs['kper']
            kwargs.pop('kper')
            try:
                kk = kk.lower()
                if kk == 'all':
                    k0 = 0
                    k1 = self.model.nper
                else:
                    k0 = 0
                    k1 = 1
            except:
                k0 = int(kk)
                k1 = k0 + 1
                # if kwargs['kper'] == 'all':
                #     kwargs.pop('kper')
                #     k0 = 0
                #     k1 = self.model.nper
                # else:
                #     k0 = int(kwargs.pop('kper'))
                #     k1 = k0 + 1
        else:
            k0 = 0
            k1 = 1

        if 'fignum' in kwargs:
            fignum = kwargs.pop('fignum')
        else:
            fignum = list(range(k0, k1))

        if 'mflay' in kwargs:
            kwargs.pop('mflay')

        axes = []
        for idx, kper in enumerate(range(k0, k1)):
            title = '{} stress period {:d}'. \
                format(self.name_base.replace('_', '').upper(),
                       kper + 1)
            if filename_base is not None:
                filename = filename_base + '_{:05d}.{}'.format(kper + 1, fext)
            else:
                filename = None
            axes.append(pu._plot_array_helper(self[kper].array, self.model,
                                              names=title, filenames=filename,
                                              fignum=fignum[idx], **kwargs))
        return axes
Plot transient 2-D model input data
dataset/ETHPy150Open modflowpy/flopy/flopy/utils/util_array.py/Transient2d.plot