<get_data_ortho
read file FILE_CELL cell
read file FILE_SYMM symm
read file FILE_REFL reflect init resol RESOL_ORT RESOL_LOW limit 0 0 0 friedel
refl key WORK_REFL sele defined end

set vari WORK_FORM global = "PR_O EN_O"

set vari MAP_FOFC global = MAP_ORT + 1
set vari MAP_2FOFC global = MAP_ORT
make map MAP_ORT conv real

ima from 0
image erase

<>utils/show.cell

!image sele hucb .a atom name C CA N end col 160 bond
!image sele hucb end col 110 bond
!<show_hucb
image sele atom name C CA N .a segm name PR_O end col 1 bond
image sele atom name C CA N .a segm name EN_O end col 160 bond
<show_other
image from erase
<re_image.cmds
image cent calc sele segm name WORK_FORM end
return