C
C     FORTRAN II FOR THE IBM 7094 FOR A HORIZONTAL DIAL
C
C     RIGHT NOW, DATA IS ENTERED INTO THE THREE STATEMENTS FOR
C                LAT
C                LNG
C                REF
      INTEGER I,J,SIZE
      REAL LAT,SLAT,HOUR,THOUR,LHA,ATNHA,DLHA
      REAL LNG,REF
C      
      WRITE(6,900)
      WRITE(6,901)
      WRITE(6,900)
      WRITE(6,902)
C     CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C     CCC     E N T E R    D A T A    B E L O W    P L E A S E     CCC
C     CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
      LAT = 32.75 
      LNG = 108.2
      REF = 105.0
C     CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C     FOR LAT 32.75 LNG 108.2 REF 105.0 THE RESULTS SHOULD BE AS BELOW
C     CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C     LAT      Hour:  -6 hour angle:  93.20  Angle: -84.10
C     32.75    Hour:  -5 hour angle:  78.20  Angle:  68.88
C              Hour:  -4 hour angle:  63.20  Angle:  46.96
C     LNG      Hour:  -3 hour angle:  48.20  Angle:  31.18
C     108.2    Hour:  -2 hour angle:  33.20  Angle:  19.49
C              Hour:  -1 hour angle:  18.20  Angle:  10.09
C     REF      noon
C     105.0    Hour:   0 hour angle:   3.20  Angle:   1.73
C              noon
C              Hour:   1 hour angle: -11.80  Angle:  -6.45
C              Hour:   2 hour angle: -26.80  Angle: -15.28
C              Hour:   3 hour angle: -41.80  Angle: -25.81
C              Hour:   4 hour angle: -56.80  Angle: -39.58
C              Hour:   5 hour angle: -71.80  Angle: -58.71
C              Hour:   6 hour angle: -86.80  Angle: -84.10
C     CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
      DO 10 I = 1,13
      J = I-7
C
C     J IS INTEGER HOUR FOR DISPLAYING
C
C     NOTE THE USE OF TRAILING .0 ETC BECAUSE OF TYPE CONVERSION RULES
C          IN FORTRAN II AS IMPLEMENTED IN THE 7079 COMPUTER SERIES
C     TAN(TIME)
C
      HOUR  = J
      DLHA  = (15.0 * HOUR) - (LNG - REF)
      LHA   = (DLHA * 2.0 * 3.1416) / 360.0
C
C     FORTRAN II ON THE IBM 7090 HERE HAS NO TAN FUNCTION HENCE SIN COS
      THOUR = SIN(LHA) / COS(LHA)
      SLAT  = SIN((LAT * 2.0 * 3.1416) / 360.0)
      ATNHA = SLAT * THOUR
      ATNHA = ATAN(ATNHA) * 360.0 / (2.0 * 3.1416)
      WRITE (6,905) J, DLHA, ATNHA
C
   10 CONTINUE
C
      WRITE(6,900)
      WRITE(6,901)
      WRITE(6,900)
      WRITE(6,906)
      WRITE(6,907) LAT, LNG, REF
      WRITE(6,900)
C
   50 STOP
C
C     NOTE - THE LENGTH BEFORE AN H IS CRITICAL - EG 45H 36H ETC
  900 FORMAT(1X, 45HXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX)
  901 FORMAT(1X, 45HXX HORIZONTAL DIAL - FORTRAN II - IBM 7094 XX)
  902 FORMAT(1X, 45HHOUR-NOON       HOUR ANGLE    HOUR LINE ANGLE)
  903 FORMAT(I4)
  904 FORMAT(F10.4)
  905 FORMAT(I5, F18.1, F15.2)
  906 FORMAT(1X, 36HLATITUDE           LONGITUDE     REF) 
  907 FORMAT(F6.2,F20.2,F15.2) 
      END

