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 TRO1
set vari FILE_MASK = mask_TRO1.xmap
read file FILE_MASK map xpl over MAP_WORK
make map MAP_WORK from MAP_FROM copy

make map MAP_WORK rescale
make map MAP_TO from MAP_WORK cell


! first copy density into mask of TRO2
set vari FILE_MASK = mask_TRO2.xmap
read file FILE_MASK map xpl over MAP_WORK
make map MAP_WORK from MAP_FROM copy

make map MAP_WORK rescale
make map MAP_TO from MAP_WORK cell


! first copy density into mask of WAT
set vari FILE_MASK = mask_WAT.xmap
read file FILE_MASK map xpl over MAP_WORK
make map MAP_WORK from MAP_FROM copy

make map MAP_WORK rescale
make map MAP_TO from MAP_WORK cell


! first copy density into mask of PPAC
set vari FILE_MASK = mask_PPAC.xmap
read file FILE_MASK map xpl over MAP_WORK
make map MAP_WORK from MAP_FROM copy

make map MAP_WORK rescale
make map MAP_TO from MAP_WORK cell


<?MAIN_UTILS:re_fft_map MAP_TO MAP_FROM
rewind file