home *** CD-ROM | disk | FTP | other *** search
/ ftp.barnyard.co.uk / 2015.02.ftp.barnyard.co.uk.tar / ftp.barnyard.co.uk / cpm / walnut-creek-CDROM / ENTERPRS / CPM / UTILS / F / F80.LBR / PL1.FOR < prev    next >
Text File  |  2000-06-30  |  14KB  |  381 lines
  1.  
  2.  
  3. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  4. C
  5. C Program PLANETF
  6. C
  7. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  8.         PROGRAM XPLANET
  9.         BYTE COMENT(25),FILNAM(25),OPTION,
  10.      1  RETRO(8),BLAH
  11.         LOGICAL HLP
  12.         DOUBLE PRECISION T,AN1(9),XJD,CONV,GMT,XL,YL,XP,YP,XN,
  13.      1  DAY,PI2,RULER(12),WEEK(7),WEEKR(7),
  14.      2  YBES,TT
  15.         DIMENSION EC0(9),EC1(9),AN0(9),P0(9),P1(9),EC(9),
  16.      1  TL(9),TB(9),TH0(9),TH1(9),XI(9),R(9),A(9),G(9),GE(9),
  17.      2  X(9),Y(9),Z(9),XX(9),YY(9),ZZ(9),IH(4,3),H(4),
  18.      3  ZOD1(12),P(120,14),ID(120,14),IM(120,14),IS(120,14),
  19.      4  GLYPH(120,14),IASP(14,14)
  20.      8  ,MONTH(12),MOND(12),XNODE(2,10),YNODE(2,10)
  21. CCCCC8CCCC
  22. C
  23. CCCCCCCCC
  24.         DATA MOND/31,0,31,30,31,30,31,31,30,31,30,31/
  25.         DATA MONTH/'JA','FE','MR','AP','MY','JN','JL','AU',
  26.      1  'SE','OC','NO','DE'/
  27.         DATA WEEK/'Monday  ','Tuesday ','Wednsday',
  28.      1  'Thursday','Friday  ','Saturday','Sunday  '/
  29.         DATA WEEKR/'Moon)   ','Mars)   ','Mercury)','Jupiter)',
  30.      1  'Venus)  ','Saturn) ','Sun)    '/
  31.         DATA ZOD/'Aries   ','Taurus  ','Gemini  ','Cancer  ',
  32.      1  'Leo     ','Virgo   ','Libra   ','Scorpio ',
  33.      2  'Sagtarus','Caprcorn','Aquarius','Pisces  '/
  34.         DATA ZOD1/'A','T','G','C','L','V','=','S','/','K','Q','P'/
  35. CCCCCCCCC
  36. C
  37. C     PLUTO ELEMENTS FROM SHARAF (1964)
  38. C
  39. CCCCCCCCC
  40.         DATA A/1.00000023,.387098599,.723331619,1.523688395,
  41.      1  5.202802875,9.53884320,19.19097811,30.0706724,39.672599/
  42.         DATA AN0/358.475833,102.279381,212.603222,319.529425,
  43.      1  225.444651,175.758444,74.313628,41.269550,231.002308/
  44.         DATA AN1/35999.049750D0,149472.515289D0,58517.803875D0,
  45.      1  19139.858500D0,3034.906654D0,1222.116782D0,428.502578D0,
  46.      2  218.466783D0,144.072477D0/
  47.         DATA EC0/.01675104,.20561421,.00682069,.09331290,
  48.      1  .04825382,.05606075,.04704433,.00853341,.24706226/
  49.         DATA EC1/-.00004180,.00002046,-.00004774,.00009206,
  50.      1  .0,.0,.0,.0,.0/
  51.         DATA P0/101.220833,75.899697,130.163833,334.218203,
  52.      1  11.907422,90.110981,169.048778,43.755611,221.592475/
  53.         DATA P1/1.719175,1.555489,1.408036,1.840758,
  54.      1  .0,.0,.0,.0,1.388888/
  55.         DATA TH0/0.,47.145944,75.779647,48.786442,
  56.      1  98.932822,112.347606,73.490250,130.678889,108.937165/
  57.         DATA TH1/0.,1.185208,.899850,.770992,
  58.      1  .0,.0,.509667,1.100972,1.358056/
  59.         DATA XI/0.,7.002881,3.393630,1.850333,
  60.      1  1.311614,2.494239,.7726658,1.779256,17.109816/
  61. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  62. C
  63. C Code
  64. C
  65. CCCCCCCCC
  66.         TAN(DUMMY)=SIN(DUMMY)/COS(DUMMY)
  67.         ARG(P,Q)=P-Q*INT(P/Q)-Q*INT(SIGN(.5,P-Q*INT(P/Q))-.5)
  68.         HLP=.FALSE.
  69.         LINEP=0
  70.         PI=ATAN(1.0)
  71.         PI2=6.283185307179586476D0
  72.         CONV=PI2/360.D0
  73.         CONVS=CONV
  74.         DO 130 I=1,120
  75.         DO 130 J=1,14
  76. 130     P(I,J)=0.
  77.         WRITE(5,10000)
  78. 10000   FORMAT('$Sample input (Gregorian date,time',
  79.      1  'conversion to GMT,latitude,longitude,OPTION)'//
  80.      2  ' 1899/12/31 23:59 +08:00 34N59 122W59) : ')
  81. CCCCCCCCC
  82. C
  83. C     HOROSCOPE LOOP STARTS HERE
  84. C
  85. CCCCCCCCC
  86.         READ(5,11000)NY,CALEND,NM,ND,NH,NMI,NCH,NCM,LATD,LATS ,
  87.      1  LATM,LONGD,LONGS,LONGM,OPTION
  88. 11000   FORMAT(I4,A1,4(I2,1X),I3,2(1X,I2),A1,I2,1X,I3,A1,I2,1X,1A)
  89. 6011    CONTINUE
  90.         WRITE(6,11005)
  91. 11005   FORMAT('$Do you care for a hardcopy ? ')
  92.         READ(5, 11006) BLAH
  93. 11006   FORMAT(A1)
  94.         IF (BLAH.NE.'Y'.AND.BLAH.NE.'y') GO TO 11035
  95.         CONTINUE
  96.         CALL OPEN(6,'PLANET.DAT')
  97.         IOUT=6
  98.         GO TO 11008
  99. 11035   IOUT=5
  100. 11008   LINEP=1
  101.         WRITE(IOUT,11500)NY,CALEND,NM,ND,NH,NMI,NCH,NCM,LATD,
  102.      1  LATS,LATM,LONGD,LONGS,LONGM
  103. 11500   FORMAT('1'//' ---Date--- -Time -Zone- -Lat- -Long- ',
  104.      1  'Comments--------'/' ',I4,A1,I2,'/',I2,I3,':',I2,
  105.      1  I4,':',I2,I3,A1,I2,I4,A1,I2)
  106. 144     GMT=((NH+NCH-12)*60+NMI+SIGN(NCM*1.,NCH*1.))/1440.
  107.         XJD=(ND-32075+1461*(NY+4800+(NM-14)/12)/4+367*(NM-2-
  108.      1  (NM-14)/12*12)/12-3*((NY+4900+(NM-14)/12)/100)/4)+GMT
  109.         IF(CALEND.NE.'B')GO TO 150
  110.         NY=1-NY
  111. 150     IF(CALEND.NE.'J'.AND.CALEND.NE.'B')GO TO 160
  112.         XJD=(ND-32075+1461*(NY+4800+(NM-14)/12)/4
  113.      1  +367*(NM-2-(NM-14)/12*12)/12-38)+GMT
  114. 160     T=(XJD-2415020.D0)/36525.D0
  115. C
  116.         JFINAL=101
  117. CCCCCCCCC
  118. C
  119. C     Planetary Ephemerides comxzted below (Newcomb/Hill,1898)
  120. C
  121. CCCCCCC
  122. 170     DO 270 J=1,JFINAL
  123.         TT=T+(J-1)/36525.D0
  124.         TP=TT+18262.D0/36525.D0
  125.         W=TT
  126.         DO 174 I=1,9
  127.         EC(I)=EC0(I)+EC1(I)*W
  128.         N=AN1(I)*TT/360.D0
  129. 174     G(I)=(AN1(I)*TT-N*360.D0+DBLE(AN0(I)))*CONV
  130.         G5=G(5)
  131.         G6=G(6)
  132.         G7=(220.169542+428.49311*TP)*CONVS
  133.         G(5)=G(5)+(.6506*SIN(2*G6-2*G5+336.9*CONVS)
  134.      1  +(3.9987-.002213*36525./4332.58*TP)
  135.      2  *SIN(5*G6-2*G5+(67.15-8197.0/3600.*TP)*CONVS)
  136.      3  +.5380*SIN(5*G6-3*G5+176.5*CONVS)
  137.      4  +.4112*SIN(2*G6-G5+1.4*CONVS)
  138.      4  +.0399278*36525./4332.58*TP*SIN(-G5+227.46*CONVS)
  139.      5  +.2763*SIN(3*G6-2*G5+127.4*CONVS)
  140.      6  +.2669*SIN(G6-G5+79.2*CONVS)  )*299.12837/3600.*CONVS
  141. CCCCCCCCC
  142. C
  143. CCCCCCCCC
  144.         G(6)=G(6)+(24.153*SIN(5*G6-2*G5+(247.11-2.277*TP)*CONVS)
  145.      1  +5.679*SIN(4*G6-2*G5+277.39*CONVS)
  146.      2  +3.505*SIN(2*G6-G5+181.43*CONVS)
  147.      3  +.657765*36525./10759.20*TP*SIN(G6+238.0*CONVS)
  148.      4  +.278*SIN(3*G6-G5+121.2*CONVS)
  149.      5  +.266*SIN(2*G6-2*G5+157.0*CONVS)
  150.      6  +.238*SIN(6*G6-2*G5-3*G7+6.9*CONVS)
  151.      7  +.223*SIN(10*G6-4*G5+(133.6-14814.5/3600.*W)*CONVS)
  152.      8  +.234*SIN(3*G7-G6+321.7*CONVS)  )*120.455/3600.*CONVS
  153. CCCCCCCCC
  154. C
  155. C
  156. CCCCCCCCC
  157.         G(7)=G(7)+33.086*W*W/3600.*CONVS
  158.         G(8)=G(8)-22.401/3600.*W*W*CONVS
  159.         DO 177 I=1,9
  160.         GE(I)=G(I)+EC(I)*SIN(G(I)+EC(I)*SIN(G(I)+EC(I)*SIN(
  161.      1  G(I)+EC(I)*SIN(G(I)+EC(I)*SIN(G(I))))))
  162. 177     X(I)=2.*ATAN(SQRT((1.+EC(I))/(1.-EC(I)))*TAN(.5*GE(I)))
  163.      1  +(P0(I)+P1(I)*W)*CONVS + (I/7-I/9)*W*5025.3/3600.*CONVS
  164.         GT=(358.415+35998.928*W)*CONVS
  165.         GJ=(225.209+3034.462*W+.332*SIN((134.4+38.5*W)*CONVS))
  166.      1  *CONVS
  167.         X(4)=X(4)+(25.384*COS(GJ-G(4)-48.9*CONVS)
  168.      1  +52.490*SIN((47.48+19.771*W)*CONVS)-37.05-13.50*W
  169.      2  +21.869*COS(2*GJ-G(4)-188.3*CONVS)
  170.      3  +16.035*COS(2*GJ-2*G(4)-191.9*CONVS)
  171.      4  +13.966*COS(-GT+2*G(4)-20.5*CONVS)
  172.      5  +8.559*COS(-GT+G(4)-35.1*CONVS)  )/3600.*CONVS
  173. CCCCCCCCC
  174. C
  175. C
  176. CCCCCCCCC
  177.         GN=(225.417+3034.904*W)*CONVS
  178.         GG=(175.753+1222.113*W)*CONVS
  179.         G1=(74.412+428.498*W)*CONVS
  180.         GP=(74.320+428.498*W)*CONVS
  181.         GPP=(41.339+218.467*W)*CONVS
  182.         X(7)=X(7)+(142.938*SIN(GG-2*G1)+19.508*COS(GG-2*G1)
  183.      1  +75.70*COS(3*G1-GG)-102.30*SIN(3*G1-GG)
  184.      2  -48.623*SIN(GN-G1)-21.320*COS(GN-G1)
  185.      3  -27.871*COS(GP-GPP)+19.869*SIN(GP-GPP)
  186.      4  +28.793*COS(2*GP-2*GPP)+10.035*SIN(2*GP-2*GPP)
  187.      5  +(18.37*SIN(3*GP-3*GPP)+8.91*COS(3*GP-3*GPP))*COS(G(7))
  188.      6  +(8.35*SIN(3*GP-3*GPP)-16.44*COS(3*GP-3*GPP))*SIN(G(7))
  189.      7  -18.585*COS(GG-G1)+12.603*SIN(GG-G1)
  190.      8  +4.327*COS(3*GP-3*GPP)+14.280*SIN(3*GP-3*GPP)
  191.      9  )/3600.*CONVS
  192. CCCCCCCCC
  193. C
  194. C
  195. CCCCCCCCC
  196.         X(7)=X(7)+((112.317*W-1.551*W*W-.516*W*W*W)*SIN(G(7))
  197.      1  -(68.339*W+9.721*W*W)*COS(G(7))+6.605*W*SIN(2*G(7))
  198.      2  -(29.44-.410*W)*SIN((20.45-22.61*W)*CONVS) )/3600.*CONVS
  199. CCCCCCCCC
  200. C
  201. C     JUP/SAT TERM & QUADRATURES OF LONG PERIOD UR`WNEP TERMS
  202. C
  203. CCCCCCCCC
  204.         X(8)=X(8)+(  33.972*W*COS(G(8))
  205.      1  +18.553*SIN((180.966+1004.034*W+.1403*W*W)*CONVS)
  206.      2  +34.138*SIN((153.267+2816.296*W-.0573*W*W)*CONVS)
  207.      3  +(26.50*W+3.92*W*W)*SIN(G(8))   )/3600.*CONVS
  208. CCCCCCCCC
  209. C
  210. C
  211. CCCCCCCCC
  212.         DO 185 I=1,9
  213.         TH=(TH0(I)+TH1(I)*W)*CONVS
  214.         DO 180 K=1,2
  215. CCCCCCCCC
  216. C
  217. C     MOTION OF JUP/SAT NODES FROM LEVERRIER
  218. C
  219. CCCCCCCCC
  220.         TH2=TH+(I/5*3636.6-I/6*493.1-I/7*3143.5)/3600.*TP*CONVS
  221.         R(I)=A(I)*(1.-EC(I)*COS(2.*ATAN(TAN((TH2+(K-1)*PI-(P0(I)
  222.      1  +P1(I)*W)*CONVS)*.5)*SQRT((1.-EC(I))/(1.+EC(I))))))
  223.         XNODE(K,I)=R(I)*COS(TH2+(K-1)*PI)
  224. 180     YNODE(K,I)=R(I)*SIN(TH2+(K-1)*PI)
  225.         TL(I)=TH+ATAN2(SIN(X(I)-TH)*COS(XI(I)*CONVS),
  226.      1  COS(X(I)-TH))+(I/5-I/7)*TP*5026.1/3600.*CONVS
  227.         TB(I)=ASIN(SIN(XI(I)*CONVS)*SIN(X(I)-TH))
  228. 185     R(I)=A(I)*(1.-EC(I)*COS(GE(I)))
  229.         TB(5)=TB(5)+4.37431*36525./4332.58*SIN(X(5)+23.62*CONVS)
  230.      1  *TP/3600.*CONVS
  231.         TB(6)=TB(6)+24.266*36525./10759.2*SIN(X(6)-13.05*CONVS)
  232.      1  *TP/3600.*CONVS
  233. CCCCCCCCC
  234. C
  235. C     TB(5) EFFECTS LESS THAN 35"T, TB(6) LESS THAN 84"T
  236. C
  237. CCCCCCCCC
  238.         R(5)=R(5)*10.**(.0002303*COS(2*G6-2*G5+336.9*CONVS)
  239.      1  +.0001679*COS(5*G6-3*G5+176.4*CONVS)
  240.      2  +.0000125634*36525./4332.58*TP*COS(-G5+227.4*CONVS)  )
  241. CCCCCCCCC
  242. C
  243. C       22",16",10"T,(7",5",3")
  244. C
  245. CCCCCCCCC
  246.         R(6)=R(6)*10.**(.0007005*COS(4*G6-2*G5+277.3*CONVS)
  247.      1  +.0003783*COS(G6-G5+79.8*CONVS)
  248.      2  +.000083491*36525./10759.20*TP*COS(G6+58.0*CONVS)
  249.      3  +.0002443*COS(2*G6-G5+176.0*CONVS)  )
  250. CCCCCCCCC
  251. C
  252. C
  253. CCCCCCCCC
  254.         DO 190 I=1,9
  255.         X(I)=R(I)*COS(TB(I))*COS(TL(I))
  256.         Y(I)=R(I)*COS(TB(I))*SIN(TL(I))
  257. 190     Z(I)=R(I)*SIN(TB(I))
  258.         DO 210 I=2,9
  259.         DO 200 K=1,2
  260.         XNODE(K,I)=XNODE(K,I)-X(1)
  261.         YNODE(K,I)=YNODE(K,I)-Y(1)
  262. 200     P(K+17,I)=ARG(ATAN2(YNODE(K,I),XNODE(K,I))/6.2831853,1.)
  263.         XX(I)=X(I)-X(1)
  264.         YY(I)=Y(I)-Y(1)
  265. 210     ZZ(I)=Z(I)-Z(1)
  266.         XX(1)=-X(1)
  267.         YY(1)=-Y(1)
  268.         ZZ(1)=-Z(1)
  269.         XNUT=-17.2327/1296000.*SIN((259.18-1934.142*W)*CONVS)
  270.         DO 230 I=1,9
  271.         P(J+19,I)=ARG(ATAN2(YY(I),XX(I))/6.28318531+XNUT,1.)
  272.         IF(J-1)230,220,230
  273. 220     DIST=SQRT(XX(I)*XX(I)+YY(I)*YY(I)+ZZ(I)*ZZ(I))
  274.         P(1,I)=ASIN(ZZ(I)/DIST)/CONVS
  275.         E=(23.4522944-.0130125*W+.002558*COS((259.18-1934.142
  276.      1  *W)*CONVS))*CONVS
  277.         P(2,I)=ASIN(COS(P(1,I)*CONVS)*SIN(P(20,I)*6.28318531)
  278.      1  *SIN(E)+SIN(P(1,I)*CONVS)*COS(E))/CONVS
  279.         P(6,I)=ARG(TL(I)/6.28318531,1.)
  280.         P(3,I)=TB(I)/CONVS
  281.         P(4,I)=ASIN(COS(TB(I))*SIN(TL(I))*SIN(E)
  282.      1  +SIN(TB(I))*COS(E))/CONVS
  283. 230     CONTINUE
  284. CCCCCCCCC
  285. C
  286. C     LUNAR EPHEMERIDES COMPUTED TO WITHIN 1' (BROWN/I.L.E.)
  287. C
  288. CCCCCCCCC
  289.         DAY=XJD-2415020.D0+(J-1.D0)
  290.         XL=.751206D0+DAY*.0366011014634D0
  291.         YL=.776935D0+DAY*.0027379092649D0
  292.         XP=.928693D0+DAY*.0003094557786D0
  293.         YP=.781169D0+DAY*.0000001307457D0
  294.         XN=.719954D0-DAY*.0001470942283D0
  295.         AL=(XL-IDINT(XL))*PI2
  296.         BL=(YL-IDINT(YL))*PI2
  297.         AP=(XP-IDINT(XP))*PI2
  298.         BP=(YP-IDINT(YP))*PI2
  299.         AN=(XN-IDINT(XN))*PI2
  300.         U=AL-AP
  301.         V=BL-BP
  302.         F=AL-AN
  303.         D=AL-BL
  304.         DL=22639*SIN(U)-4586*SIN(U-D-D)+2370*SIN(D+D)+769*SIN(U+U)
  305.      1  -668*SIN(V)-412*SIN(F+F)-212*SIN(U+U-D-D)-206*SIN(U+V-D-D)
  306.      2  +192*SIN(U+D+D)-165*SIN(V-D-D)+148*SIN(U-V)-125*SIN(D)
  307.      3  -110*SIN(U+V)-55*SIN(F+F-D-D)-45*SIN(U+F+F)+40*SIN(U-F-F)
  308.      4  -38*SIN(U-4*D)+36*SIN(3*U)-31*SIN(U+U-4*D)+28*SIN(U-V-D-D)
  309.      5  -24*SIN(V+D+D)+19*Sig(U-D)+18*SIN(V+D)+15*SIN(U-V+D+D)
  310.      6  +14*SIN(4*D)+14*SIN(U+U+D+D)-13*SIN(3*U-D-D)
  311.         P(J+19,10)=ARG(XL-IDINT(XL)+DL/1296000.+XNUT,1.)
  312.         P(J+19,11)=ARG(AN/6.28318531,1.)
  313.         IF(J-1)270,260,270
  314. 260     P(1,10)=(18461*SIN(F)+1010*SIN(U+F)-1000*SIN(F-U)
  315.      1  -624*SIN(F-2*D)-167*SIN(U+F-D-D)+199*SIN(F-U+D+D)
  316.      2  +117*SIN(F+D+D)+62*SIN(U+U+F)-33*SIN(F-U-D-D)
  317.      3  -32*SIN(F-U-U)-30*SIN(V+F-D-D)-16*SIN(U+U+F-D-D)
  318.      4  +15*SIN(U+F+D+D)+12*SIN(F-V-D-D)+9*SIN(F-U-V+D+D))/3600.
  319.         P(2,10)=ASIN(COS(P(1,10)*CONVS)*SIN(P(20,10)*6.28318531)
  320.      1  *SIN(E)+SIN(P(1,10)*CONVS)*COS(E))/CONVS
  321. 270     CONTINUE
  322. CCCCCCCCC
  323. C
  324. CCCCCCCCC
  325.         E=(23.4522944-.0130125*T+.002558*DCOS((259.18-1934.142
  326.      1  *T)*CONV))*CONVS
  327.         DO 277 I=2,9
  328.         RETRO(I-1)=' DIR '
  329.         IF(ARG(P(21,I)-P(20,I),1.)-.5)277,277,275
  330. 275     RETRO(I-1)='RETRO'
  331. 277     CONTINUE
  332. CCCCCCCCC
  333. C
  334. CCCCCCCCC
  335.         XLW=(LONGD+LONGM/60.)/360.
  336.         XLN=(LATD+LATM/60.)*CONVS
  337.         IF(LATS.EQ.'N')GO TO 280
  338.         XLN=-XLN
  339. 280     IF(LONGS.EQ.'W')GO TO 290
  340.         XLW=-XLW
  341. CCCCCCCCC
  342. C
  343. C     H(1)=GMT  H(2)=GST  H(3)=LMT  H(4)=LST
  344. C
  345. CCCCCCCCC
  346. 290     H(1)=ARG(SNGL(GMT)+.5,1.)
  347.         W=100.00213590D0*T-IDINT(100.D0*T)
  348.         H(2)=ARG(H(1)+.27691940+W,1.)
  349.         H(3)=ARG(H(1)-XLW,1.)
  350.         H(4)=ARG(H(2)-XLW,1.)
  351.         DO 300 I=1,4
  352.         W=24.*H(I)+.00013888888
  353.         IH(I,1)=W
  354.         W=60.*(W-IH(I,1))
  355.         IH(I,2)=W
  356. 300     IH(I,3)=60.*(W-IH(I,2))
  357.         WRITE(IOUT,12000),((IH(I,J),J=1,3),I=1,4)
  358. 12000   FORMAT(/' GMT=',I2,2(':',I2),'    GST=',I2,2(':',I2),
  359.      1  '    LMT=',I2,2(':',I2),'    LST=',I2,2(':',I2))
  360.         JD1=XJD
  361.         XJD1=XJD-JD1
  362.         I=JD1+INT(XJD1+1.5-XLW)-(JD1-1+INT(XJD1+1.5-XLW))/7*7
  363.         YBES=1900.D0+(XJD-2415020.31351528D0)/365.24219878125D0
  364.      1  -.00000107523D0*T*T
  365.         IYBES=YBES
  366.         YBES1=YBES-IYBES
  367.         WRITE(IOUT,13000)WEEK(I),WEEKR(I),JD1,XJD1,IYBES,YBES1
  368. 13000   FORMAT(' Day of the week (from LMT) IS ',A9,
  369.      1  ' (ruled by ',A8/' JD=',I7,F7.6,7X,' Besselian year = ',
  370.      2  I5,F10.9)
  371. CCCCCCCCC
  372. C
  373. C     HOUSE CUSP COMPUTATIONS (ACCURATE TO LAST DIGIT)
  374. C
  375. CCCCCCCCC
  376. 310   ITER=1+LATD/45
  377. CCCCCCCCC
  378. C
  379. C     "ITER" IS THE NUMBER OF PLACIDEAN ITERATIONS
  380. C@@@
  381.