rendered paste body! 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