You have to activate the pointwise card instead of the groupe wise card.
I tested with this input: thermalneu.inp and it’s working nicely as you can see in the following plots. However, you have to use the FLUKAFIX card to avoid artefacts
no FLUKAFIX
FLUKAFIX
I let you try with and without the card to see the difference and find the appropriated parameters. Here are some useful resources to understand the inputs I sent: