19 |
include 'commontracker.f' !tracker general common |
include 'commontracker.f' !tracker general common |
20 |
include 'common_mini_2.f' !common for the tracking procedure |
include 'common_mini_2.f' !common for the tracking procedure |
21 |
|
|
22 |
logical DEBUG |
c logical DEBUG |
23 |
common/dbg/DEBUG |
c common/dbg/DEBUG |
24 |
|
|
25 |
parameter (inf=1.e8) !just a huge number... |
parameter (inf=1.e8) !just a huge number... |
26 |
c------------------------------------------------------------------------ |
c------------------------------------------------------------------------ |
454 |
CALL GRKUTA(CHARGE,STEP,VECT,VOUT) |
CALL GRKUTA(CHARGE,STEP,VECT,VOUT) |
455 |
IF(VOUT(3).GT.VECT(3)) THEN |
IF(VOUT(3).GT.VECT(3)) THEN |
456 |
IFAIL=1 |
IFAIL=1 |
457 |
PRINT *,'posxy (grkuta): WARNING ===> backward track!!' |
if(WARNING) |
458 |
print*,'charge',charge |
$ PRINT *,'posxy (grkuta): WARNING ===> backward track!!' |
459 |
print*,'vect',vect |
if(WARNING)print*,'charge',charge |
460 |
print*,'vout',vout |
if(WARNING)print*,'vect',vect |
461 |
print*,'step',step |
if(WARNING)print*,'vout',vout |
462 |
|
if(WARNING)print*,'step',step |
463 |
RETURN |
RETURN |
464 |
ENDIF |
ENDIF |
465 |
Z=VOUT(3) |
Z=VOUT(3) |