c
c       Berechnung von g nach Raabe/Schoenfeldt
c   
        real k1,k2
        pi=acos(-1.)
1       write(*,*) "Eingabe geogr. Breite (Grad), Hoehe (km)"
        write(*,*) "Ende=100,100"
        read(*,*) phi,z
        if(phi.gt.95) stop
        phi=phi*pi/180.

        g0=9.8063
        k1=0.0026373
        k2=0.0000059
        g=g0*(1-pi*(z/10000))*(1-k1*cos(2.*phi)+
     +  k2*(cos(2.*phi))**2.)
        write(*,*) g
        goto 1
        end