Ignore:
Timestamp:
07/12/10 07:55:34 (15 years ago)
Author:
seroussi
Message:

steadystate running in serial

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/src/m/solutions/steadystate_core.m

    r4505 r4516  
    2525                femmodel=diagnostic_core(femmodel);
    2626
    27                 displaystring(verbose,'\n%s',['checking temperature, velocity and pressure convergence...']);
    2827                if step>1,
    29                         if steadystateconvergence(femmodel),
     28                        displaystring(verbose,'\n%s',['checking temperature, velocity and pressure convergence...']);
     29                        converged=steadystateconvergence(femmodel);
     30                        if converged,
    3031                                break;
    3132                        end
     
    3637                femmodel.elements=InputDuplicate(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VyEnum,VyOldEnum);
    3738                femmodel.elements=InputDuplicate(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VzEnum,VzOldEnum);
    38                 femmodel.elements=InputDuplicate(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,PressureEnum,PressureOlddnum);
     39                femmodel.elements=InputDuplicate(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,PressureEnum,PressureOldEnum);
    3940                femmodel.elements=InputDuplicate(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,TemperatureEnum,TemperatureOldEnum);
    4041
    4142                %Increase counter
    4243                step=step+1;
    43 
    4444        end
    4545
     
    4949                femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VyEnum);
    5050                femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VzEnum);
     51                femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VelEnum);
    5152                femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,PressureEnum);
    5253                femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,TemperatureEnum);
     54                femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,MeltingRateEnum);
    5355        end
    5456
Note: See TracChangeset for help on using the changeset viewer.