Changeset 9255

Timestamp:

08/10/11 15:19:39 (14 years ago)

Author:

bkhakbaz

Message:

added a function called GetHydrologyWaterVelocity

File:

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

issm/trunk/src/c/objects/Elements/Tria.cpp

@@ -2659 +2659 @@
 }
 /*}}}*/
+/*FUNCTION Tria::GetHydrologyWaterVelocity {{{1*/
+void Tria::GetHydrologyWaterVelocity(double* v,GaussTria* gauss){
+        /*material parameters: */
+        double mu_water=0.001787;  //unit= [N s/m2]
+        double w;
+        double rho_ice, rho_water, g;
+        double dsdx,dsdy,dbdx,dbdy;
+
+        /*Retrieve all inputs and parameters*/
+        rho_ice=matpar->GetRhoIce();
+        rho_water=matpar->GetRhoWater();
+        g=matpar->GetG();
+        Input* surfaceslopex_input=inputs->GetInput(SurfaceSlopeXEnum); _assert_(surfaceslopex_input);
+        Input* surfaceslopey_input=inputs->GetInput(SurfaceSlopeYEnum); _assert_(surfaceslopey_input);
+        Input* bedslopex_input=inputs->GetInput(BedSlopeXEnum);         _assert_(bedslopex_input);
+        Input* bedslopey_input=inputs->GetInput(BedSlopeYEnum);         _assert_(bedslopey_input);
+        Input* watercolumn_input=inputs->GetInput(WatercolumnEnum);     _assert_(watercolumn_input);
+
+        /*recover slopes: */
+        surfaceslopex_input->GetParameterValue(&dsdx,gauss);
+        surfaceslopey_input->GetParameterValue(&dsdy,gauss);
+        bedslopex_input->GetParameterValue(&dbdx,gauss);
+        bedslopey_input->GetParameterValue(&dbdy,gauss);
+
+        /*recover water column: */
+        watercolumn_input->GetParameterValue(&w,gauss);
+
+        /* Water velocity x and y components */
+        v[0]=pow(w,2)/(12 * mu_water)*(rho_ice*g*dsdx+(rho_water-rho_ice)*g*dbdx);
+        v[1]=pow(w,2)/(12 * mu_water)*(rho_ice*g*dsdy+(rho_water-rho_ice)*g*dbdy);
+}
+/*}}}*/
 /*FUNCTION Tria::Gradj {{{1*/
 void  Tria::Gradj(Vec gradient,int control_type){