c  
c 
      program  test
**************************************variables********************************
      implicit none
      real omega0,lambda0,h0,omegab,a,deltcrit
      real nspec,gamma,sigma8,sig,alpha,m,sigmacdm
      integer igwave,itrans
      common /cosmology/ omega0,lambda0,h0,omegab
      common /power_spectrum/ nspec,gamma,sigma8,igwave,itrans
      external deltcrit,sigmacdm
*******************************************************************************
      omega0=0.3
      lambda0=0.7
      h0=0.7
      omegab=0.05

      a=1.0

      write(0,*) 'delta_c= ',deltcrit(a)

      nspec=1.0
      gamma=omega0*h0
      sigma8=0.9
      igwave=0
      itrans=1

      m=1.0e+12
      sig=sigmacdm(m,alpha)
      write(0,*) 'sigma(m)= ',sig,' alpha=',alpha

      stop
      end
c------------------------------------------------------------------------------