; b04.rec
;create four fold symmetric volume
;input parameters:
;1) input file (protein positive)
;2) symmetry factor
;3) output file
x12=1.0
if (x12.eq.1)goto lb5
@rav3dal
proj3/rec/3danticlockfqinv001
(4)
proj3/rec/3danticlockfqinvsym001
;create a 3D mask:
; input parameters:
;1) input volume (protein positive, good if already low-pass filtered)
;2) first threshold in sigma units
;3) spherical mask radius
;4) number of slices to be cut from top and bottom of volume
;5) low-pass filter
;6) next threshold
;7) final threshold
;8) mask file
;9) control file (3D multiplied with mask)
lb5
@cmask3d3
proj3/rec/3danticlockwi001
(.7)
(39)
(15)
(.08)
(.5)
(.1)
proj3/mask3danticlock001
proj3/contro3danticlock001
; low pass filter mask for use in resolution test:
fq np
proj3/mask3danticlock001
proj3/mask3danticlockfq001
(5) 
(0.03) 
(0.02) 
en


; padd to larger size if final reconstruction is larger
; than the intermediate ones
;pd
;mask3dleft001
;mask3dleftpad001
;(100,100)
;(100)
;n
;(0)
;(11,11)
;(11)
; low pass filter padded mask (low pass filter should be done
; after padding to prevent an edge.
fq np
mask3dleft001
mask3dleftfq001
(5) 
(0.03) 
(0.02) 
en