+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),nterm_kcc_Tb(j,i)
+C read (itorkcc,*,end=121,err=121) nterm_kcc_Tb(j,i)
+ do k=1,nterm_kcc(j,i)
+ do l=1,nterm_kcc_Tb(j,i)
+ read (itorkcc,*,end=121,err=121)
+ & v1_chyb(l,k,j,i)
+ enddo
+ do l=1,nterm_kcc_Tb(j,i)
+ read (itorkcc,*,end=121,err=121)
+ & v2_chyb(l,k,j,i)
+ enddo
+ read (itorkcc,*,end=121,err=121)
+ & v1_kcc(k,j,i)
+ read (itorkcc,*,end=121,err=121)
+ & v2_kcc(k,j,i)
+ enddo
+ enddo
+ enddo
+C here will be the apropriate recalibrating for D-aminoacid
+C read (ithetkcc,*,end=121,err=121) nkcctyp
+ do i=1,nkcctyp
+ read (ithetkcc,*,end=121,err=121) nbend_kcc_Tb(i)
+ do j=1,nbend_kcc_Tb(i)
+ read (ithetkcc,*,end=121,err=121)
+ & v1bend_chyb(j,i)
+ enddo
+ enddo