COSMOS v7.655  COSMOSv7655
(AirShowerMC)
cprintMagF.f
Go to the documentation of this file.
1  subroutine cprintmagf(b)
2 ! print magnetic field for debug purpose
3 ! b: /magfield/ input.
4 !
5  implicit none
6 
7 #include "Zmagfield.h"
8  type(magfield)::b
9  character*100 msg
10 #ifdef UNIONMAP
11 !
12  if(b%sys .eq. 'xyz') then
13  write(*, *) 'mag_field: x=',b%x, ' y=',b%y,
14  * ' z=',b%z
15  elseif(b%sys .eq. 'hva') then
16  write(*,*) ' mag_field: horizontal=',b%h,
17  * ' vertical=', b%v, ' angle=', b%a,
18  * ' deg'
19  elseif(b%sys .eq. 'ned') then
20  write(*,*) ' mag_field: north=',b%n,
21  * ' east=',b%e, ' down=',b%d
22  else
23  write(msg,*) ' mag coord type=',b%sys,
24  * ' unkown to cprint_magfield',
25  * ' componen=', b%x, b%y, b%z
26  call cerrormsg(msg, 0)
27  endif
28 #else
29  if(b%sys .eq. 'xyz') then
30  write(*, *) 'mag_field: x=',b%x, ' y=',b%y,
31  * ' z=',b%z
32  elseif(b%sys .eq. 'hva') then
33  write(*,*) ' mag_field: horizontal=',b%x,
34  * ' vertical=', b%y, ' angle=', b%z,
35  * ' deg'
36  elseif(b%sys .eq. 'ned') then
37  write(*,*) ' mag_field: north=',b%x,
38  * ' east=',b%y, ' down=',b%z
39  else
40  write(msg,*) ' mag coord type=',b%sys,
41  * ' unkown to cprint_magfield',
42  * ' componen=', b%x, b%y, b%z
43  call cerrormsg(msg, 0)
44  endif
45 #endif
46  write(*,*) ' Unit is T. Multiply 10**4 to Gausss'
47  end
subroutine cerrormsg(msg, needrtn)
Definition: cerrorMsg.f:4
subroutine cprintmagf(b)
Definition: cprintMagF.f:2