subroutine int THE_MAP  real RFOBS real RFCALC
make map THE_MAP zero
make map THE_MAP conv complex
if ( SIGMAA_MAP .eq. "ON" .and. RFOBS .ge. 1.9 ) then
! echo "SIGMA A 2FO-FC" SIGMAA_MAP RFOBS
refl
sigma_a 2
fill-map THE_MAP select WORK_REFL .o TEST end
exit
end_if
if ( SIGMAA_MAP .eq. "ON" .and. RFOBS .le. 1.1 ) then
! echo "SIGMA A FO-FC" SIGMAA_MAP RFOBS
refl
sigma_a 1
fill-map THE_MAP select WORK_REFL .o TEST end
exit
end_if

if ( SIGMAA_MAP .eq. "OFF" ) then
refl 
weigh const
set ampl phase fwork = fcalc * 1.
set ampl scale fcalc     fwork = fobs * RFOBS - fcalc * RFCALC
fill-map THE_MAP select  WORK_REFL .o TEST end
exit
end_if

four map THE_MAP back 
make map THE_MAP rescale 
return