Changeset 26716 for issm/trunk-jpl/src/jl/solve/analyses/stressbalanceanalysis.jl

Timestamp:

12/08/21 06:09:27 (3 years ago)

Author:

Mathieu Morlighem

Message:

CHG: as usual

File:

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

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

@@ -177 +177 @@
 
         if(IsIcefront(element))
-                A = zeros(6)
 
                 #Get additional parameters and inputs
@@ -187 +186 @@
                 GetIcefrontCoordinates!(element, xyz_list_front, xyz_list, MaskIceLevelsetEnum)
                 nx, ny = NormalSection(element, xyz_list_front)
-                #println("nx ",nx," ny:",ny)
 
                 gauss = GaussTria(element, xyz_list, xyz_list_front, 3)
@@ -202 +200 @@
 
                         for i in 1:numnodes
-                                #println("1: ",gauss.weights[ig]*Jdet*term*nx*basis[i]," 2: ",gauss.weights[ig]*Jdet*term*ny*basis[i])
                                 pe.values[2*i-1] += gauss.weights[ig]*Jdet*term*nx*basis[i]
                                 pe.values[2*i  ] += gauss.weights[ig]*Jdet*term*ny*basis[i]
-                                A[2*i-1] += gauss.weights[ig]*Jdet*term*nx*basis[i]
-                                A[2*i  ] += gauss.weights[ig]*Jdet*term*ny*basis[i]
                         end
                 end