I've been reviewing the code for CSE and haven't been able to find an easy answer to the following question -
How does California Simulation Engine avoid numerical instability while using euler forward difference surface models with fixed time steps (usually 2 mins)? Thin exterior layers that see solar are particularly problematic.
Are there adaptive time steps, or are exterior layers treated as massless?