partsim with first diff analysis
add water levels init state