anchor.diff:  	use params
   	use gputils, only: gp_n_planet_pars
   	use utils1, only: logSumExp
  -	use Nested, only: MPI_COMM_WORLD
  +	!use Nested, only: MPI_COMM_WORLD
   
   	!use lib_matrix, only: inverse, determinant
   	implicit none
  
positive.code: 		allocate(times(N), rvs(N), errors(N))
  		allocate(train_var(N))
  		allocate(vel(N), dist(N), sigma(N))
  		allocate(covmat(N,N), inv_covmat(N,N))
  	end subroutine likelihood_init
  
  
negative.code: 		deallocate(observ)
  		deallocate(times, rvs, errors)
  		deallocate(vel, dist, sigma)
  		deallocate(last_vel_oversampled, vel_oversampled, times_oversampled)
  		deallocate(covmat, inv_covmat)
  	end subroutine likelihood_finish		
  
  