subroutine char SEGMENTS real KICK int MAP_REAL int MAP_IDEAL

delete atom sele segm name #* end

symm sele segm name SEGMENTS end matr 1 tran 0 0 0
copy weight sele segm name SEGMENTS end sele segm name #* end
copy temp sele segm name SEGMENTS end sele segm name #* end
show temp sele segm name SEGMENTS .a weight 0.1 1.1 end

! if you don't want to use average B values comment out the next line

set temp sele segm name #* end = RESULT_0

<?MAIN_CMDS:re_phase.cmds "#*" MAP_IDEAL

make map MAP_IDEAL zero 
make map MAP_IDEAL conv complex
refl set ampl phase scale fcalc  fwork = fcalc * 1.
refl fill-map MAP_IDEAL select defined end
four map MAP_IDEAL back 
make map MAP_IDEAL rescale

copy temp sele segm name SEGMENTS end sele segm name #* end
set coor sele segm name #* end kick KICK

anal resi init
key xx sele segm name SEGMENTS end
show key xx

set vari MAP_INDEX = nmaps + 1

set vari ISTEP = IRESULT_3 - IRESULT_2 / 12
set vari IWIND = ISTEP / 4
<?MAIN_UTILS:anal_r-omit_map_loop.com  SEGMENTS #*  MAP_REAL MAP_IDEAL \
       MAP_INDEX IRESULT_2 IRESULT_3  ISTEP IWIND

dele map MAP_INDEX

delete atom sele segm name #* end
return