Dear developers and users,
thanks for your fast and appropriate.
Once again I have a doubt. How can I apply an initial condition (for
example on saturation profile) which depends on the x and y spatial
coordinates of the domain?
//! Evaluates the initial value for a control volume. For this
//! method, the 'values' parameter stores primary variables.
void initial(PrimaryVariables &values,
const Element &element,
const FVElementGeometry &fvGeometry,
int scvIdx) const
{
values[Indices::pwIdx] =200 kPa = 2 bar
values[Indices::snIdx] = *should go here *...
}
Kind regards,
Lorenzo
_______________________________________________
Dumux mailing list
Dumux@listserv.uni-stuttgart.de
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
