SUBROUTINE KSTPRZ(T0,T1,ZET,TETA,ZETA)                            00002720
C                                                                       00002730
C     BESTIMMUNG DER PRAEZESSINSWINKEL NACH DEN ANGABEN                 00002740
C     BEI LIESKE, ASTRON.&ASTROPHYS., 73, 282-284,(1979).               00002750
C                                                                       00002760
C     EINGABE: T0 = JULIANISCHES DATUM DER AUSGANGSEPOCHE               00002770
C              T1 = JULIANISCHES DATUM DER ENDEPOCHE                    00002780
C     AUSGABE: PRAEZESSIONSWINKEL IM BOGENMASS                          00002790
C                                                                       00002800
C                                                                       00002810
      IMPLICIT REAL*8 (A-H,O-Z)                                         00002820
C                                                                       00002830
      PI = 3.141592653589793D0                                          00002840
C                                                                       00002850
C     EPOCHENDIFFERENZEN IN JULIANISCHE JAHRHUNDERTE UMRECHNEN          00002860
C                                                                       00002870
      XJ2000 = 2451545.D0                                               00002880
      DT0    = (T0-XJ2000)/36525.D0                                     00002890
      DT1    = (T1-T0)/36525.D0                                         00002900
C                                                                       00002910
      ZETA = ( 2306.2181D0 + 1.39656D0*DT0 - 0.000139D0*DT0*DT0)*DT1 +  00002920
     1       ( 0.30188D0 - 0.000344D0*DT0)*DT1*DT1 +                    00002930
     2         0.017998D0*DT1*DT1*DT1                                   00002940
      ZET  = ( 2306.2181D0 + 1.39656D0*DT0 - 0.000139D0*DT0*DT0)*DT1 +  00002950
     1       ( 1.09468D0 + 0.000066D0*DT0)*DT1*DT1 +                    00002960
     2         0.018203D0*DT1*DT1*DT1                                   00002970
      TETA = ( 2004.3109D0 - 0.85330D0*DT0 - 0.000217D0*DT0*DT0)*DT1 +  00002980
     1       (-0.42665D0 - 0.000217D0*DT0)*DT1*DT1 -                    00002990
     2         0.041833D0*DT1*DT1*DT1                                   00003000
C                                                                       00003010
C     UMWANDLUNG INS BOGENMASS                                          00003020
C                                                                       00003030
      X = (180.D0*3600.D0)                                              00003040
      ZETA = (ZETA*PI)/X                                                00003050
      TETA = (TETA*PI)/X                                                00003060
      ZET  = (ZET *PI)/X                                                00003070
      RETURN                                                            00003080
      END                                                               00003090