Bug in get_3dfldr.F

Bug reports, work arounds and fixes

Moderators: arango, robertson

Post Reply
Message
Author
austinctodd

Bug in get_3dfldr.F

#1 Unread post by austinctodd »

I noticed through some extensive debugging that there is a sign error in get_3dfldr.F where it checks if data that has a cycle_length needs to be read. Lines 149-154 say the following:

Code: Select all

        IF (Liocycle) THEN
          Trec=MOD(Trec,Nrec)-1
          IF (Trec.le.0) Trec=Nrec-Trec
        ELSE
          Trec=Trec-1
        END IF
The issue seems to be

Code: Select all

Trec=Nrec-Trec
which makes Trec greater than nrec instead of cycling the integer. It should actually be

Code: Select all

Trec=Nrec+Trec
to be consistent with get_2dfldr.F

Post Reply