--- DarthVader/TrackerLevel2/src/F77/grkuta.f 2006/06/01 09:03:09 1.2 +++ DarthVader/TrackerLevel2/src/F77/grkuta.f 2007/02/16 14:56:02 1.4 @@ -521,15 +521,17 @@ real*8 vv(3),ff(3) !inter_B.f works in double precision - do i=1,3 - vv(i)=v(i)/100. !inter_B.f works in meters - enddo -c inter_B: coordinates in m, B field in Tesla - call inter_B(vv(1),vv(2),vv(3),ff) - do i=1,3 !change back the field in kGauss - f(i)=ff(i)*10. - enddo - + do i=1,3 + vv(i)=v(i)/100. !inter_B.f works in meters + enddo +c inter_B: coordinates in m, B field in Tesla +c$$$ print*,'GUFLD: v ',v + call inter_B(vv(1),vv(2),vv(3),ff) + do i=1,3 !change back the field in kGauss + f(i)=ff(i)*10. + enddo +c$$$ print*,'GUFLD: b ',f + return end