// physical and mathematical constants

const double c_pi    = 3.141594; // Pi
const double c_light =   3.0E+5; // km/s
const double hbarc   = 197.3E-3; // (eV^2/GeV)km

// //-----------------------------------------------------
//-     vacuum oscillations parameters (19.April.2012)
double theta_12 = 34.0;    // degrees
double theta_13 = 8.93;    // degrees
double theta_23 = 45.0;    // degrees
double dmsq21   = 7.59E-5; // eV^2
double dmsq31   = 2.35E-3; // eV^2
double deltaCP  = 0.0;     // degrees

// //-----------------------------------------------------
 //-           model of the Earth's interior
double R_Earth = 6378.137.;   // km

const int Nshells = 2;

double   R_shell[Nshells] = {R_Earth/2,//km         CORE
			     R_Earth}; //km 

double   V_shell[Nshells] = {4.54323E-4, //eV^2/GeV CORE
			     1.89301E-4};//eV^2/GeV MANT

//-----------------------------------------------------
//   auxiliary stuff

double tf_max = 1.0*R_Earth/c_light;
double Enu_max = 550; // GeV

//-----------------------------------------------------
// a few global variables

double l_ss[Nshells];
double tcross[Nshells];
double times[2*Nshells];
double tlohi[2];
bool first;
double theta_13_bkp;
bool first13;