Plot r(t) as (S-1) is varied

(S-1) is the supersaturation value
Add supersaturation analysis plot