NyeCO2.cpp File Reference

#include "../io/io.h"
#include <math.h>
#include "../Numerics/types.h"

Go to the source code of this file.

Functions
IssmDouble	NyeCO2 (IssmDouble temperature)

Function Documentation

NyeCO2()

IssmDouble NyeCO2

(

IssmDouble

temperature

)

Definition at line 11 of file NyeCO2.cpp.

                                         {
 
   /*Coefficients*/
   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 */
 
   /*Beyond-melting-point cases*/
   if((temperature>200.)&&(temperature<220.)) _printf0_("CO2 ICE - POSSIBLE MELTING. Some temperature values are between 200K and 220K.\n");
   else if(temperature>=220.) _printf0_("CO2 ICE - GUARANTEED MELTING. Some temperature values are beyond 220K.\n");
 
   /*Return output*/
   return B; 
}