X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2Fsrc_MD-M%2Fparmread.F;h=e156c9782347f2f34548dcbf734020721519b5a3;hb=aedd2f72960d9cec69fcc7f0b6fd6555c7e22312;hp=c881425f9d6fd9d765060de2f24f5f7fbabd7812;hpb=2095e643aafa2fb192e066e69e5ca3b1a59c2727;p=unres.git

diff --git a/source/unres/src_MD-M/parmread.F b/source/unres/src_MD-M/parmread.F
index c881425..e156c97 100644
--- a/source/unres/src_MD-M/parmread.F
+++ b/source/unres/src_MD-M/parmread.F
@@ -775,6 +775,31 @@ C Martix of D parameters for two dimesional fourier series
       enddo
       endif
 #endif
+C read Czybyshev torsional parameters
+        read (itorkcc,*,end=121,err=121) nkcctyp
+        read (itorkcc,*,end=121,err=121) (itortyp_kcc(i),i=1,ntyp)
+        do i=-ntyp,-1
+        itortyp_kcc(i)=-itortyp_kcc(-i)
+        enddo
+        do i=1,nkcctyp
+         do j=1,nkcctyp
+C first we read the cos and sin gamma parameters
+          read (itorkcc,*,end=121,err=121) nterm_kcc(j,i)
+           do k=1,nterm_kcc(j,i)
+             read (itorkcc,*,end=121,err=121) 
+     &        v1_kcc(k,j,i),v2_kcc(k,j,i)
+           enddo
+C now Czybyshev parameters
+           read (itorkcc,*,end=121,err=121) nterm_kcc_Tb(j,i)
+           do k=1,nterm_kcc_Tb(j,i)
+             read (itorkcc,*,end=121,err=121) 
+     &        v1_chyb(k,j,i),v2_chyb(k,j,i)
+           enddo
+          enddo
+         enddo
+C here will be the apropriate recalibrating for D-aminoacid
+
+
 C Read of Side-chain backbone correlation parameters
 C Modified 11 May 2012 by Adasko
 CCC
@@ -1130,11 +1155,12 @@ C---------------------- GB or BP potential -----------------------------
 C now we start reading lipid
       do i=1,ntyp
        read (isidep,*,end=1161,err=1161)(epslip(i,j),j=i,ntyp)
-       print *,"WARNING!!"
-       do j=1,ntyp
-       epslip(i,j)=epslip(i,j)+0.05d0
-       enddo
+C       print *,"WARNING!!"
+C       do j=1,ntyp
+C       epslip(i,j)=epslip(i,j)+0.05d0
+C       enddo
       enddo
+      write(iout,*) epslip(1,1),"OK?"
 C For the GB potential convert sigma'**2 into chi'
       if (ipot.eq.4) then
 	do i=1,ntyp
@@ -1375,6 +1401,8 @@ C      endif
   118 write (iout,*) "Error reading SCp interaction parameters."
       goto 999
   119 write (iout,*) "Error reading SCCOR parameters"
+      go to 999
+  121 write (iout,*) "Error in Czybyshev parameters"
   999 continue
 #ifdef MPI
       call MPI_Finalize(Ierror)