implicit real*8 (a-h,o-z)
      character*1 sign
      pi = 4.d0*atan(1.d0)
      dtr = 180.d0/pi
C  Polaris:        $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
       Jahr = 2012
       Monat =   11                                                             Ausser Sternnummer 907 einziger
       itag  =  2                                                              variabler Eingabewert
       ihour  = 0
C Nachfolgend Astrometrie und Orbitelemente Polaris: Wielen et al. Table 1  und Table 4 retrograde
C
       alpha0 =             +037.95352509  d0                                   Grad     RA  (J2000),
       delta0 =              +89.26410349  d0                                   Grad     Dec (J2000),
       rmuafk =              +40.81        d0                                   mas/yr,  enthaelt noch cos(delta) !!!     My*-RA  (J2000),
C      rmudfk =              -15.23        d0                                   mas/yr   My -Dec (J2000),    ==> bisheriger Wert
       rmudfk =              -15.22        d0                                   mas/yr   My -Dec (J2000),
       parfk  =                7.56        d0                                   mas/yr   Parallaxe
C      rvfk   =             - 16.4         d0                                   km/s     Rad.Vel.            ==> bisheriger Wert
       rvfk   =             - 16.42        d0                                   km/s     Rad.Vel.
C
       period =               29.59        d0                                   Period
       ecc    =                0.608       d0                                   Eccent
       tper   =             1987.66        d0                                   T-Perihel
       xnode  =              167.100       d0                                   Node
       xomega =              303.010       d0                                   Omega
       xincl  =              130.2         d0                                   Incl
       axis   =               21.91        d0                                   Major Axis
C   Umrechnungen vor SBR  cappplc Aufruf
       parfk = parfk/1000.d0
C+++++++++++++++++++++++++
      call       cappplc(alpha0, delta0, rmuafk, rmudfk,
     * parfk, rvfk, Jahr, Monat, itag, ihour, racioout, raequout,
     * decout)
C
      call orbitdelta (period,tper,ecc,axis,xincl,xnode,xomega,
     &              alpha0,delta0,rmuafk,rmudfk,parfk,rvfk,
     &              jahr, monat, itag, ihour,
     &              dalf,ddel)
      alpha = raequout + dalf
      delta = decout   + ddel
C
C Ab hier Ausgabeteil zur Kontrolle
C dalf bzw. ddel sind zu recioout,raequout bzw. decout zu addieren
C
C Ausgabe
      alpha = alpha / 15.d0
      write(6,*) alpha
C
      write(6,*) delta
      idega = int(alpha)
      imina= int(alpha*60.d0 - dble(idega)*60.d0)
      seca = abs(alpha)*3600.d0
     & -dble(iabs(idega))*3600.d0-dble(imina)*60.d0
      write(6,999) idega , imina, seca
999   format(' alpha = ',I4,2x,i2,2x,f10.4,'    Equinox-Method')
      idegd = int(delta)
      imind = int(abs(delta*60.d0) - dble(iabs(idegd))*60.d0)
      secd  = abs(delta)*3600.d0
     &   -dble(iabs(idegd))*3600.d0-dble(imind)*60.d0
      sign = ' '
      if(delta.lt.0.d0) sign='-'
      write(6,998) sign,iabs(idegd) , imind, secd
998   format(' delta = ',1x,a1,i2,2x,i2,2x,f10.4,'    Equinox-Method')
C
c Step 5b  (CIO-Method)
      alpha = racioout + dalf                                                   27.06.11
      delta = decout   + ddel                                                   27.06.11
      alpha = alpha / 15.d0
      write(6,*) alpha

      write(6,*) delta
      idega = int(alpha)
      imina= int(alpha*60.d0 - dble(idega)*60.d0)
      seca = abs(alpha)*3600.d0
     &      -dble(iabs(idega))*3600.d0-dble(imina)*60.d0
      write(6,888) idega , imina, seca
      imind = int(abs(delta*60.d0) - dble(iabs(idegd))*60.d0)
      secd  = abs(delta)*3600.d0
     &      -dble(iabs(idegd))*3600.d0-dble(imind)*60.d0
      sign = ' '
      if(delta.lt.0.d0) sign='-'
      write(6,887) sign,iabs(idegd) , imind, secd
888   format(' alpha = ',I4,2x,i2,2x,f10.4,'        CIO-Method')
887   format(' delta = ',1x,a1,i2,2x,i2,2x,f10.4,'        CIO-Method')
      stop
      end
      include 'apfscompute.f'
      include 'iau_subroutines.f'
      include 'sbr_pleph.f'
      include 'orbitcompute.f'
      include 'kstprz.f'
      include 'praeznew.f'