image erase
key image sele .not all end

set vari MACRO = MY_SLIDE

image hist init
image hist atom name off
image hist pick 4638 
image hist pick 4629 
image hist pick 4605 

!save over file MACRO + .view view
< MACRO + .view

set vari PLOT = ON

if ( PLOT .eq. ON ) then

 open uni 24 file MACRO + _head.r3d over
 image plot ( raster3d unit 24 init on )
 image plot off 
 clo uni 24

end_if

!  insert code for generation of the partial images

! the ribbon plot

if ( PLOT .eq. ON ) then
 open uni 24 file MACRO + _ribbon.r3d over
 image plot ( raster3d unit 24 init on )
 rewind unit 24
end_if

key display sele segm name IA end
<?MAIN_CMDS:image_color_set.cmds 
<?MAIN_CMDS:image_col2_by_resi.cmds 
<?MAIN_CMDS:image_do_rib_stick.cmds 

if ( PLOT .eq. ON ) then
image plot off 
clo uni 24
end_if


key active sele seq 247 : 255 .a segm name IA end

! balls end sticks

if ( PLOT .eq. ON ) then
 open uni 24 file MACRO + _ball_stick.r3d over
 image plot ( raster3d unit 24 init on )
 rewind unit 24
end_if

set stick radi = 0.15
image sele active .a .not atom name H* N C O end set bond stick

set radi sele active end = 0.15
image sele active .a .not atom name H* N C O end set atom ball

if ( PLOT .eq. ON ) then
image plot off 
clo uni 24
end_if


! atom balls 

if ( PLOT .eq. ON ) then
 open uni 24 file MACRO + _ball.r3d over
 image plot ( raster3d unit 24 init on )
 rewind unit 24
end_if

set radi sele atom name CL end = 2.0
image sele atomname CL end col 198 atom ball

if ( PLOT .eq. ON ) then
image plot off 
clo uni 24
end_if

! r3d surface

if ( PLOT .eq. ON ) then
 open uni 24 file MACRO + _surf.r3d over
 image plot ( raster3d unit 24 init on )
 rewind unit 24
end_if

image col 220 sele segm name IA end poin 3angle 

if ( PLOT .eq. ON ) then
image plot off 
clo uni 24
end_if



! draw a map mask around map the active map "map_to"

if ( PLOT .eq. ON ) then
open uni 24 file MACRO + _map.r3d over
image plot ( raster3d unit 24 init on )
rewind unit 24
end_if

<MAIN_CMDS:image_map_mask.cmds 1
image map 3 stick rad 0.06
image map 3 dens map_cont stick on

if ( PLOT .eq. ON ) then
image plot off 
clo uni 24
end_if



return
end