12 if(b%sys .eq.
'xyz')
then 13 write(*, *)
'mag_field: x=',b%x,
' y=',b%y,
15 elseif(b%sys .eq.
'hva')
then 16 write(*,*)
' mag_field: horizontal=',b%h,
17 *
' vertical=', b%v,
' angle=', b%a,
19 elseif(b%sys .eq.
'ned')
then 20 write(*,*)
' mag_field: north=',b%n,
21 *
' east=',b%e,
' down=',b%d
23 write(msg,*)
' mag coord type=',b%sys,
24 *
' unkown to cprint_magfield',
25 *
' componen=', b%x, b%y, b%z
29 if(b%sys .eq.
'xyz')
then 30 write(*, *)
'mag_field: x=',b%x,
' y=',b%y,
32 elseif(b%sys .eq.
'hva')
then 33 write(*,*)
' mag_field: horizontal=',b%x,
34 *
' vertical=', b%y,
' angle=', b%z,
36 elseif(b%sys .eq.
'ned')
then 37 write(*,*)
' mag_field: north=',b%x,
38 *
' east=',b%y,
' down=',b%z
40 write(msg,*)
' mag coord type=',b%sys,
41 *
' unkown to cprint_magfield',
42 *
' componen=', b%x, b%y, b%z
46 write(*,*)
' Unit is T. Multiply 10**4 to Gausss' subroutine cerrormsg(msg, needrtn)