- read (ithep_pdb,*,err=111,end=111) a0thet(i),(athet(j,i),j=1,2),
- & (bthet(j,i),j=1,2)
+ read (ithep_pdb,*,err=111,end=111) a0thet(i),
+ & (athet(j,i,1,1),j=1,2),
+ & (bthet(j,i,1,1),j=1,2)
read (ithep_pdb,*,err=111,end=111) (polthet(j,i),j=0,3)
read (ithep_pdb,*,err=111,end=111) (gthet(j,i),j=1,3)
read (ithep_pdb,*,err=111,end=111) theta0(i),sig0(i),sigc0(i)
read (ithep_pdb,*,err=111,end=111) (polthet(j,i),j=0,3)
read (ithep_pdb,*,err=111,end=111) (gthet(j,i),j=1,3)
read (ithep_pdb,*,err=111,end=111) theta0(i),sig0(i),sigc0(i)
read (itorp,*,end=113,err=113) nterm(i,j,iblock),
read (itorp,*,end=113,err=113) nterm(i,j,iblock),
nterm(-i,-j,iblock)=nterm(i,j,iblock)
nterm(-i,-j,iblock)=nterm(i,j,iblock)
+c write (iout,*) i,j,nterm(i,j,iblock),nlor(i,j,iblock)
read (itorp,*,end=113,err=113) kk,v1(k,i,j,iblock),
& v2(k,i,j,iblock)
v1(k,-i,-j,iblock)=v1(k,i,j,iblock)
read (itorp,*,end=113,err=113) kk,v1(k,i,j,iblock),
& v2(k,i,j,iblock)
v1(k,-i,-j,iblock)=v1(k,i,j,iblock)
do j=1,ntortyp
write (iout,*) 'ityp',i,' jtyp',j
write (iout,*) 'Fourier constants'
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)
write (iout,'(3(1pe15.5))')
& vlor1(k,i,j),vlor2(k,i,j),vlor3(k,i,j)
enddo
write (iout,'(3(1pe15.5))')
& vlor1(k,i,j),vlor2(k,i,j),vlor3(k,i,j)
enddo
endif
read (itordp,*,end=114,err=114) ntermd_1(i,j,k,iblock),
& ntermd_2(i,j,k,iblock)
endif
read (itordp,*,end=114,err=114) ntermd_1(i,j,k,iblock),
& ntermd_2(i,j,k,iblock)
+ write (iout,*) ntermd_1(i,j,k,iblock),
+ & ntermd_2(i,j,k,iblock)
ntermd_1(-i,-j,-k,iblock)=ntermd_1(i,j,k,iblock)
ntermd_2(-i,-j,-k,iblock)=ntermd_2(i,j,k,iblock)
read (itordp,*,end=114,err=114) (v1c(1,l,i,j,k,iblock),l=1,
ntermd_1(-i,-j,-k,iblock)=ntermd_1(i,j,k,iblock)
ntermd_2(-i,-j,-k,iblock)=ntermd_2(i,j,k,iblock)
read (itordp,*,end=114,err=114) (v1c(1,l,i,j,k,iblock),l=1,
write (iout,'(3x,20i10)') (l,l=1,ntermd_2(i,j,k,iblock))
do l=1,ntermd_2(i,j,k,iblock)
write (iout,'(i5,20f10.5)')
write (iout,'(3x,20i10)') (l,l=1,ntermd_2(i,j,k,iblock))
do l=1,ntermd_2(i,j,k,iblock)
write (iout,'(i5,20f10.5)')
- & l,(v2s(l,m,i,j,k,iblock),m=1,ntermd_2(i,j,k,iblock))
+ & l,(v2s(l,m,i,j,k,iblock),m=1,ntermd_2(i,j,k,iblock)),
& (v2s(l,m,-i,-j,-k,iblock),m=1,ntermd_2(i,j,k,iblock))
enddo
write (iout,*)
& (v2s(l,m,-i,-j,-k,iblock),m=1,ntermd_2(i,j,k,iblock))
enddo
write (iout,*)
write (iout,*) "Coefficients of the cumulants"
endif
read (ifourier,*) nloctyp
write (iout,*) "Coefficients of the cumulants"
endif
read (ifourier,*) nloctyp
do i=0,nloctyp-1
read (ifourier,*,end=115,err=115)
read (ifourier,*,end=115,err=115) (b(ii),ii=1,13)
do i=0,nloctyp-1
read (ifourier,*,end=115,err=115)
read (ifourier,*,end=115,err=115) (b(ii),ii=1,13)