working torsional

[unres4.git] / source / unres / io.f90

diff --git a/source/unres/io.f90 b/source/unres/io.f90
index 8505d11..9388f15 100644 (file)
--- a/source/unres/io.f90
+++ b/source/unres/io.f90
@@ -720,7 +720,7 @@
 !      integer :: rescode
 !      double precision x(maxvar)
       character(len=256) :: pdbfile
-      character(len=320) :: weightcard
+      character(len=480) :: weightcard
       character(len=80) :: weightcard_t!,ucase
 !      integer,dimension(:),allocatable :: itype_pdb   !(maxres)
 !      common /pizda/ itype_pdb
@@ -742,6 +742,7 @@
       allocate(c(3,2*maxres+2)) !(3,maxres2+2) maxres2=2*maxres
       allocate(dc(3,0:2*maxres)) !(3,0:maxres2)
       allocate(itype(maxres,5)) !(maxres)
+      allocate(istype(maxres))
 !
 ! Zero out tables.
 !
@@ -783,6 +784,7 @@
       call reada(weightcard,'DELT_CORR',delt_corr,0.5d0)
       call reada(weightcard,'TEMP0',temp0,300.0d0)
       if (index(weightcard,'SOFT').gt.0) ipot=6
+      call reada(weightcard,'WBOND_NUCL',wbond_nucl,1.0D0)
 ! 12/1/95 Added weight for the multi-body term WCORR
       call reada(weightcard,'WCORRH',wcorr,1.0D0)
       if (wcorr4.gt.0.0d0) wcorr=wcorr4
@@ -996,7 +998,7 @@
 
         do i=1,nres_molec(molec)
           istype(i)=sugarcode(sequence(i,molec)(1:1),i)
-          itype(i,1)=rescode(i,sequence(i,molec)(2:4),iscode,molec)
+          itype(i,1)=rescode(i,sequence(i,molec)(2:2),iscode,molec)
         enddo
        endif