set lev 0
set vari SOLV_LEV = 0.0
set vari RAD_COMB = 2.4
set vari RAD_SOLV = 3.6
key act_mask sele by resi weight 0.9 1.1 .a work_set end
key pas_mask sele final end

<re_phase_comb.cmds
<calc_fc_map.cmds

make map MAP_MASK from MAP_FO init -9999 around 6 \
                  sele act_mask end
make map MAP_MASK atom mask RAD_COMB sele act_mask end
make map MAP_MASK from MAP_FCALC copy

make map MAP_2FOFC zero
make map MAP_2FOFC conv real
make map MAP_2FOFC set -1 1 -9999.
make map MAP_2FOFC from MAP_MASK cell

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

make map MAP_2FOFC from MAP_FO rescale
! make map MAP_2FOFC rescale  ! worse 
make map MAP_FOFC from MAP_FO copy
make map MAP_FOFC from MAP_2FOFC merge  higher ! weight 0.5 0.5 !

! solvent flattening mit combinierte dichte
make map MAP_MASK from MAP_FO init -9999 around 6 sele pas_mask end
make map MAP_MASK atom radii mask sele pas_mask end
!make map MAP_MASK + 1 from MAP_MASK init 9999 copy
!make map MAP_MASK from MAP_MASK + 1 fill 4

make map MAP_MASK from MAP_FOFC copy
make map MAP_FOFC zero
make map MAP_FOFC set -1 1 -9999.
make map MAP_FOFC from MAP_MASK cell
make map MAP_FOFC set -10000 -100 0.0
!make map MAP_FOFC set -10000 -100 9999.
!make map MAP_FOFC atom mask RAD_COMB sele segm name DERI .o SYMM_WAT end 
!make map MAP_FOFC from MAP_FO invert
!make map MAP_FOFC set -10000 -100 SOLV_LEV

four map MAP_FOFC 
reflection shells 10 r-values 
<phase_test
<calc_fo_phase_fc_map.cmds
!<calc_2fofc_map.cmds
<calc_fofc_map.cmds
set lev -2
return