move (correct) Jacobian logic into correct place: read_corr_matrix()
change read_corr_L() so that it can be used independently