Hello, IMHO in dynamics_module.f90 line 1339 IF ( mod( nrefresh, istep ) == 0 ) inv_hess(:,:) = identity( dim ) should rather be IF ( mod( istep, nrefresh ) == 0 ) inv_hess(:,:) = identity( dim ) Please fix... All the best, ciao, Nick