Changeset 27304

Timestamp:

10/10/22 07:54:47 (2 years ago)

Author:

Mathieu Morlighem

Message:

BUG: fixing AD

File:

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

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

@@ -170 +170 @@
                 }else{
                         for(int i=0;i<numvertices;i++){
-                                //if(onsurface[i]){
+                                //if(onsurface[i]>.5){
                                 slope=fabs(slopex[i]);
                                 if((atan(slope)*rad2deg)>25.){
@@ -736 +736 @@
                                                 int kk=0;
                                                 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);
@@ -837 +837 @@
                                                 test=0;
                                                 for(k=0; k<numnodes;k++){
-                                                        if(onsurface[k]){
+                                                        if(onsurface[k]>.5){
                                                                 slope=pow(slopex[k]*slopex[k],0.5);
                                                                 if(dim==2) slope=pow(pow(slopex[k],2)+pow(slopey[k],2),0.5);
@@ -859 +859 @@
                                                 IssmDouble min=1e14;
                                                 for(k=0; k<numnodes;k++){
-                                                        if(onsurface[k]){
+                                                        if(onsurface[k]>.5){
                                                                 if(surface[k]<min){
                                                                         index=k;
@@ -867 +867 @@
                                                 }
                                                 for(k=0; k<numnodes;k++){
-                                                        if(onsurface[k]){
+                                                        if(onsurface[k]>.5){
                                                                 if(k==index){
                                                                         debristhickness[k]=debristhickness[k]+debrissum;