Changeset 19725 for issm/trunk-jpl/src/c/classes/Inputs/TriaInput.cpp

Timestamp:

11/16/15 16:29:17 (9 years ago)

Author:

Mathieu Morlighem

Message:

CHG: working on the RHS of Hydro Aleah's model

File:

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

issm/trunk-jpl/src/c/classes/Inputs/TriaInput.cpp

@@ -192 +192 @@
 }
 /*}}}*/
+void TriaInput::GetInputDerivativeAverageValue(IssmDouble* derivativevalues, IssmDouble* xyz_list){/*{{{*/
+
+        int        numnodes  = this->NumberofNodes(this->interpolation_type);
+        IssmDouble numnodesd = reCast<int,IssmDouble>(numnodes);
+        IssmDouble dvalue[3];
+
+        derivativevalues[0] = 0.;
+        derivativevalues[1] = 0.;
+
+        GaussTria* gauss=new GaussTria();
+        for(int iv=0;iv<numnodes;iv++){
+                gauss->GaussNode(this->interpolation_type,iv);
+                this->GetInputDerivativeValue(&dvalue[0],xyz_list,gauss);
+
+                derivativevalues[0] += dvalue[0]/numnodesd;
+                derivativevalues[1] += dvalue[1]/numnodesd;
+        }
+        delete gauss;
+}
+/*}}}*/
 void TriaInput::GetInputAllTimeAverages(IssmDouble** pvalues,IssmDouble** ptimes, int* pnumtimes){/*{{{*/