Changeset 17267

Timestamp:

02/12/14 15:34:57 (11 years ago)

Author:

jbondzio

Message:

ADD: Penta::IsIcefront()

Location:

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

Files:

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

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

@@ -3175 +3175 @@
 }
 /*}}}*/
-
+bool  Penta::IsIcefront(void){/*{{{*/
+
+        bool isicefront;
+        int i,nrice;
+   IssmDouble ls[NUMVERTICES];
+
+        /*Retrieve all inputs and parameters*/
+        GetInputListOnVertices(&ls[0],MaskIceLevelsetEnum);
+
+        /* If only one vertex has ice, there is an ice front here */
+        isicefront=false;
+        if(IsIceInElement()){
+                nrice=0;       
+                for(i=0;i<NUMVERTICES2D;i++)
+                        if(ls[i]<0.) nrice++;
+                if(nrice==1) isicefront= true;
+        }
+        return isicefront;
+}/*}}}*/
 /*FUNCTION Penta::AverageOntoPartition {{{*/
 void  Penta::AverageOntoPartition(Vector<IssmDouble>* partition_contributions,Vector<IssmDouble>* partition_areas,IssmDouble* vertex_response,IssmDouble* qmu_part){

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

@@ -114 +114 @@
                 int    PressureInterpolation();
                 bool   IsZeroLevelset(int levelset_enum);
-                bool   IsIcefront(void){_error_("not implemented yet");};
+                bool   IsIcefront(void);
                 void   ZeroLevelsetCoordinates(IssmDouble** pxyz_zero,IssmDouble* xyz_list,int levelsetenum);
                 void   GetIcefrontCoordinates(IssmDouble** pxyz_front,IssmDouble* xyz_list,int levelsetenum){_error_("not implemented yet");};