model_ki_dist="model"//kic2//".dist"
model_ki_angle="model"//kic2//".angle"
- open (ientin,file=model_ki_dist,status='old')
+ open (imol2,file=model_ki_dist,status='old')
do irec=1,maxdim !petla do czytania wiezow na odleglosc
- read (ientin,*,end=1401) i, j, odl_temp(i+nnt-1,j+nnt-1,ki),
+ read (imol2,*,end=1401) i, j, odl_temp(i+nnt-1,j+nnt-1,ki),
& sigma_odl_temp(i+nnt-1,j+nnt-1,ki)
odl_temp(j+nnt-1,i+nnt-1,ki)=odl_temp(i+nnt-1,j+nnt-1,ki)
sigma_odl_temp(j+nnt-1,i+nnt-1,ki)=
& sigma_odl_temp(i+nnt-1,j+nnt-1,ki)
enddo
1401 continue
- close (ientin)
- open (ientin,file=model_ki_angle,status='old')
+ close (imol2)
+ open (imol2,file=model_ki_angle,status='old')
do irec=1,maxres-3 !petla do czytania wiezow na katy torsyjne
- read (ientin,*,end=1402) i, j, k,l,dih(ki,i+nnt-1),
+ read (imol2,*,end=1402) i, j, k,l,dih(ki,i+nnt-1),
& sigma_dih(ki,i+nnt-1)
write (iout,*) i, j, k,l,dih(ki,i+nnt-1),
& sigma_dih(ki,i+nnt-1)
sigma_dih(ki,i+nnt-1)=1.0d0/sigma_dih(ki,i+nnt-1)**2
enddo
1402 continue
- close (ientin)
+ close (imol2)
enddo
ii=0
write (iout,*) "nnt",nnt," nct",nct