program orbtrans6 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 c 110708: changed some method codes (EYE, LBI, CCD, PTG) 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='PTG' if ((c75 .eq. 'O') .or. (c75 .eq. 'L')) type='PTG' c if (c75 .eq. 'J') type='EYE' c if ((c75 .eq. 'E') .or. (c75 .eq. 'F')) type='CCD' c if (c75 .eq. 'K') type='LBI' c if ((c75 .eq. 'P') .or. (c75 .eq. 'Q')) type='SPE' if ((c75 .eq. 'S') .or. (c75 .eq. 'U')) type='SPE' 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