read file ?MAIN_UTILS:targ_5res.pdb coor pdb
calc bond neigh sele segm name TARG end


 set vari ITEST = 3882 !  4540 ! 4290 ! 4534 ! 2673 ! 3102
set vari ITEST int atom sele atom numb $1 end

make map MAP_TARG + 1 from MAP_TARG init 0 integ
  dele atom sele segm name NEW end
  find sec search 2 sele segm name HELP end lev 40 map MAP_ACT \
      sele segm name TARG end test ITEST contour 0.8
  set vari WORK_SEGM char global = NEW

dele atom sele segm name TARG end


<?MAIN_CMDS:anal_group.cmds RESIDUE 
<?MAIN_CMDS:anal_temp.cmds WORK_SEGM ANAL_GRP 
<?MAIN_CMDS:anal_3d_image.cmds WORK_SEGM ANAL_GRP 

menu block TRACE item TRC_DENS button off
menu block TRACE item TRC_GRAD button off

  return


ima from 5
ima era

 find sec merge sele segm name NEW end map 1 work 4 test ITEST cont MAP_CONT
 ima col 220 sele segm name NEW .a weigh 0.1 1.1 end bond stick
 ima col 230 sele segm name NEW .a weigh 1.1 2.1 end bond stick
 !dele atom sele segm name NEW .a weigh -0.1 0.1 end
 set vari WORK_SEGM char global = NEW
 ima from era

!  make bond sele segm name NEW end from topo
!  key 5 sele seq %%%5 .a segm name NEW end
!  key 6 sele seq %%%6 .a segm name NEW end
!  dele bond sele 5 end sele 6 end

<?MAIN_CMDS:anal_group.cmds RESIDUE 
<?MAIN_CMDS:anal_temp.cmds WORK_SEGM ANAL_GRP 
<?MAIN_CMDS:anal_3d_image.cmds WORK_SEGM ANAL_GRP 
!ima rem 12 24
ima cent atom ITEST
key image sele segm name WORK_SEGM .a atom name N CA C end
image sele neigh ITEST -1 end col 110 bond stick

return