subroutine int MAP_DENS int MAP_SOLV \
	real SOLV_RAD 

! MAP_SCORE is the map used for the score map
! MAP_FROM contains electron density to be analysed - unit cell map
! SOLV_RAD is the radius for the solvent statistics

! Copy the MIR density into the new MAP_SOLV, set all negative density
! values to zero and Fourier transform it.  R-value calculation is

show vari MAP_DENS
show vari MAP_SOLV
show vari SOLV_RAD


make map MAP_SOLV from MAP_DENS init 9999 cell real
make map MAP_SOLV from MAP_DENS copy
!make map MAP_SOLV from MAP_DENS mult 
make map MAP_SOLV set -99999 0 0.0

fourier map MAP_SOLV

reflect sele WORK_REFL end shells 10 r-values 

set vari NEW_MAP = nmaps + 1
make map NEW_MAP from MAP_DENS init -9999  cell real
make atom from coor 0.0 0.0 0.0
make map NEW_MAP sele atom numb natom end atom dist SOLV_RAD  linear
make map NEW_MAP set -99999 -100 0.0
dele atom sele atom numb natom end
fourier map NEW_MAP

make map MAP_SOLV from NEW_MAP mult
dele map NEW_MAP

make map MAP_SOLV conj
four map MAP_SOLV back 
make map MAP_SOLV prescale 

return