X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Fwham%2Fsrc-HCD-5D%2Fread_constr_homology.F;h=dece50f58c97fd4479a1ceb558980b3b6d423632;hb=2bc0c5f06a45a86998d13adb237a1c84e46e12e8;hp=7884fd5a52f3766d96da2a12662c6092da2a3e46;hpb=a30bd29e64da2aa47b84963fdd0bf4192ead2738;p=unres.git

diff --git a/source/wham/src-HCD-5D/read_constr_homology.F b/source/wham/src-HCD-5D/read_constr_homology.F
index 7884fd5..dece50f 100644
--- a/source/wham/src-HCD-5D/read_constr_homology.F
+++ b/source/wham/src-HCD-5D/read_constr_homology.F
@@ -145,8 +145,10 @@ c
         unres_pdb=.false.
         if (read2sigma) then
           call readpdb_template(k)
+          close(ipdbin)
         else
           call readpdb
+          close(ipdbin)
         endif
 
 c        call readpdb
@@ -163,7 +165,6 @@ c        call readpdb
         write (iout,*) "read_constr_homology: after reading pdb file"
         call flush(iout)
 #endif
-
 c
 c     Distance restraints
 c
@@ -184,14 +185,13 @@ c         write(iout,*) "tpl_k_rescore - ",tpl_k_rescore
           do irec=nnt,nct ! loop for reading res sim 
             if (read2sigma) then
              read (ientin,*,end=1401) i_tmp,rescore2_tmp,rescore_tmp,
-     &                                idomain_tmp,
      &                                rescore3_tmp,idomain_tmp
              i_tmp=i_tmp+nnt-1
              idomain(k,i_tmp)=idomain_tmp
              rescore(k,i_tmp)=rescore_tmp
              rescore2(k,i_tmp)=rescore2_tmp
              rescore3(k,i_tmp)=rescore3_tmp
-             write(iout,'(a7,i5,2f10.5,i5)') "rescore",
+             write(iout,'(a7,i5,3f10.5,i5)') "rescore",
      &                      i_tmp,rescore2_tmp,rescore_tmp,
      &                                rescore3_tmp,idomain_tmp
             else