+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=0,nkcctyp
+ do j=0,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) v11_chyb(l,k,j,i)
+ enddo
+ do l=1,nterm_kcc_Tb(j,i)
+ read (itorkcc,*,end=121,err=121) v21_chyb(l,k,j,i)
+ enddo
+ do l=1,nterm_kcc_Tb(j,i)
+ read (itorkcc,*,end=121,err=121) v12_chyb(l,k,j,i)
+ enddo
+ do l=1,nterm_kcc_Tb(j,i)
+ read (itorkcc,*,end=121,err=121) v22_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
+ if (lprint) then
+c Print valence-torsional parameters
+ write (iout,'(a)')
+ & "Parameters of the valence-torsional potentials"
+ do i=0,nkcctyp
+ do j=0,nkcctyp
+ write (iout,'(3a)') "Type ",toronelet(i),toronelet(j)
+ write (iout,'(2a20,a15)') "v_kcc","v1_chyb","v2_chyb"
+ do k=1,nterm_kcc(j,i)
+ write (iout,'(i5,f15.10,i5,2f15.10)')
+ & k,v1_kcc(k,j,i),1,v11_chyb(1,k,j,i),v21_chyb(1,k,j,i)
+ do l=2,nterm_kcc_Tb(j,i)
+ write (iout,'(20x,i5,2f15.10)')
+ & l,v11_chyb(l,k,j,i),v21_chyb(l,k,j,i)
+ enddo
+ write (iout,'(i5,f15.10,i5,2f15.10)')
+ & k,v2_kcc(k,j,i),1,v12_chyb(1,k,j,i),v22_chyb(1,k,j,i)
+ do l=2,nterm_kcc_Tb(j,i)
+ write (iout,'(20x,i5,2f15.10)')
+ & l,v12_chyb(l,k,j,i),v22_chyb(l,k,j,i)
+ enddo
+ write (iout,'(a)')
+ enddo
+ enddo
+ enddo
+ endif
+C here will be the apropriate recalibrating for D-aminoacid
+C read (ithetkcc,*,end=121,err=121) nkcctyp
+ do i=0,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
+ if (lprint) then
+ write (iout,'(a)')
+ & "Parameters of the valence-only potentials"
+ do i=0,nkcctyp
+ write (iout,'(2a)') "Type ",toronelet(i)
+ do k=1,nbend_kcc_Tb(i)
+ write(iout,'(i5,f15.10)') k,v1bend_chyb(k,i)
+ enddo
+ enddo
+ endif