set the minimul Bq per model-particle to a very low number

[heikoklein]

Avoid numerical subnormals not being removed from the model by setting a slightly higher limit of Bq/model-particle to 1000.*1e-38(=positive_subnormal).