ground/src/read_B.for

*************************************************************************
*     
*     Subroutine read_B.f
*     
*     it calls the subroutines which read the magnetic field maps for
*     the PAMELA spectrometer
*     
*     needs:
*     - ./read_B_inner.f inner map reading subroutine
*     - ./read_B_outer.f outer map reading subroutine
*     
*     to be called before ./inter_B.f (interpolation subroutine)
*     
*************************************************************************

      function readb2maps(idata_dir)
c      subroutine read_B(idata_dir)
      CHARACTER*250 idata_dir
      integer lastnb2
c      integer lung
      parameter(NWPAWC=8500000)
      common/pawc/hmem(NWPAWC)
      save/pawc/
      common/quest/iquest(100)
      save/quest/
      call hlimit(nwpawc)
c
      lastnb2 = length(idata_dir)+1
c     call the subroutine which reads the maps of the measurements taken
c     inside the magnetic cavity
      call read_B_inner(idata_dir(1:lastnb2))

c     call the subroutine which reads the maps of the measurements taken
c     outside the magnetic cavity
cc      call read_B_outer(idata_dir(1:lastnb2))        

      return
      end

      include './read_B_inner.for'
c      include './read_B_outer.for' 

      
      integer function length(str)
c     return the string length without the blanks characters

      implicit integer (k-l)
      character *(*) str    
      
      lmax=len(str)      
       
c     search the last non blank character
      do i = 1, lmax
         if( str(i:i).eq.'.' )then
            length=i
            return
         endif
      enddo

      length=lmax

      return
      end     
      
1	*************************************************************************
2	*
3	* Subroutine read_B.f
4	*
5	* it calls the subroutines which read the magnetic field maps for
6	* the PAMELA spectrometer
7	*
8	* needs:
9	* - ./read_B_inner.f inner map reading subroutine
10	* - ./read_B_outer.f outer map reading subroutine
11	*
12	* to be called before ./inter_B.f (interpolation subroutine)
13	*
14	*************************************************************************
15
16	function readb2maps(idata_dir)
17	c subroutine read_B(idata_dir)
18	CHARACTER*250 idata_dir
19	integer lastnb2
20	c integer lung
21	parameter(NWPAWC=8500000)
22	common/pawc/hmem(NWPAWC)
23	save/pawc/
24	common/quest/iquest(100)
25	save/quest/
26	call hlimit(nwpawc)
27	c
28	lastnb2 = length(idata_dir)+1
29	c call the subroutine which reads the maps of the measurements taken
30	c inside the magnetic cavity
31	call read_B_inner(idata_dir(1:lastnb2))
32
33	c call the subroutine which reads the maps of the measurements taken
34	c outside the magnetic cavity
35	cc call read_B_outer(idata_dir(1:lastnb2))
36
37	return
38	end
39
40	include './read_B_inner.for'
41	c include './read_B_outer.for'
42
43
44	integer function length(str)
45	c return the string length without the blanks characters
46
47	implicit integer (k-l)
48	character () str
49
50	lmax=len(str)
51
52	c search the last non blank character
53	do i = 1, lmax
54	if( str(i:i).eq.'.' )then
55	length=i
56	return
57	endif
58	enddo
59
60	length=lmax
61
62	return
63	end
64