set vari MAP_MASK = 4
set vari RAD_COMB = 2.4
set vari RAD_SOLV = 6.

<re_phase.cmds
make map MAP_FOFC zero
make map MAP_FOFC conv complex
refl set ampl phase fwork = fcalc * 1.
!refl set ampl phase scale fcalc fwork = fcalc * 1.
refl fill-map MAP_FOFC resol 3. 8.  defined 
four map MAP_FOFC back 
make map MAP_FOFC rescale 


make map MAP_MASK from MAP_MIR init -9999 around 6 sele segm name S* end
make map MAP_MASK atom mask RAD_COMB sele segm name S* .a weight 0.9 1.1 end
make map MAP_MASK from MAP_FOFC 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_MIR rescale
make map MAP_FOFC from MAP_MIR 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_MIR init -9999 around 6 sele segm name S* end
make map MAP_MASK atom mask RAD_SOLV sele segm name S* end
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


four map MAP_FOFC 
reflection shells 10 r-values 
<calc_fo_map.cmds
return