Index: trunk/testsuite/ZH_implementation/makefile
===================================================================
--- trunk/testsuite/ZH_implementation/makefile	(revision 594)
+++ trunk/testsuite/ZH_implementation/makefile	(revision 595)
@@ -1,16 +1,16 @@
 F90C = gfortran
 F90FLAGS = -fbounds-check -ffixed-line-length-none
 
 HBLIBS = -L../../HiggsBounds-5/ -lHB
 HBINCLUDE = -I../../HiggsBounds-5/
 
 HSLIBS = -L../../HiggsSignals-2/ -lHS
 HSINCLUDE = -I../../HiggsSignals-2/
 
 
 
-test_ZH:
+test_ZH: test_ZH.f90
 	$(F90C) $(F90FLAGS) $(HBINCLUDE) test_ZH.f90 -o test_ZH $(HBLIBS)
 
 test_HS:
 	$(F90C) $(F90FLAGS) $(HBINCLUDE) $(HSINCLUDE) test_HS.f90 -o test_HS $(HBLIBS) $(HSLIBS)
Index: trunk/testsuite/ZH_implementation/test_ZH.f90
===================================================================
--- trunk/testsuite/ZH_implementation/test_ZH.f90	(revision 594)
+++ trunk/testsuite/ZH_implementation/test_ZH.f90	(revision 595)
@@ -1,69 +1,69 @@
 program test_ZH
 
  use theory_XS_SM_functions
  implicit none
  
  integer :: i
  double precision :: ghz,ght,ghb,gat,gab,mh,nnlo_sum,nnlo_incl
  
  call initialize_HiggsBounds(1,0,'LandH')
  
- ghz = 0.0D0
- ght = 0.0D0
- ghb = 0.0D0
+ ghz = 1.0D0
+ ght = 1.0D0
+ ghb = 1.0D0
  gat = 0.0D0
- gab = 10.0D0
+ gab = 1.0D0
  
 !  write(*,*) "Mass      NNLO,incl. k(NNLO/NLO,incl) NNLO(gg,qq,bb) k(NNLO/NLO,gg,qq,bb) ",&
 !  &          "k(NNLO/NLO,gg) k(NNLO/NLO,qq) k(NNLO/NLO,bb)"
  
 !   write(*,*) "Tevatron"
 !   do i = 1,10
 !   mh = 10.0D0 + (i-1)*25.0D0
 !   nnlo_sum = ZH_cpmix_nnlo_ggqqbb(mh,'TEV  ',ghz,ght,ghb,gat,gab)
 !   nnlo_incl = ZH_cpmix_nnlo(mh,'TEV  ',ghz,ght,ghb,gat,gab)  
 !   write(*,*) mh, nnlo_sum, nnlo_sum/nnlo_incl, XS_tev_HZ_SM(mh) 
 !   enddo  
 
   
 !   write(*,*) "LHC7"
 !   do i = 1,10
 !   mh = 10.0D0 + (i-1)*25.0D0
 !   nnlo_sum = ZH_cpmix_nnlo_ggqqbb(mh,'LHC7 ',ghz,ght,ghb,gat,gab)
 !   nnlo_incl = ZH_cpmix_nnlo(mh,'LHC7 ',ghz,ght,ghb,gat,gab)  
 !   write(*,*) mh, nnlo_sum, nnlo_sum/nnlo_incl, XS_lhc7_HZ_SM(mh,.True.,.False.), XS_lhc7_HZ_SM(mh,.True.,.True.)
 !   enddo  
 
 !   write(*,*) "LHC8"
 !   do i = 1,10
 !   mh = 10.0D0 + (i-1)*25.0D0
 !   nnlo_sum = ZH_cpmix_nnlo_ggqqbb(mh,'LHC8 ',ghz,ght,ghb,gat,gab)
 !   nnlo_incl = ZH_cpmix_nnlo(mh,'LHC8 ',ghz,ght,ghb,gat,gab)  
 !   write(*,*) mh, nnlo_sum, nnlo_sum/nnlo_incl, XS_lhc8_HZ_SM(mh,.True.,.False.), XS_lhc8_HZ_SM(mh,.True.,.True.)
 !   enddo  
 
   write(*,*) "LHC13"
   do i = 1,10
   mh = 10.0D0 + (i-1)*25.0D0
   nnlo_sum = ZH_cpmix_nnlo_ggqqbb(mh,'LHC13',ghz,ght,ghb,gat,gab)
   nnlo_incl = ZH_cpmix_nnlo(mh,'LHC13',ghz,ght,ghb,gat,gab)  
   write(*,*) mh, nnlo_sum, nnlo_sum/nnlo_incl, XS_lhc13_HZ_SM(mh,.True.,.False.), XS_lhc13_HZ_SM(mh,.True.,.True.),&
   & ZH_cpmix_nnlo_gg(mh,'LHC13',ghz,ght,ghb,gat,gab), ZH_cpmix_nnlo_bb(mh,'LHC13',ghb,gab)
   enddo  
 
 !   write(*,*) "Tevatron WH"
 !   do i = 1,10
 !   mh = 10.0D0 + (i-1)*25.0D0
 !   write(*,*) mh, WH_nnlo(mh,'TEV  ',ghz,ght,ghb), WH_nnlo_SM(mh,'TEV  '), XS_tev_HW_SM(mh,.True.)
 !   enddo 
 ! 
 !   write(*,*) "LHC 13"
 !   do i = 1,10
 !   mh = 10.0D0 + (i-1)*25.0D0
 !   write(*,*) mh, WH_nnlo(mh,'LHC13',ghz,ght,ghb), WH_nnlo_SM(mh,'LHC13'), XS_lhc13_HW_SM(mh,.True.,.False.)
 !   enddo 
 
   
 ! '(11F12.5)'
  call finish_HiggsBounds
 end program test_ZH
\ No newline at end of file