Support the leakage scheme for CCSN

[hfhsieh]

• Support the neutrino leakage scheme for CCSN simulations

  • Use NEUTRINO_SCHEME to control which neutrino scheme is used (LIGHTBULB/LEAKAGE)
  • The background profiles of leakage rays are constructed by Aux_ComputeRay()
    • Supports four configurations: r, (r, theta), (r, phi), and (r, theta, phi)
  • Add the runtime parameter OPT__OUTPUT_LEAKAGE to control whether to dump the heating rate (dEdt_Nu and dYedt_Nu) with the actual signs
    • Temporarily disable the checks on these two passive fields to avoid runtime errors in restart runs (since they can be negative in the snapshots when enabling OPT__OUTPUT_LEAKAGE)

• Miscellaneous

  • Remove the obsolete IsInit_dEdt_Nu
  • Update the output format in Record__CentralQuant
  • Update Detect_Shock() to include different measures of mean shock radius (using weight of volume and inverse volume)

Note:

• If the bad_alloc() runtime error in Aux_ComputeRay(), uncomment the line #define LEAKAGE_SEQUENTIAL_MAPPING in CPU_Src_Leakage.cpp to reduce memory usage