Update to FDTD/MoM example.

tests/numerical/vs_MoM_MATLAB/antenna_bowtie_fs/antenna_bowtie_fs.in

@@ -5,32 +5,27 @@ from gprMax.input_cmd_funcs import *
 title = 'antenna_bowtie_fs'
 print('#title: {}'.format(title))
 
-domain = domain(0.200, 0.120, 0.120)
+#domain = domain(0.400, 0.300, 0.300)
+domain = domain(0.200, 0.200, 0.100)
 dxdydz = dx_dy_dz(0.001, 0.001, 0.001)
 time_window = time_window(30e-9)
 bowtie_dims = (0.050, 0.100) # Length, height
-tx_pos = (domain[0]/2, domain[1]/2, 0.020)
+tx_pos = (domain[0]/2, domain[1]/2, domain[2]/2)
 
 # Source excitation and type
 print('#waveform: gaussian 1 1.5e9 mypulse')
 print('#transmission_line: x {:g} {:g} {:g} 50 mypulse'.format(tx_pos[0], tx_pos[1], tx_pos[2]))
 
-# Output point
-print('#rx: {:g} {:g} {:g}'.format(tx_pos[0], tx_pos[1], tx_pos[2] + 0.020))
-
 # Bowtie - upper x half
-triangle(tx_pos[0] + dxdydz[0], tx_pos[1] - dxdydz[1], tx_pos[2], tx_pos[0] + bowtie_dims[0] + dxdydz[0], tx_pos[1] - bowtie_dims[1]/2 - dxdydz[1], tx_pos[2], tx_pos[0] + bowtie_dims[0] + dxdydz[0], tx_pos[1] + bowtie_dims[1]/2 - dxdydz[1], tx_pos[2], 0, 'pec')
+triangle(tx_pos[0] + dxdydz[0], tx_pos[1], tx_pos[2], tx_pos[0] + bowtie_dims[0] + dxdydz[0], tx_pos[1] - bowtie_dims[1]/2, tx_pos[2], tx_pos[0] + bowtie_dims[0] + dxdydz[0], tx_pos[1] + bowtie_dims[1]/2, tx_pos[2], 0, 'pec')
 
 # Bowtie - lower x half
-triangle(tx_pos[0], tx_pos[1], tx_pos[2], tx_pos[0] - bowtie_dims[0], tx_pos[1] - bowtie_dims[1]/2, tx_pos[2], tx_pos[0] - bowtie_dims[0], tx_pos[1] + bowtie_dims[1]/2, tx_pos[2], 0, 'pec')
+triangle(tx_pos[0] + dxdydz[0], tx_pos[1], tx_pos[2], tx_pos[0] - bowtie_dims[0], tx_pos[1] - bowtie_dims[1]/2, tx_pos[2], tx_pos[0] - bowtie_dims[0], tx_pos[1] + bowtie_dims[1]/2, tx_pos[2], 0, 'pec')
 
 # Uncomment to temporarily check location of source
 #edge(tx_pos[0], tx_pos[1], tx_pos[2], tx_pos[0] + dxdydz[0], tx_pos[1], tx_pos[2], 'pec')
 
-# Geometry view of entire domain
+# Detailed geometry view around bowtie and feed position
-geometry_view(0, 0, 0, domain[0], domain[1], domain[2], dxdydz[0], dxdydz[1], dxdydz[2], title)
-
-# Detailed geometry view around feed position
 geometry_view(tx_pos[0] - bowtie_dims[0] - 2*dxdydz[0], tx_pos[1] - bowtie_dims[1]/2  - 2*dxdydz[1], tx_pos[2] - 2*dxdydz[2], tx_pos[0] + bowtie_dims[0] + 2*dxdydz[0], tx_pos[1] + bowtie_dims[1]/2  + 2*dxdydz[1], tx_pos[2] + 2*dxdydz[2], dxdydz[0], dxdydz[1], dxdydz[2], title + '_pcb', type='f')
 
 #end_python: