Changeset 24718

Timestamp:

04/21/20 21:14:41 (5 years ago)

Author:

ibsmith

Message:

Updated CO2 parameters based on new study.

File:

• issm/trunk-jpl/src/c/shared/Elements/NyeCO2.cpp (modified) (1 diff)

issm/trunk-jpl/src/c/shared/Elements/NyeCO2.cpp

@@ -12 +12 @@
 
         /*Coefficients*/
-        const IssmPDouble Rg      = 8.3144598;     /* J mol^-1 K^-1 */ 
-        const IssmPDouble A_const = pow(10.,10.8); /*s^-1 MPa       */ 
-        const IssmPDouble Q       = 63000.;        /*J mol^-1       */ 
-        const IssmPDouble n       = 7.;            /*Glen's exponent*/
+        const IssmPDouble Rg      = 8.3144598;     /* J mol^-1 K^-1   */ 
+        const IssmPDouble A_const = pow(10.,13.0); /* s^-1 MPa        */ 
+        const IssmPDouble Q       = 66900.;        /* J mol^-1        */ 
+        const IssmPDouble n       = 8.;            /* Glen's exponent */
 
         /*Arrhenius Law*/
-        IssmDouble A = A_const *exp(-Q/(temperature*Rg));  /*s^-1 MPa   */
-        IssmDouble B = 1e6*pow(A,-1/n);                    /*s^(1/n) Pa */
+        IssmDouble A = A_const *exp(-Q/(temperature*Rg));  /* s^-1 MPa   */
+        IssmDouble B = 1e6*pow(A,-1/n);                    /* s^(1/n) Pa */
 
         /*Beyond-melting-point cases*/