Changeset 26730 for issm/trunk-jpl/src/jl/core/analyses/stressbalanceanalysis.jl

Timestamp:

12/12/21 14:16:13 (3 years ago)

Author:

Mathieu Morlighem

Message:

CHG: adding mass transport, not quite working yet

File:

• issm/trunk-jpl/src/jl/core/analyses/stressbalanceanalysis.jl (modified) (3 diffs)

issm/trunk-jpl/src/jl/core/analyses/stressbalanceanalysis.jl

@@ -33 +33 @@
 
         #Provide node indices to element
-        for i in 1:md.mesh.numberofvertices
+        for i in 1:md.mesh.numberofelements
                 Update(elements[i],inputs,i,md,P1Enum)
         end
@@ -72 +72 @@
 function Core(analysis::StressbalanceAnalysis,femmodel::FemModel)# {{{
 
+        #Set current analysis to Stressnalance
+        SetCurrentConfiguration!(femmodel, analysis)
+
         #Fetch parameters relevant to solution sequence
         maxiter = FindParam(femmodel.parameters,StressbalanceMaxiterEnum)
@@ -78 +81 @@
         abstol  = FindParam(femmodel.parameters,StressbalanceAbstolEnum)
 
+        #Call solution sequence to compute new speeds
         solutionsequence_nonlinear(femmodel,analysis,maxiter,restol,reltol,abstol)
+
+        #Save output
+        RequestedOutputsx(femmodel, [VxEnum,VyEnum,VelEnum])
 
 end #}}}