/[PAMELA software]/DarthVader/TrackerLevel2/src/F77/grkuta.f
ViewVC logotype

Diff of /DarthVader/TrackerLevel2/src/F77/grkuta.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.3 by pam-fi, Wed Nov 8 16:42:28 2006 UTC revision 1.5 by pam-fi, Fri Apr 27 11:19:51 2007 UTC
# Line 500  c      print*,'- ',VOUT(3),z,VOUT(1),x,V Line 500  c      print*,'- ',VOUT(3),z,VOUT(1),x,V
500       $     +(VOUT(IY)-VECT(IY))**2       $     +(VOUT(IY)-VECT(IY))**2
501       $     +(VOUT(IZ)-VECT(IZ))**2       $     +(VOUT(IZ)-VECT(IZ))**2
502       $     )       $     )
503        print*,'WARNING: GRKUTA2 --> '  c      print*,'WARNING: GRKUTA2 --> '
504       $     ,'helix :-( ... length evaluated with straight line'  c     $     ,'helix :-( ... length evaluated with straight line'
505    
506  *  *
507    999 END    999 END
# Line 521  c      print*,'- ',VOUT(3),z,VOUT(1),x,V Line 521  c      print*,'- ',VOUT(3),z,VOUT(1),x,V
521        real*8 vv(3),ff(3)        !inter_B.f works in double precision        real*8 vv(3),ff(3)        !inter_B.f works in double precision
522    
523    
524          do i=1,3        do i=1,3
525            vv(i)=v(i)/100.       !inter_B.f works in meters           vv(i)=v(i)/100.        !inter_B.f works in meters
526          enddo        enddo
527  c       inter_B: coordinates in m, B field in Tesla  c     inter_B: coordinates in m, B field in Tesla
528  c        print*,'GUFLD: v ',v  c$$$      print*,'GUFLD: v ',v
529          call inter_B(vv(1),vv(2),vv(3),ff)        call inter_B(vv(1),vv(2),vv(3),ff)
530          do i=1,3                !change back the field in kGauss        do i=1,3                  !change back the field in kGauss
531            f(i)=ff(i)*10.           f(i)=ff(i)*10.
532          enddo        enddo
533  c        print*,'GUFLD: b ',f  c$$$      print*,'GUFLD: b ',f
534  c        print*,vv,ff        
   
535        return        return
536        end        end
537    
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.23