Changeset 27861

Timestamp:

07/27/23 10:32:21 (2 years ago)

Author:

vverjans

Message:

CHG: adding case for levelset==0 at a vertex in GetIcefrontArea()

File:

• issm/trunk-jpl/src/c/classes/Elements/Tria.cpp (modified) (1 diff)

issm/trunk-jpl/src/c/classes/Elements/Tria.cpp

@@ -2854 +2854 @@
                                 H[3]=Haux[1]+s[1]*(Haux[2]-Haux[1]);
                                 Haverage=(H[2]+H[3])/2;
+                                break;
+                        case 3:
+                                if(counter==1) distance = 0; //front has 0 width on this element because levelset is 0 at a single vertex
+                                else if(counter==2){ //two vertices with levelset=0: averaging ice front depth over both
+                                        Haverage = 0;
+                                        for(int i=0;i<NUMVERTICES;i++){
+                                                if(lsf[indices[i]]==0.) Haverage -= Haux[indices[i]]/2;
+                                                if(Haverage<Haux[indices[i]]/2-1e-3) break; //done with the two vertices
+                                        }
+                                }
                                 break;
                         default: