Fixed bug in SBS apparent density plot variable.

SteveMaas1978 · SteveMaas1978 · commit a8dea2ca0492 · 2026-03-02T16:16:52.000-07:00
diff --git a/FEBioChem/FEBioChemPlot.cpp b/FEBioChem/FEBioChemPlot.cpp
@@ -396,12 +396,13 @@ bool FEChemPlotSBSApparentDensity::Save(FEDomain &dom, FEDataStream& a)
 	FEChemReactionDomain& rdom = static_cast<FEChemReactionDomain&>(dom);
 
 	FEChemReactionDiffusionMaterial* mat = dynamic_cast<FEChemReactionDiffusionMaterial*>(rdom.GetMaterial());
+	int nsol = mat->Species();
 	int nsbm = m_sbmName.size();
 	vector<int> lid(nsbm, -1);
 	for (int i = 0; i < nsbm; ++i)
 	{
 		FEChemReactiveSpeciesBase* rs = mat->FindSpecies(m_sbmName[i]);
-		if (rs) lid[i] = rs->GetLocalID();
+		if (rs) lid[i] = rs->GetLocalID() - nsol;
 	}
 
 	if ((nsbm == 1) && (lid[0] == -1)) return false;

Original file line number	Diff line number	Diff line change
`@@ -396,12 +396,13 @@ bool FEChemPlotSBSApparentDensity::Save(FEDomain &dom, FEDataStream& a)`
`396`	`396`	`FEChemReactionDomain& rdom = static_cast<FEChemReactionDomain&>(dom);`
`397`	`397`
`398`	`398`	`FEChemReactionDiffusionMaterial* mat = dynamic_cast<FEChemReactionDiffusionMaterial*>(rdom.GetMaterial());`
	`399`	`+ int nsol = mat->Species();`
`399`	`400`	`int nsbm = m_sbmName.size();`
`400`	`401`	`vector<int> lid(nsbm, -1);`
`401`	`402`	`for (int i = 0; i < nsbm; ++i)`
`402`	`403`	`{`
`403`	`404`	`FEChemReactiveSpeciesBase* rs = mat->FindSpecies(m_sbmName[i]);`
`404`		`- if (rs) lid[i] = rs->GetLocalID();`
	`405`	`+ if (rs) lid[i] = rs->GetLocalID() - nsol;`
`405`	`406`	`}`
`406`	`407`
`407`	`408`	`if ((nsbm == 1) && (lid[0] == -1)) return false;`