-
Notifications
You must be signed in to change notification settings - Fork 1
/
make_frealign_param.spi
107 lines (73 loc) · 1.34 KB
/
make_frealign_param.spi
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
; b20_get_sel_angles_sel_shifts_final.spi
;
;
; This batch gets sel_angles and sel_shifts_final and puts them in 1 doc file in
; preparation for frealign.
;
;;;;;;;;;;;;;
FR
?Input angular file? [angular]
FR
?Input shifts file? [shifts]
FR
?Input CTF-parameter file? [ctf]
RR [mag]
?Input actual magnification?
;;;;;;;;;;;;;
UD N x12
[shifts]
;x12=12921;#of particles
;x13=2 ;round of refinement # (directory #)
x19=360
x20=180
x21=-1
x70=0.0
; magnification
x71=[mag]
DO LB1 x11=1,x12
;
; Get angular file
;
UD IC,x11,x14,x15,x16
[angular]
;x14=x19-x14 ;make psi negative
;x16=x19-x16 ;make phi negative
;
; Get shifts
;
UD IC,x11,x17,x18
[shifts]
;x17=2*x17 ;multiply x shift by 2 to account for decimation factor
;x18=2*x18 ;multiply y shift by 2 to account for decimation factor
if(x14.LT.0.0)x14=x14+360.
if(x14.LT.0.0)x14=x14+360.
if(x14.GT.360.0)x14=x14-360.
if(x14.GT.360.0)x14=x14-360.
SD x11,x14,x15,x16,x17,x18
angles_and_shifts
;
; get ctf parameters, add mag. and "phase residuals"
;
UD IC,x11,x21,x22,x23,x24,x25
[ctf]
sd x71,x21,x22,x23,x24,x70,x70
ctf_param
LB1
SD E
angles_and_shifts
UD ICE
[angular]
UD ICE
[shifts]
UD ICE
[ctf]
sd e
ctf_param
;
; Format the parameters for frealign.
;
VM
fre_format_param.sh angles_and_shifts.spi ctf_param.spi parameters_in_01.spi
VM
rm -f angles_and_shifts.spi ctf_param.spi
EN D