Changeset 21898

Timestamp:

07/30/17 20:13:36 (8 years ago)

Author:

Mathieu Morlighem

Message:

CHG: added Max function for Seg

Location:

issm/trunk-jpl/src/c/classes/Elements

Files:

• Seg.cpp (modified) (1 diff)
• Seg.h (modified) (1 diff)

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

@@ -168 +168 @@
         *pxyz_front=xyz_front;
 }/*}}}*/
+IssmDouble Seg::GetGroundedPortion(IssmDouble* xyz_list){/*{{{*/
+        /*Computeportion of the element that is grounded*/ 
+
+        bool              mainlyfloating = true;
+        const IssmPDouble epsilon        = 1.e-15;
+        IssmDouble        phi;
+        IssmDouble        gl[NUMVERTICES];
+
+        /*Recover parameters and values*/
+        parameters->FindParam(&domaintype,DomainTypeEnum);
+        GetInputListOnVertices(&gl[0],MaskGroundediceLevelsetEnum);
+
+        /*Be sure that values are not zero*/
+        if(gl[0]==0.) gl[0]=gl[0]+epsilon;
+        if(gl[1]==0.) gl[1]=gl[1]+epsilon;
+
+        if(gl[0]>0 && gl[1]>0) phi=1; // All grounded
+        else if(gl[0]<0 && gl[1]<0) phi=0; // All floating
+        else if(gl[0]<0 && gl[1]>0){ //1 grounded
+                phi=1./(1.-gl[0]/gl[1]);
+        }
+        else if(gl[1]<0 && gl[0]>0){ //0 grounded
+                phi=1./(1.-gl[1]/gl[0]);
+        }
+
+        if(phi>1 || phi<0) _error_("Error. Problem with portion of grounded element: value should be between 0 and 1");
+
+        return phi;
+}
+/*}}}*/
 int        Seg::GetNumberOfNodes(void){/*{{{*/
         return this->NumberofNodes(this->element_type);

issm/trunk-jpl/src/c/classes/Elements/Seg.h

@@ -62 +62 @@
                 int         GetElementType(void){_error_("not implemented yet");};
                 void        GetGroundedPart(int* point1,IssmDouble* fraction1, IssmDouble* fraction2,bool* mainlyfloating){_error_("not implemented yet");};
-                IssmDouble  GetGroundedPortion(IssmDouble* xyz_list){_error_("not implemented yet");};
+                IssmDouble  GetGroundedPortion(IssmDouble* xyz_list);
                 void               GetIcefrontCoordinates(IssmDouble** pxyz_front,IssmDouble* xyz_list,int levelsetenum);
                 void        GetInputValue(IssmDouble* pvalue,Node* node,int enumtype){_error_("not implemented yet");};