source: issm/oecreview/Archive/23390-24306/ISSM-24103-24104.diff

Last change on this file was 24307, checked in by Mathieu Morlighem, 5 years ago

NEW: adding Archive/23390-24306
File size: 1.9 KB
RevLine 
[24307]1Index: ../trunk-jpl/src/c/shared/Elements/NyeCO2.cpp
2===================================================================
3--- ../trunk-jpl/src/c/shared/Elements/NyeCO2.cpp (revision 24103)
4+++ ../trunk-jpl/src/c/shared/Elements/NyeCO2.cpp (revision 24104)
5@@ -4,15 +4,17 @@
6 * where rigidigty (in s^(1/n)Pa) is the flow law paramter in the flow law sigma=B*e(1/n) (Nye, p2000).
7 */
8
9-#include "../io/io.h" #include <math.h> #include "../Numerics/types.h"
10+#include "../io/io.h"
11+#include <math.h>
12+#include "../Numerics/types.h"
13
14 IssmDouble NyeCO2(IssmDouble temperature){
15
16 /*Coefficients*/
17- const IssmPDouble Rg = 8.3144598; /* J mol^-1 K^-1 */ const
18- IssmPDouble A_const = pow(10.,10.8); /*s^-1 MPa */ const IssmPDouble
19- Q = 63000.; /*J mol^-1 */ const IssmPDouble n =
20- 7.; /*Glen's exponent*/
21+ const IssmPDouble Rg = 8.3144598; /* J mol^-1 K^-1 */
22+ const IssmPDouble A_const = pow(10.,10.8); /*s^-1 MPa */
23+ const IssmPDouble Q = 63000.; /*J mol^-1 */
24+ const IssmPDouble n = 7.; /*Glen's exponent*/
25
26 /*Arrhenius Law*/
27 IssmDouble A = A_const *exp(-Q/(temperature*Rg)); /*s^-1 MPa */
28@@ -19,10 +21,8 @@
29 IssmDouble B = 1e6*pow(A,-1/n); /*s^(1/n) Pa */
30
31 /*Beyond-melting-point cases*/
32- if((temperature>200.)&&(temperature<220.)) printf("CO2 ICE - POSSIBLE
33- MELTING. Some temperature values are between 200K and 220K.\n");
34- else if(temperature>=220.) _printf0_("CO2 ICE - GUARANTEED MELTING. Some
35- temperature values are beyond 220K.\n");
36+ if((temperature>200.)&&(temperature<220.)) _printf0_("CO2 ICE - POSSIBLE MELTING. Some temperature values are between 200K and 220K.\n");
37+ else if(temperature>=220.) _printf0_("CO2 ICE - GUARANTEED MELTING. Some temperature values are beyond 220K.\n");
38
39 /*Return output*/
40 return B;
Note: See TracBrowser for help on using the repository browser.