+ sint1t2n1=sint1t2n
+ sint1t2n=sint1t2n*sint1t2
+ sumth1tyb1=tschebyshev(1,nval,v11_chyb(1,j,itori,itori1),
+ & costheti12)
+ gradth1tyb1=-0.5d0*sinthet1*gradtschebyshev(0,nval-1,
+ & v11_chyb(1,j,itori,itori1),costheti12)
+c write (iout,*) "v11",(v11_chyb(k,j,itori,itori1),k=1,nval),
+c & " sumth1tyb1",sumth1tyb1," gradth1tyb1",gradth1tyb1
+ sumth2tyb1=tschebyshev(1,nval,v21_chyb(1,j,itori,itori1),
+ & costheti22)
+ gradth2tyb1=-0.5d0*sinthet2*gradtschebyshev(0,nval-1,
+ & v21_chyb(1,j,itori,itori1),costheti22)
+c write (iout,*) "v21",(v21_chyb(k,j,itori,itori1),k=1,nval),
+c & " sumth2tyb1",sumth2tyb1," gradth2tyb1",gradth2tyb1
+ sumth1tyb2=tschebyshev(1,nval,v12_chyb(1,j,itori,itori1),
+ & costheti12)
+ gradth1tyb2=-0.5d0*sinthet1*gradtschebyshev(0,nval-1,
+ & v12_chyb(1,j,itori,itori1),costheti12)
+c write (iout,*) "v12",(v12_chyb(k,j,itori,itori1),k=1,nval),
+c & " sumth1tyb2",sumth1tyb2," gradth1tyb2",gradth1tyb2
+ sumth2tyb2=tschebyshev(1,nval,v22_chyb(1,j,itori,itori1),
+ & costheti22)
+ gradth2tyb2=-0.5d0*sinthet2*gradtschebyshev(0,nval-1,
+ & v22_chyb(1,j,itori,itori1),costheti22)
+c write (iout,*) "v22",(v22_chyb(k,j,itori,itori1),k=1,nval),
+c & " sumth2tyb2",sumth2tyb2," gradth2tyb2",gradth2tyb2