anchor.diff:  
   	use params
   	use utils1, only: logSumExp
  +	use Nested, only: MPI_COMM_WORLD
   
   	!use lib_matrix, only: inverse, determinant
   	implicit none
  
positive.code:  		tau = 1.d0
  
   		if (using_gp) then
  			slhood = gp1%get_lnlikelihood(times, rvs, Cube)
  			return
   		end if
  
  
negative.code: 		real(kind=8), intent(out) :: slhood
  		integer, intent(in) :: context
  		real(kind=8) :: lhood, lhood_test(1,1), det, jitter
  		integer :: i, n, ierr
  
  		! times, rvs and errors are defined in params and initialized/read in main
  		! Cube(1:nest_nPar) = P, K, ecc, omega, t0, Vsys
  