subroutine int MAP_IN int MAP_OUT 
!The zero flat
make map MAP_IN set -10000 -9900 0.0

fourier map MAP_IN 
reflect sele WORK_REFL end shell 10  r-values
 
! ! ! At this point you can call the MAIN atomic probability function
! ! if ( HIST_TEMP .gt. 2.0 ) then
! ! refl select WORK_REFL end set ampl phase fwork = fcalc * 1
! ! <>utils/hist_match MAP_IN HIST_TEMP
! ! end_if

<?MAIN_CMDS:calc_fofc_diff_map.cmds MAP_OUT 2.0 1.0
 return