太感激了

首先感谢老师您在百忙中，在凌晨1点给我指导。谢谢您！
  根据我的要求，我又思考了下，应该是每个p对应一个conr(p)，而不是每个i对应一个conr(i)。根据您的指导，再结合我的程序我做了修改：
q=1
  do i=1,4*m*n-1
   conduc=0.0
   k=0
   do p=q,450000 
     if((eigenvalue(i)<=Ef(p)).and.(eigenvalue(i+1)>Ef(p))) then ！如果p满足条件，就执行下面的操作，再看下满足这个条件的p有几个
        q=p+1
        k=k+1
         do temp=1,i
         do j=i+1,4*m*n
          ea=matmul(traneigv(temp,:),hx)
          ec=matmul(traneigv(temp,:),hy)
          eb=matmul(traneigv(j,:),hx)
          ed=matmul(traneigv(j,:),hy)
          conduc=conduc+(dot_product(ea,eigenvector(:,j))*dot_product(ed,eigenvector(:,temp))    
        end do       
       end do  
  conr(p)=2*pi*aimag(conduc)/S
       end if    
    end do
     
  end do

[quote]首先感谢老师您在百忙中，在凌晨1点给我指导。谢谢您！
  根据我的要求，我又思考了下，应该是每个p对应一个conr(p)，而不是每个i对应一个conr(i)。根据您的指导，再结合我的程序我做了修改：
q=1
  do i=1,4*m*n-1
   conduc=0.0
   k=0
   do p=q,450000 
     if((eigenvalue(i)<=Ef(p)).and.(eigenvalue(i+1)>Ef(p))) then ！如果p满足条件，就执行下面的操作，再看下满足这个条件的p有几个
        q=p+1
        k=k+1
         do temp=1,i
         do j=i+1,4*m*n
          ea=matmul(traneigv(temp,:),hx)
          ec=matmul(traneigv(temp,:),hy)
          eb=matmul(traneigv(j,:),hx)
          ed=matmul(traneigv(j,:),hy)
          conduc=conduc+(dot_product(ea,eigenvector(:,j))*dot_product(ed,eigenvector(:,temp))    
        end do       
       end do  
  conr(p)=2*pi*aimag(conduc)/S
       end if    
    end do
     
  end do[/quote]

    你一会儿每个p对应一个conr(p),一会儿每个i对应一个conr(i)?!每个p对应一个conr(p),又回到了我第17楼的贴子,用那段程序即可!
    每个i找满足的p,就先对i循环,每个p找对应的i,就先对p循环,这样程序可读性强一些!
    关键就是到底是每个p对应一个conr(p),还是i对应一个conr(i)?从一开始我就揣摩你到底要哪一个,你自己似乎对这个都捉摸不定!

哈哈，是每个p对应一个conr(p)

用17楼的帖子吧。哈哈