/[PAMELA software]/DarthVader/TrackerLevel2/src/F77/readB.f
ViewVC logotype

Diff of /DarthVader/TrackerLevel2/src/F77/readB.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.1 by mocchiut, Fri May 19 13:15:56 2006 UTC revision 1.5 by pam-fi, Thu Nov 2 15:31:06 2006 UTC
# Line 15  Line 15 
15    
16        subroutine readB(cpath)        subroutine readB(cpath)
17        include 'common_c2f.f'        include 'common_c2f.f'
18        character*80 cpath  
19        character*80 ppath        LOGICAL DEBUG
20          LOGICAL VERBOSE
21          LOGICAL WARNING
22          COMMON/DBG/DEBUG,VERBOSE,WARNING
23          SAVE/DBG/
24    
25          character*256 cpath
26          character*256 ppath
27    
28                
29        b_error=0        b_error=0
30          ppath=''
31    
32        la=80        la=256
33        do i=1,80        do i=1,256
34           if(cpath(i:i).eq.'/')la=i           if(cpath(i:i).eq.'/')la=i
35        enddo        enddo
36        ppath=cpath(1:la)        ppath=cpath(1:la)
# Line 30  Line 38 
38        b_path    = ppath        b_path    = ppath
39        b_pathlen = la        b_pathlen = la
40    
41  c     FUNZIONAAAAAAAAAAAAAAA!!!!!!!!!!!!!!!!!!        if(DEBUG)print*,'Field loaded: ',b_loaded
   
       if(b_debug.eq.1)print*,'Field loaded: ',b_loaded  
42        if(b_loaded.eq.0)then        if(b_loaded.eq.0)then
43    
44  c     call the subroutine which reads the maps of the measurements taken  c     call the subroutine which reads the maps of the measurements taken
# Line 79  c     outside the magnetic cavity Line 85  c     outside the magnetic cavity
85  c     implicit double precision (a-h,o-z)  c     implicit double precision (a-h,o-z)
86        include 'common_B.f'        include 'common_B.f'
87        include 'common_c2f.f'        include 'common_c2f.f'
88        character*80 path  
89          LOGICAL DEBUG
90          LOGICAL VERBOSE
91          LOGICAL WARNING
92          COMMON/DBG/DEBUG,VERBOSE,WARNING
93          SAVE/DBG/
94    
95          character*256 path
96                
97  C  C
98        REAL hmemor(10000000)        REAL hmemor(10000000)
# Line 137  c     print*,' ' Line 150  c     print*,' '
150        Bmap_file='measure_n3_290302.rz'        Bmap_file='measure_n3_290302.rz'
151                
152  c     opens magnetic field map first file  c     opens magnetic field map first file
153        if(b_debug.eq.1)print *                  c      if(b_debug.eq.1)print *                
154          if(VERBOSE)print *                
155       $     ,path(1:LNBLNK(path))//Bmap_file       $     ,path(1:LNBLNK(path))//Bmap_file
156        call HROPEN        call HROPEN
157       $     (lun_Bmap_file,'Bmap'       $     (lun_Bmap_file,'Bmap'
# Line 261  c     print*,' ' Line 275  c     print*,' '
275    
276        Bmap_file='measure_n4_110402_corrected.rz'        Bmap_file='measure_n4_110402_corrected.rz'
277  c     opens magnetic field map first file  c     opens magnetic field map first file
278        if(b_debug.eq.1)print *                   !,'Opening file: '  c      if(b_debug.eq.1)print *                   !,'Opening file: '
279          if(VERBOSE)print *                   !,'Opening file: '
280       $     ,path(1:LNBLNK(path))//Bmap_file       $     ,path(1:LNBLNK(path))//Bmap_file
281        call HROPEN        call HROPEN
282       $     (lun_Bmap_file,'Bmap'       $     (lun_Bmap_file,'Bmap'
# Line 381  c--------------------------------------- Line 396  c---------------------------------------
396   21   continue   21   continue
397                
398        b_error = 1        b_error = 1
399        if(b_debug.eq.1)  c      if(b_debug.eq.1)
400          if(DEBUG)
401       $     print*       $     print*
402       $     ,'read_B_inner: ERROR OPENING MAGNETIC FIELD MAP FILE: '       $     ,'read_B_inner: ERROR OPENING MAGNETIC FIELD MAP FILE: '
403       $     ,Bmap_file       $     ,Bmap_file
# Line 398  c--------------------------------------- Line 414  c---------------------------------------
414   22   continue   22   continue
415                
416        b_error = 1        b_error = 1
417        if(b_debug.eq.1)  c      if(b_debug.eq.1)
418          if(DEBUG)
419       $     print*,'read_B_inner: ERROR WHILE READING NTUPLE, at entry       $     print*,'read_B_inner: ERROR WHILE READING NTUPLE, at entry
420       $     : ',iev       $     : ',iev
421                
# Line 441  c     implicit double precision (a-h,o-z Line 458  c     implicit double precision (a-h,o-z
458        include 'common_B.f'        include 'common_B.f'
459        include 'common_c2f.f'        include 'common_c2f.f'
460  C  C
461        character*80 path        LOGICAL DEBUG
462          LOGICAL VERBOSE
463          LOGICAL WARNING
464          COMMON/DBG/DEBUG,VERBOSE,WARNING
465          SAVE/DBG/
466    
467          character*256 path
468    
469        REAL hmemor(10000000)        REAL hmemor(10000000)
470        integer Iquest(100)            integer Iquest(100)    
# Line 500  c     print*,' ' Line 523  c     print*,' '
523        Bmap_file='External_top_map_n4_150402.rz'        Bmap_file='External_top_map_n4_150402.rz'
524                
525  c     opens magnetic field map first file  c     opens magnetic field map first file
526        if(b_debug.eq.1)print *  c      if(b_debug.eq.1)print *
527          if(VERBOSE)print *
528       $     ,path(1:LNBLNK(path))//Bmap_file       $     ,path(1:LNBLNK(path))//Bmap_file
529        call HROPEN        call HROPEN
530       $     (lun_Bmap_file,'Bmap'       $     (lun_Bmap_file,'Bmap'
# Line 621  c--------------------------------------- Line 645  c---------------------------------------
645   21   continue   21   continue
646                
647        b_error = 1        b_error = 1
648        if(b_debug.eq.1)  c      if(b_debug.eq.1)
649          if(DEBUG)
650       $     print*       $     print*
651       $     ,'read_B_inner: ERROR OPENING MAGNETIC FIELD MAP FILE: '       $     ,'read_B_inner: ERROR OPENING MAGNETIC FIELD MAP FILE: '
652       $     ,Bmap_file       $     ,Bmap_file
# Line 638  c--------------------------------------- Line 663  c---------------------------------------
663   22   continue   22   continue
664                
665        b_error = 1        b_error = 1
666        if(b_debug.eq.1)  c      if(b_debug.eq.1)
667          if(DEBUG)
668       $     print*       $     print*
669       $     ,'read_B_inner: ERROR WHILE READING NTUPLE, at event       $     ,'read_B_inner: ERROR WHILE READING NTUPLE, at event
670       $     : ',iev       $     : ',iev
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.23