call reada(controlcard,'SCALSCP',scalscp,1.0d0)
call reada(controlcard,'CUTOFF',cutoff_corr,7.0d0)
call reada(controlcard,'DELT_CORR',delt_corr,0.5d0)
+ call reada(controlcard,'LIPSCALE',lipscale,1.0D0)
+c write (iout,*) "lipscale",lipscale
r0_corr=cutoff_corr-delt_corr
write (iout,*) "iparm",iparm," myparm",myparm
V2SS = 7.61d0
V3SS = 13.7d0
- do i=1,maxres-1
- do j=i+1,maxres
+ do i=1,max_cyst-1
+ do j=i+1,max_cyst
dyn_ssbond_ij(i,j)=1.0d300
enddo
enddo
EEold(2,2,-i)=-b(10,i)+b(11,i)
EEold(2,1,-i)=-b(12,i)+b(13,i)
EEold(1,2,-i)=-b(12,i)-b(13,i)
-c write(iout,*) "TU DOCHODZE"
-c print *,"JESTEM"
+ write(iout,*) "TU DOCHODZE"
+ print *,"JESTEM"
c ee(1,1,i)=1.0d0
c ee(2,2,i)=1.0d0
c ee(2,1,i)=0.0d0