Linux下fortran报错

spring627

请教各位，有个fortran程序，编译通过了，但是调试的时候出现下面的情况

1. 
   
2. [hujinhe@console ~]$ gdb ./mdord.exe
   
3. GNU gdb Red Hat Linux (6.3.0.0-1.96rh)
   
4. Copyright 2004 Free Software Foundation, Inc.
   
5. GDB is free software, covered by the GNU General Public License, and you are
   
6. welcome to change it and/or distribute copies of it under certain conditions.
   
7. Type "show copying" to see the conditions.
   
8. There is absolutely no warranty for GDB.  Type "show warranty" for details.
   
9. This GDB was configured as "x86_64-redhat-linux-gnu"...Using host libthread_db library "/lib64/tls/libthread_db.so.1".
   
10. 
    
11. (gdb) run <in.dip> out.out
    
12. Starting program: /export/home/hujinhe/mdord.exe <in.dip> out.out
    
13. [Thread debugging using libthread_db enabled]
    
14. [New Thread 182894224864 (LWP 13692)]
    
15. 
    
16. Program received signal SIGSEGV, Segmentation fault.
    
17. [Switching to Thread 182894224864 (LWP 13692)]
    
18. 0x000000000041e998 in rand_ ()
    
19. (gdb) backtrace
    
20. #0  0x000000000041e998 in rand_ ()
    
21. #1  0x0000000000405c6a in start () at mdrdipord.f:559
    
22. #2  0x0000000000403074 in main$mdrdipord_$BLK () at mdrdipord.f:36
    
23. #3  0x000000000040302a in main ()
    
24. (gdb)
    

复制代码

刚接触编程，有些迷糊。请问是不是第36和559行有问题？

1. 
   
2. (gdb) l 36
   
3. 31      C
   
4. 32      C     READ IN AND INITIALISE DATA.
   
5. 33            CALL SYSDEF
   
6. 34      C
   
7. 35      C     SET UP STARTING CONFIGURATION
   
8. 36            CALL START
   
9. 37      C
   
10. 38      C     PRINT OUT SYSTEM PARAMETERS
    
11. 39            CALL OUTPUT(1)
    
12. 40      C***********************************************************************
    

复制代码

第36行是调用子程序START，第559行是START子程序里的一部分，如下：

1. 
   
2. (gdb) l 559
   
3. 554           COMMON/CMVELO/VMOM(3,NSP),ANGM(3,NSP)
   
4. 555           COMMON/CMFRCS/FXYZ(3,NSP),TXYZ(3,NSP)
   
5. 556           DIMENSION XB(4),YB(4),ZB(4)
   
6. 557           DATA XB/0.5,-.5,0.5,-.5/,YB/0.5,-.5,-.5,0.5/,ZB/0.5,0.5,-.5,-.5/
   
7. 558     C     XRANF(IDUM)=2.0*RANF()-1.0
   
8. 559           XRANF(IDUM)=2.0*RAND()-1.0
   
9. 560           IDUM=0
   
10. 561           IF(.NOT.NLRES)GO TO 10
    
11. 562           CALL AGAIN(1)
    
12. 563           GO TO 32
    

复制代码

可我还是看不出来是哪有问题，所以要请教各位大侠，不胜感激