loglik function

Wrapper to the Fortran subroutines computing the log-likelihood