| 15 |
|
|
| 16 |
subroutine readB(cpath) |
subroutine readB(cpath) |
| 17 |
include 'common_c2f.f' |
include 'common_c2f.f' |
| 18 |
|
|
| 19 |
|
LOGICAL DEBUG |
| 20 |
|
LOGICAL VERBOSE |
| 21 |
|
LOGICAL WARNING |
| 22 |
|
COMMON/DBG/DEBUG,VERBOSE,WARNING |
| 23 |
|
SAVE/DBG/ |
| 24 |
|
|
| 25 |
character*80 cpath |
character*80 cpath |
| 26 |
character*80 ppath |
character*80 ppath |
| 27 |
|
|
| 28 |
|
|
| 29 |
b_error=0 |
b_error=0 |
| 30 |
|
ppath='' |
| 31 |
|
|
| 32 |
la=80 |
la=80 |
| 33 |
do i=1,80 |
do i=1,80 |
| 40 |
|
|
| 41 |
c FUNZIONAAAAAAAAAAAAAAA!!!!!!!!!!!!!!!!!! |
c FUNZIONAAAAAAAAAAAAAAA!!!!!!!!!!!!!!!!!! |
| 42 |
|
|
| 43 |
if(b_debug.eq.1)print*,'Field loaded: ',b_loaded |
c if(b_debug.eq.1)print*,'Field loaded: ',b_loaded |
| 44 |
|
if(DEBUG)print*,'Field loaded: ',b_loaded |
| 45 |
if(b_loaded.eq.0)then |
if(b_loaded.eq.0)then |
| 46 |
|
|
| 47 |
c call the subroutine which reads the maps of the measurements taken |
c call the subroutine which reads the maps of the measurements taken |
| 88 |
c implicit double precision (a-h,o-z) |
c implicit double precision (a-h,o-z) |
| 89 |
include 'common_B.f' |
include 'common_B.f' |
| 90 |
include 'common_c2f.f' |
include 'common_c2f.f' |
| 91 |
|
|
| 92 |
|
LOGICAL DEBUG |
| 93 |
|
LOGICAL VERBOSE |
| 94 |
|
LOGICAL WARNING |
| 95 |
|
COMMON/DBG/DEBUG,VERBOSE,WARNING |
| 96 |
|
SAVE/DBG/ |
| 97 |
|
|
| 98 |
character*80 path |
character*80 path |
| 99 |
|
|
| 100 |
C |
C |
| 153 |
Bmap_file='measure_n3_290302.rz' |
Bmap_file='measure_n3_290302.rz' |
| 154 |
|
|
| 155 |
c opens magnetic field map first file |
c opens magnetic field map first file |
| 156 |
if(b_debug.eq.1)print * |
c if(b_debug.eq.1)print * |
| 157 |
|
if(VERBOSE)print * |
| 158 |
$ ,path(1:LNBLNK(path))//Bmap_file |
$ ,path(1:LNBLNK(path))//Bmap_file |
| 159 |
call HROPEN |
call HROPEN |
| 160 |
$ (lun_Bmap_file,'Bmap' |
$ (lun_Bmap_file,'Bmap' |
| 278 |
|
|
| 279 |
Bmap_file='measure_n4_110402_corrected.rz' |
Bmap_file='measure_n4_110402_corrected.rz' |
| 280 |
c opens magnetic field map first file |
c opens magnetic field map first file |
| 281 |
if(b_debug.eq.1)print * !,'Opening file: ' |
c if(b_debug.eq.1)print * !,'Opening file: ' |
| 282 |
|
if(VERBOSE)print * !,'Opening file: ' |
| 283 |
$ ,path(1:LNBLNK(path))//Bmap_file |
$ ,path(1:LNBLNK(path))//Bmap_file |
| 284 |
call HROPEN |
call HROPEN |
| 285 |
$ (lun_Bmap_file,'Bmap' |
$ (lun_Bmap_file,'Bmap' |
| 399 |
21 continue |
21 continue |
| 400 |
|
|
| 401 |
b_error = 1 |
b_error = 1 |
| 402 |
if(b_debug.eq.1) |
c if(b_debug.eq.1) |
| 403 |
|
if(DEBUG) |
| 404 |
$ print* |
$ print* |
| 405 |
$ ,'read_B_inner: ERROR OPENING MAGNETIC FIELD MAP FILE: ' |
$ ,'read_B_inner: ERROR OPENING MAGNETIC FIELD MAP FILE: ' |
| 406 |
$ ,Bmap_file |
$ ,Bmap_file |
| 417 |
22 continue |
22 continue |
| 418 |
|
|
| 419 |
b_error = 1 |
b_error = 1 |
| 420 |
if(b_debug.eq.1) |
c if(b_debug.eq.1) |
| 421 |
|
if(DEBUG) |
| 422 |
$ print*,'read_B_inner: ERROR WHILE READING NTUPLE, at entry |
$ print*,'read_B_inner: ERROR WHILE READING NTUPLE, at entry |
| 423 |
$ : ',iev |
$ : ',iev |
| 424 |
|
|
| 461 |
include 'common_B.f' |
include 'common_B.f' |
| 462 |
include 'common_c2f.f' |
include 'common_c2f.f' |
| 463 |
C |
C |
| 464 |
|
LOGICAL DEBUG |
| 465 |
|
LOGICAL VERBOSE |
| 466 |
|
LOGICAL WARNING |
| 467 |
|
COMMON/DBG/DEBUG,VERBOSE,WARNING |
| 468 |
|
SAVE/DBG/ |
| 469 |
|
|
| 470 |
character*80 path |
character*80 path |
| 471 |
|
|
| 472 |
REAL hmemor(10000000) |
REAL hmemor(10000000) |
| 526 |
Bmap_file='External_top_map_n4_150402.rz' |
Bmap_file='External_top_map_n4_150402.rz' |
| 527 |
|
|
| 528 |
c opens magnetic field map first file |
c opens magnetic field map first file |
| 529 |
if(b_debug.eq.1)print * |
c if(b_debug.eq.1)print * |
| 530 |
|
if(VERBOSE)print * |
| 531 |
$ ,path(1:LNBLNK(path))//Bmap_file |
$ ,path(1:LNBLNK(path))//Bmap_file |
| 532 |
call HROPEN |
call HROPEN |
| 533 |
$ (lun_Bmap_file,'Bmap' |
$ (lun_Bmap_file,'Bmap' |
| 648 |
21 continue |
21 continue |
| 649 |
|
|
| 650 |
b_error = 1 |
b_error = 1 |
| 651 |
if(b_debug.eq.1) |
c if(b_debug.eq.1) |
| 652 |
|
if(DEBUG) |
| 653 |
$ print* |
$ print* |
| 654 |
$ ,'read_B_inner: ERROR OPENING MAGNETIC FIELD MAP FILE: ' |
$ ,'read_B_inner: ERROR OPENING MAGNETIC FIELD MAP FILE: ' |
| 655 |
$ ,Bmap_file |
$ ,Bmap_file |
| 666 |
22 continue |
22 continue |
| 667 |
|
|
| 668 |
b_error = 1 |
b_error = 1 |
| 669 |
if(b_debug.eq.1) |
c if(b_debug.eq.1) |
| 670 |
|
if(DEBUG) |
| 671 |
$ print* |
$ print* |
| 672 |
$ ,'read_B_inner: ERROR WHILE READING NTUPLE, at event |
$ ,'read_B_inner: ERROR WHILE READING NTUPLE, at event |
| 673 |
$ : ',iev |
$ : ',iev |