subroutine char SEGMENTS int MAP_REAL int MAP_IDEAL

<?MAIN_CMDS:re_phase.cmds SEGMENTS MAP_IDEAL
<?MAIN_CMDS:calc_2fofc_map.cmds MAP_REAL

symm sele segm name SEGMENTS end matr 1 tran 0 0 0
copy weight sele segm name SEGMENTS end sele segm name #* end
show temp sele segm name SEGMENTS .a weight 0.1 1.1 end
set temp sele segm name #* end = RESULT_0
<?MAIN_CMDS:re_phase.cmds "#*" MAP_IDEAL
<?MAIN_CMDS:calc_fc_map.cmds MAP_IDEAL

set vari MAP_INDEX = nmaps + 1
make map MAP_INDEX from MAP_IDEAL ini 0 inte \
         around 8 sele segm name SEGMENTS end
make map MAP_INDEX from MAP_IDEAL atom index 0.0 1000. \
                  sele segm name SEGMENTS .a weigh 0.1 1.1 end

anal resi init
anal dens ( correl base-map = MAP_REAL ideal-map = MAP_IDEAL \
            index = MAP_INDEX ) \
     sele weight 0.1 1.1 .and. segm name SEGMENTS  end
dele map MAP_INDEX
return