source/magnet/main_B.f

      program main_B


      parameter (NWPAWC=6500000)
      common/PAWC/HMEM(NWPAWC)


      real*8 vv(3),ff(3)

      call HLIMIT(NWPAWC)


      call read_B
      print*,' '
      print*,'FINITA LETTURA CAMPO MAGNETICO'
      print*,' '

      print*,' '
      print*,'x,y,z? (cm)'
      do i=1,3
         read(5,*) vv(i)         !in cm
      enddo
      print*,vv(1),vv(2),vv(3)
      print*,' '
      
      do i=1,3                  !trasformo le coordinate in m per
        vv(i)=vv(i)*0.01        !inter_B.f
      enddo
      
      call inter_B(vv(1),vv(2),vv(3),ff) !coordinate in m, campo in tesla

      print*,' '
      print*,' point (m)',vv
      print*,' B (T)',ff
      print*,' '

      print*,' '
      print*,'   --- THE END ---'
      print*,' '
      print*,' '
      
      stop
      end

      include './read_B.f' !routine che legge il campo misurato
      include './inter_B.f'     !routine che calcola il campo
1	program main_B
2
3
4	parameter (NWPAWC=6500000)
5	common/PAWC/HMEM(NWPAWC)
6
7
8	real*8 vv(3),ff(3)
9
10	call HLIMIT(NWPAWC)
11
12
13
14	call read_B
15	print*,' '
16	print*,'FINITA LETTURA CAMPO MAGNETICO'
17	print*,' '
18
19	print*,' '
20	print*,'x,y,z? (cm)'
21	do i=1,3
22	read(5,*) vv(i) !in cm
23	enddo
24	print*,vv(1),vv(2),vv(3)
25	print*,' '
26
27	do i=1,3 !trasformo le coordinate in m per
28	vv(i)=vv(i)*0.01 !inter_B.f
29	enddo
30
31	call inter_B(vv(1),vv(2),vv(3),ff) !coordinate in m, campo in tesla
32
33	print*,' '
34	print*,' point (m)',vv
35	print*,' B (T)',ff
36	print*,' '
37
38	print*,' '
39	print*,' --- THE END ---'
40	print*,' '
41	print*,' '
42
43	stop
44	end
45
46	include './read_B.f' !routine che legge il campo misurato
47	include './inter_B.f' !routine che calcola il campo