Changeset 3180 for issm/trunk/src/c/objects/Penta.cpp

Timestamp:

03/04/10 09:31:08 (15 years ago)

Author:

Mathieu Morlighem

Message:

Added Misfit 3

File:

• issm/trunk/src/c/objects/Penta.cpp (modified) (1 diff)

issm/trunk/src/c/objects/Penta.cpp

@@ -4209 +4209 @@
 }
 /*}}}*/
+/*FUNCTION SurfaceArea {{{1*/
+#undef __FUNCT__ 
+#define __FUNCT__ "Penta::SurfaceArea"
+double Penta::SurfaceArea(void* inputs,int analysis_type,int sub_analysis_type){
+
+        double S;
+        Tria* tria=NULL;
+
+        /*If on water, return 0: */
+        if(onwater)return 0;
+
+        /*Bail out if this element if:
+         * -> Non collapsed and not on the surface
+         * -> collapsed (2d model) and not on bed) */
+        if ((!collapse && !onsurface) || (collapse && !onbed)){
+                return 0;
+        }
+        else if (collapse){
+
+                /*This element should be collapsed into a tria element at its base. Create this tria element, 
+                 * and compute SurfaceArea*/
+                tria=(Tria*)SpawnTria(0,1,2); //grids 0, 1 and 2 make the new tria (lower face).
+                S=tria->SurfaceArea(inputs,analysis_type,sub_analysis_type);
+                delete tria;
+                return S;
+        }
+        else{
+
+                tria=(Tria*)SpawnTria(3,4,5); //grids 3, 4 and 5 make the new tria (upper face).
+                S=tria->SurfaceArea(inputs,analysis_type,sub_analysis_type);
+                delete tria;
+                return S;
+        }
+}
+/*}}}*/
 /*FUNCTION SurfaceNormal {{{1*/
 #undef __FUNCT__