Print frequency of max power for waveform.

gprMax/grid.py

@@ -199,6 +199,7 @@ def dispersion_analysis(G):
 
                 # Get frequency for max power
                 freqmaxpower = np.where(power[1::] == np.amax(power[1::]))[0][0]
+                print(freqs[freqmaxpower])
 
                 # Set maximum frequency to a threshold drop from maximum power, ignoring DC value
                 freq = np.where((np.amax(power[freqmaxpower::]) - power[freqmaxpower::]) > G.highestfreqthres)[0][0] + 1