18   integer, 
parameter :: b2gd_ntmax = 101
 
   19   integer, 
parameter :: b2gd_nvars = 4
 
   24   real :: b2gd_treal(b2gd_ntmax)
 
   25   real :: b2gd_xreal(b2gd_ntmax)
 
   29   real(s2_dp) :: b2gd_Omega_matter
 
   30   real(s2_dp) :: b2gd_Omega_Lambda
 
   31   real(s2_dp) :: b2gd_alpha, b2gd_Bianchi_h
 
   32   real(s2_dp) :: b2gd_RH_start
 
   33   real(s2_dp) :: b2gd_ze
 
   34   real(s2_dp) :: b2gd_theta_0, b2gd_phi_0     
 
   36   real(s2_dp) :: b2gd_xarr(b2gd_nvars, b2gd_ntmax)
 
   37   real(s2_dp) :: b2gd_tarr(b2gd_ntmax)
 
   41   real(s2_dp) :: b2gd_tstop
 
   42   real(s2_dp) :: b2gd_deltat