| 查看: 1091 | 回复: 0 | ||
[求助]
fortran运行提示错误：program exception—Access violation
|
本人新手，运行出错了，我不知道该怎么改，求助，求助高手，大牛帮忙看看吧 program main ! ! implicit none ! integer(4) :: i,j,k,n,pnblk,nbl integer(4), dimension(  , allocatable :: pim,pjm,pkminteger(4), dimension( , allocatable :: idim,jdim,kdimreal(8), dimension(:,:, , pointer :: x1,y1,z1real(8), dimension(:,:, , pointer :: x,y,ztype thblk real(8), dimension(:,:, , pointer :: x,y,zend type thblk ! type twblk real(8), dimension(:,:, , pointer :: x1,y1,z1end 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 |
回复此楼» 猜你喜欢
人气不行了
已经有4人回复
-
看来今天不会放榜了？
已经有5人回复
-
我面上完蛋了
已经有8人回复
-
只有每年这种时候来逛逛小木虫
已经有25人回复
-
什么时候开奖？
已经有7人回复
-
时间戳又变了
已经有18人回复
-
时间戳又变了8-15
已经有27人回复
-
今天基金会出结果吗？20260819
已经有16人回复
-
科研孤儿太难了
已经有16人回复
-
应该是下周三26日公布了吧？
已经有4人回复
40