RICHARDS_FLOW or DEFORMATION? In my opinion, it works well with DEFORMATION only.
*************************************************************************
GeoSys-FEM Function:
06/2006 YD Implementation
02/2008 JOD removed
03/2008 HS/KG activated for adaptive time step
Reload primary variable
**************************************************************************/
void CRFProcess::PrimaryVariableReload()
{
// char pcsT;
// pcsT = _pcs_type_name[0];
// switch(pcsT){
// case 'L':
// break;
// case 'U':
// break;
// case 'G':
// break;
// case 'T':
// break;
// case 'C':
// break;
// case 'M':
// PrimaryVariableReloadTransport();
// break;
// case 'R': // Richards flow
// PrimaryVariableReloadRichards();
// break;
// }
switch (this->getProcessType())
{
case FiniteElement::MASS_TRANSPORT:
PrimaryVariableReloadTransport();
break;
case FiniteElement::RICHARDS_FLOW:
PrimaryVariableReloadRichards();
break;
default:
break;
}
}
Moreover,I am not sure what do you mean by
$RELOAD
3 100
if (reload == 1 || reload == 3)
$RELOAD
0 Do Nothing
1 Read
2 Read what has already written
3 ?
Regards,
Mohammad