#include <cmath>
#include <iostream>
#include <fstream>
#include <vector>
#include "earthmodel.hh"
Go to the source code of this file.
|
const double | FIRNDEPTH =-150. |
|
const double | scale_factor =0.97276901289 |
|
const double | ellipsoid_inv_f = 298.257223563 |
|
const double | ellipsoid_b = EarthModel::R_EARTH*(1-(1/ellipsoid_inv_f)) |
|
const double | eccentricity = sqrt((1/ellipsoid_inv_f)*(2-(1/ellipsoid_inv_f))) |
|
const double | bedmap_a_bar = pow(eccentricity,2)/2 + 5*pow(eccentricity,4)/24 + pow(eccentricity,6)/12 + 13*pow(eccentricity,8)/360 |
|
const double | bedmap_b_bar = 7*pow(eccentricity,4)/48 + 29*pow(eccentricity,6)/240 + 811*pow(eccentricity,8)/11520 |
|
const double | bedmap_c_bar = 7*pow(eccentricity,6)/120 + 81*pow(eccentricity,8)/1120 |
|
const double | bedmap_d_bar = 4279*pow(eccentricity,8)/161280 |
|
const double | bedmap_c_0 = (2*EarthModel::R_EARTH / sqrt(1-pow(eccentricity,2))) * pow(( (1-eccentricity) / (1+eccentricity) ),eccentricity/2) |
|
const string | crust20_in ="data/outcr" |
|
const string | crust20_out ="altitudes.txt" |
|
const string crust20_in ="data/outcr" |
const string crust20_out ="altitudes.txt" |
const double ellipsoid_inv_f = 298.257223563 |
const double FIRNDEPTH =-150. |
const double scale_factor =0.97276901289 |