!make pair sele ca .a numb segm 5 end sele ca .a numb segm 4 .a seq A* end
!make pair sele ca .a numb segm 6 end sele ca .a numb segm 4 .a seq B* end
<>utils/get_top_par_19_csd
set vari MAP_HEX global = 1
set vari MAP_ORT global = 3
set vari RESOL_LOW global = 10.
set vari RESOL_HEX global = 3.3
set vari RESOL_ORT global = 3.2

set vari FILE_ATOM = SAVE_FILE.PDB

read file FILE_ATOM atom pdb
read file SAVE_FILE.CTAB ctab
key MOL_HEX sele segm name PR_H EN_H end
key MOL_ORT sele segm name PR_O EN_O  end
set vari WORK_SEGM global = "PR* EN*"

<get_data_hex
read file FILE_CELL cell
read file FILE_SYMM symm

make map MAP_HEX from 0 init -9999 grid  AGRID BGRID CGRID cell real
make map MAP_HEX + 1 from 0 init -9999 grid  AGRID BGRID CGRID cell real


<get_data_ortho
read file FILE_CELL cell
read file FILE_SYMM symm

make map MAP_ORT from 0 init -9999 grid  AGRID BGRID CGRID cell real
make map MAP_ORT + 1 from 0 init -9999 grid  AGRID BGRID CGRID cell real


! <>utils/protein.bond

 <DEF_ALL  "*"

<load_switch
image init
<work_ort
<re_image.cmds
return