do i=-ntyp,-1
itortyp(i)=-itortyp(-i)
enddo
-c write (iout,*) 'ntortyp',ntortyp
+ write (iout,*) 'ntortyp',ntortyp
do i=0,ntortyp-1
- do j=-ntortyp,ntortyp
+ do j=-ntortyp+1,ntortyp-1
read (itorp,*,end=113,err=113) nterm(i,j,iblock),
& nlor(i,j,iblock)
nterm(-i,-j,iblock)=nterm(i,j,iblock)
& v2(k,i,j,iblock)
v1(k,-i,-j,iblock)=v1(k,i,j,iblock)
v2(k,-i,-j,iblock)=-v2(k,i,j,iblock)
- v0ij=v0ij+si*v1(k,i,j)
+ v0ij=v0ij+si*v1(k,i,j,iblock)
si=-si
+ write(iout,*) i,j,k,iblock,nterm(i,j,iblock),v1(k,-i,-j,iblock)
enddo
do k=1,nlor(i,j,iblock)
read (itorp,*,end=113,err=113) kk,vlor1(k,i,j),
do j=1,ntortyp
write (iout,*) 'ityp',i,' jtyp',j
write (iout,*) 'Fourier constants'
- do k=1,nterm(i,j)
- write (iout,'(2(1pe15.5))') v1(k,i,j),v2(k,i,j)
+ do k=1,nterm(i,j,iblock)
+ write (iout,'(2(1pe15.5))') v1(k,i,j,iblock),
+ & v2(k,i,j,iblock)
enddo
write (iout,*) 'Lorenz constants'
- do k=1,nlor(i,j)
+ do k=1,nlor(i,j,iblock)
write (iout,'(3(1pe15.5))')
& vlor1(k,i,j),vlor2(k,i,j),vlor3(k,i,j)
enddo