read atom
*
SGX

set weight sele numb atom natom end 1.0
set temp sele numb atom natom end HIST_TEMP
make map MAP_SOLV zero
make map MAP_SOLV conv real
make map MAP_SOLV atom dens sele numb atom natom end 
delete atom sele numb atom natom end 
four map MAP_SOLV 
reflect set ampl phase fwork = fwork * 1.0  *  fcalc * 1.0

return