subroutine int MAP_OF_ORI
! the map MAP_ENVEL will be created for the map presentations (as the last map)
! CUSHION is the number of grids around the box of your "active" atoms
! CUT_OFF tells how low can a density be still 
!	considered to belong to the model
! MAP_OF_ORIG tells from which map you want to extract the density 
!	The variable "map_to" is set with the MAP_ACT? item 


set vari MAP_ENVEL global int ?= nmaps + 1
set vari CUSHION = 5
set vari CUT_OFF = 0.0

! don't touch below if you are not sure what you are doing

if ( MAP_OF_ORI .eq. MAP_ENVEL ) return
make map MAP_ENVEL from MAP_OF_ORI around CUSHION sele active end \
     init 9999 real
make map MAP_ENVEL + 1 from MAP_OF_ORI around CUSHION sele active end \
     init integer 
make map MAP_ENVEL + 1 from MAP_OF_ORI atom ( index  CUT_OFF 1000. )
make map MAP_ENVEL from MAP_ENVEL + 1 copy
make map MAP_ENVEL set 1 1000000 9999.
make map MAP_ENVEL set -100000 0 -9999.
make map MAP_ENVEL from MAP_OF_ORI copy
delete map MAP_ENVEL + 1
image map MAP_ENVEL dens 1.0 15 15 15
return