simple FORTRAN code for GEMPAK that fails on 64-bit Linux machines

• To: gembud@xxxxxxxxxxxxxxxx
• Subject: simple FORTRAN code for GEMPAK that fails on 64-bit Linux machines
• From: David Ovens <ovens@xxxxxxxxxxxxxxxxxxxx>
• Date: Tue, 24 Apr 2007 22:09:07 -0700

Gembuds,

I have this incredibly simple FORTAN code, test.f, that works fine on
32-bit Linux machines with any version of GEMPAK, but that fails
miserably on all of my 64-bit Linux machines with any version of
GEMPAK (I've tried 5.9.4 and 5.8.3ag3 specially built for 64-bit).  I
am curious if anyone out there can successfully run this on a 64-bit
Linux machine.  And, if you can, please share your
$NAWIPS/config/Makefile.linux.

Code test.f:
     program test
     call in_bdta(iret)
     u = 5
     v = 4
     spd = pr_sped(u,v)
     print *,'spd = ',spd
     stop
     end

Compile and run:
 ifort test.f $GEMLIB/gemlib.a -lg2c; a.out

Correct answer:

64-bit Linux box answer:
spd =  -1.0842022E-19


Cheers,

David
--
David Ovens              e-mail: ovens@xxxxxxxxxxxxxxxxxxxx
Research Meteorologist    phone: (206) 685-8108
Dept of Atm. Sciences      plan: Real-time MM5 forecasting for the
Box 351640                        Pacific Northwest
University of Washington          http://www.atmos.washington.edu/mm5rt
Seattle, WA  98195               Weather Graphics and Loops
                                 http://www.atmos.washington.edu/~ovens/loops

• Follow-Ups:
  • Re: simple FORTRAN code for GEMPAK that fails on 64-bit Linux machines
    • From: David Ovens