program noorbtrans 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 [ orbits.2 --> orbits.3 ] c c 110708: changed some method codes (EYE, LBI, CCD, PTG) c 120709: updated for new WDS format c character*10 date,wds character*9 rho character*7 theta character*1 tcode,rcode,tech1 character*3 type character*2 tech2 character*35 hcode character*100 title c c enter file names for i/o c open(10,file='noorb.2',status='UNKNOWN') open(11,file='noorb.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,tcode,theta,rcode,rho, $ hcode,tech1,tech2 908 format(a10,t18,a10,t29,a1,a7,t45,a1,a9,t93,a35,t122,a1,t122,a2) if (wds .eq. ' ') go to 500 if (date .eq. ' . ') go to 400 if (tcode .eq. 'V') go to 400 if (rcode .eq. '<') go to 400 if (rcode .eq. '>') go to 400 if (rcode .eq. 'U') go to 400 if (rcode .eq. 'E') go to 400 if (theta .eq. ' . ') go to 400 if (rho .eq. ' . ') go to 400 c type='VIS' c if (tech1 .eq. 'A') type='SPE' if (tech1 .eq. 'C') type='CCD' if (tech1 .eq. 'E') type='CCD' if (tech1 .eq. 'H') type='HIP' if (tech1 .eq. 'I') type='SPE' if (tech1 .eq. 'J') type='EYE' if (tech1 .eq. 'K') type='LBI' if (tech1 .eq. 'O') type='OCC' if (tech1 .eq. 'P') type='PTG' if (tech1 .eq. 'S') type='SPE' c weight=1.0 c write(11,909) date,tcode,theta,rcode,rho,type,weight,hcode 909 format(6x,a10,t19,a1,a7,t32,a1,a9,t50,a3,f5.1,5x,a35) go to 400 c 500 write(11,910) 910 format(95x,' ') go to 200 999 stop 'data formatting program finished' end