====resolve_map.com=========
#!/bin/csh -v
#
fft:
# Calculate maps from density modified map.

fft hklin resolve MAPOUT resolve.map << +
TITLE Map from solvent flattening and Hist. Matching
!AXIS  X Y Z
!limit 0 0.5 0 0.5 0 0.25
reso 20.0 2.5
LABIN F1=FP PHI=PHIM W=FOMM
END
+

mapmask mapin resolve.map mapout resolve.map_X << END
AXIS   X Y Z
XYZLIM 0.0 1.0   0.0 1.0  0.0 1.0
END

mapman -b << EOF
re m1 resolve.map_X CCP4
li m1
no m1
bones skel m1 1.5 1 10
BOnes CO resolve.bones skel 5
BOnes PRu bones.pdb 5 5.
br m1 resolve.brix
quit
EOF
====resolve_map.com========