Using octave and a modified example from Julius Piso I got the following output. It does correlate well with the beam pattern from the real device with the 3 lobes. next step will be to extract the beam pattern from the simulation and plot it alongside the measured beam pattern.