* \param [in] int : int corresponding to the enum CELLS, NODES or FACES
* \param [out] void
* */
- void setInitialField(string fileName, string fieldName, int timeStepNumber, int field_support_type);
+ void setInitialField(string fileName, string fieldName, int timeStepNumber, int order, int meshLevel, int field_support_type);
/** \fn setInitialField
* \brief sets the initial field from a field in a med file
* \param [in] EntityType : CELLS, NODES or FACES
* \param [out] void
* */
- void setInitialField(string fileName, string fieldName, int timeStepNumber, EntityType typeField = CELLS);
+ void setInitialField(string fileName, string fieldName, int timeStepNumber, int order = 0, int meshLevel=0, EntityType typeField = CELLS);
/** \fn setInitialFieldConstant
* \brief sets a constant initial field on a mesh stored in a med file
}
//Function needed because swig of enum EntityType fails
-void ProblemCoreFlows::setInitialField(string fileName, string fieldName, int timeStepNumber, int field_support_type)
+void ProblemCoreFlows::setInitialField(string fileName, string fieldName, int timeStepNumber, int order, int meshLevel, int field_support_type)
{
if(_FECalculation && field_support_type!= NODES)
cout<<"Warning : finite element simulation should have initial field on nodes!!!"<<endl;
switch(field_support_type)
{
case CELLS:
- VV = Field(fileName, CELLS, fieldName, timeStepNumber, 0);
+ VV = Field(fileName, CELLS, fieldName, timeStepNumber, order, meshLevel);
break;
case NODES:
- VV = Field(fileName, NODES, fieldName, timeStepNumber, 0);
+ VV = Field(fileName, NODES, fieldName, timeStepNumber, order, meshLevel);
break;
case FACES:
- VV = Field(fileName, FACES, fieldName, timeStepNumber, 0);
+ VV = Field(fileName, FACES, fieldName, timeStepNumber, order, meshLevel);
break;
default:
std::ostringstream message;
setInitialFieldConstant( nDim, Vconstant, xmin, xmax, nx, leftSide, rightSide, ymin, ymax, ny, backSide, frontSide, zmin, zmax, nz, bottomSide, topSide, typeField);
}
-void ProblemCoreFlows::setInitialField(string fileName, string fieldName, int timeStepNumber, EntityType typeField)
+void ProblemCoreFlows::setInitialField(string fileName, string fieldName, int timeStepNumber, int order, int meshLevel, EntityType typeField)
{
if(_FECalculation && typeField!= NODES)
cout<<"Warning : finite element simulation should have initial field on nodes!!!"<<endl;
- Field VV(fileName, typeField, fieldName, timeStepNumber, 0);
+ Field VV(fileName, typeField, fieldName, timeStepNumber, order, meshLevel);
setInitialField(VV);
}