Changeset 3975 for issm/trunk

Timestamp:

05/27/10 15:23:58 (15 years ago)

Author:

seroussi

Message:

modifications in balanced thickness (serial)

File:

• issm/trunk/src/m/solutions/jpl/balancedthickness_core.m (modified) (1 diff)

issm/trunk/src/m/solutions/jpl/balancedthickness_core.m

@@ -6 +6 @@
 
         %get FE model
-        m=models.bt;
+        verbose=models.bt.parameters.Verbose;
         results.time=0;
         results.step=1;
 
-        displaystring(m.parameters.Verbose,'\n%s',['call computational core:']);
-        results.h_g=diagnostic_core_linear(m,analysis_type,sub_analysis_type);
+        displaystring(verbose,'\n%s',['depth averaging velocity...']);
+        [models.bt.elements,models.bt.nodes,models.bt.vertices,models.bt.loads,models.bt.materials,models.bt.parameters]=DepthAverageInput(models.bt.elements,models.bt.nodes,models.bt.vertices,models.bt.loads,models.bt.materials,models.bt.parameters,VxEnum);
+        [models.bt.elements,models.bt.nodes,models.bt.vertices,models.bt.loads,models.bt.materials,models.bt.parameters]=DepthAverageInput(models.bt.elements,models.bt.nodes,models.bt.vertices,models.bt.loads,models.bt.materials,models.bt.parameters,VyEnum);
 
-        displaystring(m.parameters.Verbose,'\n%s',['extrude computed thickness on all layers:']);
-        results.h_g=FieldExtrude(m.elements,m.nodes,m.vertices,m.loads,m.materials,m.parameters,results.h_g,'thickness',0);
+        displaystring(verbose,'\n%s',['call computational core:']);
+        h_g=diagnostic_core_linear(models.bt,analysis_type,sub_analysis_type);
 
+        %Update
+        models.bt.elements=UpdateInputsFromSolution(models.bt.elements,models.bt.nodes,models.bt.vertices,models.bt.loads,models.bt.materials,models.bt.parameters,h_g,BalancedthicknessAnalysisEnum,NoneAnalysisEnum);
+
+        displaystring(verbose,'\n%s',['extrude computed thickness on all layers:']);
+        models.bt.elements=ExtrudeInput(models.bt.elements,models.bt.nodes,models.bt.vertices,models.bt.loads,models.bt.materials,models.bt.parameters,ThicknessEnum);
+
+        %NEED TO BE CLEANED
+        results.h_g=h_g;
 end %end function