@@ -517,8 +517,7 @@ def _radiation_vectors(self, theta: ArrayLikeN2F, phi: ArrayLikeN2F):
517
517
518
518
return N_theta , N_phi , L_theta , L_phi
519
519
520
- def fields_spherical (
521
- self , r : float , theta : ArrayLikeN2F , phi : ArrayLikeN2F ) -> xr .Dataset :
520
+ def fields_spherical (self , r : float , theta : ArrayLikeN2F , phi : ArrayLikeN2F ) -> xr .Dataset :
522
521
"""Get fields at a point relative to monitor center in spherical coordinates.
523
522
524
523
Parameters
@@ -546,8 +545,13 @@ def fields_spherical(
546
545
k = self .k
547
546
eta = self .eta
548
547
549
- scalar_proj_r = (- self .phasor_positive_sign * 1j * k * np .exp (
550
- self .phasor_positive_sign * 1j * k * r ) / (4 * np .pi * r ))
548
+ scalar_proj_r = (
549
+ - self .phasor_positive_sign
550
+ * 1j
551
+ * k
552
+ * np .exp (self .phasor_positive_sign * 1j * k * r )
553
+ / (4 * np .pi * r )
554
+ )
551
555
552
556
# assemble E felds
553
557
Et_array = - scalar_proj_r * (L_phi + eta * N_theta )
@@ -571,12 +575,12 @@ def fields_spherical(
571
575
Hp = xr .DataArray (data = Hp_array [None , ...], coords = coords , dims = dims )
572
576
573
577
field_data = xr .Dataset (
574
- {"E_r" : Er , "E_theta" : Et , "E_phi" : Ep , "H_r" : Hr , "H_theta" : Ht , "H_phi" : Hp })
578
+ {"E_r" : Er , "E_theta" : Et , "E_phi" : Ep , "H_r" : Hr , "H_theta" : Ht , "H_phi" : Hp }
579
+ )
575
580
576
581
return field_data
577
582
578
- def fields_cartesian (
579
- self , x : ArrayLikeN2F , y : ArrayLikeN2F , z : ArrayLikeN2F ) -> xr .Dataset :
583
+ def fields_cartesian (self , x : ArrayLikeN2F , y : ArrayLikeN2F , z : ArrayLikeN2F ) -> xr .Dataset :
580
584
"""Get fields at a point relative to monitor center in cartesian coordinates.
581
585
582
586
Parameters
@@ -614,10 +618,12 @@ def fields_cartesian(
614
618
r , theta , phi = self ._car_2_sph (_x , _y , _z )
615
619
_field_data = self .fields_spherical (r , theta , phi )
616
620
617
- Er , Etheta , Ephi = [_field_data [comp ].values \
618
- for comp in ["E_r" , "E_theta" , "E_phi" ]]
619
- Hr , Htheta , Hphi = [_field_data [comp ].values \
620
- for comp in ["H_r" , "H_theta" , "H_phi" ]]
621
+ Er , Etheta , Ephi = [
622
+ _field_data [comp ].values for comp in ["E_r" , "E_theta" , "E_phi" ]
623
+ ]
624
+ Hr , Htheta , Hphi = [
625
+ _field_data [comp ].values for comp in ["H_r" , "H_theta" , "H_phi" ]
626
+ ]
621
627
622
628
Ex_data [i , j , k ], Ey_data [i , j , k ], Ez_data [i , j , k ] = self ._sph_2_car_field (
623
629
Er , Etheta , Ephi , theta , phi
@@ -637,8 +643,7 @@ def fields_cartesian(
637
643
Hy = xr .DataArray (data = Hy_data , coords = coords , dims = dims )
638
644
Hz = xr .DataArray (data = Hz_data , coords = coords , dims = dims )
639
645
640
- field_data = xr .Dataset (
641
- {"Ex" : Ex , "Ey" : Ey , "Ez" : Ez , "Hx" : Hx , "Hy" : Hy , "Hz" : Hz })
646
+ field_data = xr .Dataset ({"Ex" : Ex , "Ey" : Ey , "Ez" : Ez , "Hx" : Hx , "Hy" : Hy , "Hz" : Hz })
642
647
643
648
return field_data
644
649
0 commit comments