program Tester

        use ModelData
        use CAMBinterface
        implicit none
        integer l,Maxl
        real*8 Maxk
	character(LEN=80) :: BesselsFile

        type(CAMBparams)  P

        Maxl=1500
        Maxk=3000

	BesselsFile='spherbes.dat'
        
        call Initialize(BesselsFile,Maxl,Maxk)
      
        
        P%Out_scalar = .true.
        P%Out_tensor = .false.

        P%Om_baryon = .045
        P%Om_cdm = 0.355
        P%Om_lambda = 0.6
        P%Om_nu = 0.0
        P%Hubble_0=65
       
        P%T_CMB=2.726
        P%Y_He=0.24
        P%Num_Nu_massless=3.04
        P%Num_Nu_massive=0

        P%Scalar_spectral_index=1
        P%Scalar_initial_condition=initial_adiabatic

        P%Tensor_spectral_index=1
        
        P%RECFAST_recombination =.true.
        P%Reionization = .false.
     
        P%Max_l=1500
        P%Max_k=3000
        
        call ZeroCls !clear the Cl arrays (in modules.f90, ModelData)    
  
        call GetCls(P)
       
        do l=2,1500
            print *, l,clts(l,1) !print out temp Cls
        end do
        

        end program Tester