This one works for sure, tested against old Fortran code.