program orb6orbits c c This program reads data from orb6.master and c generates html files containing the "keeper" orbits. c character*1 junk,keep,note,ref2 character*20 rad,dd,ads,hd,hip,v1,v2,p,pe,a,ae,ii,iie,nn, $ nne,t,te,e,ee,o,oe,eqnx,last character*7 ref1 character*4 ads1 character*10 wds character*18 png character*168 header character*232 data character*264 header2 c open(10,file='../Catalog/orb6.master', $ status='UNKNOWN') open(11,file='orb6orbits.header',status='UNKNOWN') open(12,file='orb6orbits.html',status='UNKNOWN') open(13,file='orb6orbits.body',status='UNKNOWN') open(14,file='orb6orbits.txtheader',status='UNKNOWN') open(15,file='orb6orbits.txt',status='UNKNOWN') c c read header lines from orbit catalog c do 100 n=1,4 100 read(10,901) junk 901 format(a1) c do 120 n=1,7 read(14,922) header2 922 format(a264) write(15,922) header2 120 continue c norbits=0 c c write header lines for orb6orbits.body c write(13,902) 902 format(''/''/''/ $ ''/'

'/'

')
c
c	read and write header lines for html file
c
  200	read(11,903,end=300) header
  903	format(a168)
  	write(12,903) header
	go to 200
c
c	read star name and orbital elements. Loop back for another 
c	if orbit isn't a "keeper"
c
  300	read(10,904,end=700) rad,wds,dd,ads,hd,hip,v1,v2,p,pe,a,ae,
     $     ii,iie,nn,nne,t,te,e,ee,o,oe,eqnx,last,grade,note,keep,
     $     ref1,ref2,png,data
  904	format(a18,1x,a10,1x,a15,a6,2a7,2a7,1x,a13,a12,a11,a9,
     $     2a9,a10,a9,a13,a12,2a9,2a9,2a5,10x,f3.1,1x,a1,1x,a1,1x,
     $     a7,a1,1x,a18,t1,a232)
	if (wds .eq. '          ') go to 300
	if (keep .eq. 'Y') keep='y'
	if (note .eq. 'N') note='n'
	if (keep .ne. 'y') go to 300
	igrade=nint(grade)
	if (grade .lt. 1.49) igrade=1
	if (grade .gt. 4.51) igrade=5
	if (grade .gt. 5.51) igrade=6
	if (grade .gt. 6.51) igrade=7
	if (grade .gt. 7.51) igrade=8
	if (grade .gt. 8.51) igrade=9
c
c	write lines for orb6orbits.html
c
        ads1='ADS '
        if (ads .eq. '.     ') ads='      '
        if (ads .eq. '      ') ads1='    '
c
	if ((note .ne. 'n') .and. (igrade .ne. 6) .and. (igrade .ne. 7))
     $     write(12,905) wds,rad,wds,dd,hd,v1,p,a,ii,
     $     nn,t,e,o,eqnx,igrade,note,png,wds,ref1,ref2,ref1,ref2
  905   format('',a18,1x,a10,1x,a15,a7,1x,a7,1x,
     $     a13,2x,a11,1x,a9,1x,a10,1x,a13,2x,a9,1x,a9,1x,a5,i2,1x,a1,
     $     ' P',
     $     ' E',
     $     '  ',a7,a1,'')
c
	if ((note .ne. 'n') .and. (igrade .eq. 6))
     $     write(12,915) wds,rad,wds,dd,hd,v1,p,a,ii,
     $     nn,t,e,o,eqnx,igrade,note,wds,ref1,ref2,ref1,ref2
  915   format('',a18,1x,a10,1x,a15,a7,1x,a7,1x,
     $     a13,2x,a11,1x,a9,1x,a10,1x,a13,2x,a9,1x,a9,1x,a5,i2,1x,a1,2x,
     $     ' E',
     $     '  ',a7,a1,'')
c
	if ((note .ne. 'n') .and. (igrade .eq. 7))
     $     write(12,915) wds,rad,wds,dd,hd,v1,p,a,ii,
     $     nn,t,e,o,eqnx,igrade,note,wds,ref1,ref2,ref1,ref2
c
	if ((note .eq. 'n') .and. (igrade .ne. 6) .and. (igrade .ne. 7)) 
     $     write(12,906) wds,rad,wds,dd,hd,v1,p,a,ii,
     $     nn,t,e,o,eqnx,igrade,wds,png,wds,ref1,ref2,ref1,ref2
  906   format('',a18,1x,a10,1x,a15,a7,1x,a7,1x,
     $     a13,2x,a11,1x,a9,1x,a10,1x,a13,2x,a9,1x,a9,1x,a5,i2,
     $     ' N',
     $     ' P',
     $     ' E',
     $     '  ',a7,a1,'')
c
	if ((note .eq. 'n') .and. (igrade .eq. 6)) 
     $     write(12,916) wds,rad,wds,dd,hd,v1,p,a,ii,
     $     nn,t,e,o,eqnx,igrade,wds,wds,ref1,ref2,ref1,ref2
  916   format('',a18,1x,a10,1x,a15,a7,1x,a7,1x,
     $     a13,2x,a11,1x,a9,1x,a10,1x,a13,2x,a9,1x,a9,1x,a5,i2,
     $     ' N',2x,
     $     ' E',
     $     '  ',a7,a1,'')
c
	if ((note .eq. 'n') .and. (igrade .eq. 7)) 
     $     write(12,916) wds,rad,wds,dd,hd,v1,p,a,ii,
     $     nn,t,e,o,eqnx,igrade,wds,wds,ref1,ref2,ref1,ref2
c
        write(12,907) ads1,ads,hip,v2,pe,ae,iie,nne,te,ee,oe,last
  907   format(30x,a4,a6,5x,a7,1x,a7,2x,a12,3x,a9,2x,a9,1x,a9,3x,
     $     a12,2x,a9,1x,a9,1x,a5/)
c
c	write lines for orb6orbits.body
c
	if ((note .ne. 'n') .and. (ref2 .ne. ' ') .and. 
     $     (igrade .ne. 6) .and. (igrade .ne. 7))
     $     write(13,908) rad,wds,dd,hd,v1,p,a,ii,nn,t,e,o,eqnx,
     $     igrade,note,png,wds,ref1,ref2,ref1,ref2
  908	format(a18,1x,a10,1x,a15,a7,1x,a7,1x,
     $     a13,2x,a11,1x,a9,1x,a10,1x,a13,2x,a9,1x,a9,1x,a5,i2,1x,a1,
     $     ' P',
     $     ' E ',
     $     ' ',a7,a1,'')
c
	if ((note .ne. 'n') .and. (ref2 .ne. ' ') .and. (igrade .eq. 6))
     $     write(13,918) rad,wds,dd,hd,v1,p,a,ii,nn,t,e,o,eqnx,
     $     igrade,note,wds,ref1,ref2,ref1,ref2
  918	format(a18,1x,a10,1x,a15,a7,1x,a7,1x,
     $     a13,2x,a11,1x,a9,1x,a10,1x,a13,2x,a9,1x,a9,1x,a5,i2,1x,a1,2x,
     $     ' E ',
     $     ' ',a7,a1,'')
c
	if ((note .ne. 'n') .and. (ref2 .ne. ' ') .and. (igrade .eq. 7))
     $     write(13,918) rad,wds,dd,hd,v1,p,a,ii,nn,t,e,o,eqnx,
     $     igrade,note,wds,ref1,ref2,ref1,ref2
c
  	if ((note .eq. 'n') .and. (ref2 .ne. ' ') .and. 
     $     (igrade .ne. 6) .and. (igrade .ne. 7))
     $     write(13,909) rad,wds,dd,hd,v1,p,a,ii,nn,t,e,o,eqnx,
     $     igrade,wds,png,wds,ref1,ref2,ref1,ref2
  909	format(a18,1x,a10,1x,a15,a7,1x,a7,1x,a13,2x,a11,1x,a9,1x,a10,
     $     1x,a13,2x,a9,1x,a9,1x,a5,i2,
     $     ' N',
     $     ' P',
     $     ' E ',
     $     ' ',a7,a1,'')
c
  	if ((note .eq. 'n') .and. (ref2 .ne. ' ') .and. (igrade .eq. 6))
     $     write(13,919) rad,wds,dd,hd,v1,p,a,ii,nn,t,e,o,eqnx,
     $     igrade,wds,wds,ref1,ref2,ref1,ref2
  919	format(a18,1x,a10,1x,a15,a7,1x,a7,1x,a13,2x,a11,1x,a9,1x,a10,
     $     1x,a13,2x,a9,1x,a9,1x,a5,i2,
     $     ' N',2x,
     $     ' E ',
     $     ' ',a7,a1,'')
c
  	if ((note .eq. 'n') .and. (ref2 .ne. ' ') .and. (igrade .eq. 7))
     $     write(13,919) rad,wds,dd,hd,v1,p,a,ii,nn,t,e,o,eqnx,
     $     igrade,wds,wds,ref1,ref2,ref1,ref2
c
	if ((note .ne. 'n') .and. (ref2 .eq. ' ') .and. 
     $     (igrade .ne. 6) .and. (igrade .ne. 7))
     $     write(13,910) rad,wds,dd,hd,v1,p,a,ii,nn,t,e,o,eqnx,
     $     igrade,note,png,wds,ref1,ref1
  910	format(a18,1x,a10,1x,a15,a7,1x,a7,1x,
     $     a13,2x,a11,1x,a9,1x,a10,1x,a13,2x,a9,1x,a9,1x,a5,i2,1x,a1,
     $     ' P',
     $     ' E ',
     $     ' ',a7,'')
c
	if ((note .ne. 'n') .and. (ref2 .eq. ' ') .and. (igrade .eq. 6))
     $     write(13,920) rad,wds,dd,hd,v1,p,a,ii,nn,t,e,o,eqnx,
     $     igrade,note,wds,ref1,ref1
  920	format(a18,1x,a10,1x,a15,a7,1x,a7,1x,
     $     a13,2x,a11,1x,a9,1x,a10,1x,a13,2x,a9,1x,a9,1x,a5,i2,1x,a1,2x,
     $     ' E ',
     $     ' ',a7,'')
c
	if ((note .ne. 'n') .and. (ref2 .eq. ' ') .and. (igrade .eq. 7))
     $     write(13,920) rad,wds,dd,hd,v1,p,a,ii,nn,t,e,o,eqnx,
     $     igrade,note,wds,ref1,ref1
c
  	if ((note .eq. 'n') .and. (ref2 .eq. ' ') .and. 
     $     (igrade .ne. 6) .and. (igrade .ne. 7)) 
     $     write(13,911) rad,wds,dd,hd,v1,p,a,ii,nn,t,e,o,eqnx,
     $     igrade,wds,png,wds,ref1,ref1
  911	format(a18,1x,a10,1x,a15,a7,1x,a7,1x,a13,2x,a11,1x,a9,1x,a10,
     $     1x,a13,2x,a9,1x,a9,1x,a5,i2,
     $     ' N',
     $     ' P',
     $     ' E ',
     $     ' ',a7,'')
c
  	if ((note .eq. 'n') .and. (ref2 .eq. ' ') .and. (igrade .eq. 6))
     $     write(13,921) rad,wds,dd,hd,v1,p,a,ii,nn,t,e,o,eqnx,
     $     igrade,wds,wds,ref1,ref1
  921	format(a18,1x,a10,1x,a15,a7,1x,a7,1x,a13,2x,a11,1x,a9,1x,a10,
     $     1x,a13,2x,a9,1x,a9,1x,a5,i2,
     $     ' N',2x,
     $     ' E ',
     $     ' ',a7,'')
c
  	if ((note .eq. 'n') .and. (ref2 .eq. ' ') .and. (igrade .eq. 7))
     $     write(13,921) rad,wds,dd,hd,v1,p,a,ii,nn,t,e,o,eqnx,
     $     igrade,wds,wds,ref1,ref1
c
        write(13,907) ads1,ads,hip,v2,pe,ae,iie,nne,te,ee,oe,last
c
c       write lines for orb6orbits2.txt
c
        write(15,912) data,igrade,note,ref1,ref2,png
  912   format(a232,i2,1x,a1,1x,a7,a1,1x,a18)
c
	go to 300
c
  700	write(12,913)
  913	format('

') write(13,914) 914 format(''/'') c 999 stop end