Common function for calculation of atmospheric averages:
func=1: muon production height
func=2: 1/cos~1/rho
func=3: muon track length
all others: normalization.
Class contains all parametrization table io routines, history/debugging info/secondary
output definitions, and a fast double precision variable formatting routine