# Parameter Sweep Configuration for Bursting Bubble (Newtonian)
# This file defines a parameter sweep for batch simulations

# ============================================================
# Base Configuration
# ============================================================
# Start from this parameter file and override sweep variables
BASE_CONFIG=default.params

# ============================================================
# Case Number Range
# ============================================================
# Starting case number (4-digit: 1000-9999)
# Each parameter combination gets auto-incremented CaseNo
CASE_START=1000

# Ending case number (inclusive)
# Number of combinations should equal CASE_END - CASE_START + 1
# For this sweep: 16 Oh × 1 MAXlevel × 1 zWall = 16 combinations, so CASE_END = 1015
CASE_END=1015

# ============================================================
# Sweep Variables
# ============================================================
# Format: SWEEP_<variable>=value1,value2,value3,...
# The sweep will generate all combinations (cartesian product)

# Sweep Ohnesorge number (16 logarithmically spaced points from 0.01 to 0.1)
SWEEP_Oh=0.0100,0.0117,0.0136,0.0158,0.0185,0.0215,0.0251,0.0293,0.0341,0.0398,0.0464,0.0541,0.0631,0.0736,0.0858,0.1000

# Sweep mesh refinement level
# SWEEP_MAXlevel=12

# Sweep wall distance (uncomment to use)
# zWall: distance from bubble south pole to bottom wall
# Domain size: Ldomain = min(zWall + 6.0, 16.0)
# SWEEP_zWall=0.025

# NOTE: Bond is fixed at 1e-3 (only Bo0.0010.dat geometry available)
# To sweep Bond values, you would need to generate additional Bo*.dat files

# ============================================================
# Resolution sweep examples (commented out)
# ============================================================
# The new adaptive space/time knobs can be swept exactly like Oh, using the
# same SWEEP_<key> mechanism. Uncomment ONE block (and comment out SWEEP_Oh
# above) to run a focused convergence/sensitivity study.
#
# Mesh-convergence study (3 cases): vary the max refinement level.
# SWEEP_MAXlevel=10,11,12
#
# Timestep-ceiling sensitivity (2 cases): vary the adaptive dt ceiling.
# SWEEP_dtmax=1e-2,5e-3
#
# CFL sensitivity (2 cases): vary the CFL number.
# SWEEP_CFL=0.1,0.2
#
# IMPORTANT: CASE_END must equal CASE_START + (Cartesian product size) - 1,
# i.e. the product of the lengths of ALL ACTIVE (uncommented) SWEEP_ variables.
# Adding a second active sweep variable MULTIPLIES the case count: e.g.
# SWEEP_MAXlevel=10,11,12 alone => 3 cases (CASE_END = CASE_START + 2), but
# SWEEP_MAXlevel=10,11,12 together with SWEEP_dtmax=1e-2,5e-3 => 3 x 2 = 6
# cases (CASE_END = CASE_START + 5). Update CASE_END accordingly before running.

# ============================================================
# Output Configuration
# ============================================================
# Output folders are created in simulationCases/<CaseNo>/
# CaseNo auto-increments from CASE_START for each parameter combination
# Cases run sequentially (one at a time)