set vari SOLV_METHOD = SOLOMON

set vari MAP_ORIG = 1  ! usual place for FOBS or 2FO-FC map
set vari MAP_DENS = 2  ! working map
set vari MAP_SOLV = 3
set vari SOLV_CUT = 0.50
set vari SOLV_RAD = 6.0
set vari HOLE_SIZE = 6
set vari HIST_TEMP = 0.0

! ========== end of file header

! does the map modifications
make map MAP_DENS conv real
make map MAP_DENS zero
make map MAP_DENS set -1 1 9999.
make map MAP_DENS from MAP_ORIG copy

make map MAP_SOLV from MAP_DENS init -9999

if ( SOLV_METHOD .eq. "SOLOMON" ) then
< ?MAIN_UTILS:solv_flat_solomon MAP_DENS MAP_SOLV MAP_SOLV + 1 \
	SOLV_RAD SOLV_CUT 
end_if

if ( SOLV_METHOD .eq. "WANG" ) then
! classical Wang algorithm
< ?MAIN_UTILS:solv_flat_wang MAP_DENS MAP_SOLV MAP_SOLV + 1 \
	SOLV_RAD SOLV_CUT 
end_if

if ( SOLV_METHOD .eq. "FAST" ) then
< ?MAIN_UTILS:solv_flat_fast MAP_DENS MAP_SOLV MAP_SOLV + 1 \
	SOLV_RAD SOLV_CUT 
end_if

< ?MAIN_UTILS:re_fft_map MAP_SOLV MAP_DENS 10.0

return