CMR Header

CCMC Metadata Registry (CMR)

Menu Page

Go to:

CMR: View Simulation Run Info

CORHEL-CME/CME/1.0/ROR/Weihang_Zhang_042424_SH_2

Run Description

CCMC ROR run for CORHEL-CME/CME 1.0 with RORID Weihang_Zhang_042424_SH_2

Simulation Time:

Simulation Start Time: 2009-02-13T00:00:00Z
Simulation End Time: 2009-02-13T00:00:00Z
Time Step in second:
Simulation Time Description:
Simulation Time Caveats:

Temporal Dependence (whether the outputs of the simulation code are time-dependent?): true

Regions:

Heliosphere.Inner
Sun.Corona

Simulation Domain :

Contacts :

Weihang.Zhang, ModelUser

List of Run Input Parameter Grouping and their associated properties :

Name

Description

Caveats

Properties

1. Magnetic Map Selection

A set of parameters for selecting the date, time, and magnetic field map in the CORHEL-CME web interface.

Name	Description	Caveats	Units	Default Value or Actual Value
br_map.cr_no	Carrington rotation number during the period of interest.			2080
br_map.cr_ut	A specific universal date and time of interest. The magnetic field map to be downloaded and/or used is for the Carrington rotation period in which the central meridian longitude of this date and time falls in.			2009-02-13 00:00
br_map.obs	Data source (observatory).	Will be set to 'user' if the user provided a tar file with a previously saved session.		user
br_map.orig_obs	Original data source (observatory)	This can be different from br_map.obs in case of user-provided data files.		gong
br_map.cpal	Color map of the Br map plot.			blue_red
br_map.cmin	Minimum scale for the Br map plot.		Gauss	-206.76000000000000
br_map.cmax	Maximum scale for the Br map plot.		Gauss	206.76000000000000
br_map.cr_in	Indicates whether to pick the event from the date field or the Carrington rotation number field.			cr_ut

2.1 Source Region Selection

A set of parameters for the Active Region (AR) selection in the RBSL Flux Rope Designer.

Name	Description	Caveats	Units	Default Value or Actual Value
ar0.lon0	Top left longitude of the active region.		deg	270.50000000000000000
ar0.lat0	Top left latitude of the active region.		deg	0.50000000000000000
ar0.lon1	Bottom right longitude of the active region.		deg	281.77000000000000000
ar0.lat1	Bottom right latitude of the active region.		deg	-9.78000000000000000
ar0.min_ds	Minimum delta of the mesh resolution.	Set this value to the minimum resolution that is used inside the selected region. Changing this value affects the number of points required to resolve the region as well as the smoothness of the magnetic field.		0.004000
ar0.max_npts	Maximum points of the mesh resolution.	Set this value to change the maximum number of points used to resolve the smaller side of the selected region. This ensures that for every size of the box you draw, it will be resolved appropriately. Changing this value affects the number of points required to resolve the region as well as the smoothness of the magnetic field.		70
ar0.ds	Mesh resolution delta.	This is the mesh resolution (in radians) that is used to resolve the selected region. It is determined by the size of the selected region and the two mesh parameters (minimum delta and maxim number of points).	radians	0.004000
ar0.cpal	Color map of the Active region plot.			blue_red
ar0.cmin	Minimum scaling of the Active region map plot.		Gauss	-206.76000000000000
ar0.cmax	Maximum scaling of the Active region map image.		Gauss	206.76000000000000
ar0.br_r_scale	PIL Radial Height	Radial height at which the PIL is computed. Greater radial height value gives simpler PILs. A typical choice ranges from 1.00 to 1.02 Rs.	Rs	1.010
ar0.name	User-specified name of the active region.			Active Region 1

2.2 Polarity Inversion Line

A set of parameters for the PIL selection in the RBSL Flux Rope Designer.

Name	Description	Units	Default Value or Actual Value
nl0.lon0	Starting longitude of the Polarity Inversion Line.	deg	276.64
nl0.lat0	Starting latitude of the Polarity Inversion Line.	deg	-3.09
nl0.lon1	Ending longitude of the Polarity Inversion Line.	deg	275.72
nl0.lat1	Ending latitude of the Polarity Inversion Line.	deg	-6.95
nl0.bp_cmin	Minimum scaling of the External Field (Bp) plot.	Gauss	-28.42646101812790
nl0.bp_cmax	Maximum scaling of the External Field (Bp) plot.	Gauss	7.10661525453196

2.3 RBSL Flux Rope Designer

A set of parameters for the flux rope selection in the RBSL Flux Rope Designer.

Name	Description	Caveats	Units	Default Value or Actual Value
fr0.fp_pos_lon	Positive longitude of the foot point of the flux rope.		deg	277.33
fr0.fp_pos_lat	Positive latitude of the foot point of the flux rope.		deg	-3.99
fr0.fp_neg_lon	Negative longitude of the foot point of the flux rope.		deg	274.79
fr0.fp_neg_lat	Negative latitude of the foot point of the flux rope.		deg	-6.11
fr0.apex_height	Apex height of the flux rope.	Height of the default axis peak above the photosphere. The default axis peak is located at the flux-rope center.	Rs	0.0250
fr0.fract_width	Fractional width of the flux rope.	Ratio of the flux rope minor radius (thickness) to the apex height. Thin flux ropes are likely to become kink-unstable. To suppress the instability, a value of 0.7 or larger should be chosen.	Rs	0.7000
fr0.rbsl_frac_curr	Fraction of the flux-rope current with respect to the 'optimized' current.	A factor of 1 would minimize the total Lorentz forces along the body of the initial flux-rope geometry. In practice, the exact threshold for eruption depends on the details of the configuration and needs to be established through experimentation. Typically, numbers lower than the range 1.2 to 1.5 yield a stable configuration, while numbers above that lead to an eruption. The default value (1.5) should produce an eruption in most cases.		1.500
fr0.axial_current	Axial current of the flux rope converted to code units.	Because of how we choose to define the RBSL path, the sign signifies whether the current flows from the negative footpoint to the positive footpoint (+) or the positive footpoint to the negative footpoint (-). This will match what is in the MAS input file.		2.768552
fr0.insertion_delay_mins	Insertion delay time of the flux rope in the final CME simulation.	Use this only if multiple CMEs are modeled. The maximum delay time is 240 minutes. Note: This only applies to the CME simulation in Step 4, not to the zero-beta runs in Step 1.	min	0.00
fr0.rbsl_minor_r	Minor radius (half diameter) of the flux rope.		Rs	0.0175
fr0.rbsl_hmod_pos	The position of the peak of the flux-rope axis. 0 is the left-most point and 1 is the right-most point along the axis.			0.50
fr0.rbsl_hmod_width	The width of the axis peak.			0.50
fr0.rbsl_hmod_amp	The amplitude of the axis peak.			1.00

MAS Settings - CME Eruption

High-resolution thermodynamic CME simulation in the corona, where a flux rope is inserted and erupts depending on the user-defined parameter settings in the interface.

Name	Description	Caveats	Units	Valid Min	Valid Max
cme_er.mas.nr	cme_er.mas.nr
cme_er.mas.nt	cme_er.mas.nt
cme_er.mas.np	cme_er.mas.np
cme_er.mas.r0	cme_er.mas.r0
cme_er.mas.r1	cme_er.mas.r1
cme_er.mas.rfrac	cme_er.mas.rfrac
cme_er.mas.nfrmesh	cme_er.mas.nfrmesh
cme_er.mas.tfrac	cme_er.mas.tfrac
cme_er.mas.dtratio	cme_er.mas.dtratio
cme_er.mas.nftmesh	cme_er.mas.nftmesh
cme_er.mas.pfrac	cme_er.mas.pfrac
cme_er.mas.dpratio	cme_er.mas.dpratio
cme_er.mas.nfpmesh	cme_er.mas.nfpmesh
cme_er.mas.phishift	cme_er.mas.phishift
cme_er.mas.rs_final	cme_er.mas.rs_final
cme_er.mas.long_sequence_numbers	cme_er.mas.long_sequence_numbers
cme_er.mas.ipltxint	cme_er.mas.ipltxint
cme_er.mas.plotlist	List of output plot file separated by comma with no space in between.	One or more of: all, vr, vt, vp, br, bt, bp, rho
cme_er.mas.tmax	Maximum time of simulation in code units (CU), where 1 CU = 24 minutes		CU
cme_er.mas.slice_txint	cme_er.mas.slice_txint
cme_er.mas.slice_plotlist	cme_er.mas.slice_plotlist
cme_er.mas.slice_tp_radii	cme_er.mas.slice_tp_radii
cme_er.mas.ntmax	Maximum number of time steps
cme_er.mas.dtmax	Maximum time step in code unit (CU), where 1 CU = 24 minutes		CU
cme_er.mas.dtmin	Minimum time step in code units (CU), where 1 CU = 24 minutes
cme_er.mas.dt_init	cme_er.mas.dt_init
cme_er.mas.use_dt_init_on_restart	cme_er.mas.use_dt_init_on_restart
cme_er.mas.dt_max_wave_cfl	cme_er.mas.dt_max_wave_cfl
cme_er.mas.slund	Lundquist number
cme_er.mas.eta_profile	Resistivity profile
cme_er.mas.eta_profile_tp_file	cme_er.mas.eta_profile_tp_file
cme_er.mas.visc	Viscosity
cme_er.mas.visc_profile	Kind of viscosity (uniform or custom profile)
cme_er.mas.visc_profile_tp_file	cme_er.mas.visc_profile_tp_file
cme_er.mas.calculation_frame	cme_er.mas.calculation_frame
cme_er.mas.restart_calculation_frame	cme_er.mas.restart_calculation_frame
cme_er.mas.initial_field	Initial state for magnetic field (see routine initialize_magnetic_field).
cme_er.mas.bnfile	File containing the magnetogram data
cme_er.mas.potential_field_bc	cme_er.mas.potential_field_bc
cme_er.mas.b_in_gauss	Specify magnetic field in Gauss or code units (1CU = 2.205 G).	False - in code units, True - in Gauss.
cme_er.mas.time_at_start	cme_er.mas.time_at_start
cme_er.mas.use_sts_visc	cme_er.mas.use_sts_visc
cme_er.mas.use_sts_tc	cme_er.mas.use_sts_tc
cme_er.mas.ifprec_t	cme_er.mas.ifprec_t
cme_er.mas.ifprec_v	cme_er.mas.ifprec_v
cme_er.mas.ifprec_pot2d	cme_er.mas.ifprec_pot2d
cme_er.mas.ifprec_divb	cme_er.mas.ifprec_divb
cme_er.mas.epscg_a	Convergence criterion for matrix inversion for resistivity
cme_er.mas.epscg_v	cme_er.mas.epscg_v
cme_er.mas.epscg_potfld	cme_er.mas.epscg_potfld
cme_er.mas.upwind_a	Flag to set upwinding in the advection of the vector potential (1 - upwind, 0 - do not upwind).
cme_er.mas.upwind_v	Flag to set upwinding in the advection of the velocity (1 - upwind, 0 - do not upwind).
cme_er.mas.use_old_vdgv	Use old formulation to calculate advection of the velocity.
cme_er.mas.ncgmax	Maximum number of iterations in matrix inversion
cme_er.mas.tbc0	cme_er.mas.tbc0
cme_er.mas.rho0	Base density in code units (1 CU = 10^8 cm^−3)		CU
cme_er.mas.char_bc0	cme_er.mas.char_bc0
cme_er.mas.char_bc1	cme_er.mas.char_bc1
cme_er.mas.ubzero	cme_er.mas.ubzero
cme_er.mas.cfl	In fully explict runs, fraction of the Courant-Friedrichs-Lewy limit			0	1
cme_er.mas.advance_tc	cme_er.mas.advance_tc
cme_er.mas.tcond	cme_er.mas.tcond
cme_er.mas.t_cutoff1	cme_er.mas.t_cutoff1
cme_er.mas.p_cutoff1	cme_er.mas.p_cutoff1
cme_er.mas.nfilt_kappa	cme_er.mas.nfilt_kappa
cme_er.mas.radloss	cme_er.mas.radloss
cme_er.mas.ifimplrad	cme_er.mas.ifimplrad
cme_er.mas.rad_law	cme_er.mas.rad_law
cme_er.mas.tc_r	cme_er.mas.tc_r
cme_er.mas.tc_dr	cme_er.mas.tc_dr
cme_er.mas.tc_fac0	cme_er.mas.tc_fac0
cme_er.mas.tc_fac1	cme_er.mas.tc_fac1
cme_er.mas.print_matrix_t	cme_er.mas.print_matrix_t
cme_er.mas.alpha_nocoll	cme_er.mas.alpha_nocoll
cme_er.mas.tc_nocoll_r	cme_er.mas.tc_nocoll_r
cme_er.mas.tc_nocoll_dr	cme_er.mas.tc_nocoll_dr
cme_er.mas.pole_filter_rho	cme_er.mas.pole_filter_rho
cme_er.mas.pole_filter_t	cme_er.mas.pole_filter_t
cme_er.mas.pole_filter_pw	cme_er.mas.pole_filter_pw
cme_er.mas.iftfloor	cme_er.mas.iftfloor
cme_er.mas.tfloor	cme_er.mas.tfloor
cme_er.mas.ifcheck0temp	cme_er.mas.ifcheck0temp
cme_er.mas.ifrholimit	cme_er.mas.ifrholimit
cme_er.mas.rho_limit_factor	cme_er.mas.rho_limit_factor
cme_er.mas.rho_limit_file	cme_er.mas.rho_limit_file
cme_er.mas.max_alfven_speed_rho_mod	cme_er.mas.max_alfven_speed_rho_mod
cme_er.mas.advance_pw	Flag to turn on Alfven wave pressure
cme_er.mas.pw0	cme_er.mas.pw0
cme_er.mas.rho_aw	cme_er.mas.rho_aw
cme_er.mas.use_pw_theta_profile	cme_er.mas.use_pw_theta_profile
cme_er.mas.awthprof_th	cme_er.mas.awthprof_th
cme_er.mas.awthprof2_th	cme_er.mas.awthprof2_th
cme_er.mas.awthprof_dth	cme_er.mas.awthprof_dth
cme_er.mas.fluxrope_preserve_br0	cme_er.mas.fluxrope_preserve_br0
cme_er.mas.rbsl_fluxrope_1.add	cme_er.mas.rbsl_fluxrope(1)%add
cme_er.mas.rbsl_fluxrope_1.path_file	cme_er.mas.rbsl_fluxrope(1)%path_file
cme_er.mas.rbsl_fluxrope_1.axial_current	cme_er.mas.rbsl_fluxrope(1)%axial_current
cme_er.mas.rbsl_fluxrope_1.minor_radius	cme_er.mas.rbsl_fluxrope(1)%minor_radius
cme_er.mas.rbsl_fluxrope_1.handedness	cme_er.mas.rbsl_fluxrope(1)%handedness
cme_er.mas.rbsl_fluxrope_1.mirror_poloidal	cme_er.mas.rbsl_fluxrope(1)%mirror_poloidal
cme_er.mas.rbsl_fluxrope_1.insertion_time	cme_er.mas.rbsl_fluxrope(1)%insertion_time
cme_er.mas.rbsl_fluxrope_2.add	cme_er.mas.rbsl_fluxrope(2)%add
cme_er.mas.rbsl_fluxrope_2.path_file	cme_er.mas.rbsl_fluxrope(2)%path_file
cme_er.mas.rbsl_fluxrope_2.axial_current	cme_er.mas.rbsl_fluxrope(2)%axial_current
cme_er.mas.rbsl_fluxrope_2.minor_radius	cme_er.mas.rbsl_fluxrope(2)%minor_radius
cme_er.mas.rbsl_fluxrope_2.handedness	cme_er.mas.rbsl_fluxrope(2)%handedness
cme_er.mas.rbsl_fluxrope_2.mirror_poloidal	cme_er.mas.rbsl_fluxrope(2)%mirror_poloidal
cme_er.mas.rbsl_fluxrope_2.insertion_time	cme_er.mas.rbsl_fluxrope(2)%insertion_time
cme_er.mas.trace_particles	cme_er.mas.trace_particles
cme_er.mas.trace_seed_file	cme_er.mas.trace_seed_file
cme_er.mas.trace_ixint	cme_er.mas.trace_ixint

MAS Settings - CME Propagation

The continuation of the thermodynamic CME simulation, where the high-resolution results from the first part are remeshed into a lower resolution to propagate the CME out into the rest of the corona

Name	Description	Caveats	Units	Valid Min	Valid Max
cme_pr.mas.nr	cme_pr.mas.nr
cme_pr.mas.nt	cme_pr.mas.nt
cme_pr.mas.np	cme_pr.mas.np
cme_pr.mas.r0	cme_pr.mas.r0
cme_pr.mas.r1	cme_pr.mas.r1
cme_pr.mas.rfrac	cme_pr.mas.rfrac
cme_pr.mas.drratio	cme_pr.mas.drratio
cme_pr.mas.nfrmesh	cme_pr.mas.nfrmesh
cme_pr.mas.tfrac	cme_pr.mas.tfrac
cme_pr.mas.dtratio	cme_pr.mas.dtratio
cme_pr.mas.nftmesh	cme_pr.mas.nftmesh
cme_pr.mas.pfrac	cme_pr.mas.pfrac
cme_pr.mas.dpratio	cme_pr.mas.dpratio
cme_pr.mas.nfpmesh	cme_pr.mas.nfpmesh
cme_pr.mas.phishift	cme_pr.mas.phishift
cme_pr.mas.rs_final	cme_pr.mas.rs_final
cme_pr.mas.long_sequence_numbers	cme_pr.mas.long_sequence_numbers
cme_pr.mas.ipltxint	cme_pr.mas.ipltxint
cme_pr.mas.plotlist	List of output plot file separated by comma with no space in between.	One or more of: all, vr, vt, vp, br, bt, bp, rho
cme_pr.mas.tmax	Maximum time of simulation in code units (CU), where 1 CU = 24 minutes		CU
cme_pr.mas.slice_txint	cme_pr.mas.slice_txint
cme_pr.mas.slice_plotlist	cme_pr.mas.slice_plotlist
cme_pr.mas.slice_tp_radii	cme_pr.mas.slice_tp_radii
cme_pr.mas.ntmax	Maximum number of time steps
cme_pr.mas.dtmax	Maximum time step in code unit (CU), where 1 CU = 24 minutes		CU
cme_pr.mas.dtmin	Minimum time step in code units (CU), where 1 CU = 24 minutes
cme_pr.mas.dt_init	cme_pr.mas.dt_init
cme_pr.mas.use_dt_init_on_restart	cme_pr.mas.use_dt_init_on_restart
cme_pr.mas.dt_max_wave_cfl	cme_pr.mas.dt_max_wave_cfl
cme_pr.mas.slund	Lundquist number
cme_pr.mas.eta_profile	Resistivity profile
cme_pr.mas.eta_profile_tp_file	cme_pr.mas.eta_profile_tp_file
cme_pr.mas.visc	Viscosity
cme_pr.mas.visc_profile	Kind of viscosity (uniform or custom profile)
cme_pr.mas.visc_profile_tp_file	cme_pr.mas.visc_profile_tp_file
cme_pr.mas.calculation_frame	cme_pr.mas.calculation_frame
cme_pr.mas.restart_calculation_frame	cme_pr.mas.restart_calculation_frame
cme_pr.mas.initial_field	Initial state for magnetic field (see routine initialize_magnetic_field).
cme_pr.mas.bnfile	File containing the magnetogram data
cme_pr.mas.potential_field_bc	cme_pr.mas.potential_field_bc
cme_pr.mas.b_in_gauss	Specify magnetic field in Gauss or code units (1CU = 2.205 G).	False - in code units, True - in Gauss.
cme_pr.mas.time_at_start	cme_pr.mas.time_at_start
cme_pr.mas.use_sts_visc	cme_pr.mas.use_sts_visc
cme_pr.mas.use_sts_tc	cme_pr.mas.use_sts_tc
cme_pr.mas.ifprec_t	cme_pr.mas.ifprec_t
cme_pr.mas.ifprec_v	cme_pr.mas.ifprec_v
cme_pr.mas.ifprec_pot2d	cme_pr.mas.ifprec_pot2d
cme_pr.mas.ifprec_divb	cme_pr.mas.ifprec_divb
cme_pr.mas.epscg_a	Convergence criterion for matrix inversion for resistivity
cme_pr.mas.epscg_v	cme_pr.mas.epscg_v
cme_pr.mas.epscg_potfld	cme_pr.mas.epscg_potfld
cme_pr.mas.upwind_a	Flag to set upwinding in the advection of the vector potential (1 - upwind, 0 - do not upwind).
cme_pr.mas.upwind_v	Flag to set upwinding in the advection of the velocity (1 - upwind, 0 - do not upwind).
cme_pr.mas.use_old_vdgv	Use old formulation to calculate advection of the velocity.
cme_pr.mas.ncgmax	Maximum number of iterations in matrix inversion
cme_pr.mas.tbc0	cme_pr.mas.tbc0
cme_pr.mas.rho0	Base density in code units (1 CU = 10^8 cm^−3)		CU
cme_pr.mas.char_bc0	cme_pr.mas.char_bc0
cme_pr.mas.char_bc1	cme_pr.mas.char_bc1
cme_pr.mas.ubzero	cme_pr.mas.ubzero
cme_pr.mas.cfl	In fully explict runs, fraction of the Courant-Friedrichs-Lewy limit			0	1
cme_pr.mas.advance_tc	cme_pr.mas.advance_tc
cme_pr.mas.tcond	cme_pr.mas.tcond
cme_pr.mas.t_cutoff1	cme_pr.mas.t_cutoff1
cme_pr.mas.p_cutoff1	cme_pr.mas.p_cutoff1
cme_pr.mas.nfilt_kappa	cme_pr.mas.nfilt_kappa
cme_pr.mas.radloss	cme_pr.mas.radloss
cme_pr.mas.ifimplrad	cme_pr.mas.ifimplrad
cme_pr.mas.rad_law	cme_pr.mas.rad_law
cme_pr.mas.tc_r	cme_pr.mas.tc_r
cme_pr.mas.tc_dr	cme_pr.mas.tc_dr
cme_pr.mas.tc_fac0	cme_pr.mas.tc_fac0
cme_pr.mas.tc_fac1	cme_pr.mas.tc_fac1
cme_pr.mas.print_matrix_t	cme_pr.mas.print_matrix_t
cme_pr.mas.alpha_nocoll	cme_pr.mas.alpha_nocoll
cme_pr.mas.tc_nocoll_r	cme_pr.mas.tc_nocoll_r
cme_pr.mas.tc_nocoll_dr	cme_pr.mas.tc_nocoll_dr
cme_pr.mas.pole_filter_rho	cme_pr.mas.pole_filter_rho
cme_pr.mas.pole_filter_t	cme_pr.mas.pole_filter_t
cme_pr.mas.pole_filter_pw	cme_pr.mas.pole_filter_pw
cme_pr.mas.iftfloor	cme_pr.mas.iftfloor
cme_pr.mas.tfloor	cme_pr.mas.tfloor
cme_pr.mas.ifcheck0temp	cme_pr.mas.ifcheck0temp
cme_pr.mas.ifrholimit	cme_pr.mas.ifrholimit
cme_pr.mas.rho_limit_factor	cme_pr.mas.rho_limit_factor
cme_pr.mas.rho_limit_file	cme_pr.mas.rho_limit_file
cme_pr.mas.max_alfven_speed_rho_mod	cme_pr.mas.max_alfven_speed_rho_mod
cme_pr.mas.advance_pw	Flag to turn on Alfven wave pressure
cme_pr.mas.pw0	cme_pr.mas.pw0
cme_pr.mas.rho_aw	cme_pr.mas.rho_aw
cme_pr.mas.use_pw_theta_profile	cme_pr.mas.use_pw_theta_profile
cme_pr.mas.awthprof_th	cme_pr.mas.awthprof_th
cme_pr.mas.awthprof2_th	cme_pr.mas.awthprof2_th
cme_pr.mas.awthprof_dth	cme_pr.mas.awthprof_dth
cme_pr.mas.trace_particles	cme_pr.mas.trace_particles
cme_pr.mas.trace_seed_file	cme_pr.mas.trace_seed_file
cme_pr.mas.trace_ixint	cme_pr.mas.trace_ixint

2.1 Source Region Selection (AR#2)

A set of parameters for an additional AR in the RBSL Flux Rope Designer. The CORHEL-CME code allows modeling of up to two events concurrently, with the second one being optional.

Name	Description	Caveats	Units	Default Value or Actual Value
ar1.lon0	Top left longitude of the second active region.		deg
ar1.lat0	Top left latitude of the second active region.		deg
ar1.lon1	Bottom right longitude of the second active region.		deg
ar1.lat1	Bottom right latitude of the second active region.		deg
ar1.min_ds	Minimum delta of the mesh resolution.	Set this value to the minimum resolution that is used inside the selected region. Changing this value affects the number of points required to resolve the region as well as the smoothness of the magnetic field.
ar1.max_npts	Maximum points of the mesh resolution.	Set this value to change the maximum number of points used to resolve the smaller side of the selected region. This ensures that for every size of the box you draw, it will be resolved appropriately. Changing this value affects the number of points required to resolve the region as well as the smoothness of the magnetic field.
ar1.ds	Mesh resolution delta.	This is the mesh resolution (in radians) that is used to resolve the selected region. It is determined by the size of the selected region and the two mesh parameters (minimum delta and maxim number of points).	radians
ar1.cpal	Color map of the Active region plot.			blue_red
ar1.cmin	Minimum scaling of the Active region map plot.		Gauss
ar1.cmax	Maximum scaling of the Active region map image.		Gauss
ar1.br_r_scale	PIL Radial Height.	Radial height at which the PIL is computed. Greater radial height value gives simpler PILs. A typical choice ranges from 1.00 to 1.02 Rs.
ar1.name	User-specified name of the second active region.

2.2 Polarity Inversion Line (PIL#2)

A set of parameters for an additional PIL in the RBSL Flux Rope Designer. The CORHEL-CME code allows modeling of up to two events concurrently, with the second one being optional.

Name	Description	Units
nl1.lon0	Starting longitude of the second Polarity Inversion Line.	deg
nl1.lat0	Starting latitude of the second Polarity Inversion Line.	deg
nl1.lon1	Ending longitude of the second Polarity Inversion Line.	deg
nl1.lat1	Ending latitude of the second Polarity Inversion Line.	deg
nl1.bp_cmin	Minimum scaling of the External Field (Bp) plot.	Gauss
nl1.bp_cmax	Maximum scaling of the External Field (Bp) plot.	Gauss

2.3 RBSL Flux Rope Designer (#2)

A set of parameters for an additional flux rope in the RBSL Flux Rope Designer. The CORHEL-CME code allows modeling of up to two events concurrently, with the second one being optional.

Name	Description	Caveats	Units	Default Value or Actual Value
fr1.fp_pos_lon	Positive longitude of the foot point of the second flux rope.
fr1.fp_pos_lat	Positive latitude of the foot point of the second flux rope.		deg
fr1.fp_neg_lon	Negative longitude of the foot point of the second flux rope.		deg
fr1.fp_neg_lat	Negative latitude of the foot point of the second flux rope.		deg
fr1.apex_height	Apex height of the second flux rope.	Height of the default axis peak above the photosphere. The default axis peak is located at the flux-rope center.	Rs
fr1.fract_width	Fractional width of the second flux rope.	Ratio of the flux rope minor radius (thickness) to the apex height. Thin flux ropes are likely to become kink-unstable. To suppress the instability, a value of 0.7 or larger should be chosen.	Rs
fr1.rbsl_frac_curr	Fraction of the flux-rope current with respect to the 'optimized' current.	A factor of 1 would minimize the total Lorentz forces along the body of the initial flux-rope geometry. In practice, the exact threshold for eruption depends on the details of the configuration and needs to be established through experimentation. Typically, numbers lower than the range 1.2 to 1.5 yield a stable configuration, while numbers above that lead to an eruption. The default value (1.5) should produce an eruption in most cases.		1.5
fr1.axial_current	Axial current of the flux rope converted to code units.	Because of how we choose to define the RBSL path, the sign signifies whether the current flows from the negative footpoint to the positive footpoint (+) or the positive footpoint to the negative footpoint (-). This will match what is in the MAS input file.
fr1.insertion_delay_mins	Insertion delay time of the second flux rope in the final CME simulation.	Use this only if multiple CMEs are modeled. The maximum delay time is 240 minutes. Note: This only applies to the CME simulation in Step 4, not to the zero-beta runs in Step 1.	min	0
fr1.rbsl_minor_r	Minor radius (half diameter) of the second flux rope.		Rs
fr1.rbsl_hmod_pos	The position of the peak of the flux-rope axis. 0 is the left-most point and 1 is the right-most point along the axis.
fr1.rbsl_hmod_width	The width of the axis peak.
fr1.rbsl_hmod_amp	The amplitude of the axis peak.

MAS Settings - CME Propagation Inner Heliosphere

The CME propagation from close to the Sun through the inner heliosphere.

Name	Description	Caveats	Units
cme_hel.mas.nr	cme_hel.mas.nr
cme_hel.mas.nt	cme_hel.mas.nt
cme_hel.mas.np	cme_hel.mas.np
cme_hel.mas.r0	cme_hel.mas.r0
cme_hel.mas.r1	cme_hel.mas.r1
cme_hel.mas.rfrac	cme_hel.mas.rfrac
cme_hel.mas.drratio	cme_hel.mas.drratio
cme_hel.mas.nfrmesh	cme_hel.mas.nfrmesh
cme_hel.mas.tfrac	cme_hel.mas.tfrac
cme_hel.mas.dtratio	cme_hel.mas.dtratio
cme_hel.mas.nftmesh	cme_hel.mas.nftmesh
cme_hel.mas.pfrac	cme_hel.mas.pfrac
cme_hel.mas.dpratio	cme_hel.mas.dpratio
cme_hel.mas.nfpmesh	cme_hel.mas.nfpmesh
cme_hel.mas.phishift	cme_hel.mas.phishift
cme_hel.mas.rs_final	cme_hel.mas.rs_final
cme_hel.mas.plotlist	List of output plot file separated by comma with no space in between.	One or more of: all, vr, vt, vp, br, bt, bp, rho
cme_hel.mas.long_sequence_numbers	cme_hel.mas.long_sequence_numbers
cme_hel.mas.tmax	Maximum time of simulation in code units (CU), where 1 CU = 24 minutes		CU
cme_hel.mas.tpltxint	Plot field diagnostics every tpltxint Alfven times, in code units (CU), where 1 CU = 24 minutes.		CU
cme_hel.mas.ntmax	Maximum number of time steps
cme_hel.mas.dtmax	Maximum time step in code unit (CU), where 1 CU = 24 minutes		CU
cme_hel.mas.dtmin	Minimum time step in code units (CU), where 1 CU = 24 minutes		CU
cme_hel.mas.dt_init	cme_hel.mas.dt_init
cme_hel.mas.use_sts_visc	cme_hel.mas.use_sts_visc
cme_hel.mas.use_sts_tc	cme_hel.mas.use_sts_tc
cme_hel.mas.ifprec_v	cme_hel.mas.ifprec_v
cme_hel.mas.ifprec_t	cme_hel.mas.ifprec_t
cme_hel.mas.ifprec_pot2d	cme_hel.mas.ifprec_pot2d
cme_hel.mas.ifprec_divb	cme_hel.mas.ifprec_divb
cme_hel.mas.epscg_a	Convergence criterion for matrix inversion for resistivity
cme_hel.mas.epscg_v	cme_hel.mas.epscg_v
cme_hel.mas.epscg_potfld	cme_hel.mas.epscg_potfld
cme_hel.mas.ncgmax	Maximum number of iterations in matrix inversion
cme_hel.mas.pole_filter_rho	cme_hel.mas.pole_filter_rho
cme_hel.mas.pole_filter_t	cme_hel.mas.pole_filter_t
cme_hel.mas.slice_txint	cme_hel.mas.slice_txint
cme_hel.mas.slice_plotlist	cme_hel.mas.slice_plotlist
cme_hel.mas.slice_tp_radii	cme_hel.mas.slice_tp_radii
cme_hel.mas.slund	Lundquist number
cme_hel.mas.visc	Viscosity
cme_hel.mas.char_bc0	cme_hel.mas.char_bc0
cme_hel.mas.char_bc1	cme_hel.mas.char_bc1
cme_hel.mas.tbc0	cme_hel.mas.tbc0
cme_hel.mas.upwind_a	Flag to set upwinding in the advection of the vector potential (1 - upwind, 0 - do not upwind).
cme_hel.mas.upwind_v	Flag to set upwinding in the advection of the velocity (1 - upwind, 0 - do not upwind).
cme_hel.mas.use_old_vdgv	Use old formulation to calculate advection of the velocity.
cme_hel.mas.ncgmax	Maximum number of iterations in matrix inversion
cme_hel.mas.betapc_a	cme_hel.mas.betapc_a
cme_hel.mas.betapc_v_flow	cme_hel.mas.betapc_v_flow
cme_hel.mas.betapc_v_wave	cme_hel.mas.betapc_v_wave
cme_hel.mas.betapc_t_flow	cme_hel.mas.betapc_t_flow
cme_hel.mas.betapc_t_wave	cme_hel.mas.betapc_t_wave
cme_hel.mas.betapc_rho_flow	cme_hel.mas.betapc_rho_flow
cme_hel.mas.betapc_rho_wave	cme_hel.mas.betapc_rho_wave
cme_hel.mas.use_old_jxb	cme_hel.mas.use_old_jxb
cme_hel.mas.gamma	Ratio of specific heats
cme_hel.mas.ubzero	cme_hel.mas.ubzero
cme_hel.mas.n_ssdip	cme_hel.mas.n_ssdip
cme_hel.mas.advance_tc	cme_hel.mas.advance_tc
cme_hel.mas.tcond	cme_hel.mas.tcond
cme_hel.mas.ifimplrad	cme_hel.mas.ifimplrad
cme_hel.mas.iftfloor	cme_hel.mas.iftfloor
cme_hel.mas.omega_corotate	cme_hel.mas.omega_corotate
cme_hel.mas.calculation_frame	cme_hel.mas.calculation_frame
cme_hel.mas.restart_calculation_frame	cme_hel.mas.restart_calculation_frame
cme_hel.mas.long_sequence_numbers_input	cme_hel.mas.long_sequence_numbers_input

Run Access Information :

Keyword(s) :

CORHEL

Linked to Parent Spase Resource(s) (example: SimulationModel) :

spase://CCMC/SimulationModel/CORHEL-CME/1

List of Publications Referencing this Run :

CMR Footer

Curator: Chiu Wiegand | NASA Official: Dr. Masha Kuznetsova | Privacy and Security Notices | Accessibility | CCMC Data Collection Consent Agreement