set vari FILE_ATOM = "amy.pdb"
set vari FILE_CTAB = 

set vari FILE_SYMM = ">symm/p21_21_21.symm"
set vari FILE_FOBS = "amylase_mir.fobs"
set vari FILE_CELL = "cell.dat"

set vari RESOL_MIN global = 8.0
set vari RESOL_MAX global = 3.0
set vari AGRID int = 90
set vari BGRID int = 60
set vari CGRID int = 40
set vari H_LIMIT int = 0
set vari K_LIMIT int = 0
set vari L_LIMIT int = 0

set vari MAP_2FOFC global = 1
set vari MAP_FOFC global = 2
set vari MAP_WORK global = 2
set vari MAP_FO global = 1

<>utils/get_top_par_19_csd
<>cmds/load_depp_page

!read file FILE_ATOM coor xpl
! set vari WORK_SEGM char global segm sele all end 

set vari WORK_SEGM global char = "MO*"

read file FILE_CTAB ctable
show bonds
if ( IRESULT_0 .le. 0 ) <>utils/calculate.bonds WORK_SEGM

!<xplor_2_main.com

!<DEF_ALL WORK_SEGM

define init group 
set charge sele .not atom name C CA N H O .a \
    resi name ASP GLU LYS ARG end 0.0

read file FILE_CELL cell
read file FILE_SYMM symm

read file FILE_FOBS refl init limit H_LIMIT K_LIMIT L_LIMIT friedel \
     reso RESOL_MIN RESOL_MAX


make map 1 from 0 init -9999 grid AGRID BGRID CGRID cell real
make map 2 from 1 init -9999 
copy map 2 3

ener all on dens on dens map MAP_2FOFC dens scal 25.

!<?MAIN_CMDS:re_phase.cmds WORK_SEGM MAP_FOFC
!<?MAIN_CMDS:calc_2fofc_map.cmds MAP_2FOFC

ima init cent calc
<model.view
!<>utils/rama
<>utils/show.cell

ima sele atom name CA C N end col 1 bond

ima from erase
 ima col 245 map 1 dens 0.8  12 12 12

<?MAIN_CMDS:re_image.cmds WORK_SEGM

return