Changeset 4516 for issm/trunk/src/m/solutions/steadystate_core.m

Timestamp:

07/12/10 07:55:34 (15 years ago)

Author:

seroussi

Message:

steadystate running in serial

File:

• issm/trunk/src/m/solutions/steadystate_core.m (modified) (3 diffs)

issm/trunk/src/m/solutions/steadystate_core.m

@@ -25 +25 @@
                 femmodel=diagnostic_core(femmodel); 
 
-                displaystring(verbose,'\n%s',['checking temperature, velocity and pressure convergence...']);
                 if step>1,
-                        if steadystateconvergence(femmodel),
+                        displaystring(verbose,'\n%s',['checking temperature, velocity and pressure convergence...']);
+                        converged=steadystateconvergence(femmodel);
+                        if converged,
                                 break;
                         end
@@ -36 +37 @@
                 femmodel.elements=InputDuplicate(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VyEnum,VyOldEnum);
                 femmodel.elements=InputDuplicate(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VzEnum,VzOldEnum);
-                femmodel.elements=InputDuplicate(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,PressureEnum,PressureOlddnum);
+                femmodel.elements=InputDuplicate(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,PressureEnum,PressureOldEnum);
                 femmodel.elements=InputDuplicate(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,TemperatureEnum,TemperatureOldEnum);
 
                 %Increase counter
                 step=step+1;
-
         end
 
@@ -49 +49 @@
                 femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VyEnum);
                 femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VzEnum);
+                femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VelEnum);
                 femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,PressureEnum);
                 femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,TemperatureEnum);
+                femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,MeltingRateEnum);
         end