Here are some thoughts. Simple models of a cylindrical magnet are 1) a uniform volume magnetization, w) an equivalent surface current K that replaces the volume with an effective solenoid or 3) as equivalent magnetic charges on the pole faces. All of these are strictly incorrect--but they still can be used to estimate the field everywhere outside of the magnet. (A dipole model, on the other hand, is useful only at distances that are large compared with the magnet dimensions.) The reason it may not pay to construct a high fidelity model, which includes non-uniform magnetization resulting from the demagnetizing field, is because of the large uncertainties you have about the actual magnet. For instance, your link lists a field strength but no detail about how it was obtained. Is that the highest saturated value of B inside the material? An average value? Was it measured outside? Where? (If a Hall effect probe was used, it likely represents some sort of average across the pole face at a distance of some mm--which is not small compared to the magnet dimensions).
In short I'm agreeing with "dipole" above; a high fidelity model is not justified by your low fidelity data, so you might as well use a lower fidelity model. You will get guesstimates with rather large error bars but this may still be useful to you.