Commit ca17e5d4 authored by jouke's avatar jouke
Browse files

Making sure doussan write initial values to tables, avoiding 

extrapolation warnings
parent 4af5f8ae

OpenSimRoot/src/modules/WaterModule/DoussanModel.cpp

+5 −2
Original line number Diff line number Diff line
@@ -504,8 +504,11 @@ void WaterUptakeDoussanModel::writeValues(const Time &t) { 
		if (wuList[i]->evaluateTime(t)) {

			wuList[i]->set(t, jf);

			//std::cout<<std::endl<<t<<' '<<i<<" doussan "<<wuList[i]->getPath();

			if (xpList[i])

				xpList[i]->set(t, psi_x[i]-gravity[i]);

			if (xpList[i]){

				const double p=psi_x[i]-gravity[i];

				xpList[i]->setInit(p);//will only do so if table is empty

				xpList[i]->set(t, p);

			}

		} else {

			if (fabs(jf) > 1e-10)

				msg::warning(