X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2Fio.F90;h=ffce3d7df38b702d46b416a9a93c9f62a4f1da81;hb=603baff0aabf21ffced01092a8e60ebce34f1adc;hp=b2ebc22463690f46358d9915555f12ebe0683752;hpb=6838b92b256fbeb9953f1e77eed0a0ec68b12192;p=unres4.git

diff --git a/source/unres/io.F90 b/source/unres/io.F90
index b2ebc22..ffce3d7 100644
--- a/source/unres/io.F90
+++ b/source/unres/io.F90
@@ -1099,13 +1099,24 @@
         endif
         if(.not.allocated(vbld_inv)) allocate(vbld_inv(2*nres))
         do i=2,nres
+          if (molnum(i).eq.1) then
           vbld(i)=vbl
           vbld_inv(i)=vblinv
+
+          else
+          vbld(i)=7.0
+          vbld_inv(i)=1.0/7.0
+          endif
         enddo
         do i=2,nres-1
+           if (molnum(i).eq.1) then
 !          print *, "molnum",molnum(i),itype(i,molnum(i)),nres,i 
           vbld(i+nres)=dsc(iabs(itype(i,molnum(i))))
           vbld_inv(i+nres)=dsc_inv(iabs(itype(i,molnum(i))))
+           else
+          vbld(i+nres)=vbldsc0_nucl(1,iabs(itype(i,molnum(i))))
+          vbld_inv(i+nres)=1.0/vbldsc0_nucl(1,iabs(itype(i,molnum(i))))
+           endif
 !          write (iout,*) "i",i," itype",itype(i,1),
 !     &      " dsc",dsc(itype(i,1))," vbld",vbld(i),vbld(i+nres)
         enddo