read in parameters for CME shock ellipsoid in Kwon's model
=======if using data or observation at 1au to calculate Va and Vs======= bm densw,tempsw could be replaced with those from observation
=======if using data at rc1 to calculate Va and Vs directly===== magnetosonic speed (Vms): Vms2 = 0.5 * {Va2 + Vs2 + [(Va2 + Vs2)2 - 4 * Va2 * Vs2 * (cos (theta))2]0.5} ==================================================== ----------choose smaller one-----------------------------------
if (tauc2_0Arguments
None
Common Blocks
read_param
(subroutine)
cofm
(subroutine)
solarwind1
(function)
f0mod
(subroutine)
split
(subroutine)
vfunc
(subroutine)
solarwind
(subroutine)
drvbmag
(subroutine)
inorout
(subroutine)
walk3d
(subroutine)
sim3d
(program)
walk3d
(subroutine)
sim3d_em
(program)
walk3d
(subroutine)
sksurface
(subroutine)
shockfront
(program)
">
common /bpark/
Type
Attributes
Name
Initial
real(kind=real64)
::
b1au
real(kind=real64)
::
densw0
real(kind=real64)
::
facip
real(kind=real64)
::
k4ok2
real(kind=real64)
::
k6ok2
real(kind=real64)
::
omega
real(kind=real64)
::
vom
real(kind=real64)
::
vsw
Type
Attributes
Name
Initial
real(kind=real64)
::
acsk
integer
::
nsk
real(kind=real64)
::
pexsk(0:144,8)
real(kind=real64)
::
pska(20,8)
real(kind=real64)
::
tska(20)
real(kind=real64)
::
vskf0
read_dmumuDG
(subroutine)
read_dmumuAH
(subroutine)
read_dmumuBK
(subroutine)
preparefb
(subroutine)
read_param
(subroutine)
read_dxx
(subroutine)
prepareptcl
(subroutine)
sim3d
(program)
sim3d_em
(program)
shockfront
(program)
">
common /dir/
Type
Attributes
Name
Initial
character(len=256)
::
dir
prepareptcl
(subroutine)
split
(subroutine)
inorout
(subroutine)
walk3d
(subroutine)
sim3d
(program)
walk3d
(subroutine)
sim3d_em
(program)
walk3d
(subroutine)
sksurface
(subroutine)
shockfront
(program)
">
common /tmprm/
Type
Attributes
Name
Initial
real(kind=real64)
::
dmapjul
real(kind=real64)
::
t0org
real(kind=real64)
::
tcme0
real(kind=real64)
::
tdl
real(kind=real64)
::
te