1 jones 1.1         program conback
 2           
 3           
 4                   implicit none
 5           

 6 jones 1.3         include 'db.inc'

 7 jones 1.1         include 'conv.inc'
 8           
 9                   character*256 line

10 jones 1.2         character*120 infile

11 jones 1.4         character*120 outfile

12 jones 1.1         integer  l_num
13           
14            
15           	record /ofp_targ/ var(TARG_VAR)	! data to use 
16                   integer i,j,k,lines				! i/o unit to use
17           	character*8    type			! type of ME to operate on
18           
19           
20                         
21                   type = 'set'        
22                   l_num=0
23                   lines=0
24           

25 jones 1.2         
26                   write(*,*) 'Optimize Matrix element coefficient filename'
27                   read(*,*) infile

28 jones 1.4         write(*,*) 'Output filename'
29                   read(*,*) outfile

30 jones 1.2         open(unit=lun,status='old',name=infile,err=92)

31 jones 1.1 
32           10	continue
33           	read (lun, '(a)', end=1000) line	! read in a line
34           	l_num = l_num + 1		! present line #
35           
36           	call me_db_read_parse(line, l_num)! parse the line
37           
38           	go to 10				! loop over the file
39           
40           1000	continue
41                   write(6,*)'opening file'

42 jones 1.4         open(unit=chan,status='unknown',name=outfile,err=93)

43 jones 1.1 c	write(chan,*)'! focal plane rotation coeffs(4)'
44           c	write(chan,*)'! detector offsets'
45           c	write(chan,*)'! Z position of true focus'
46           c	write(chan,*)'! recostruction matrix elements'
47           c	write(chan,*)'! <phi y theta delta | nmpq>;(x**n xp**m y**p yp**q)'
48           c	write(chan,*)'h_ang_slope_x  = 0.0'
49           c	write(chan,*)'h_ang_slope_y  = 0.0'
50           c	write(chan,*)'h_ang_offset_x = 0.0'
51           c	write(chan,*)'h_ang_offset_y = 0.0'
52           c	write(chan,*)'h_det_offset_x = 0.0'
53           c	write(chan,*)'h_det_offset_y = 0.0'
54           c	write(chan,*)'h_z_true_focus = 0.0'
55                   write(chan,*)'---------------------------------------------------------
56                >------------'
57           
58                   call printme(ofp_targ,lines)
59           
60                   do k=1,lines

61 jones 1.4          write(chan,1200)(coeff(i,k),i=1,4),(expon(j,k),j=1,5)

62 jones 1.1         enddo
63           
64                   write(chan,*)"--------------------------------------------------------
65                >-------------"
66           
67                   goto 100        

68 jones 1.2  92     write(*,*)'could not open ',infile

69 jones 1.1         goto 100

70 jones 1.4  93     write(*,*)'could not open ',outfile

71 jones 1.1         goto 100
72           

73 jones 1.4  1200   format(1x,4(G16.9),1x,5i1)

74 jones 1.1 
75            100    close(unit=lun)
76                   close(unit=chan)

77 jones 1.4         write(*,*) ' new elements in ',outfile

78 jones 1.1         end  
79           
80           
81           
82