{"action":"opened","issue":{"url":"https://api.github.com/repos/modflowpy/flopy/issues/48","labels_url":"https://api.github.com/repos/modflowpy/flopy/issues/48/labels{/name}","comments_url":"https://api.github.com/repos/modflowpy/flopy/issues/48/comments","events_url":"https://api.github.com/repos/modflowpy/flopy/issues/48/events","html_url":"https://github.com/modflowpy/flopy/issues/48","id":104289199,"number":48,"title":"SEAWAT; Itype==-1 and itype==-15 not counted as sink/sources","user":{"login":"bdestombe","id":11296133,"avatar_url":"https://avatars.githubusercontent.com/u/11296133?v=3","gravatar_id":"","url":"https://api.github.com/users/bdestombe","html_url":"https://github.com/bdestombe","followers_url":"https://api.github.com/users/bdestombe/followers","following_url":"https://api.github.com/users/bdestombe/following{/other_user}","gists_url":"https://api.github.com/users/bdestombe/gists{/gist_id}","starred_url":"https://api.github.com/users/bdestombe/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/bdestombe/subscriptions","organizations_url":"https://api.github.com/users/bdestombe/orgs","repos_url":"https://api.github.com/users/bdestombe/repos","events_url":"https://api.github.com/users/bdestombe/events{/privacy}","received_events_url":"https://api.github.com/users/bdestombe/received_events","type":"User","site_admin":false},"labels":[],"state":"open","locked":false,"assignee":null,"milestone":null,"comments":0,"created_at":"2015-09-01T15:14:59Z","updated_at":"2015-09-01T15:14:59Z","closed_at":null,"body":"The cells with itype==-1 and itype==15 are counted by the SEAWAT routine as sink/sources. Not by flopy. \r\n\r\nin flopy\\mt3d\\mtssm.py lines: 50 - 59, to:\r\n```python\r\nself.stress_period_data = mflist(self, model=self.parent.mf,\r\n                                         data=stress_period_data)\r\n\r\nself.__maxssm = np.sum(self.stress_period_data.data[0].itype == -1)\r\nself.__maxssm += np.sum(self.stress_period_data.data[0].itype == -15)\r\n\r\n#        self.__maxssm = 0\r\n#if (self.parent.btn.icbund != None):\r\n\r\nif isinstance(self.parent.btn.icbund, np.ndarray):\r\n\tself.maxssm += (self.parent.btn.icbund < 0).sum()\r\nfor p in self.__SsmPackages:\r\n\tprint 'label', p.label\r\n\tprint 'instance', p.instance\r\n\tif ((p.label == 'BAS6') and (p.instance != None)):\r\n\t\tself.__maxssm += (p.instance.ibound < 0).sum()\r\n\telif p.instance != None:\r\n\t\tself.__maxssm += p.instance.ncells()\r\n```\r\n\r\nAlso moved the function call that assigns the stress_period_data variable upwards. That garantees that the stress_period_data is in mflist format."}}