subroutine char SEGMENTS integer MAP_WORK
! pay attention to the  "main_2_sigmaa_2_main.ccp4" macro which should
! be present in your working directory.  You can find a sample file in
! $MAIN/doc/dens_mod/phase_comb/ directory.

<?MAIN_CMDS:re_phase.cmds SEGMENTS MAP_WORK

write over file main_fcalc.phase refl polar select WORK_REFL end
shell main_2_sigmaa_2_main.ccp4  > sigmma.out
read file sigmaa.phase refl over resol RESOL_MIN RESOL_MAX friedel
<MAIN_CMDS:calc_fo_map.cmds MAP_WORK
return