subroutine char SEGMENTS int THE_MAP 

make map THE_MAP conv real
make map THE_MAP zero
make map THE_MAP set -1 1 -9999.

set vari ATOM_RAD = 6.0 
set vari MAP_WORK = 3

<?MAIN_UTILS:make_atom_mask  SEGMENTS THE_MAP ATOM_RAD MAP_WORK  mask.xmap
make map 3 set 9000 99999 1.0

make map THE_MAP from MAP_WORK cell

make map THE_MAP set -99990. -9990 0.0

fourier map THE_MAP 

refl set ampl phase temp 50.0 const 0.35 scale fcalc fpart = fcalc * 1.0

< ?MAIN_CMDS:re_phase.cmds SEGMENTS THE_MAP

refl set ampl phase fcalc = fcalc * 1. - fpart * 1.0  

< MAIN_CMDS:calc_fofc_diff_map.com THE_MAP 2.0 1.0 

return