1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37
| do I=1,120
do J=1,120
do i=1,30
do j=1,30
H6(i,j)=(0.d0)
c write(*,*) 'H6(i,j)=',H6(i,j)
if(i.eq.j) then
call sub1(A,i,y1,y2)
H1(i,j)=-(y1+y2)
c write(*,*) 'H1(i,j)=',H1(i,j)
call sub5(A,i,y7)
H4(i,j)=(y7)
c write(*,*) 'H4(i,j)=',H4(i,j)
call sub3(A,i,y5)
H3(i,j)=(y5)
c write(*,*) 'H3(i,j)=',H3(i,j)
call sub11(A,i,y3,y4)
H2(i,j)=-(y3+y4)
c write(*,*) 'H2(i,j)=',H2(i,j)
else
call sub2(A,i,y11,y22)
H1(i,j)=-(y11+y22)
c write(*,*) 'H1(i,j)=',H1(i,j)
call sub4(A,i,y6)
H3(i,j)=(y6)
c write(*,*) 'H3(i,j)=',H3(i,j)
call sub6(A,i,y8)
H4(i,j)=(y8)
c write(*,*) 'H4(i,j)=',H4(i,j)
call sub22(A,i,y33,y44)
H2(i,j)=-(y33+y44)
c write(*,*) 'H2(i,j)=',H2(i,j)
endif
enddo
enddo
Hbv(i,j)=
enddo |
Partager