Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/Balancethickness/CreateConstraintsBalancethickness.cpp =================================================================== --- /issm/trunk-jpl/src/c/modules/ModelProcessorx/Balancethickness/CreateConstraintsBalancethickness.cpp (revision 15464) +++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/Balancethickness/CreateConstraintsBalancethickness.cpp (revision 15465) @@ -14,7 +14,4 @@ Constraints* constraints=*pconstraints; - /*Create constraints if they do not exist yet*/ - if(!constraints) constraints = new Constraints(); - /*Do not add constraints in DG*/ if(stabilization!=3){ Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/Balancethickness/CreateLoadsBalancethickness.cpp =================================================================== --- /issm/trunk-jpl/src/c/modules/ModelProcessorx/Balancethickness/CreateLoadsBalancethickness.cpp (revision 15464) +++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/Balancethickness/CreateLoadsBalancethickness.cpp (revision 15465) @@ -18,7 +18,4 @@ /*Recover pointer: */ Loads* loads=*ploads; - - /*Create loads if they do not exist yet*/ - if(!loads) loads = new Loads(); /*Loads only in DG*/ Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/BedSlope/CreateConstraintsBedSlope.cpp =================================================================== --- /issm/trunk-jpl/src/c/modules/ModelProcessorx/BedSlope/CreateConstraintsBedSlope.cpp (revision 15464) +++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/BedSlope/CreateConstraintsBedSlope.cpp (revision 15465) @@ -10,11 +10,5 @@ void CreateConstraintsBedSlope(Constraints** pconstraints, IoModel* iomodel){ - /*Recover pointer: */ - Constraints* constraints=*pconstraints; + /*No constraints*/ - /*Create constraints if they do not exist yet*/ - if(!constraints) constraints = new Constraints(); - - /*Assign output pointer: */ - *pconstraints=constraints; } Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/BedSlope/CreateLoadsBedSlope.cpp =================================================================== --- /issm/trunk-jpl/src/c/modules/ModelProcessorx/BedSlope/CreateLoadsBedSlope.cpp (revision 15464) +++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/BedSlope/CreateLoadsBedSlope.cpp (revision 15465) @@ -9,11 +9,5 @@ void CreateLoadsBedSlope(Loads** ploads, IoModel* iomodel){ - /*Recover pointer: */ - Loads* loads=*ploads; + /*No loads*/ - /*Create loads if they do not exist yet*/ - if(!loads) loads = new Loads(); - - /*Assign output pointer: */ - *ploads=loads; } Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/CreateDataSets.cpp =================================================================== --- /issm/trunk-jpl/src/c/modules/ModelProcessorx/CreateDataSets.cpp (revision 15464) +++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/CreateDataSets.cpp (revision 15465) @@ -32,6 +32,8 @@ materials = *pmaterials; - /*Creates Nodes dataset if empty*/ - if(!*pnodes) *pnodes = new Nodes(); + /*Creates Nodes and constraints datasets if empty*/ + if(!*pnodes) *pnodes = new Nodes(); + if(!*pconstraints) *pconstraints = new Constraints(); + if(!*ploads) *ploads = new Loads(); /*Now, branch onto analysis dependent model generation: */ Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/DiagnosticHoriz/CreateConstraintsDiagnosticHoriz.cpp =================================================================== --- /issm/trunk-jpl/src/c/modules/ModelProcessorx/DiagnosticHoriz/CreateConstraintsDiagnosticHoriz.cpp (revision 15464) +++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/DiagnosticHoriz/CreateConstraintsDiagnosticHoriz.cpp (revision 15465) @@ -56,7 +56,4 @@ constraints=*pconstraints; - /*Create constraints if they do not exist yet*/ - if(!constraints) constraints = new Constraints(); - /*Now, is the flag macayaealpattyn on? otherwise, do nothing: */ if(!ismacayealpattyn & !isstokes & !isl1l2){ Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/DiagnosticHoriz/CreateLoadsDiagnosticHoriz.cpp =================================================================== --- /issm/trunk-jpl/src/c/modules/ModelProcessorx/DiagnosticHoriz/CreateLoadsDiagnosticHoriz.cpp (revision 15464) +++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/DiagnosticHoriz/CreateLoadsDiagnosticHoriz.cpp (revision 15465) @@ -40,12 +40,6 @@ Loads* loads=*ploads; - /*Create loads if they do not exist yet*/ - if(!loads) loads = new Loads(); - /*Now, is the flag macayaealpattyn on? otherwise, do nothing: */ - if(!ismacayealpattyn & !isstokes & !isl1l2){ - *ploads=loads; - return; - } + if(!ismacayealpattyn & !isstokes & !isl1l2) return; /*Create pressure loads as boundary conditions. Pay attention to the partitioning if we are running in parallel (the nodes Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/DiagnosticHutter/CreateConstraintsDiagnosticHutter.cpp =================================================================== --- /issm/trunk-jpl/src/c/modules/ModelProcessorx/DiagnosticHutter/CreateConstraintsDiagnosticHutter.cpp (revision 15464) +++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/DiagnosticHutter/CreateConstraintsDiagnosticHutter.cpp (revision 15465) @@ -26,12 +26,6 @@ iomodel->Constant(&ishutter,FlowequationIshutterEnum); - /*Create constraints if they do not exist yet*/ - if(!constraints) constraints = new Constraints(); - /*Now, is the flag ishutter on? otherwise, do nothing: */ - if (!ishutter){ - *pconstraints=constraints; - return; - } + if (!ishutter) return; /*Fetch data: */ Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/DiagnosticHutter/CreateLoadsDiagnosticHutter.cpp =================================================================== --- /issm/trunk-jpl/src/c/modules/ModelProcessorx/DiagnosticHutter/CreateLoadsDiagnosticHutter.cpp (revision 15464) +++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/DiagnosticHutter/CreateLoadsDiagnosticHutter.cpp (revision 15465) @@ -9,11 +9,4 @@ void CreateLoadsDiagnosticHutter(Loads** ploads, IoModel* iomodel){ - /*Recover pointer: */ - Loads* loads=*ploads; - - /*Create loads if they do not exist yet*/ - if(!loads) loads = new Loads(); - - /*Assign output pointer: */ - *ploads=loads; + /*No loads*/ } Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/DiagnosticVert/CreateConstraintsDiagnosticVert.cpp =================================================================== --- /issm/trunk-jpl/src/c/modules/ModelProcessorx/DiagnosticVert/CreateConstraintsDiagnosticVert.cpp (revision 15464) +++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/DiagnosticVert/CreateConstraintsDiagnosticVert.cpp (revision 15465) @@ -20,12 +20,6 @@ Constraints* constraints=*pconstraints; - /*Create constraints if they do not exist yet*/ - if(!constraints) constraints = new Constraints(); - /*return if 2d mesh*/ - if(iomodel->dim==2){ - *pconstraints=constraints; - return; - } + if(iomodel->dim==2) *pconstraints=constraints; /*Fetch data: */ Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/DiagnosticVert/CreateLoadsDiagnosticVert.cpp =================================================================== --- /issm/trunk-jpl/src/c/modules/ModelProcessorx/DiagnosticVert/CreateLoadsDiagnosticVert.cpp (revision 15464) +++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/DiagnosticVert/CreateLoadsDiagnosticVert.cpp (revision 15465) @@ -9,11 +9,5 @@ void CreateLoadsDiagnosticVert(Loads** ploads, IoModel* iomodel){ - /*Recover pointer: */ - Loads* loads=*ploads; + /*No loads*/ - /*Create loads if they do not exist yet*/ - if(!loads) loads = new Loads(); - - /*Assign output pointer: */ - *ploads=loads; } Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/Enthalpy/CreateConstraintsEnthalpy.cpp =================================================================== --- /issm/trunk-jpl/src/c/modules/ModelProcessorx/Enthalpy/CreateConstraintsEnthalpy.cpp (revision 15464) +++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/Enthalpy/CreateConstraintsEnthalpy.cpp (revision 15465) @@ -30,12 +30,6 @@ Constraints* constraints=*pconstraints; - /*Create constraints if they do not exist yet*/ - if(!constraints) constraints = new Constraints(); - /*return if 2d mesh*/ - if(iomodel->dim==2){ - *pconstraints=constraints; - return; - } + if(iomodel->dim==2) return; /*Fetch data: */ Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/Enthalpy/CreateLoadsEnthalpy.cpp =================================================================== --- /issm/trunk-jpl/src/c/modules/ModelProcessorx/Enthalpy/CreateLoadsEnthalpy.cpp (revision 15464) +++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/Enthalpy/CreateLoadsEnthalpy.cpp (revision 15465) @@ -9,11 +9,4 @@ void CreateLoadsEnthalpy(Loads** ploads, IoModel* iomodel){ - /*Recover pointer: */ - Loads* loads=*ploads; - - /*Create loads if they do not exist yet*/ - if(!loads) loads = new Loads(); - - /*Assign output pointer: */ - *ploads=loads; + /*No loads */ } Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/Gia/CreateConstraintsGia.cpp =================================================================== --- /issm/trunk-jpl/src/c/modules/ModelProcessorx/Gia/CreateConstraintsGia.cpp (revision 15464) +++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/Gia/CreateConstraintsGia.cpp (revision 15465) @@ -10,11 +10,5 @@ void CreateConstraintsGia(Constraints** pconstraints, IoModel* iomodel){ - /*Recover pointer: */ - Constraints* constraints=*pconstraints; + /*No constraints*/ - /*Create constraints if they do not exist yet*/ - if(!constraints) constraints = new Constraints(); - - /*Assign output pointer: */ - *pconstraints=constraints; } Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/Gia/CreateLoadsGia.cpp =================================================================== --- /issm/trunk-jpl/src/c/modules/ModelProcessorx/Gia/CreateLoadsGia.cpp (revision 15464) +++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/Gia/CreateLoadsGia.cpp (revision 15465) @@ -9,11 +9,5 @@ void CreateLoadsGia(Loads** ploads, IoModel* iomodel){ - /*Recover pointer: */ - Loads* loads=*ploads; + /*No loads*/ - /*Create loads if they do not exist yet*/ - if(!loads) loads = new Loads(); - - /*Assign output pointer: */ - *ploads=loads; } Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/HydrologyDCEfficient/CreateConstraintsHydrologyDCEfficient.cpp =================================================================== --- /issm/trunk-jpl/src/c/modules/ModelProcessorx/HydrologyDCEfficient/CreateConstraintsHydrologyDCEfficient.cpp (revision 15464) +++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/HydrologyDCEfficient/CreateConstraintsHydrologyDCEfficient.cpp (revision 15465) @@ -12,24 +12,15 @@ /*Recover pointer: */ - bool isefficientlayer; - int hydrology_model; Constraints* constraints=*pconstraints; - /*Create constraints if they do not exist yet*/ - if(!constraints) constraints = new Constraints(); - /*Do we really want DC?*/ + int hydrology_model; iomodel->Constant(&hydrology_model,HydrologyModelEnum); - if(hydrology_model!=HydrologydcEnum){ - *pconstraints=constraints; - return; - } + if(hydrology_model!=HydrologydcEnum) return; /*Do we want an efficient layer*/ + bool isefficientlayer; iomodel->Constant(&isefficientlayer,HydrologydcIsefficientlayerEnum); - if(!isefficientlayer){ - *pconstraints=constraints; - return; - } + if(!isefficientlayer) return; IoModelToConstraintsx(constraints,iomodel,HydrologydcSpceplHeadEnum,HydrologyDCEfficientAnalysisEnum); Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/HydrologyDCEfficient/CreateLoadsHydrologyDCEfficient.cpp =================================================================== --- /issm/trunk-jpl/src/c/modules/ModelProcessorx/HydrologyDCEfficient/CreateLoadsHydrologyDCEfficient.cpp (revision 15464) +++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/HydrologyDCEfficient/CreateLoadsHydrologyDCEfficient.cpp (revision 15465) @@ -9,27 +9,16 @@ void CreateLoadsHydrologyDCEfficient(Loads** ploads, IoModel* iomodel){ - /*Intermediary*/ - bool isefficientlayer; - int hydrology_model; - /*Recover pointer: */ Loads* loads=*ploads; - /*Create loads if they do not exist yet*/ - if(!loads) loads = new Loads(); - /*Do we really want DC?*/ + int hydrology_model; iomodel->Constant(&hydrology_model,HydrologyModelEnum); - if(hydrology_model!=HydrologydcEnum){ - *ploads=loads; - return; - } + if(hydrology_model!=HydrologydcEnum) return; /*Do we want an efficient layer*/ + bool isefficientlayer; iomodel->Constant(&isefficientlayer,HydrologydcIsefficientlayerEnum); - if(!isefficientlayer){ - *ploads=loads; - return; - } + if(!isefficientlayer) return; /*Nothing for now*/ Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/HydrologyDCInefficient/CreateConstraintsHydrologyDCInefficient.cpp =================================================================== --- /issm/trunk-jpl/src/c/modules/ModelProcessorx/HydrologyDCInefficient/CreateConstraintsHydrologyDCInefficient.cpp (revision 15464) +++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/HydrologyDCInefficient/CreateConstraintsHydrologyDCInefficient.cpp (revision 15465) @@ -12,17 +12,10 @@ /*Recover pointer: */ - int hydrology_model; Constraints* constraints=*pconstraints; /*retrieve some parameters: */ + int hydrology_model; iomodel->Constant(&hydrology_model,HydrologyModelEnum); - - /*Create constraints if they do not exist yet*/ - if(!constraints) constraints = new Constraints(); - - if(hydrology_model!=HydrologydcEnum){ - *pconstraints=constraints; - return; - } + if(hydrology_model!=HydrologydcEnum) return; IoModelToConstraintsx(constraints,iomodel,HydrologydcSpcsedimentHeadEnum,HydrologyDCInefficientAnalysisEnum); Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/HydrologyDCInefficient/CreateLoadsHydrologyDCInefficient.cpp =================================================================== --- /issm/trunk-jpl/src/c/modules/ModelProcessorx/HydrologyDCInefficient/CreateLoadsHydrologyDCInefficient.cpp (revision 15464) +++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/HydrologyDCInefficient/CreateLoadsHydrologyDCInefficient.cpp (revision 15465) @@ -9,20 +9,11 @@ void CreateLoadsHydrologyDCInefficient(Loads** ploads, IoModel* iomodel){ - /*Intermediary*/ - int hydrology_model; - /*Recover pointer: */ Loads* loads=*ploads; /*Fetch parameters: */ + int hydrology_model; iomodel->Constant(&hydrology_model,HydrologyModelEnum); - - /*Create loads if they do not exist yet*/ - if(!loads) loads = new Loads(); - - if(hydrology_model!=HydrologydcEnum){ - *ploads=loads; - return; - } + if(hydrology_model!=HydrologydcEnum) return; iomodel->FetchData(1,MeshVertexonbedEnum); Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/HydrologyShreve/CreateConstraintsHydrologyShreve.cpp =================================================================== --- /issm/trunk-jpl/src/c/modules/ModelProcessorx/HydrologyShreve/CreateConstraintsHydrologyShreve.cpp (revision 15464) +++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/HydrologyShreve/CreateConstraintsHydrologyShreve.cpp (revision 15465) @@ -12,18 +12,11 @@ /*Recover pointer: */ - int hydrology_model; - bool isefficientlayer; Constraints* constraints=*pconstraints; /*retrieve some parameters: */ + int hydrology_model; iomodel->Constant(&hydrology_model,HydrologyModelEnum); - /*Create constraints if they do not exist yet*/ - if(!constraints) constraints = new Constraints(); - - if(hydrology_model!=HydrologyshreveEnum){ - *pconstraints=constraints; - return; - } + if(hydrology_model!=HydrologyshreveEnum) return; IoModelToConstraintsx(constraints,iomodel,HydrologyshreveSpcwatercolumnEnum,HydrologyShreveAnalysisEnum); Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/HydrologyShreve/CreateLoadsHydrologyShreve.cpp =================================================================== --- /issm/trunk-jpl/src/c/modules/ModelProcessorx/HydrologyShreve/CreateLoadsHydrologyShreve.cpp (revision 15464) +++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/HydrologyShreve/CreateLoadsHydrologyShreve.cpp (revision 15465) @@ -9,15 +9,5 @@ void CreateLoadsHydrologyShreve(Loads** ploads, IoModel* iomodel){ - /*Intermediary*/ - int numberofvertices; - Pengrid *pengrid = NULL; + /*No loads*/ - /*Recover pointer: */ - Loads* loads=*ploads; - - /*Create loads if they do not exist yet*/ - if(!loads) loads = new Loads(); - - /*Assign output pointer: */ - *ploads=loads; } Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/Melting/CreateConstraintsMelting.cpp =================================================================== --- /issm/trunk-jpl/src/c/modules/ModelProcessorx/Melting/CreateConstraintsMelting.cpp (revision 15464) +++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/Melting/CreateConstraintsMelting.cpp (revision 15465) @@ -10,11 +10,5 @@ void CreateConstraintsMelting(Constraints** pconstraints, IoModel* iomodel){ - /*Recover pointer: */ - Constraints* constraints=*pconstraints; + /*No constraints*/ - /*Create constraints if they do not exist yet*/ - if(!constraints) constraints = new Constraints(); - - /*Assign output pointer: */ - *pconstraints=constraints; } Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/Melting/CreateLoadsMelting.cpp =================================================================== --- /issm/trunk-jpl/src/c/modules/ModelProcessorx/Melting/CreateLoadsMelting.cpp (revision 15464) +++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/Melting/CreateLoadsMelting.cpp (revision 15465) @@ -14,7 +14,4 @@ /*Recover pointer: */ Loads* loads=*ploads; - - /*Create loads if they do not exist yet*/ - if(!loads) loads = new Loads(); //create penalties for nodes: no node can have a temperature over the melting point Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/Prognostic/CreateConstraintsPrognostic.cpp =================================================================== --- /issm/trunk-jpl/src/c/modules/ModelProcessorx/Prognostic/CreateConstraintsPrognostic.cpp (revision 15464) +++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/Prognostic/CreateConstraintsPrognostic.cpp (revision 15465) @@ -14,7 +14,4 @@ Constraints* constraints=*pconstraints; - /*Create constraints if they do not exist yet*/ - if(!constraints) constraints = new Constraints(); - /*Do not add constraints in DG, they are weakly imposed*/ if(stabilization!=3){ Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/Prognostic/CreateLoadsPrognostic.cpp =================================================================== --- /issm/trunk-jpl/src/c/modules/ModelProcessorx/Prognostic/CreateLoadsPrognostic.cpp (revision 15464) +++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/Prognostic/CreateLoadsPrognostic.cpp (revision 15465) @@ -21,7 +21,4 @@ /*Recover pointer: */ Loads* loads=*ploads; - - /*Create loads if they do not exist yet*/ - if(!loads) loads = new Loads(); /*Loads only in DG*/ Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/SurfaceSlope/CreateConstraintsSurfaceSlope.cpp =================================================================== --- /issm/trunk-jpl/src/c/modules/ModelProcessorx/SurfaceSlope/CreateConstraintsSurfaceSlope.cpp (revision 15464) +++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/SurfaceSlope/CreateConstraintsSurfaceSlope.cpp (revision 15465) @@ -10,11 +10,4 @@ void CreateConstraintsSurfaceSlope(Constraints** pconstraints, IoModel* iomodel){ - /*Recover pointer: */ - Constraints* constraints=*pconstraints; - - /*Create constraints if they do not exist yet*/ - if(!constraints) constraints = new Constraints(); - - /*Assign output pointer: */ - *pconstraints=constraints; + /*No constraints*/ } Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/SurfaceSlope/CreateLoadsSurfaceSlope.cpp =================================================================== --- /issm/trunk-jpl/src/c/modules/ModelProcessorx/SurfaceSlope/CreateLoadsSurfaceSlope.cpp (revision 15464) +++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/SurfaceSlope/CreateLoadsSurfaceSlope.cpp (revision 15465) @@ -9,11 +9,5 @@ void CreateLoadsSurfaceSlope(Loads** ploads, IoModel* iomodel){ - /*Recover pointer: */ - Loads* loads=*ploads; + /*No loads*/ - /*Create loads if they do not exist yet*/ - if(!loads) loads = new Loads(); - - /*Assign output pointer: */ - *ploads=loads; } Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/Thermal/CreateConstraintsThermal.cpp =================================================================== --- /issm/trunk-jpl/src/c/modules/ModelProcessorx/Thermal/CreateConstraintsThermal.cpp (revision 15464) +++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/Thermal/CreateConstraintsThermal.cpp (revision 15465) @@ -14,7 +14,4 @@ Constraints* constraints=*pconstraints; - /*Create constraints if they do not exist yet*/ - if(!constraints) constraints = new Constraints(); - /*Only 3d mesh supported*/ if(iomodel->dim==3){ Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/Thermal/CreateLoadsThermal.cpp =================================================================== --- /issm/trunk-jpl/src/c/modules/ModelProcessorx/Thermal/CreateLoadsThermal.cpp (revision 15464) +++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/Thermal/CreateLoadsThermal.cpp (revision 15465) @@ -14,8 +14,5 @@ if(iomodel->dim==2) _error_("2d meshes not supported yet"); - /*Create loads if they do not exist yet*/ - if(!loads) loads = new Loads(); - - //create penalties for nodes: no node can have a temperature over the melting point + /*create penalties for nodes: no node can have a temperature over the melting point*/ iomodel->FetchData(1,ThermalSpctemperatureEnum); CreateSingleNodeToElementConnectivity(iomodel);