program orbtrans5 c c This program "translates" the WDS-format input file c to the format used in the orbit and plotting programs. c Weights are set equal to 1.0 - they'll be determined c in second program. This one adds the decimal points, c etc. so the second program can read the data properly. c [ orbgrades.2* --> orbgrades.3* ] c character*9 date character*6 theta character*7 rho character*3 dd1,type,tel character*1 c75,c76,c77 character*20 hcode character*100 title character*10 wds c c enter file names for i/o c open(10,file='orbits.2',status='UNKNOWN') open(11,file='orbits.3',status='UNKNOWN') c c read and write wds title and orbital elements lines c 200 read(10,907,end=999,err=999) title 907 format(a100) write(11,907) title read(10,907) title write(11,907) title c c read wds data, determine method, and write in orbit file format c 400 read(10,908,end=999,err=999) wds,date,theta,rho,hcode, $ dd1,tel,c75,c76,c77,idate 908 format(a10,t24,a9,1x,a6,1x,a7,t61,a20,t64,a3,t72,a3,3a1,t24,i4) if (wds .eq. ' ') go to 500 if (date .eq. ' . ') go to 400 if (theta .eq. ' . ') go to 400 if (rho .eq. ' . ') go to 400 c type='VIS' c if ((c75 .eq. 'G') .or. (c75 .eq. 'H')) type='VPH' if ((c75 .eq. 'O') .or. (c75 .eq. 'L')) type='VPH' c if (c75 .eq. 'J') type='XVI' c if (c75 .eq. 'K') type='XSP' if ((c75 .eq. 'P') .or. (c75 .eq. 'Q')) type='XSP' if ((c75 .eq. 'S') .or. (c75 .eq. 'U')) type='XSP' if ((dd1 .eq. 'MKT') .or. (dd1 .eq. 'NOI')) type='XSP' c if ((dd1 .eq. 'McA') .or. (dd1 .eq. 'Hrt')) type='CSP' if ((dd1 .eq. 'Msn') .or. (dd1 .eq. 'Fek')) type='CSP' if ((dd1 .eq. 'Trn') .or. (dd1 .eq. 'TtB')) type='CSP' c if (dd1 .eq. 'WSI') type='WSI' if ((dd1 .eq. 'Msn') .and. (tel .eq. ' 26')) type='WSI' if ((dd1 .eq. 'Msn') .and. (tel .eq. ' 61')) type='WSI' if ((dd1 .eq. 'Msn') .and. (tel .eq. ' 82')) type='WSI' if ((dd1 .eq. 'Msn') .and. (tel .eq. ' 83')) type='WSI' if ((dd1 .eq. 'Msn') .and. (tel .eq. ' 84')) type='WSI' if ((dd1 .eq. 'Msn') .and. (tel .eq. '150') $ .and. (idate .ge. 1999)) type='WSI' if ((dd1 .eq. 'Msn') .and. (tel .eq. '157') $ .and. (idate .ge. 1999)) type='WSI' if ((dd1 .eq. 'Msn') .and. (tel .eq. '158') $ .and. (idate .ge. 1999)) type='WSI' c if (dd1 .eq. 'HIP') type='HIP' if (dd1 .eq. 'TYC') type='TYC' c weight=1.0 c if(theta .eq. ' . ') weight=0.0 if(rho .eq. ' . ') weight=0.0 write(11,909) date,theta,rho,type,weight,hcode 909 format(6x,a9,4x,a6,7x,a7,10x,a3,f5.1,5x,a20) go to 400 c 500 write(11,910) 910 format(80x,' ') go to 200 999 stop 'data formatting program finished' end