This part of the MODSIM custom code needs some flexibility in the case where we don't want to have convergence checked on certain segments (i.e., segments with no GW-SW interactions because streambed conductance is set to 0):
foreach (DataRow m_Row in m_SyncTblSEG.Rows)// i = 0; i < MF_Acc_Dep.Length; i++)
{
//MF_Acc_Dep_Identifier[i] = (double) m_Row["iseg"];
if (i != (int)((double)m_Row["iseg"] - 1)) throw new Exception("Iseg doesn't match the index of the array");
MS_Flows[i] = 0;
IDivert[i] = (int) (double) m_Row["Diversion"];
IRelease[i] = (int)(double)m_Row["ResRelease"];
m_Link = myModel.FindLink((string)m_Row["Link Name"]); // Use .FindLink() instead
MS_Links[i] = m_Link;
i += 1;
}
Perhaps we can schedule an online mini-sprint to address this, though without the use of webex this might be more difficult than in the past.