!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

<get_data_hex
read file FILE_CELL cell
read file FILE_SYMM symm
read file FILE_ATOM atom pdb

make map MAP_HEX from 0 init -9999 auto_gridcell real
make map MAP_HEX + 1 from MAP_HEX init -9999 cell real

key MOL_HEX sele all end

<get_data_ortho
read file FILE_CELL cell
read file FILE_SYMM symm
read file FILE_ATOM atom pdb

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

key MOL_ORT sele .not MOL_HEX end

! <>utils/protein.bond
read file model_both.ctab ctab
set vari WORK_SEGM char global segm sele all end 

 <DEF_ALL WORK_SEGM

<load_switch
image init
<work_ortho
<re_image.cmds 
return