subroutine char SEGMENTS int THE_MAP 

! renames new waters to prevent deletions of already accepted solvents
rename segm WAT sele segm name WAT2 end  
set vari II int segm sele segm name WAT end
if ( II .gt. 0 ) then
set vari WORK_SEGM char global - 
    segm sele segm name WORK_SEGM .o segm name WAT end
end_if

set weight sele atom name H* .a segm name SEGMENTS end = 0.0 
set temp sele atom name H* .a segm name SEGMENTS end = 0.0 
set temp sele temp -0.1 0.1 .a .not atom name H* .a \ 
    segm name SEGMENTS end = 25.0 ! resets Bs for rebuilt atoms


define sele segm name SEGMENTS end init group single

ener b-val group = -1 \ 
     b-val bond = 0.004 \ 
     b-val angle = 0.003 \ 
     b-val mini = 5 \ 
     b-val maxi = 60

! no kick applied 

show temp sele weight 0.9 1.1 end 
show temp sele weight 0.9 1.1 end bond

if ( REF_METHOD .eq.  "REF_DIFF" ) then
ener dens on \ 
     dens map 3 \ 
     b-val scal = 5000

make map 3 from 1 init cell real

 mini b-val individua  \ 
      sele segm name SEGMENTS end \ 
     phase sele WORK_REFL end \ 
      step = 20 \ 
      write = 1 \ 
      gain = 100.
dele map 3
end_if

if ( REF_METHOD .eq. "REF_XLSQ" ) then
 ener x-targ on \ 
      x-targ scale = 800 \ 
      x-targ map 3
 make map 3 from 1 init cell real 
 mini b-val individua \ 
      sele segm name = SEGMENTS end \  
      step = 20 \ 
     phase sele WORK_REFL end \ 
      writ = 1
 delete map 3
 ener dens on
end_if


show temp sele weight 0.9 1.1 end 
show temp sele weight 0.9 1.1 end bond

mini posit writ 10 step 80
return