c---------------------------------
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
       PARAMETER ( pi =3.1415926535d0)
       PARAMETER ( sq2 = 1.414213562d0)
       PARAMETER ( tpisq = 2.d0*pi**2)
       Common /pars/ aexpn,hubble,Om0,OmL,sig8
       Common /Auxil/ xout,eps1,eps2, Rtophat,Anorm

