read (ithep,*,err=111,end=111) a0thet(i),(athet(j,i,1,1),j=1,2),
& (bthet(j,i,1,1),j=1,2)
read (ithep,*,err=111,end=111) (polthet(j,i),j=0,3)
read (ithep,*,err=111,end=111) a0thet(i),(athet(j,i,1,1),j=1,2),
& (bthet(j,i,1,1),j=1,2)
read (ithep,*,err=111,end=111) (polthet(j,i),j=0,3)
- read (ithep,*,err=111,end=111) (gthet(j,i),j=1,3)
- read (ithep,*,err=111,end=111) theta0(i),sig0(i),sigc0(i)
- sigc0(i)=sigc0(i)**2
+ read (ithep,*,err=111,end=111) (gthet(j,i),j=1,3)
+ read (ithep,*,err=111,end=111) theta0(i),sig0(i),sigc0(i)
+ sigc0(i)=sigc0(i)**2
enddo
do i=1,ntyp
athet(1,i,1,-1)=athet(1,i,1,1)
enddo
do i=1,ntyp
athet(1,i,1,-1)=athet(1,i,1,1)
read (ithep,*,err=111,end=111) nthetyp,ntheterm,ntheterm2,
& ntheterm3,nsingle,ndouble
nntheterm=max0(ntheterm,ntheterm2,ntheterm3)
read (ithep,*,err=111,end=111) nthetyp,ntheterm,ntheterm2,
& ntheterm3,nsingle,ndouble
nntheterm=max0(ntheterm,ntheterm2,ntheterm3)
read (ithep,*,err=111,end=111) (ithetyp(i),i=1,ntyp1)
read (ithep,*,err=111,end=111) (ithetyp(i),i=1,ntyp1)
bsc(1,i)=0.0D0
read(irotam,*,end=112,err=112)(censc(k,1,i),k=1,3),
& ((blower(k,l,1),l=1,k),k=1,3)
bsc(1,i)=0.0D0
read(irotam,*,end=112,err=112)(censc(k,1,i),k=1,3),
& ((blower(k,l,1),l=1,k),k=1,3)
- censc(1,1,-i)=censc(1,1,i)
- censc(2,1,-i)=censc(2,1,i)
- censc(3,1,-i)=-censc(3,1,i)
-
do j=2,nlob(i)
read (irotam,*,end=112,err=112) bsc(j,i)
read (irotam,*,end=112,err=112) (censc(k,j,i),k=1,3),
& ((blower(k,l,j),l=1,k),k=1,3)
do j=2,nlob(i)
read (irotam,*,end=112,err=112) bsc(j,i)
read (irotam,*,end=112,err=112) (censc(k,j,i),k=1,3),
& ((blower(k,l,j),l=1,k),k=1,3)
- censc(1,j,-i)=censc(1,j,i)
- censc(2,j,-i)=censc(2,j,i)
- censc(3,j,-i)=-censc(3,j,i)
-C BSC is amplitude of Gaussian
- do j=-ntortyp+1,ntortyp-1
- read (itorp,*,end=113,err=113) nterm(i,j,iblock),
+ 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)
nlor(-i,-j,iblock)=nlor(i,j,iblock)
v0ij=0.0d0
si=-1.0d0
& nlor(i,j,iblock)
nterm(-i,-j,iblock)=nterm(i,j,iblock)
nlor(-i,-j,iblock)=nlor(i,j,iblock)
v0ij=0.0d0
si=-1.0d0
- do k=1,nterm(i,j,iblock)
- read (itorp,*,end=113,err=113) kk,v1(k,i,j,iblock),
- & v2(k,i,j,iblock)
+ do k=1,nterm(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)
v2(k,-i,-j,iblock)=-v2(k,i,j,iblock)
v0ij=v0ij+si*v1(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,iblock)
c write(iout,*) v1(k,-i,-j,iblock),v1(k,i,j,iblock),
c &v2(k,-i,-j,iblock),v2(k,i,j,iblock)
enddo
c write(iout,*) v1(k,-i,-j,iblock),v1(k,i,j,iblock),
c &v2(k,-i,-j,iblock),v2(k,i,j,iblock)
enddo
- do k=1,nlor(i,j,iblock)
+ do k=1,nlor(i,j,iblock)
read (itorp,*,end=113,err=113) kk,vlor1(k,i,j),
read (itorp,*,end=113,err=113) kk,vlor1(k,i,j),
- & vlor2(k,i,j),vlor3(k,i,j)
+ & vlor2(k,i,j),vlor3(k,i,j)
v0ij=v0ij+vlor1(k,i,j)/(1+vlor3(k,i,j)**2)
enddo
v0(i,j,iblock)=v0ij
v0ij=v0ij+vlor1(k,i,j)/(1+vlor3(k,i,j)**2)
enddo
v0(i,j,iblock)=v0ij
write (iout,*) 'ityp',i,' jtyp',j
write (iout,*) 'Fourier constants'
do k=1,nterm(i,j,iblock)
write (iout,*) 'ityp',i,' jtyp',j
write (iout,*) 'Fourier constants'
do k=1,nterm(i,j,iblock)
- write (iout,'(2(1pe15.5))') v1(k,i,j,iblock),
+ write (iout,'(2(1pe15.5))') v1(k,i,j,iblock),
& vlor1(k,i,j),vlor2(k,i,j),vlor3(k,i,j)
enddo
enddo
& vlor1(k,i,j),vlor2(k,i,j),vlor3(k,i,j)
enddo
enddo
- read (itordp,*,end=114,err=114) ntermd_1(i,j,k,iblock),
+ read (itordp,*,end=114,err=114) 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_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))
- read (itordp,*,end=114,err=114) (v1s(1,l,i,j,k,iblock),l=1,
+ read (itordp,*,end=114,err=114)(v1s(1,l,i,j,k,iblock),l=1,
& ntermd_1(i,j,k,iblock))
& ntermd_1(i,j,k,iblock))
- read (itordp,*,end=114,err=114) (v1c(2,l,i,j,k,iblock),l=1,
+ read (itordp,*,end=114,err=114)(v1c(2,l,i,j,k,iblock),l=1,
& ntermd_1(i,j,k,iblock))
& ntermd_1(i,j,k,iblock))
- read (itordp,*,end=114,err=114) (v1s(2,l,i,j,k,iblock),l=1,
+ read (itordp,*,end=114,err=114)(v1s(2,l,i,j,k,iblock),l=1,
& ntermd_1(i,j,k,iblock))
C Martix of D parameters for one dimesional foureir series
do l=1,ntermd_1(i,j,k,iblock)
& ntermd_1(i,j,k,iblock))
C Martix of D parameters for one dimesional foureir series
do l=1,ntermd_1(i,j,k,iblock)
write (iout,*) 'Pairs of angles:'
write (iout,'(3x,20i10)') (l,l=1,ntermd_2(i,j,k,iblock))
do l=1,ntermd_2(i,j,k,iblock)
write (iout,*) 'Pairs of angles:'
write (iout,'(3x,20i10)') (l,l=1,ntermd_2(i,j,k,iblock))
do l=1,ntermd_2(i,j,k,iblock)
& l,(v2c(l,m,i,j,k,iblock),m=1,ntermd_2(i,j,k,iblock))
enddo
write (iout,*)
write (iout,'(3x,20i10)') (l,l=1,ntermd_2(i,j,k,iblock))
do l=1,ntermd_2(i,j,k,iblock)
& l,(v2c(l,m,i,j,k,iblock),m=1,ntermd_2(i,j,k,iblock))
enddo
write (iout,*)
write (iout,'(3x,20i10)') (l,l=1,ntermd_2(i,j,k,iblock))
do l=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
& 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
- read (isccor,*,end=119,err=119) nsccortyp
+ read (isccor,*,end=1113,err=1113) nsccortyp
- read (isccor,*,end=119,err=119) (isccortyp(i),i=1,ntyp)
+ write (iout,*) "Tu wchodze"
+ read (isccor,*,end=1113,err=1113) (isccortyp(i),i=1,ntyp)
- do j=1,nsccortyp
- read (isccor,*,end=119,err=119) nterm_sccor(i,j),nlor_sccor(i,j)
+ do j=1,nsccortyp
+ read (isccor,*,end=1113,err=1113) nterm_sccor(i,j),
+ & nlor_sccor(i,j)
nterm_sccor(-i,j)=nterm_sccor(i,j)
nterm_sccor(-i,-j)=nterm_sccor(i,j)
nterm_sccor(i,-j)=nterm_sccor(i,j)
nterm_sccor(-i,j)=nterm_sccor(i,j)
nterm_sccor(-i,-j)=nterm_sccor(i,j)
nterm_sccor(i,-j)=nterm_sccor(i,j)
- do k=1,nterm_sccor(i,j)
- read (isccor,*,end=119,err=119) kk,v1sccor(k,l,i,j)
+ do k=1,nterm_sccor(i,j)
+ read (isccor,*,end=1113,err=1113) kk,v1sccor(k,l,i,j)
if (j.eq.iscprol) then
if (i.eq.isccortyp(10)) then
v1sccor(k,l,i,-j)=v1sccor(k,l,i,j)
if (j.eq.iscprol) then
if (i.eq.isccortyp(10)) then
v1sccor(k,l,i,-j)=v1sccor(k,l,i,j)
v0ijsccor=v0ijsccor+si*v1sccor(k,l,i,j)
v0ijsccor1=v0ijsccor+si*v1sccor(k,l,-i,j)
v0ijsccor2=v0ijsccor+si*v1sccor(k,l,i,-j)
v0ijsccor3=v0ijsccor+si*v1sccor(k,l,-i,-j)
si=-si
enddo
v0ijsccor=v0ijsccor+si*v1sccor(k,l,i,j)
v0ijsccor1=v0ijsccor+si*v1sccor(k,l,-i,j)
v0ijsccor2=v0ijsccor+si*v1sccor(k,l,i,-j)
v0ijsccor3=v0ijsccor+si*v1sccor(k,l,-i,-j)
si=-si
enddo
- do k=1,nlor_sccor(i,j)
- read (isccor,*,end=119,err=119) kk,vlor1sccor(k,i,j),
+ do k=1,nlor_sccor(i,j)
+ read (isccor,*,end=1113,err=1113) kk,vlor1sccor(k,i,j),
& vlor2sccor(k,i,j),vlor3sccor(k,i,j)
v0ijsccor=v0ijsccor+vlor1sccor(k,i,j)/
&(1+vlor3sccor(k,i,j)**2)
& vlor2sccor(k,i,j),vlor3sccor(k,i,j)
v0ijsccor=v0ijsccor+vlor1sccor(k,i,j)/
&(1+vlor3sccor(k,i,j)**2)
read (isccor,*,end=119,err=119) (isccortyp(i),i=1,ntyp)
c write (iout,*) 'ntortyp',ntortyp
maxinter=3
read (isccor,*,end=119,err=119) (isccortyp(i),i=1,ntyp)
c write (iout,*) 'ntortyp',ntortyp
maxinter=3
write (iout,*) "Coefficients of the cumulants"
endif
read (ifourier,*) nloctyp
write (iout,*) "Coefficients of the cumulants"
endif
read (ifourier,*) nloctyp
read (ifourier,*,end=115,err=115)
read (ifourier,*,end=115,err=115) (b(ii),ii=1,13)
if (lprint) then
read (ifourier,*,end=115,err=115)
read (ifourier,*,end=115,err=115) (b(ii),ii=1,13)
if (lprint) then
- EE(1,1,-i)= b(10)+b(11)
- EE(2,2,-i)=-b(10)+b(11)
- EE(2,1,-i)=-b(12)+b(13)
- EE(1,2,-i)=-b(12)-b(13)
-