Changeset 16137 for issm/trunk/src/c/shared/Elements/Arrhenius.cpp
- Timestamp:
- 09/16/13 09:43:55 (12 years ago)
- Location:
- issm/trunk
- Files:
-
- 5 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk
- Property svn:ignore
-
old new 1 nightlylog 2 configure.sh 1 3 par 2 4 ad
-
- Property svn:mergeinfo changed
/issm/trunk-jpl merged: 15397-15401,15403-15487,15489-15701,15704-15735,15737-16076,16082-16133
- Property svn:ignore
-
issm/trunk/src
- Property svn:mergeinfo changed
-
issm/trunk/src/c
- Property svn:ignore
-
old new 14 14 probe.results 15 15 stXXXX* 16 16 .deps 17 .dirstamp
-
- Property svn:ignore
-
issm/trunk/src/c/shared/Elements
-
Property svn:ignore
set to
.deps
.dirstamp
-
Property svn:ignore
set to
-
issm/trunk/src/c/shared/Elements/Arrhenius.cpp
r15396 r16137 26 26 27 27 /*Some physical constants (Payne2000)*/ 28 IssmDouble beta=8.66* pow(10.,-4.);28 IssmDouble beta=8.66*1.e-4; 29 29 IssmDouble R=8.314; 30 30 … … 33 33 34 34 /*convert temperature to absolute temperature*/ 35 _assert_(depth );35 _assert_(depth>0); 36 36 Tstar=temperature-beta*depth; 37 _assert_(Tstar );37 _assert_(Tstar>0); 38 38 39 39 /*Get A*/ 40 40 if(Tstar<263.15){ 41 A=3.61 *pow(10.,-13.) * exp( -6.*pow(10.,4.)/(R*Tstar));41 A=3.61e-13*exp( -6.e+4/(R*Tstar)); 42 42 } 43 43 else{ 44 A=1.73 *pow(10., 3.) * exp(-13.9*pow(10.,4.)/(R*Tstar));44 A=1.73e+3 *exp(-13.9e+4/(R*Tstar)); 45 45 } 46 46
Note:
See TracChangeset
for help on using the changeset viewer.