<get_data_hex
read file FILE_CELL cell
read file FILE_SYMM symm
read file FILE_REFL reflect init resol RESOL_HEX RESOL_LOW limit 0 0 0 friedel
refl key WORK_REFL sele defined end
set vari WORK_FORM global = "PR_H EN_H"
set vari MAP_FOFC global = MAP_HEX + 1
set vari MAP_2FOFC global = MAP_HEX 
make map MAP_HEX conv real

ima from 0
image erase

<>utils/show.cell
image sele segm name PR_H EN_H .a atom name C CA N end col 1 bond
<show_other
image from erase
<re_image.cmds
image cent calc sele segm name WORK_FORM end
return