Changeset 27850

Timestamp:

07/24/23 07:30:19 (20 months ago)

Author:

rueckamp

Message:

CHG: exclude ComputeMeanEla which is corrupt

Location:

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

Files:

• Element.h (modified) (1 diff)
• Penta.cpp (modified) (1 diff)
• Penta.h (modified) (1 diff)
• Seg.h (modified) (1 diff)
• Tetra.h (modified) (1 diff)
• Tria.h (modified) (1 diff)

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

@@ -257 +257 @@
                 virtual void       ComputeStressTensor(void)=0;
                 virtual void       ComputeEsaStrainAndVorticity(void)=0;
-                virtual void       ComputeMeanEla(IssmDouble* paltitude,int* pcounter)=0;
+                //virtual void       ComputeMeanEla(IssmDouble* paltitude,int* pcounter)=0;
                 virtual void       Configure(Elements* elements,Loads* loads,Nodes* nodes,Vertices* vertices,Materials* materials,Parameters* parameters,Inputs* inputsin)=0;
                 virtual void       ControlInputSetGradient(IssmDouble* gradient,int enum_type,int control_index,int offset,int M,int N,int interp)=0;

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

@@ -5283 +5283 @@
 /*}}}*/
 #endif
-void       Penta::ComputeMeanEla(IssmDouble* paltitude, int* pcounter){/*{{{*/
-
-        if(!this->IsOnSurface() & !this->IsIceInElement()) return;
-        int     ipos=0,ineg=0, counter=0;
-        IssmDouble mean_element_surface=0;
+//void       Penta::ComputeMeanEla(IssmDouble* paltitude, int* pcounter){/*{{{*/
+//
+//      if(!this->IsOnSurface() & !this->IsIceInElement()) return;
+//        int     ipos=0,ineg=0, counter=0;
+//        IssmDouble mean_element_surface=0;
         
         /* input */
-        IssmDouble smb[NUMVERTICES],surface[NUMVERTICES];
-        Element::GetInputListOnVertices(&smb[0],SmbMassBalanceEnum);
-        Element::GetInputListOnVertices(&surface[0],SurfaceEnum);
-
-        for(int iv=0;iv<NUMVERTICES;iv++){
-                if(smb[iv]>0) ipos=ipos+1;
-                if(smb[iv]<0) ineg=ineg+1;
-        }  
+//        IssmDouble smb[NUMVERTICES],surface[NUMVERTICES];
+//      Element::GetInputListOnVertices(&smb[0],SmbMassBalanceEnum);
+//      Element::GetInputListOnVertices(&surface[0],SurfaceEnum);
+
+//      for(int iv=0;iv<NUMVERTICES;iv++){
+//              if(smb[iv]>0) ipos=ipos+1;
+//              if(smb[iv]<0) ineg=ineg+1;
+//      }  
 
         /* we define ELA if an element has pos and neg smb */
-        if(ipos>0 & ineg>0){
-                for(int iv=0;iv<NUMVERTICES;iv++) mean_element_surface=mean_element_surface+surface[iv]/double(NUMVERTICES);
-                *paltitude+=mean_element_surface;
-                *pcounter+=counter+1;
-        }
-}
-/*}}}*/
+//      if(ipos>0 & ineg>0){
+//              for(int iv=0;iv<NUMVERTICES;iv++) mean_element_surface=mean_element_surface+surface[iv]/double(NUMVERTICES);
+//              *paltitude+=mean_element_surface;
+//              *pcounter+=counter+1;
+//      }
+//}
+/*}}}*/

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

@@ -64 +64 @@
                 void           ComputeSigmaNN(){_error_("not implemented yet");};
                 void           ComputeStressTensor();
+                //void           ComputeMeanEla(IssmDouble* paltitude, int* pcounter);
                 void           Configure(Elements* elements,Loads* loads,Nodes* nodes,Vertices* vertices,Materials* materials,Parameters* parameters,Inputs* inputsin);
                 void           ControlInputSetGradient(IssmDouble* gradient,int enum_type,int control_index,int offset,int M,int N,int interp);

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

@@ -50 +50 @@
                 void        ComputeSigmaNN(){_error_("not implemented yet");};
                 void        ComputeStressTensor(){_error_("not implemented yet");};
-                void        ComputeMeanEla(IssmDouble* paltitude, int* pcounter){_error_("not implemented yet");};
+                //void        ComputeMeanEla(IssmDouble* paltitude, int* pcounter){_error_("not implemented yet");};
                 void        Configure(Elements* elements,Loads* loads,Nodes* nodesin,Vertices* verticesin,Materials* materials,Parameters* parameters,Inputs* inputsin){_error_("not implemented yet");};
                 void        ControlInputSetGradient(IssmDouble* gradient,int enum_type,int control_index,int offset,int M,int N,int interp){_error_("not implemented yet");};

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

@@ -46 +46 @@
                 void        ComputeSigmaNN(){_error_("not implemented yet");};
                 void        ComputeStressTensor(){_error_("not implemented yet");};
-                void        ComputeMeanEla(IssmDouble* paltitude, int* pcounter){_error_("not implemented yet");};
+                //void        ComputeMeanEla(IssmDouble* paltitude, int* pcounter){_error_("not implemented yet");};
                 void        ComputeDeviatoricStressTensor(){_error_("not implemented yet");};
                 void        ComputeEsaStrainAndVorticity(){_error_("not implemented yet!");};

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

@@ -72 +72 @@
                 void        ComputeStressTensor();
                 void        ComputeSurfaceNormalVelocity();
-                void        ComputeMeanEla(IssmDouble* paltitude, int* pcounter){_error_("not implemented yet");};
+                //void        ComputeMeanEla(IssmDouble* paltitude, int* pcounter){_error_("not implemented yet");};
                 void        Configure(Elements* elements,Loads* loads,Nodes* nodesin,Vertices* verticesin,Materials* materials,Parameters* parameters,Inputs* inputsin);
                 void        ControlInputSetGradient(IssmDouble* gradient,int enum_type,int control_index,int offset,int M,int N,int interp);