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	pam-fi	1.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