/[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.2 by pam-fi, Thu Jun 1 09:03:09 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          call inter_B(vv(1),vv(2),vv(3),ff)  c$$$      print*,'GUFLD: v ',v
529          do i=1,3                !change back the field in kGauss        call inter_B(vv(1),vv(2),vv(3),ff)
530            f(i)=ff(i)*10.        do i=1,3                  !change back the field in kGauss
531          enddo           f(i)=ff(i)*10.
532          enddo
533    c$$$      print*,'GUFLD: b ',f
534          
535        return        return
536        end        end
537    
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.23