rr0(i)=0.0D0
a0thet(i)=0.0D0
do j=1,2
rr0(i)=0.0D0
a0thet(i)=0.0D0
do j=1,2
- athet(j,i)=0.0D0
- bthet(j,i)=0.0D0
+ do k=-1,1
+ do kk=-1,1
+ athet(j,i,k,kk)=0.0D0
+ bthet(j,i,k,kk)=0.0D0
+ enddo
+ enddo
enddo
do j=0,3
polthet(j,i)=0.0D0
enddo
do j=0,3
polthet(j,i)=0.0D0
enddo
nlob(ntyp1)=0
dsc(ntyp1)=0.0D0
enddo
nlob(ntyp1)=0
dsc(ntyp1)=0.0D0
do i=1,maxtor
itortyp(i)=0
do j=1,maxtor
do k=1,maxterm
do i=1,maxtor
itortyp(i)=0
do j=1,maxtor
do k=1,maxterm
- v1(k,j,i)=0.0D0
- v2(k,j,i)=0.0D0
+ v1(k,j,i,iblock)=0.0D0
+ v2(k,j,i,iblock)=0.0D0
do i=1,maxres
itype(i)=0
itel(i)=0
do i=1,maxres
itype(i)=0
itel(i)=0