CMR Header

CCMC Metadata Registry (CMR)

Menu Page

Go to:

 CMR: View Simulation Run Info

CORHEL-CME/ZB/1.0/ROR/Abril_Sahade_050124_SH_1

Run Description 

CCMC ROR run for CORHEL-CME/ZB 1.0 with RORID Abril_Sahade_050124_SH_1

Simulation Time:

Simulation Start Time: 2021-11-02T01:30:00Z
Simulation End Time: 2021-11-02T01:30: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 :

abril.sahade, 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 Valid Min Valid Max
br_map.cr_no 
Carrington rotation number during the period of interest.
 2250
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.
 2021-11-02 01:30
br_map.obs 
Data source (observatory).
 Will be set to 'user' if the user provided a tar file with a previously saved session. hmi
br_map.orig_obs 
Original data source (observatory)
 This can be different from br_map.obs in case of user-provided data files. None
br_map.cpal 
Color map of the Br map plot.
 blue_red
br_map.cmin 
Minimum scale for the Br map plot.
 Gauss -500.00000000000000
br_map.cmax 
Maximum scale for the Br map plot.
 Gauss 500.00000000000000
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 Valid Min Valid Max
ar0.lon0 
Top left longitude of the active region.
 deg 186.33000000000000000
ar0.lat0 
Top left latitude of the active region.
 deg 38.53000000000000000
ar0.lon1 
Bottom right longitude of the active region.
 deg 216.95000000000000000
ar0.lat1 
Bottom right latitude of the active region.
 deg 9.77000000000000000
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.007173
ar0.cpal 
Color map of the Active region plot.
 blue_red
ar0.cmin 
Minimum scaling of the Active region map plot.
Gauss -100.00000000000000
ar0.cmax 
Maximum scaling of the Active region map image.
 Gauss 100.00000000000000
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.040
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 Caveats Units Default Value or Actual Value Valid Min Valid Max
nl0.lon0 
Starting longitude of the Polarity Inversion Line.
 deg 202.60
nl0.lat0 
Starting latitude of the Polarity Inversion Line.
 deg 25.59
nl0.lon1 
Ending longitude of the Polarity Inversion Line.
 deg 209.25
nl0.lat1 
Ending latitude of the Polarity Inversion Line.
 deg 20.20
nl0.bp_cmin 
Minimum scaling of the External Field (Bp) plot.
 Gauss -91.44442045883760
nl0.bp_cmax 
Maximum scaling of the External Field (Bp) plot.
 Gauss 22.86110511470940
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 Valid Min Valid Max
fr0.fp_pos_lon 
Positive longitude of the foot point of the flux rope.
 deg 211.21
fr0.fp_pos_lat 
Positive latitude of the foot point of the flux rope.
 deg 15.59
fr0.fp_neg_lon 
Negative longitude of the foot point of the flux rope.
 deg 197.94
fr0.fp_neg_lat 
Negative latitude of the foot point of the flux rope.
 deg 29.15
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. 8.246449
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 Zero-Beta Simulation Run 
Settings for the MAS zero-beta simulation run.
Name Description Caveats Units Default Value or Actual Value Valid Min Valid Max
zb.mas.tpltxint 
Plot field diagnostics every tpltxint Alfven times, in code units (CU), where 1 CU = 24 minutes
 CU 4.821256e-02 0
zb.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 vr,vt,vp,br,bt,bp,jr,jt,jp,rho
zb.mas.tmax 
Maximum time of simulation in code units (CU), where 1 CU = 24 minutes
 CU 0.578551 0
zb.mas.ntmax 
Maximum number of time steps
 10000000 0
zb.mas.dtmax 
Maximum time step in code unit (CU), where 1 CU = 24 minutes
 CU 0.010719 0
zb.mas.dtmin 
Minimum time step in code units (CU), where 1 CU = 24 minutes
 CU 2.143896e-06 0
zb.mas.dt_init 
zb.mas.dt_init
 1.071948e-04
zb.mas.slund 
Lundquist number
 1.e6 0
zb.mas.eta_profile 
Resistivity profile
 uniform
zb.mas.visc 
Viscosity
 .01 0
zb.mas.visc_profile 
Kind of viscosity (uniform or custom profile)
 CUSTOM_PROFILE
zb.mas.visc_profile_file 
zb.mas.visc_profile_file
 dsprof_lin_ds01_limited.hdf
zb.mas.char_bc 
Use characteristics boundary conditions
 True
zb.mas.cfl 
In fully explict runs, fraction of the Courant-Friedrichs-Lewy limit
 .6 0 1
zb.mas.epscg_a 
Convergence criterion for matrix inversion for resistivity
 1.e-11
zb.mas.ncgmax 
Maximum number of iterations in matrix inversion
 1000000
zb.mas.upwind_a 
Flag to set upwinding in the advection of the vector potential (1 - upwind, 0 - do not upwind).
 1.
zb.mas.upwind_v 
Flag to set upwinding in the advection of the velocity (1 - upwind, 0 - do not upwind).
 1.
zb.mas.use_old_vdgv 
Use old formulation to calculate advection of the velocity.
 True
zb.mas.advance_t 
Flag to advance the temperature
 False
zb.mas.advance_rho 
Flag to advance the density
 False
zb.mas.advance_tc 
zb.mas.advance_tc
 False
zb.mas.advance_pw 
Flag to turn on Alfven wave pressure
 False
zb.mas.zero_v_parallel 
zb.mas.zero_v_parallel
 True
zb.mas.gamma 
Ratio of specific heats
 1.05 1 5/3
zb.mas.rho0 
Base density in code units (1 CU = 10^8 cm^−3)
 CU 2.
zb.mas.g0 
Gravity at the inner radial boundary in code units (CU)
 CU 0.
zb.mas.initial_field 
Initial state for magnetic field (see routine initialize_magnetic_field).
 One of these: ALFVENWAVE1, ALFVENWAVE2, POTENTIALFIELD, DIPOLE. POTENTIAL_FIELD
zb.mas.bnfile 
File containing the magnetogram data
 br_zb_cme_tp.hdf
zb.mas.b_in_gauss 
Specify magnetic field in Gauss or code units (1CU = 2.205 G).
 False - in code units, True - in Gauss. True
zb.mas.initial_plasma 
Initial plasma equilibrium.
 One of these: HYDROSTATIC, ZERO-BETA, 1DFILE, STREAMER. ZERO-BETA
zb.mas.zb_rho_profile 
zb.mas.zb_rho_profile
 READ_FROM_FILE
zb.mas.zb_rho_file 
zb.mas.zb_rho_file
 rho_zb.hdf
zb.mas.trace_particles 
zb.mas.trace_particles
 True
zb.mas.trace_seed_file 
zb.mas.trace_seed_file
 tracer_input.hdf
zb.mas.trace_txint 
zb.mas.trace_txint
 4.821256e-02
zb.mas.fluxrope_preserve_br0 
zb.mas.fluxrope_preserve_br0
 False
zb.mas.fluxrope_1.add 
zb.mas.fluxrope_1.add
zb.mas.fluxrope_1.path_file 
zb.mas.fluxrope_1.path_file
zb.mas.fluxrope_1.axial_current 
zb.mas.fluxrope_1.axial_current
zb.mas.fluxrope_1.minor_radius 
zb.mas.fluxrope_1.minor_radius
zb.mas.fluxrope_1.handedness 
zb.mas.fluxrope_1.handedness
zb.mas.fluxrope_1.mirror_poloidal 
zb.mas.fluxrope_1.mirror_poloidal
zb.mas.fluxrope_2.add 
zb.mas.fluxrope_2.add
zb.mas.fluxrope_2.path_file 
zb.mas.fluxrope_2.path_file
zb.mas.fluxrope_2.axial_current 
zb.mas.fluxrope_2.axial_current
zb.mas.fluxrope_2.minor_radius 
zb.mas.fluxrope_2.minor_radius
zb.mas.fluxrope_2.handedness 
zb.mas.fluxrope_2.handedness
zb.mas.fluxrope_2.mirror_poloidal 
zb.mas.fluxrope_2.mirror_poloidal
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 Valid Min Valid Max
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 Caveats Units Default Value or Actual Value Valid Min Valid Max
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 Valid Min Valid Max
fr1.fp_pos_lon 
Positive longitude of the foot point of the second flux rope.
deg
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.

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