[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

feff8.2 rmax vs ratx check



Hi all,

I upgraded from feff800 to feff82 and found some strange behaviour (with the precompiled DOS version as well as with the version I generated myself from the downloaded source code on an AIX system):
During my EXAFS calculations I always got the warning message

   WARNING:  rmax > distance to most distant atom.
             Some paths may be missing.
             rmax, ratx   1.00000E+01  0.00000E+00

although the most distant atoms in the input file were at 15A.
I compared the source code of both feff versions and found the line

         if (rtmp .gt. ratx)  ratx = rtmp

(subroutine paths, line 18933 in feff800.f) but not in feff82.f. After inserting this line after line 19032 into the respective feff82.f subroutine, the behaviour of the warning message appears correct.

As I don't know fortran well and since I'm not familiar with the detailed structure of feff, I'd like to know whether my "bugfix" is complete or if the omitted setting of ratx affects other parts of the program as well.

Thanks,
Peter

PS: Somewhere I read that reports on portability problems are welcome. Here is a minor one: To compile feff on our AIX 4.3 machine I had to rename the subroutine "trap". Otherwise I got an error message complaining about the wrong number of parameters for the INTRINSIC subroutine trap.
--
--------------------------------------------------------------
Peter Pfalzer
Universitaet Augsburg                    Tel: +49-821-598-3215
Lehrstuhl fuer Experimentalphysik II     Fax: +49-821-598-3411
Universitaetsstr. 1
D-86135 Augsburg
Germany                   Peter.Pfalzer@physik.uni-augsburg.de
--------------------------------------------------------------