Changeset 27302

Timestamp:

10/10/22 05:12:40 (2 years ago)

Author:

Mathieu Morlighem

Message:

CHG: minor AD

File:

• issm/trunk-jpl/src/c/analyses/DebrisAnalysis.cpp (modified) (2 diffs)

issm/trunk-jpl/src/c/analyses/DebrisAnalysis.cpp

@@ -604 +604 @@
 
         int numnodes = element->GetNumberOfNodes();
-        IssmDouble* thickness    = xNew<IssmDouble>(numnodes);
-        IssmDouble* thicknessold    = xNew<IssmDouble>(numnodes);
-        IssmDouble* newthickness = xNew<IssmDouble>(numnodes);
-        IssmDouble* icethickness    = xNew<IssmDouble>(numnodes);
-        IssmDouble* bedslopex    = xNew<IssmDouble>(numnodes);
-        IssmDouble* surfaceslopex    = xNew<IssmDouble>(numnodes);
+        IssmDouble* thickness     = xNew<IssmDouble>(numnodes);
+        IssmDouble* thicknessold  = xNew<IssmDouble>(numnodes);
+        IssmDouble* newthickness  = xNew<IssmDouble>(numnodes);
+        IssmDouble* icethickness  = xNew<IssmDouble>(numnodes);
+        IssmDouble* bedslopex     = xNew<IssmDouble>(numnodes);
+        IssmDouble* surfaceslopex = xNew<IssmDouble>(numnodes);
 
         /*Use the dof list to index into the solution vector: */
@@ -701 +701 @@
                                         if(isminthicknessinelement){
                                                 for(k=0; k<numnodes;k++){
-                                                        if(onsurface[k]){
+                                                        if(onsurface[k]>0.5){
                                                                 slope=fabs(slopex[k]);
                                                                 if(dim==2) slope=pow(pow(slopex[k],2)+pow(slopey[k],2),0.5);