#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 |