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