@@ -323,24 +323,24 @@ def ephemeris(time, location, pressure=101325, temperature=12):
323
323
324
324
while np .max (abs (EccenAnom - E )) > 0.0001 :
325
325
E = EccenAnom
326
- EccenAnom = MeanAnom + np .degrees (Eccen ) * ( np .sin (np .radians (E ) ))
326
+ EccenAnom = MeanAnom + np .degrees (Eccen )* np .sin (np .radians (E ))
327
327
328
328
TrueAnom = (
329
329
2 * np .mod (np .degrees (np .arctan2 (((1 + Eccen ) / (1 - Eccen )) ** 0.5 *
330
330
np .tan (np .radians (EccenAnom ) / 2. ), 1 )), 360 ))
331
331
EcLon = np .mod (MlPerigee + TrueAnom , 360 ) - Abber
332
332
EcLonR = np .radians (EcLon )
333
- DecR = np .arcsin (np .sin (ObliquityR )* ( np .sin (EcLonR ) ))
333
+ DecR = np .arcsin (np .sin (ObliquityR )* np .sin (EcLonR ))
334
334
335
- RtAscen = np .degrees (np .arctan2 (np .cos (ObliquityR ) * (( np .sin (EcLonR )) ),
336
- np .cos (EcLonR )))
335
+ RtAscen = np .degrees (np .arctan2 (np .cos (ObliquityR )* np .sin (EcLonR ),
336
+ np .cos (EcLonR ) ))
337
337
338
338
HrAngle = LocAST - RtAscen
339
339
HrAngleR = np .radians (HrAngle )
340
340
HrAngle = HrAngle - (360 * ((abs (HrAngle ) > 180 )))
341
341
342
- SunAz = np .degrees (np .arctan2 (- 1 * np .sin (HrAngleR ), np . cos ( LatR ) *
343
- ( np .tan (DecR )) - np .sin (LatR )* ( np .cos (HrAngleR )) ))
342
+ SunAz = np .degrees (np .arctan2 (- np .sin (HrAngleR ),
343
+ np . cos ( LatR ) * np .tan (DecR ) - np .sin (LatR )* np .cos (HrAngleR ) ))
344
344
SunAz [SunAz < 0 ] += 360
345
345
346
346
SunEl = np .degrees (np .arcsin (
@@ -355,7 +355,7 @@ def ephemeris(time, location, pressure=101325, temperature=12):
355
355
Refract = pd .Series (0 , index = time_utc )
356
356
357
357
Refract [(Elevation > 5 ) & (Elevation <= 85 )] = (
358
- 58.1 / TanEl - 0.07 / (TanEl ** 3 ) + 8.6e-05 / (TanEl ** 5 ))
358
+ 58.1 / TanEl - 0.07 / (TanEl ** 3 ) + 8.6e-05 / (TanEl ** 5 ) )
359
359
360
360
Refract [(Elevation > - 0.575 ) & (Elevation <= 5 )] = ( Elevation *
361
361
(- 518.2 + Elevation * (103.4 + Elevation * (- 12.79 + Elevation * 0.711 ))) +
0 commit comments