%                            pz90con.m
%  Scope:   This MATLAB macro sets the PZ-90 (earth's parameters 1990) most 
%           used constants as global variables.
%  Usage:   pz90con
%  Description of global parameters:
%           a_smaxis  -  output, PZ-90 earth semi-major axis in meters
%           b_smaxis  -  output, PZ-90 earth semi-minor axis in meters
%           eccentr   -  output, PZ-90 earth eccentricity 
%           eccentr2  -  output, PZ-90 earth eccentricity squared
%           flatness  -  output, PZ-90 flatness (ellipticity)
%           gravpar   -  output, PZ-90 gravity parameter in meters**3/sec**2
%                        earth's gravitational constant (mass of earth
%                        including earth's atmosphere)
%           gravatm   -  output, PZ-90 gravitational constant of atmosphere
%                        in m**3/s**2
%           rot_rate  -  output, PZ-90 earth rotation rate in radians/second
%           c_speed   -  output, PZ-90 speed of light in vacuum in meters/second
%           g0        -  output, PZ-90 equatorial acceleration of gravity in 
%                        meters/sec**2
%           gcoratm   -  output, PZ-90 correction to the acceleration of gravity
%                        at sea level because of atmosphere, in meters/sec**2
%           j2        -  output, PZ-90 second zonal harmonic of the geopotential
%  Reference:
%           [1] Langley, R. B., GLONASS: Review and update. GPS World, vol. 8,
%               September 1997, p. 71.
%  Last update: 10/15/97
%  Copyright (C) 1997 by LL Consulting. All Rights Reserved.

global  a_smaxis b_smaxis eccentr eccentr2 flatness gravpar gravatm rot_rate
global c_speed g0 gcoratm j2

%  Initialize the constants

a_smaxis = 6378136.;   
b_smaxis = 6356751.361745712;
eccentr  = 0.08181910652836;
eccentr2 = 0.00669436619310;     %  eccentr2 = flatness * (2. - flatness)
flatness = 0.00335280374302;     %  1. / 298.25784 = 1 - b_smaxis/a_smaxis

gravpar  = 3.9860044e+14;
gravatm  = 0.35e+9;
rot_rate = 7.292115e-5; 
c_speed  = 2.99792458e+8;
g0       = 9.780328;
gcoratm  = - 0.9e-5;
j2       = 1082.6257e-6;