Skip to content

Submits rudimentary sputtering bug fix for PR

vermarien requested to merge bugs/sputtering into master

Created by: tomasjames

At sufficiently low initial shock velocity, the value of the drift velocity (driftVel) through the C-type shock might be 0 km/s at certain time steps.

As the fractional sputtering abundance is dependent on 1/driftVel, a value of driftVel=0 causes this abundance to take on a NaN value. The error tolerances of the integrator require the exact values of the abundance, so NaN values break the code at this point.

The fix is to not call the sputtering routine from within the sublimation routine if driftVel=0.

Merge request reports