fortran编程过程中出现的错误，不知道怎么改了

lindaleehy

原来有个大程序，里面有一个三维数组level（，，）

我想加进去些程序，然后让三维数组中的数值和一个二维数组中的数值进行比较，可以直接比较吗？

if (level(ibj,jbj,kbj).ne.shzrwlb(izd))then
goto 3004

这样编完程序之后就出现了错误

the shapes of the array expressions do not conform

拜托各位帮忙解决下，在下也算是刚刚学，还不太懂，万分感谢！！！

四叶草

按照楼主的表述这个比较应该是可行的，但从你给的信息上，看不出问题，建议贴出全部代码。

lindaleehy

四叶草 发表于 2016-2-20 17:09
按照楼主的表述这个比较应该是可行的，但从你给的信息上，看不出问题，建议贴出全部代码。

全部代码实在是太多了。我自己编的倒是不多，
    do 3006 ir=1,200
         imax(ir)=2
         jmax(ir)=2
         kkmax(ir)=2

         imin(ir)=nrow-1
         jmin(ir)=nrow-1
         kkmin(ir)=nrow-1

         izhx(ir)=(nrow+1)/2
         jzhx(ir)=(nrow+1)/2
         kzhx(ir)=(nrow+1)/2

         shzrwlb(ir)=1
         ryx(ir)=0
         lyx(ir)=0
         nb(ir)=0

3006 continue、


上面是第一段


         ixzh=i
         jxzh=j
         kxzh=k

2002  do 2003 ixz=1,ibl
        ibian=level(ixzh,jxzh,kxzh)
        ibianbian=shzrwlb(ixz)
        if(ibian.ne.ibianbian)then
        goto 2002
      
         iq=(ixzh-izhx(ixz))**2
         iw=(jxzh-jzhx(ixz))**2
         ie=(kxzh-kzhx(ixz))**2

         D=(iq+iw+ie)**(.5)

          ryx(ixz)=0.34*nb(ixz)
          lyx(ixz)=L/R*ryx(ixz)
         
          xzh=ryx(ixz)-lyx(ixz)
        end if

        if(D.lt.xzh)then
        goto 2001
        end if
2003  continue

第二段

ibl=ibl+(inabrs-1)
        
3000     do 3001 ib=ibl-inabrs+2,ibl
        shzrwlb(ib)=inewlb+i-inblev
3001  continue


        inblev=ibl

        do 3002 kkk=2,nrow-1
          do 3002 jjj=2,nrow-1
             do 3002 iii=2,nrow-1

             if (u(iii,jjj,kkk).gt.0)then

               ibj=iii
               jbj=jjj
               kbj=kkk

3004           do 3003 izd=1,ibl
                 
                 ibianbianbian=level(ibj,jbj,kbj)
                 ibianbianbianbian=shzrwlb(izd)
                 if (ibianbianbian.ne.ibianbianbianbian)then
                 goto 3004

                 if(ibj.gt.imax(izd))then
                 imax(izd)=ibj
                 else if(ibj.lt.imin(izd))then
                        imin(izd)=ibj
                 end if

                 if(jbj.gt.jmax(izd))then
                 jmax(izd)=jbj
                 else if(jbj.lt.jmin(izd))then
                        jmin(izd)=jbj
                 end if

                 if(kbj.gt.kkmax(izd))then
                 kkmax(izd)=kbj
                 else if(kbj.lt.kkmin(izd))then
                        kkmin(izd)=kbj                    
                 end if

                 end if
3003    continue
              
              end if

3002  continue



             do 3005 izx=1,ibl
             izhx(izx)=(imax(izx)+imin(izx))/2
             jzhx(izx)=(jmax(izx)+jmin(izx))/2
             kzhx(izx)=(kkmax(izx)+kkmin(izx))/2

             ibc=imax(izx)-imin(izx)
             jbc=jmax(izx)-jmin(izx)
             kbc=kkmax(izx)-kkmin(izx)

             ibbjj=ibc
             jbbjj=jbc
             kbbjj=kbc

             bbjj=ibbjj

             if (jbbjj.lt.bbjj)then
             bbjj=jbbjj
             else if (kbbjj.lt.bbjj)then
                  bbjj=kbbjj
             end if

             nb(izx)=bbjj

3005        continue
第三段

上面的是插入到原来的长程序里的，原来的程序没有问题，就是代码太长了

lindaleehy

level的最开始定义

integer level(71,71,71)

中间我看也有用level直接比较的代码
比如：
if(level(ii,jj,kk).eq.max)goto 151
里面max是等于20000的数值

四叶草

lindaleehy 发表于 2016-2-21 10:54
level的最开始定义

integer level(71,71,71)

   if (ibianbianbian.ne.ibianbianbianbian)then
      goto 3004
上面的两语句和你提问的语句表达的是一个意思，你的程序里用的是哪一种？你的程序错误提示出在level()这个数组，你回过去好好看看，这个数组的定义是否正确？以前的程序定义了，是动态数组level(,,) ？你又重新定义过，integer level(71,71,71)？