Program WTPROC
c program to process Wind Tunnel Data Files into TECPLOT format.
      real D,CWF,L,RM,PM,YM,AI,ASS,DP
      real CD,CCWF,CL,CRM,CPM,CYM,V
      character*40 infil,outfil
      character*70 title
      character*80 buffer
      character*1  infilc(40),outfilc(40),bufferc(80),titlec(70)
      equivalence  (infil,infilc),(outfil,outfilc),(buffer,bufferc)
      equivalence  (title,titlec)
c
      c1=7.67345e-4
      c0=0.059862
      Sref=0.09
      Swing=1.0
      write (*,*) 'Wind Tunnel Data Processing'
      write (*,*) '---------------------------'
13    write (*,*) ' '
      write (*,*) ' Give Name of INPUT wind tunnel data file :'
      read (*,1) infil
1     format(A40)
      if (infilc(1).eq.' ') goto 999
      open(1,file=infil,status='old',err=11)
      goto 12
11    write (*,*) 'Cannot Find File. Try Again.'
      goto 13
12    write (*,*) ' '
      write (*,*) ' Give Name of OUTPUT Tecplot file : '
      read (*,1) outfil
      if (outfilc(1).eq.' ') then
       write (*,*) ' Using default OUTPUT file : OUT.TEC'
       outfil='out.tec'
      endif
      open (2,file=outfil,status='UNKNOWN')
      read (1,1001) title
1001  format(A70)
99    read (1,1002,err=88,end=88) buffer      
1002  format(A80)
      if (bufferc(22).eq.'A') then
       do 201 i=1,20
        bufferc(i)=bufferc(i+29)
201    continue
       read (buffer,*) Swing
      endif 
      if (bufferc(1).ne.'R') goto 99
      read (1,1002) buffer
      read (1,1002) buffer
      bufferc(1)='T'
      bufferc(2)='I'
      bufferc(3)='T'
      bufferc(4)='L'
      bufferc(5)='E'
      bufferc(6)=' '
      bufferc(7)='='
      bufferc(8)=' '
      bufferc(9)='"'
      do 21 i=1,65
       j=71-i
       if (titlec(j).ne.' ') goto 22
21    continue
22    do 23 i=1,j
       bufferc(i+9)=titlec(i)
23    continue
      j=j+10
      bufferc(j)='"'
      write (2,2001) (bufferc(i),i=1,j)
2001  format(80a1)
      write (2,*) 'VARIABLES = AI,CL,CD,CPM,CL2'
      write (2,*) 'ZONE F=POINT'
      cl0=-1000.0
199   read (1,1002,end=188,err=188) buffer
      read (1,1003,end=188,err=188) D,CWF,L,RM,PM,YM,AI,ASS,DP
1003  format(9F8.2)
      read (1,1004,end=188,err=188) CD,CCWF,CL,CRM,CPM,CYM,V
1004  format(6F8.4,16X,F9.4)
      CL2=CL*CL
      CDS=C0-C1*AI
      CD=CD-CDS*Sref/Swing
      write(2,1005) AI,CL,CD,CPM,CL2
1005  format(5F10.4)
      goto 199
188   write (*,*) 'Processing Complete.'
      close (1)
      close (2)
      goto 999
88    write (*,*) 'No Wind Tunnel Data found.'          
999   write (*,*) '.....Exiting'
      stop
      end