Merge branch 'main' of https://github.com/G1OJS/NECBOL

Author: G1OJS

docs/outline.html

@@ -190,7 +190,11 @@
         'E_theta_mag', 'E_theta_phase_deg', 'E_phi_mag', 'E_phi_phase_deg'
     ]
 
+    n_theta, n_phi, theta_step, phi_step = _extract_pattern_grid_info(pattern)
+    print(n_theta, n_phi, theta_step, phi_step )
+    
     with open(file_path, "w") as f:
+        f.write(f" ***** DATA CARD NO.  5   RP   0   {n_theta}    {n_phi}  1003 0 0  {theta_step:.1f}  {phi_step:.1f} 0.0 0.0\n\n")
         f.write("                                                - - - RADIATION PATTERNS - - -\n")
         f.write("\n")
         f.write("  - - ANGLES - -           - POWER GAINS -       - - - POLARIZATION - - -    - - - E(THETA) - - -    - - - E(PHI) - - -\n")
@@ -208,7 +212,7 @@
                 line += f"{val:{fmt}}"
             f.write(line.rstrip() + "\n")
 
-
+        f.write("\n RUN TIME =   0\n")
 
 
 </pre>
@@ -297,7 +301,7 @@
 
     # Total field magnitude squared
     total_power = abs(E_theta)**2 + abs(E_phi)**2
-    total_gain_dBi = 10 * math.log10(total_power) if total_power &gt; 0 else -math.inf
+    total_gain_dBi = 10 * math.log10(total_power) if total_power &gt; 0 else -999
 
     # RHCP and LHCP components
     E_rhcp = (E_theta - 1j * E_phi) / math.sqrt(2)
@@ -306,19 +310,19 @@
     rhcp_power = abs(E_rhcp)**2
     lhcp_power = abs(E_lhcp)**2
 
-    rhcp_gain_dBi = 10 * math.log10(rhcp_power) if rhcp_power &gt; 0 else -math.inf
-    lhcp_gain_dBi = 10 * math.log10(lhcp_power) if lhcp_power &gt; 0 else -math.inf
+    rhcp_gain_dBi = 10 * math.log10(rhcp_power) if rhcp_power &gt; 0 else -999
+    lhcp_gain_dBi = 10 * math.log10(lhcp_power) if lhcp_power &gt; 0 else -999
 
     # Axial ratio in dB
     max_pol_power = max(rhcp_power, lhcp_power)
     min_pol_power = min(rhcp_power, lhcp_power)
     if min_pol_power == 0:
-        axial_ratio_dB = float(&#x27;inf&#x27;)
+        axial_ratio_dB = 999
     else:
         axial_ratio_dB = 10 * math.log10(max_pol_power / min_pol_power)
 
     # Polarization sense
-    if axial_ratio_dB == float(&#x27;inf&#x27;):
+    if axial_ratio_dB == 999:
         polarization_sense = &quot;Linear&quot;
     elif rhcp_power &gt; lhcp_power:
         polarization_sense = &quot;RHCP&quot;
@@ -344,8 +348,8 @@
     vert_power = abs(E_theta)**2
     horiz_power = abs(E_phi)**2
 
-    vert_gain_dBi = 10 * math.log10(vert_power) if vert_power &gt; 0 else -math.inf
-    horiz_gain_dBi = 10 * math.log10(horiz_power) if horiz_power &gt; 0 else -math.inf
+    vert_gain_dBi = 10 * math.log10(vert_power) if vert_power &gt; 0 else -999
+    horiz_gain_dBi = 10 * math.log10(horiz_power) if horiz_power &gt; 0 else -999
 
     return {
         &#x27;vert_gain_dBi&#x27;: vert_gain_dBi,
@@ -363,6 +367,36 @@
     }
 
 
+</pre>
+</details>
+<details><summary><span class ='def signature'> def _extract_pattern_grid_info(pattern):</span></summary>
+<pre class ='def content'>    # Compute theta and phi for each entry
+    for entry in pattern:
+        entry[&#x27;theta_deg&#x27;] = 90.0 - entry[&#x27;elevation_deg&#x27;]
+        entry[&#x27;phi_deg&#x27;] = entry[&#x27;azimuth_deg&#x27;]
+
+    # Extract sorted unique theta and phi values
+    theta_vals = sorted(set(round(e[&#x27;theta_deg&#x27;], 5) for e in pattern))
+    phi_vals = sorted(set(round(e[&#x27;phi_deg&#x27;], 5) for e in pattern))
+
+    # Determine steps (assuming uniform grid)
+    if len(theta_vals) &gt; 1:
+        theta_step = round(theta_vals[1] - theta_vals[0], 5)
+    else:
+        theta_step = 0.0
+
+    if len(phi_vals) &gt; 1:
+        phi_step = round(phi_vals[1] - phi_vals[0], 5)
+    else:
+        phi_step = 0.0
+
+    n_theta = len(theta_vals)
+    n_phi = len(phi_vals)
+
+    return n_theta, n_phi, theta_step, phi_step
+
+
+
 </pre>
 </details>
 <details><summary><span class ='def signature'> def _get_available_results(model):</span></summary>
@@ -629,7 +663,7 @@
                 LDTAG = LD[&#x27;iTag&#x27;]
                 R_Ohms, L_uH , C_pF = LD[&#x27;RoLuCp&#x27;]
                 # these strings are set programatically so shouldn&#x27;t need an error trap
-                f.write(f&quot;LD {LDTYP} {LDTAG} 0 0 {R_Ohms} {L_uH * 1e-6} {C_pF * 1e-12}\n&quot;)
+                f.write(f&quot;LD {LDTYP} {LDTAG} 0 0 {R_Ohms:8.4e} {L_uH * 1e-6:8.4e} {C_pF * 1e-12:8.4e}\n&quot;)
 
             # Feed
             f.write(f&quot;EX 0 {self.EX_tag} 1 0 1 0\n&quot;)

docs/text_list.txt

@@ -11,6 +11,7 @@ def _add_wire(self, iTag, nS, x1, y1, z1, x2, y2, z2, wr):
 def _compute_full_farfield_metrics(E_theta, E_phi):
 def _cos_sin(self,angle_deg):
 def _cosine_series(s, terms):
+def _extract_pattern_grid_info(pattern):
 def _fast_perpendicular(v, threshold=100):
 def _format_params(self, params):
 def _from_suffixed_dimensions(self, params: dict, whitelist=[]) -> dict: