gpamela/gpfield/inter_B.F

*************************************************************************
*     
*     Subroutine inter_B.f (from tracker software analysis)
*     
*     it computes the magnetic field in a chosen point x,y,z inside or 
*     outside the magnetic cavity, using a trilinear interpolation of
*     B field measurements (read before by means of ./read_B.f)
*     if the point falls outside the interpolation volume, set the field to 0
*     
*     needs:
*     - ../common/common_B_inner.f common file for the inner magnetic field map
*     - ./inter_B_inner.f common file for the inner magnetic field map
*     
*     to be called after ./read_B.f (magnetic field map reading subroutine)
*     
*     input: coordinates in m
*     output: magnetic field in T
*     
*************************************************************************

      subroutine inter_B(x,y,z,res) !coordinates in m, magnetic field in T
      IMPLICIT DOUBLE PRECISION (A-H,O-Z)
#include "gpfield.inc"

c------------------------------------------------------------------------
c     
c     local variables
c     
c------------------------------------------------------------------------

      real*8 x,y,z              !point of interest
      real*8 res(3)             !interpolated B components: res = (Bx, By, Bz)


c------------------------------------------------------------------------
c     
c     set the field outside the interpolation volume to be 0
c     
c------------------------------------------------------------------------

      do ip=1,3
        res(ip)=0.
      enddo


c------------------------------------------------------------------------
c     
c     check if the point falls inside the interpolation volume
c     
c------------------------------------------------------------------------
      
      if((x.ge.edgexmin).and.(x.le.edgexmax)
     $     .and.(y.ge.edgeymin).and.(y.le.edgeymax)
     $     .and.(z.ge.edgezmin).and.(z.le.edgezmax)) then

        call inter_B_inner(x,y,z,res)
       
      endif

c     TODO: add inter_B_outer subroutine(s) for the external magnetic field map

      return
      end


1	cafagna	3.1	*************************************************************************
2			*
3			* Subroutine inter_B.f (from tracker software analysis)
4			*
5			* it computes the magnetic field in a chosen point x,y,z inside or
6			* outside the magnetic cavity, using a trilinear interpolation of
7			* B field measurements (read before by means of ./read_B.f)
8			* if the point falls outside the interpolation volume, set the field to 0
9			*
10			* needs:
11			* - ../common/common_B_inner.f common file for the inner magnetic field map
12			* - ./inter_B_inner.f common file for the inner magnetic field map
13			*
14			* to be called after ./read_B.f (magnetic field map reading subroutine)
15			*
16			* input: coordinates in m
17			* output: magnetic field in T
18			*
19			*************************************************************************
20
21			subroutine inter_B(x,y,z,res) !coordinates in m, magnetic field in T
22			IMPLICIT DOUBLE PRECISION (A-H,O-Z)
23			#include "gpfield.inc"
24
25			c------------------------------------------------------------------------
26			c
27			c local variables
28			c
29			c------------------------------------------------------------------------
30
31			real*8 x,y,z !point of interest
32			real*8 res(3) !interpolated B components: res = (Bx, By, Bz)
33
34
35			c------------------------------------------------------------------------
36			c
37			c set the field outside the interpolation volume to be 0
38			c
39			c------------------------------------------------------------------------
40
41			do ip=1,3
42			res(ip)=0.
43			enddo
44
45
46			c------------------------------------------------------------------------
47			c
48			c check if the point falls inside the interpolation volume
49			c
50			c------------------------------------------------------------------------
51
52			if((x.ge.edgexmin).and.(x.le.edgexmax)
53			$ .and.(y.ge.edgeymin).and.(y.le.edgeymax)
54			$ .and.(z.ge.edgezmin).and.(z.le.edgezmax)) then
55
56			call inter_B_inner(x,y,z,res)
57
58			endif
59
60			c TODO: add inter_B_outer subroutine(s) for the external magnetic field map
61
62			return
63			end
64
65