anchor.diff:  
                   ! Calculate P_gen for the second excitation
                   pspawn_jk = X__ * abs(hmatel_jk) / Pgen_jk
  +
                   ! Attempt spawning
  -                psuccess = genrand_real2()
  -
  -                ! Multiple offspring again...
  -                nspawn_jk = int(pspawn_jk)
  -                pspawn_jk = pspawn_jk - nspawn_jk
  -
  -                ! Stochastic offspring
  -                if (pspawn_jk > psuccess) nspawn_jk = nspawn_jk + 1
  +                nspawn_jk = spawn_from_prob(pspawn_jk)
   
                   if (nspawn_jk > 0) then
                   ! Successful spawning on jk
  
positive.code:     end subroutine fs_stochastic_death
  
  
  end module folded_spectrum_utils
  
  
  
negative.code:             
              ! Calculate P_gen for the first excitation
              pspawn_ki = Xo_ * abs(hmatel_ki) / Pgen_ki
  
              ! Attempt spawning
              nspawn_ki = spawn_from_prob(pspawn_ki)
  
              if(nspawn_ki > 0 ) then
              ! Successful spawning on ki
  