Changeset 20616

Timestamp:

05/16/16 15:59:49 (9 years ago)

Author:

Mathieu Morlighem

Message:

NEW: enabling inversions for frictiontemp

Location:

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

Files:

• Friction.cpp (modified) (3 diffs)
• Friction.h (modified) (1 diff)

issm/trunk-jpl/src/c/classes/Loads/Friction.cpp

@@ -49 +49 @@
                         GetAlphaHydroComplement(palpha_complement,gauss);
                         break;
+                case 4:
+                        GetAlphaTempComplement(palpha_complement,gauss);
+                        break;
           default:
                         _error_("not supported");
@@ -54 +57 @@
 
 }/*}}}*/
-
 void Friction::GetAlphaViscousComplement(IssmDouble* palpha_complement, Gauss* gauss){/*{{{*/
 
@@ -191 +193 @@
 }
 /*}}}*/
+void Friction::GetAlphaTempComplement(IssmDouble* palpha_complement, Gauss* gauss){/*{{{*/
+        /*Here, we want to parameterize the friction as a function of temperature
+         *
+         * alpha2 = alpha2_viscous * 1/f(T)
+         *
+         * where f(T) = exp((T-Tpmp)/gamma)
+         */
+
+        /*Intermediaries: */
+        IssmDouble  f,T,pressure,Tpmp,gamma;
+        IssmDouble  alpha_complement;
+
+        /*Get viscous part*/
+        this->GetAlphaViscousComplement(&alpha_complement,gauss);
+
+        /*Get pressure melting point (Tpmp) for local pressure and get current temperature*/
+        element->GetInputValue(&T,gauss,TemperatureEnum);
+        element->GetInputValue(&pressure,gauss,PressureEnum);
+        Tpmp = element->TMeltingPoint(pressure);
+
+        /*Compute scaling parameter*/
+        element->parameters->FindParam(&gamma,FrictionGammaEnum);
+        alpha_complement = alpha_complement/ exp((T-Tpmp)/gamma);
+
+        /*Assign output pointers:*/
+        *palpha_complement=alpha_complement;
+}/*}}}*/
 void Friction::GetAlpha2(IssmDouble* palpha2, Gauss* gauss){/*{{{*/
 

issm/trunk-jpl/src/c/classes/Loads/Friction.h

@@ -32 +32 @@
                 void  GetAlphaViscousComplement(IssmDouble* alpha_complement,Gauss* gauss);
                 void  GetAlphaHydroComplement(IssmDouble* alpha_complement,Gauss* gauss);
+                void  GetAlphaTempComplement(IssmDouble* alpha_complement,Gauss* gauss);
                 void  GetAlpha2(IssmDouble* palpha2,Gauss* gauss);
                 void  GetAlpha2Coulomb(IssmDouble* palpha2,Gauss* gauss);