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 |
|
|
39 |
|
|
40 |
c FUNZIONAAAAAAAAAAAAAAA!!!!!!!!!!!!!!!!!! |
c FUNZIONAAAAAAAAAAAAAAA!!!!!!!!!!!!!!!!!! |
41 |
|
|
42 |
if(b_debug.eq.1)print*,'Field loaded: ',b_loaded |
c if(b_debug.eq.1)print*,'Field loaded: ',b_loaded |
43 |
|
if(DEBUG)print*,'Field loaded: ',b_loaded |
44 |
if(b_loaded.eq.0)then |
if(b_loaded.eq.0)then |
45 |
|
|
46 |
c call the subroutine which reads the maps of the measurements taken |
c call the subroutine which reads the maps of the measurements taken |
87 |
c implicit double precision (a-h,o-z) |
c implicit double precision (a-h,o-z) |
88 |
include 'common_B.f' |
include 'common_B.f' |
89 |
include 'common_c2f.f' |
include 'common_c2f.f' |
90 |
|
|
91 |
|
LOGICAL DEBUG |
92 |
|
LOGICAL VERBOSE |
93 |
|
LOGICAL WARNING |
94 |
|
COMMON/DBG/DEBUG,VERBOSE,WARNING |
95 |
|
SAVE/DBG/ |
96 |
|
|
97 |
character*80 path |
character*80 path |
98 |
|
|
99 |
C |
C |
152 |
Bmap_file='measure_n3_290302.rz' |
Bmap_file='measure_n3_290302.rz' |
153 |
|
|
154 |
c opens magnetic field map first file |
c opens magnetic field map first file |
155 |
if(b_debug.eq.1)print * |
c if(b_debug.eq.1)print * |
156 |
|
if(VERBOSE)print * |
157 |
$ ,path(1:LNBLNK(path))//Bmap_file |
$ ,path(1:LNBLNK(path))//Bmap_file |
158 |
call HROPEN |
call HROPEN |
159 |
$ (lun_Bmap_file,'Bmap' |
$ (lun_Bmap_file,'Bmap' |
277 |
|
|
278 |
Bmap_file='measure_n4_110402_corrected.rz' |
Bmap_file='measure_n4_110402_corrected.rz' |
279 |
c opens magnetic field map first file |
c opens magnetic field map first file |
280 |
if(b_debug.eq.1)print * !,'Opening file: ' |
c if(b_debug.eq.1)print * !,'Opening file: ' |
281 |
|
if(VERBOSE)print * !,'Opening file: ' |
282 |
$ ,path(1:LNBLNK(path))//Bmap_file |
$ ,path(1:LNBLNK(path))//Bmap_file |
283 |
call HROPEN |
call HROPEN |
284 |
$ (lun_Bmap_file,'Bmap' |
$ (lun_Bmap_file,'Bmap' |
398 |
21 continue |
21 continue |
399 |
|
|
400 |
b_error = 1 |
b_error = 1 |
401 |
if(b_debug.eq.1) |
c if(b_debug.eq.1) |
402 |
|
if(DEBUG) |
403 |
$ print* |
$ print* |
404 |
$ ,'read_B_inner: ERROR OPENING MAGNETIC FIELD MAP FILE: ' |
$ ,'read_B_inner: ERROR OPENING MAGNETIC FIELD MAP FILE: ' |
405 |
$ ,Bmap_file |
$ ,Bmap_file |
416 |
22 continue |
22 continue |
417 |
|
|
418 |
b_error = 1 |
b_error = 1 |
419 |
if(b_debug.eq.1) |
c if(b_debug.eq.1) |
420 |
|
if(DEBUG) |
421 |
$ print*,'read_B_inner: ERROR WHILE READING NTUPLE, at entry |
$ print*,'read_B_inner: ERROR WHILE READING NTUPLE, at entry |
422 |
$ : ',iev |
$ : ',iev |
423 |
|
|
460 |
include 'common_B.f' |
include 'common_B.f' |
461 |
include 'common_c2f.f' |
include 'common_c2f.f' |
462 |
C |
C |
463 |
|
LOGICAL DEBUG |
464 |
|
LOGICAL VERBOSE |
465 |
|
LOGICAL WARNING |
466 |
|
COMMON/DBG/DEBUG,VERBOSE,WARNING |
467 |
|
SAVE/DBG/ |
468 |
|
|
469 |
character*80 path |
character*80 path |
470 |
|
|
471 |
REAL hmemor(10000000) |
REAL hmemor(10000000) |
525 |
Bmap_file='External_top_map_n4_150402.rz' |
Bmap_file='External_top_map_n4_150402.rz' |
526 |
|
|
527 |
c opens magnetic field map first file |
c opens magnetic field map first file |
528 |
if(b_debug.eq.1)print * |
c if(b_debug.eq.1)print * |
529 |
|
if(VERBOSE)print * |
530 |
$ ,path(1:LNBLNK(path))//Bmap_file |
$ ,path(1:LNBLNK(path))//Bmap_file |
531 |
call HROPEN |
call HROPEN |
532 |
$ (lun_Bmap_file,'Bmap' |
$ (lun_Bmap_file,'Bmap' |
647 |
21 continue |
21 continue |
648 |
|
|
649 |
b_error = 1 |
b_error = 1 |
650 |
if(b_debug.eq.1) |
c if(b_debug.eq.1) |
651 |
|
if(DEBUG) |
652 |
$ print* |
$ print* |
653 |
$ ,'read_B_inner: ERROR OPENING MAGNETIC FIELD MAP FILE: ' |
$ ,'read_B_inner: ERROR OPENING MAGNETIC FIELD MAP FILE: ' |
654 |
$ ,Bmap_file |
$ ,Bmap_file |
665 |
22 continue |
22 continue |
666 |
|
|
667 |
b_error = 1 |
b_error = 1 |
668 |
if(b_debug.eq.1) |
c if(b_debug.eq.1) |
669 |
|
if(DEBUG) |
670 |
$ print* |
$ print* |
671 |
$ ,'read_B_inner: ERROR WHILE READING NTUPLE, at event |
$ ,'read_B_inner: ERROR WHILE READING NTUPLE, at event |
672 |
$ : ',iev |
$ : ',iev |