Changeset 24104
- Timestamp:
- 07/22/19 16:51:26 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/shared/Elements/NyeCO2.cpp
r24103 r24104 5 5 */ 6 6 7 #include "../io/io.h" #include <math.h> #include "../Numerics/types.h" 7 #include "../io/io.h" 8 #include <math.h> 9 #include "../Numerics/types.h" 8 10 9 11 IssmDouble NyeCO2(IssmDouble temperature){ 10 12 11 13 /*Coefficients*/ 12 const IssmPDouble Rg = 8.3144598; /* J mol^-1 K^-1 */ const13 IssmPDouble A_const = pow(10.,10.8); /*s^-1 MPa */ const IssmPDouble14 Q = 63000.; /*J mol^-1 */ const IssmPDouble n =15 14 const IssmPDouble Rg = 8.3144598; /* J mol^-1 K^-1 */ 15 const IssmPDouble A_const = pow(10.,10.8); /*s^-1 MPa */ 16 const IssmPDouble Q = 63000.; /*J mol^-1 */ 17 const IssmPDouble n = 7.; /*Glen's exponent*/ 16 18 17 19 /*Arrhenius Law*/ … … 20 22 21 23 /*Beyond-melting-point cases*/ 22 if((temperature>200.)&&(temperature<220.)) printf("CO2 ICE - POSSIBLE 23 MELTING. Some temperature values are between 200K and 220K.\n"); 24 else if(temperature>=220.) _printf0_("CO2 ICE - GUARANTEED MELTING. Some 25 temperature values are beyond 220K.\n"); 24 if((temperature>200.)&&(temperature<220.)) _printf0_("CO2 ICE - POSSIBLE MELTING. Some temperature values are between 200K and 220K.\n"); 25 else if(temperature>=220.) _printf0_("CO2 ICE - GUARANTEED MELTING. Some temperature values are beyond 220K.\n"); 26 26 27 27 /*Return output*/
Note:
See TracChangeset
for help on using the changeset viewer.