请教一下，这个Fortran程序是否有问题，运行时候总是报错

齐天笑猴

请教大家一下，我的Fortran程序是否有问题，运行时，总是报错Segmentation fault (core dumped)，急求各位大神！

1. PROGRAM MAIN
   
2. IMPLICIT NONE
   
3. CHARACTER*40:: fdate,UTC_date
   
4. CHARACTER*40:: fdate_char,UTC_date_char
   
5. CHARACTER*40:: cTemp
   
6. CHARACTER*80:: plot_path
   
7. CHARACTER*40:: id
   
8. CHARACTER*80:: st_name='SURFACE / CHINA'
   
9. CHARACTER*80:: plat_form='FM-12 SYNOP'
   
10. CHARACTER*80:: source='MICAPS DIAMOND ONE'
    
11. INTEGER i,j,k,qc
    
12. 
    
13. INTEGER yr,mon,date,hr,stat_num  
    
14. INTEGER,ALLOCATABLE:: stat_id(:),sta_class(:)  
    
15. REAL,ALLOCATABLE:: stat_elev(:),temp(:),longi(:),lati(:),dew(:),win_dir(:),V(:),&
    
16. cloud(:),presea(:),pre3h(:),weather1(:),weather2(:),rain6h(:),&
    
17. lcloudz(:),lcloudl(:),lcloudh(:),visib(:),weathn(:),mcloudz(:),&
    
18. hcloudz(:),mark1(:), mark2(:), vtemp24(:), vpre24(:),psfc(:),stat_class(:)
    
19. 
    
20. REAL a(13)
    
21. INTEGER b(13)
    
22. LOGICAL is_sound,bogus,discard
    
23. CHARACTER*40:: stat_name='aaa'
    
24. !.......................INCLUDE.................
    
25.        INCLUDE 'PARMS3.EXT'
    
26.        INCLUDE 'FDESC3.EXT'
    
27.        INCLUDE 'IODECL3.EXT'
    
28.        INCLUDE 'netcdf.inc'
    
29. !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    
30. 
    
31. fdate='15010208'   !beijing time
    
32. UTC_date='15010200' !UTC
    
33. plot_path='/home/sxw/xiaowei/' !PATH
    
34. 
    
35. !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    
36. 
    
37. 
    
38. qc=0
    
39. is_sound=.false.
    
40. bogus=.false.
    
41. discard=.false.
    
42. fdate_char='20'//trim(fdate)//'0000'
    
43. UTC_date_char='20'//trim(UTC_date)//'0000'
    
44. a=-888888.0
    
45. b=0
    
46. 
    
47. OPEN (101,FILE=trim(plot_path)//'plot/'//trim(adjustl(fdate))//'.000')
    
48. READ(101,*)
    
49. READ(101,*) yr,mon,date,hr,stat_num
    
50. CLOSE(101)
    
51. 
    
52. ALLOCATE (stat_id(stat_num))
    
53. ALLOCATE (longi(stat_num))
    
54. ALLOCATE (lati(stat_num))
    
55. ALLOCATE (stat_elev(stat_num))
    
56. ALLOCATE (temp(stat_num))
    
57. ALLOCATE (dew(stat_num))
    
58. ALLOCATE (win_dir(stat_num))
    
59. ALLOCATE (V(stat_num))
    
60. ALLOCATE (sta_class(stat_num))
    
61. ALLOCATE (cloud(stat_num))
    
62. ALLOCATE (presea(stat_num))
    
63. ALLOCATE (pre3h(stat_num))
    
64. ALLOCATE (weather1(stat_num))
    
65. ALLOCATE (weather2(stat_num))
    
66. ALLOCATE (rain6h(stat_num))
    
67. ALLOCATE (lcloudz(stat_num))
    
68. ALLOCATE (lcloudl(stat_num))
    
69. ALLOCATE (lcloudh(stat_num))
    
70. ALLOCATE (visib(stat_num))
    
71. ALLOCATE (weathn(stat_num))
    
72. ALLOCATE (mcloudz(stat_num))
    
73. ALLOCATE (hcloudz(stat_num))
    
74. ALLOCATE (mark1(stat_num))
    
75. ALLOCATE (mark2(stat_num))
    
76. ALLOCATE (vtemp24(stat_num))
    
77. ALLOCATE (vpre24(stat_num))
    
78. ALLOCATE (stat_class(stat_num))
    
79.       OPEN(101,FILE=trim(plot_path)//'plot/'//trim(adjustl(fdate))//'.000')
    
80.                 READ(101,*)
    
81.                 READ(101,*)
    
82.                 DO i=1,stat_num
    
83.                     READ(101,*) stat_id(i),longi(i),lati(i),stat_elev(i),sta_class(i),cloud(i),&
    
84.                                 win_dir(i),V(i),presea(i),pre3h(i),weather1(i),weather2(i),&
    
85.                                 rain6h(i),lcloudz(i),lcloudl(i),lcloudh(i),dew(i),visib(i),&
    
86.                                 weathn(i),temp(i),mcloudz(i),hcloudz(i),mark1(i),mark2(i),vtemp24(i),vpre24(i)
    
87.                 ENDDO
    
88.       CLOSE(101)
    
89. 
    
90. OPEN(101,FILE='obs.20'//trim(adjustl(UTC_date))//'',ACCESS='append')
    
91.        DO i=1,stat_num
    
92.                
    
93.                WRITE(ctemp,'(i5)') stat_id(i)
    
94.                WRITE(id,'(i5)') stat_id(i)
    
95.                WRITE(101,'(2f20.5,4A40,f20.5,5I10,3L10,2I10,A20,13(f13.5,I7))') lati(i), longi(i),adjustl(id),&
    
96.                                                                                 st_name,plat_form,source,stat_elev(i),-888888,-888888,-888888,&
    
97.                                                                                 i,-888888,is_sound,bogus,discard,-888888,-888888,trim(UTC_date_char),&  
    
98.                                                                                 -888888.,0,((a(k),b(k)),k=1,12)
    
99.                                 IF(presea(i)<=500.) psfc(i)=1000+presea(i)/10
    
100.                                 IF(presea(i)>500.) psfc(i)=900+presea(i)/10
     
101.                                 IF(presea(i)==9999.) psfc(i)=-888888.
     
102.                                 IF(temp(i)==9999.) temp(i)=-888888.-273.16
     
103.                                 IF(dew(i)==9999.) dew=-888888.-273.16
     
104.                                 IF(win_dir(i)==9999.) win_dir(i)=-888888.
     
105.                                 IF(V(i)==9999.) V(i)=-888888.
     
106.                                 WRITE(101,'(10(F13.5,I7))') psfc(i)*100,qc,stat_elev(i),qc,(temp(i)+273.16),qc,(dew(i)+273.16),qc,&
     
107.                                                             V(i),qc,win_dir(i),qc,((a(k),b(k)),k=1,4)
     
108. 
     
109.                 WRITE(101,'(10(F13.5,I7))') -777777.,0,-777777.,0,((a(k),b(k)),k=1,8)
     
110.                 WRITE(101,*) 50,0,0
     
111.        ENDDO
     
112. CLOSE(101)
     
113. 
     
114. DEALLOCATE (stat_id)
     
115. DEALLOCATE (longi)
     
116. DEALLOCATE (lati)
     
117. DEALLOCATE (stat_elev)
     
118. DEALLOCATE (stat_class)
     
119. 
     
120. DEALLOCATE (temp)
     
121. 
     
122. DEALLOCATE (win_dir)
     
123. DEALLOCATE (V)
     
124. DEALLOCATE (cloud)
     
125. DEALLOCATE (presea)
     
126. DEALLOCATE (pre3h)
     
127. DEALLOCATE (weather1)
     
128. DEALLOCATE (weather2)
     
129. DEALLOCATE (rain6h)
     
130. DEALLOCATE (lcloudz)
     
131. DEALLOCATE (lcloudl)
     
132. DEALLOCATE (lcloudh)
     
133. DEALLOCATE (visib)
     
134. DEALLOCATE (weathn)
     
135. DEALLOCATE (mcloudz)
     
136. DEALLOCATE (hcloudz)
     
137. DEALLOCATE (mark1)
     
138. DEALLOCATE (mark2)
     
139. DEALLOCATE (vtemp24)
     
140. DEALLOCATE (vpre24)
     
141. 
     
142. 
     
143. 
     
144. END

复制代码

lqouc

报错信息不全，一般自己百度报错信息就能找到解决方案

齐天笑猴

lqouc 发表于 2018-1-10 10:28
报错信息不全，一般自己百度报错信息就能找到解决方案

我在网上查看了，但是仍然没有找到，用gdb调试查找core报错文件，显示如下：

1. Program terminated with signal 11, Segmentation fault.
   
2. #0  0x000000000040e4e4 in main () at ./little_r_micaps.f90:91
   
3. 91                     WRITE(101,'(2f20.5,4A40,f20.5,5I10,3L10,2I10,A20,13(f13.5,I7))') lati(i), longi(i),adjustl(id),&
   
4. Missing separate debuginfos, use: debuginfo-install cyrus-sasl-lib-2.1.23-13.el6_3.1.x86_64 glibc-2.12-1.132.el6.x86_64 keyutils-libs-1.4-4.el6.x86_64 krb5-libs-1.10.3-10.el6_4.6.x86_64 libcom_err-1.41.12-18.el6.x86_64 libcurl-7.19.7-37.el6_4.x86_64 libidn-1.18-2.el6.x86_64 libselinux-2.0.94-5.3.el6_4.1.x86_64 libssh2-1.4.2-1.el6.x86_64 nspr-devel-4.10.0-1.el6.x86_64 nss-3.15.1-15.el6.x86_64 nss-softokn-freebl-3.14.3-9.el6.x86_64 nss-util-3.15.1-3.el6.x86_64 openssl-1.0.1e-15.el6.x86_64 zlib-1.2.3-29.el6.x86_64
   

复制代码

但是感觉这个地方没有错误啊，不知道您知道吗？

lqouc

齐天笑猴 发表于 2018-1-10 17:06
我在网上查看了，但是仍然没有找到，用gdb调试查找core报错文件，显示如下：

但是感觉这个地方没有错 ...

看上去没问题，我也不知道该怎么修改，只能多尝试了
比如trim adjust 这些都直接在前面定义成新的变量，带着函数在输出语句里容易出错
((a(k),b(k)),k=1,12)这个最好也都拆开写

齐天笑猴

哦，好的。非常感谢。我尝试一下