24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1063  |  回复: 0

wuyifeijx

铜虫 (初入文坛)

[求助] fortran运行提示错误:program exception—Access violation

本人新手,运行出错了,我不知道该怎么改,求助,求助高手,大牛帮忙看看吧
program main
  !   
  !
  implicit none
  !
  integer(4) :: i,j,k,n,pnblk,nbl
  integer(4), dimension(, allocatable :: pim,pjm,pkm
  integer(4), dimension(, allocatable :: idim,jdim,kdim
  real(8), dimension(:,:,, pointer  ::  x1,y1,z1
  real(8), dimension(:,:,, pointer  ::  x,y,z
type thblk
    real(8), dimension(:,:,, pointer  ::  x,y,z
  end type thblk
  !
  type twblk
    real(8), dimension(:,:,, pointer  ::  x1,y1,z1
  end type twblk
  type(thblk), allocatable :: pblk3(
  type(twblk), allocatable :: pblk2(
  allocate(pblk2(1:nbl))
  allocate(pblk3(1:pnblk))
  allocate(pim(1:pnblk),pjm(1:pnblk),pkm(1:pnblk))
  allocate(idim(1:nbl),jdim(1:nbl),kdim(1:nbl))
! allocate arrays
      do n=1,nbl
            allocate(pblk2(n)%x1(0:idim(n),0:jdim(n),0:kdim(n)),                &
                 pblk2(n)%y1(0:idim(n),0:jdim(n),0:kdim(n)),                 &
                 pblk2(n)%z1(0:idim(n),0:jdim(n),0:kdim(n)))
          end do
          !   
      do n=1,pnblk
        allocate(pblk3(n)%x(0:pim(n),0:pjm(n),0:pkm(n)),              &
                 pblk3(n)%y(0:pim(n),0:pjm(n),0:pkm(n)),              &
                 pblk3(n)%z(0:pim(n),0:pjm(n),0:pkm(n)))
      end do
      !
  open(11,file='backstep5_3levdn.p3dfmt')
   rewind(11)                  
      read(11,*) nbl
      read(11,*) (idim(n),jdim(n),kdim(n),n=1,nbl)
     !
      do n=1,nbl
       read(11,*) (((x1(i,j,k),i=1,idim(n)),j=1,jdim(n)),k=1,kdim(n)),   &
                  (((y1(i,j,k),i=1,idim(n)),j=1,jdim(n)),k=1,kdim(n)),  &
                    (((z1(i,j,k),i=1,idim(n)),j=1,jdim(n)),k=1,kdim(n))   
       end do
           !
      close(11)
!
   do n=1,pnblk
    do k=0,pkm(n)
       do j=0,pjm(n)
        do i=0,pim(n)
          pblk3(n)%x=pblk2(n)%x1
          pblk3(n)%y=pblk2(n)%y1
          pblk3(n)%z=pblk2(n)%z1  
        end do
       end do
    end do
   end do
!  
  open(11,file='prep.grd',form='unformatted')
   rewind(11)
      do n=1,pnblk        
                  write(11) pblk3(n)%x,pblk3(n)%y,pblk3(n)%z
      end do
      close(11)
  !
end program main
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wuyifeijx 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自然上会要求 +5 无名者登山 2026-05-18 9/450 2026-05-18 17:50 by BlakeReary
[考博] 博士申请 +4 星…… 2026-05-18 5/250 2026-05-18 17:34 by 炎甲00
[硕博家园] 考博自荐 +5 科研狗111 2026-05-13 6/300 2026-05-18 11:22 by 糊糊涂涂好
[找工作] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 ky2p12rrjj 2026-05-15 4/200 2026-05-17 19:47 by Equinoxhua
[考博] 2026博士还有哪些学校有名额 +6 小王求读研 2026-05-15 7/350 2026-05-17 16:54 by 知音湖畔
[考博] 光量子物理方向 博士招生 1人(2026.09) +3 sandyworld 2026-05-15 4/200 2026-05-17 14:38 by sandyworld
[基金申请] 这年头没有找到涵评专家,还有中面上的可能吗 +13 dd921ww 2026-05-12 16/800 2026-05-17 12:38 by 于轩
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 v9tggjlwd0 2026-05-15 4/200 2026-05-17 08:11 by 11n4dfd8yn
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 v9tggjlwd0 2026-05-15 4/200 2026-05-17 08:06 by 11n4dfd8yn
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 l7k6xnh0yc 2026-05-14 8/400 2026-05-17 07:26 by 11n4dfd8yn
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 ky2p12rrjj 2026-05-15 4/200 2026-05-17 00:50 by ue3ir18jc3
[公派出国] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 x0mp7owy2b 2026-05-15 4/200 2026-05-17 00:35 by ue3ir18jc3
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 x0mp7owy2b 2026-05-15 4/200 2026-05-17 00:25 by ue3ir18jc3
[基金申请] 请问大佬b0816评完了吗 +4 市民华南虎 2026-05-12 8/400 2026-05-16 19:54 by Equinoxhua
[高分子] 本人最近太闲了,谁有问题可以提,每天会统一回复 +9 一切都是空工 2026-05-12 20/1000 2026-05-16 19:52 by Equinoxhua
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 x0mp7owy2b 2026-05-15 4/200 2026-05-16 17:45 by j6b2pdz07o
[有机交流] 求有机合成大神指点三硫酸乙烯酯(CAS:2793408-99-6)的合成路线 30+3 Leekmid 2026-05-13 10/500 2026-05-16 16:37 by czyzsu
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 l7k6xnh0yc 2026-05-14 6/300 2026-05-16 11:29 by h3oerqvkv9
[硕博家园] 申请博士 +3 呃?呃 2026-05-15 3/150 2026-05-16 11:01 by a4742549
[考博] 材料类只有一篇综述能申博么 +4 乐逍遥谷 2026-05-13 4/200 2026-05-14 12:05 by zhyzzh
信息提示
请填处理意见