Re: Fortran overflow
From: Leonid Gibiansky - lgibiansky@emmes.com
Subject: Re: [NMusers] Fortran overflow
Date: 2/12/2004 9:31 AM
As a first try:
With overflow error, there might be a division by zero. If you look
on your control stream, one place is
S3 = VSS-S2
K32 = Q/S3
Here S3 can be anything. You may try
IF(S3.LE.0.001) EXIT
K32 = Q/S3
Alternative is to use
S2 = THETA(4)*EXP(ETA(4))
S3 = THETA(5)*EXP(ETA(5))
VSS = S2 +S3
Leonid