From fdd7cee1f29a05ace359164e6e484514002f381f Mon Sep 17 00:00:00 2001 From: craig-warren Date: Sat, 11 Mar 2023 18:47:19 -0700 Subject: [PATCH] Tidying up patterns --- ...na_like_GSSI_1500_patterns_E_rxsorigin.txt | 3 -- examples/antenna_like_GSSI_1500_patterns_H | 48 ------------------- 2 files changed, 51 deletions(-) delete mode 100644 antenna_like_GSSI_1500_patterns_E_rxsorigin.txt delete mode 100644 examples/antenna_like_GSSI_1500_patterns_H diff --git a/antenna_like_GSSI_1500_patterns_E_rxsorigin.txt b/antenna_like_GSSI_1500_patterns_E_rxsorigin.txt deleted file mode 100644 index e16ca1d0..00000000 --- a/antenna_like_GSSI_1500_patterns_E_rxsorigin.txt +++ /dev/null @@ -1,3 +0,0 @@ -0.125000 -0.190000 -0.190000 diff --git a/examples/antenna_like_GSSI_1500_patterns_H b/examples/antenna_like_GSSI_1500_patterns_H deleted file mode 100644 index eae5bbb9..00000000 --- a/examples/antenna_like_GSSI_1500_patterns_H +++ /dev/null @@ -1,48 +0,0 @@ -#title: GSSI 1.5GHz antenna field patterns -#dx_dy_dz: 0.001 0.001 0.001 -#pml_cells: 14 - -#python: -import os -import numpy as np - -from gprMax.input_cmd_funcs import * -from user_libs.antennas.GSSI import antenna_like_GSSI_1500 - -filename = os.path.splitext(os.path.split(inputfile)[1])[0] - -timewindows = np.array([4.5e-9]) # For 0.3m max -radii = np.linspace(0.1, 0.3, 20) -theta = np.linspace(3, 357, 60) -materials = ['5 0 1 0 er5'] # Can add more to list and use selector integer to choose -selector = 0 - -fs = np.array([0.040, 0.040, 0.040]) -domain = np.array([2 * fs[0] + 2 * radii[-1], 2 * fs[1] + 0.107, 2 * fs[2] + 2 * radii[-1]]) -antennaposition = np.array([fs[0] + radii[-1], domain[1] / 2, fs[2] + radii[-1]]) -antenna_like_GSSI_1500(antennaposition[0], antennaposition[1], antennaposition[2]) - -print('#domain: {:.3f} {:.3f} {:.3f}'.format(domain[0], domain[1], domain[2])) -print('#time_window: {:.3e}'.format(timewindows[selector])) - -## Can introduce soil model -#print('#soil_peplinski: 0.5 0.5 2.0 2.66 0.001 0.25 mySoil') -#print('#fractal_box: 0 0 0 {} {} {} 1.5 1 1 1 50 mySoil mySoilBox 1'.format(domain[0], domain[1], fs[2] + radii[-1])) - -print('#material: {}'.format(materials[selector])) -print('#box: 0 0 0 {} {} {} {} n'.format(domain[0], domain[1], fs[2] + radii[-1], materials[selector].split()[-1])) - -## Save the position of the antenna to file for use when processing results -np.savetxt(os.path.join(os.path.dirname(inputfile), filename + '_rxsorigin.txt'), antennaposition, fmt="%f") - -## Generate receiver points for pattern -for radius in range(len(radii)): - ## H-plane circle (xz plane, y=0, phi=0,pi) - x = radii[radius] * np.sin(theta * np.pi /180) * np.cos(180 * np.pi / 180) - y = radii[radius] * np.sin(theta * np.pi /180) * np.sin(180 * np.pi / 180) - z = radii[radius] * np.cos(theta * np.pi /180) - for rxpt in range(len(theta)): - print('#rx: {:.3f} {:.3f} {:.3f}'.format(x[rxpt] + antennaposition[0], y[rxpt] + antennaposition[1], z[rxpt] + antennaposition[2])) - -geometry_view(0, 0, 0, domain[0], domain[1], domain[2], 0.001, 0.001, 0.001, filename, 'n') -#end_python: \ No newline at end of file