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