subroutine int MAP_FROM int MAP_TO int CYCLE int CYCLE_N
set vari CYCLE = CYCLE + 1
show vari CYCLE
if ( CYCLE .gt. CYCLE_N ) return

set vari MAP_WORK = MAP_TO + 1
set vari MAP_ADD = MAP_WORK + 1

make map MAP_TO zero
make map MAP_TO conv real 
make map MAP_TO set -1 1 -9999

! rotating and adding density maps together


! first copy density into mask of MOLA
set vari FILE_MASK = mask_MOLA.xmap
read file FILE_MASK map xpl over MAP_WORK
make map MAP_WORK from MAP_FROM copy
< mol_MOLA_to_MOLB.com
<?MAIN_UTILS:add_a_map FILE_MASK MAP_ADD MAP_FROM MAP_WORK

make map MAP_WORK rescale
make map MAP_TO from MAP_WORK cell


! first copy density into mask of MOLB
set vari FILE_MASK = mask_MOLB.xmap
read file FILE_MASK map xpl over MAP_WORK
make map MAP_WORK from MAP_FROM copy
< mol_MOLB_to_MOLA.com
<?MAIN_UTILS:add_a_map FILE_MASK MAP_ADD MAP_FROM MAP_WORK

make map MAP_WORK rescale
make map MAP_TO from MAP_WORK cell


! first copy density into mask of IA
set vari FILE_MASK = mask_IA.xmap
read file FILE_MASK map xpl over MAP_WORK
make map MAP_WORK from MAP_FROM copy
< mol_IA_to_IB.com
<?MAIN_UTILS:add_a_map FILE_MASK MAP_ADD MAP_FROM MAP_WORK

make map MAP_WORK rescale
make map MAP_TO from MAP_WORK cell


! first copy density into mask of IB
set vari FILE_MASK = mask_IB.xmap
read file FILE_MASK map xpl over MAP_WORK
make map MAP_WORK from MAP_FROM copy
< mol_IB_to_IA.com
<?MAIN_UTILS:add_a_map FILE_MASK MAP_ADD MAP_FROM MAP_WORK

make map MAP_WORK rescale
make map MAP_TO from MAP_WORK cell


<?MAIN_UTILS:re_fft_map MAP_TO MAP_FROM
rewind file