24 jones 1.2
25 write(*,*) 'Optimize Matrix element coefficient filename'
26 read(*,*) infile
27 open(unit=lun,status='old',name=infile,err=92)
|
28 jones 1.1
29 10 continue
30 read (lun, '(a)', end=1000) line ! read in a line
31 l_num = l_num + 1 ! present line #
32
33 call me_db_read_parse(line, l_num)! parse the line
34
35 go to 10 ! loop over the file
36
37 1000 continue
38 write(6,*)'opening file'
39 open(unit=chan,status='unknown',name='hms_recon_coeff.new',err=93)
40 c write(chan,*)'! focal plane rotation coeffs(4)'
41 c write(chan,*)'! detector offsets'
42 c write(chan,*)'! Z position of true focus'
43 c write(chan,*)'! recostruction matrix elements'
44 c write(chan,*)'! <phi y theta delta | nmpq>;(x**n xp**m y**p yp**q)'
45 c write(chan,*)'h_ang_slope_x = 0.0'
46 c write(chan,*)'h_ang_slope_y = 0.0'
47 c write(chan,*)'h_ang_offset_x = 0.0'
48 c write(chan,*)'h_ang_offset_y = 0.0'
49 jones 1.1 c write(chan,*)'h_det_offset_x = 0.0'
50 c write(chan,*)'h_det_offset_y = 0.0'
51 c write(chan,*)'h_z_true_focus = 0.0'
52 write(chan,*)'---------------------------------------------------------
53 >------------'
54
55 call printme(ofp_targ,lines)
56
57 do k=1,lines
58 write(chan,1200)(coeff(i,k),i=1,4),(expon(j,k),j=1,4)
59 enddo
60
61 write(chan,*)"--------------------------------------------------------
62 >-------------"
63
64 goto 100
|
66 jones 1.1 goto 100
67 93 write(*,*)'could not open hms_recon_coeff.new'
68 goto 100
69
70 1200 format(1x,4(G16.9),1x,4i1)
71
72 100 close(unit=lun)
73 close(unit=chan)
|