24小时热门版块排行榜    

查看: 4269  |  回复: 1

汤小家有只鹿

铁虫 (初入文坛)

[求助] 运行fortran出现forrtl: severe (174): SIGSEGV, segmentation fault occurred 已有1人参与

(linux)我读取两个nc文件,再建立一个新的nc文件,把前两个nc文件中读取出的内容放进新的nc文件。求助,要交作业,拜托
以下是我的程序:
      ! This program is for reading area.nc and calendar
      program link         
      implicit none
      character*299 cass_cal
      character*299 cass_har
      character*299 link_ac
      integer ierr,ncid1,ncid2,ncid3,varid1,varid2,varid3,len_file
      integer, parameter :: longitude=3600, latitude=1800
      real lon_1(longitude), lat_1(latitude)
      real lon_2(longitude), lat_2(latitude)
      real harveste_1(longitude,latitude)
      real(8) area_2(longitude,latitude)                 

      integer lon,lat          !dimension IDs
      integer area,harveste    !variable IDs
      integer areadims(2),hedims(2) !variable shapes           
                                 
      include 'netcdf.inc'
      cass_cal='/home/iga/data/xycao/calendar/output/Cassava.op.nc'
      cass_har='/home/iga/data/xycao/05area/cassava/op/cass_harvest_2005_01.nc'
      link_ac='/home/iga/data/xycao/05area/cassava/op/link_ac.nc'
      !read calendar
      len_file=len_trim(cass_cal)
      ierr=nf_open(trim(cass_cal),nf_nowrite,ncid1) !open nc file,get ID number (ncid)

      ierr=nf_inq_varid(ncid1,'longitude',varid1)    !get"longitude"
      ierr=nf_get_var_real(ncid1,varid1,lon_1)

      ierr=nf_inq_varid(ncid1,'latitude',varid1)    !get"latitude"
      ierr=nf_get_var_real(ncid1,varid1,lat_1)

      !ierr=nf_inq_varid (ncid1, 'plant.end', varid1)  !open 'plant.end',get ID number
      !ierr=nf_get_var_real (ncid1,varid1,plante_1)      

      ierr=nf_inq_varid (ncid1, 'harvest.end', varid1) !open 'harvest.end',get ID number
      ierr=nf_get_var_real (ncid1,varid1,harveste_1)

      !read area
      len_file=len_trim(cass_har)
      ierr=nf_open(trim(cass_har),nf_write,ncid2)  !open nc file,get ID number (ncid)

      ierr=nf_inq_varid(ncid2,'longitude',varid2) !get"longitude"
      ierr=nf_get_var_real(ncid2,varid2,lon_1)

      ierr=nf_inq_varid(ncid2,'latitude',varid2) !get"latitude"
      ierr=nf_get_var_real(ncid2,varid2,lat_2)

      ierr=nf_inq_varid (ncid2, 'Band1', varid2) !open 'Band1',get ID number
      ierr=nf_get_var_real (ncid2,varid2,area_2)
      
      print*,'ierr'
      !print*, plante(1145,702)
      print*, harveste_1(1145,704)                       !check
      print*, area_2(1145,702),area_2(1145,704)

      !create a new file to link area and calendar
      len_file=len_trim(link_ac)
      ierr=nf_create(trim(link_ac),nf_clobber,ncid3)
   
      ierr=nf_def_dim(ncid3,'lat_1',1800,lat)
      ierr=nf_def_dim(ncid3,'lon_1',3600,lon)
      hedims(1)=lat
      hedims(2)=lon
      ierr=nf_def_var(ncid3,'harveste_1',nf_real,2,harveste)
      areadims(1)=lat
      areadims(2)=lon
      ierr=nf_def_var(ncid3,'area_2',nf_double,2,area)

      end program link

运行时出现了以下错误:
forrtl: severe (174): SIGSEGV, segmentation fault occurred
Image              PC                Routine            Line        Source            
link.exe           000000000046F579  Unknown               Unknown  Unknown
link.exe           000000000046DEF0  Unknown               Unknown  Unknown
link.exe           000000000043FD52  Unknown               Unknown  Unknown
link.exe           0000000000424973  Unknown               Unknown  Unknown
link.exe           00000000004035FB  Unknown               Unknown  Unknown
libpthread.so.0    0000003B0EA0F710  Unknown               Unknown  Unknown
libnetcdff.so.5    00002B059A138AEC  Unknown               Unknown  Unknown
link.exe           0000000000403342  Unknown               Unknown  Unknown
link.exe           0000000000402EF6  Unknown               Unknown  Unknown
libc.so.6          0000003B0E61ED5D  Unknown               Unknown  Unknown
link.exe           0000000000402DE9  Unknown               Unknown  Unknown
回复此楼
小蜗牛,快快爬
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1014511134

专家顾问 (正式写手)

【答案】应助回帖

★
感谢参与,应助指数 +1
汤小家有只鹿: 金币+1, ★有帮助, 谢谢回复 2017-12-04 18:50:20
访问非法地址,查看一下数组是否越界什么的,比如合并后的数据大小等等,只是看过,没用过netcdf读写
2楼2017-12-04 15:05:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 汤小家有只鹿 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[职场人生] 40+中年人失业后,还能找到好工作吗? +9 xqy2009 2026-10-13 11/550 2026-10-15 16:37 by beefly
[有机交流] 用NBS溴化时溶剂除了CCl4外还有什么 +3 18834002860 2026-10-12 6/300 2026-10-15 14:30 by 长沙摩尔信息科�
[基金申请] 申请科研项目纯靠实力能成功吗? 200+7 David1982 2026-10-13 11/550 2026-10-15 10:30 by zyqchem
[硕博家园] 极限拉扯,如何抉择 +11 Isak0320 2026-10-12 15/750 2026-10-14 21:55 by 猫鱼不是鱼
[论文投稿] 二审4个多月了还是Awaiting Reviewer Scores 10+3 纳斯达克ML 2026-10-11 4/200 2026-10-14 17:30 by babero
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +8 wangyl_123 2026-10-10 13/650 2026-10-14 15:12 by zuesli
[教师之家] 课题组会帮你通过考核么? +4 好运来~9527 2026-10-10 7/350 2026-10-13 08:41 by xli1984
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
信息提示
请填处理意见