! Comp berechnen
     DO ii = 1, anz
        comp(ii) = acos(sin(delta(ii))*sin(par(2)) &
             & + cos(delta(ii))*cos(par(2))*cos(t(ii))) &
             & - par(3)*tan(z(ii)) + (2*mod(ii,2)-1)*par(4)
     END DO

! Jacobi Matrix
     DO ii = 1,anz
        t(ii) = theta0 + (1+mue)*ut(ii) + par(1) - alpha(ii)
        azi = atan(cos(delta(ii))*sin(t(ii))/ &
             & (cos(delta(ii))*sin(par(2))*cos(t(ii)) - sin(delta(ii))*cos(par(2))))
!       jacobi(ii,1) = -sin(t(ii))*sin(azi)
        jacobi(ii,1) = -cos(par(2))*sin(azi)
        jacobi(ii,2) = -cos(azi)
        jacobi(ii,3) = -tan(z(ii))
        jacobi(ii,4) = 2*mod(ii,2)-1
     END DO