def __init__(self, model, npln=1, istrat=1, iswizt=53, nprn=1,
                 toeslope=0.05, tipslope=0.05, zetamin=0.005, delzeta=0.05,
                 nu=0.025, zeta=[], ssz=[], isource=0, extension='swi',
                 fname_output='swi.zta'):
        Package.__init__(self, model)  # Call ancestor's init to set self.parent
        nrow, ncol, nlay, nper = self.parent.nrow_ncol_nlay_nper
        self.unit_number = [29, 53]
        self.extension = extension
        self.file_name = [self.parent.name + '.' + self.extension, fname_output]
        self.name = ['SWI', 'DATA(BINARY)']
        self.heading = '# Salt Water Intrusion package file for MODFLOW-2000, generated by Flopy.'
        self.npln = npln
        self.istrat = istrat
        self.iswizt = iswizt
        self.nprn = nprn
        self.toeslope = toeslope
        self.tipslope = tipslope
        self.zetamin = zetamin
        self.delzeta = delzeta
        # Create arrays so that they have the correct size
        if self.istrat == 1:
            self.nu = Util2d(model, (self.npln + 1,), np.float32, nu, name='nu')
        else:
            self.nu = Util2d(model, (self.npln + 2,), np.float32, nu, name='nu')
        self.zeta = []
        for i in range(nlay):
            # self.zeta.append( empty((nrow, ncol, self.npln)) )
            self.zeta
        for i in range(nlay):
            self.zeta.append(Util2d(model, (self.npln, nrow, ncol), np.float32, zeta[i], name='zeta_' + str(i + 1)))
        self.ssz = Util3d(model, (nlay, nrow, ncol), np.float32, ssz, name='ssz')
        self.isource = Util3d(model, (nlay, nrow, ncol), np.int, isource, name='isource')
        self.parent.add_package(self)
Package constructor.
dataset/ETHPy150Open modflowpy/flopy/flopy/modflow/mfswi.py/ModflowSwi.__init__