; FIXMRC : cut and paste the spider-generated power spectum onto the
;          model file generated by ctffind

x88 = 1   ; micrograph number

FR G
<mrc>pow{***x88}   ; input power spectrum file generated by ctffind2
                 ; (already in SPIDER format)
FR G
<spi>power/pw_avg{***x88}   ; input power spectrum file generated by spider

FR G
<out>fix{***x88}    ; output SPIDER file

; ------------------------------------------------

FI x11,x12   ; get image size
<spi>
(12,2)

; this part estimates a background and subtracts it from the power spectrum
; for better visualization

x13 = (x11/2)

FQ
<spi>
_2   
(7)              ; filter type (BW LP)
(0.002,0.025)        ; passband, stopband

SU
<spi>
_2
_3
*

MA
_3          ; input - the filtered file
_4          ; output - filt file center masked
(0.0,35)    ; radius of the central dot
D
E
(0.0)
(x13+1,x13+1)

; resize and recale the background -subtracted spider image

FI x11,x12   ; get image size
<mrc>
(12,2)

x21 = (x11/2) + 1   ; image center
x22 = (x12/2) + 1


IP
_4
_1
x11,x12

WI        ; extract left half of mrc image
<mrc>
_5
(x11/2, x12)  ; dimensions
1,1       ; upper left coordinates

FS x31,x32      ; get max & min values
_5

WI        ; extract right half of spider image
_1
_2
(x11/2, x12)  ; dimensions
x21,1       ; upper left coordinates

x40 = x31 - x32  ; range of values
x41 = x31 + (x40/2) ; rescaled max
x42 = x32 - (x40/2) ; rescaled min

AR SCA       ; rescale the spider window
_2
_3
x42,x41

CP
<mrc>
<out>

IN
_3
<out>
x21,1

EN D