Wednesday, May 28, 2003, 1:21:48 AM, you wrote:
RCØ> Dear Group
RCØ> I am doing my master thesis on finite element modeling using fortran 90.
RCØ> (NAG on a Unix machine). I have now got a problem.
RCØ> My problem is that when I turn a "Write(*,*) 'Anything'" command on, the
RCØ> result gets wrong.
RCØ> All calculations are with double precision.
RCØ> I have used "IMPLICIT NONE" in all subroutines.
RCØ> My reals are declared as ex.: REAL(KIND=8) :: stress(neqn)
RCØ> All formulas are written as ex.: dfdx=dsqrt(x)*2.d0*....
RCØ> Could you give my a idea, where to look for my error?
RCØ> I am relatively new in Fortran so could you please keep the answer in basic
RCØ> terms, please?
RCØ> Ms.Sc., Rasmus C. Østergaard,
RCØ> Solids Mechanichs Department, Technical University of Denmark.
RCØ> [log in to unmask]
Dear Rasmus,
Check subscripts of arrays in your program. May be some of them goes
outside limits in array declaration. I saw the "Write Ghost" many
years ago exactly in such situation.
Regards,
Sergei Nemnyugin
|