created PsiToV function, added read Jacobian in phase1