Changeset 5195

Timestamp:

08/12/10 12:16:59 (15 years ago)

Author:

seroussi

Message:

output_frequency matlab

Location:

issm/trunk/src/m

Files:

• classes/@model/model.m (modified) (1 diff)
• classes/@model/setdefaultparameters.m (modified) (1 diff)
• classes/public/marshall.m (modified) (1 diff)
• solutions/transient2d_core.m (modified) (2 diffs)
• solutions/transient3d_core.m (modified) (2 diffs)
• solvers/solver_linear.m (modified) (1 diff)

issm/trunk/src/m/classes/@model/model.m

@@ -233 +233 @@
 
         %Results fields
+        md.output_frequency=0;
         md.inputfilename='';
         md.outputfilename='';

issm/trunk/src/m/classes/@model/setdefaultparameters.m

@@ -231 +231 @@
 md.plot=0;
 
+%How often to save results, default is 1 so save every step
+md.output_frequency=1;
+
 %Parallelisation parameters
 

issm/trunk/src/m/classes/public/marshall.m

@@ -110 +110 @@
 WriteData(fid,md.verbose,'Integer','verbose');
 WriteData(fid,md.plot,'Integer','plot');
+WriteData(fid,md.output_frequency,'Integer','output_frequency');
 WriteData(fid,md.artificial_diffusivity,'Integer','artificial_diffusivity');
 WriteData(fid,md.prognostic_DG,'Integer','prognostic_DG');

issm/trunk/src/m/solutions/transient2d_core.m

@@ -11 +11 @@
         yts=femmodel.parameters.Yts;
         solution_type=femmodel.parameters.SolutionType;
+        output_frequency=femmodel.parameters.OutputFrequency;
 
         %Initialize
@@ -34 +35 @@
                 [femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters]=UpdateGeometry(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters);
 
-                if solution_type==Transient2DSolutionEnum,
+                if (solution_type==Transient2DSolutionEnum & mod(step,output_frequency)==0),
                         displaystring(verbose,'\n%s',['saving results...']);
                         femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VxEnum,step,time);

issm/trunk/src/m/solutions/transient3d_core.m

@@ -11 +11 @@
         yts=femmodel.parameters.Yts;
         solution_type=femmodel.parameters.SolutionType;
+        output_frequency=femmodel.parameters.OutputFrequency;
 
         %Initialize
@@ -40 +41 @@
                 [femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters]=UpdateVertexPositions(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters);
 
-                if solution_type==Transient3DSolutionEnum,
+                if (solution_type==Transient3DSolutionEnum & mod(step,output_frequency)==0),
                         displaystring(verbose,'\n%s',['saving results...']);
                         femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VxEnum,step,time);

issm/trunk/src/m/solvers/solver_linear.m

@@ -27 +27 @@
         %Update inputs using new solution
         [femmodel.elements,femmodel.materials]=InputUpdateFromSolution(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,u_g);
-