@@ -292,13 +292,6 @@ def calculate_mag_drift_unit_vectors_ecef(latitude, longitude, altitude,
292292 ecef_x , ecef_y , ecef_z = trans .geodetic_to_ecef (latitude , longitude ,
293293 altitude )
294294
295- # This shouldn't have been pushed to the repo
296- # idx, = np.where(np.isnan(ecef_x))
297- # if len(idx) > 0:
298- # print("Encountered nan starting ecef locations ", ecef_x[idx],
299- # ecef_y[idx], ecef_z[idx], latitude[idx], longitude[idx],
300- # altitude[idx])
301-
302295 # Begin method calculation.
303296
304297 # Magnetic field at root location
@@ -358,13 +351,6 @@ def calculate_mag_drift_unit_vectors_ecef(latitude, longitude, altitude,
358351 ecef_input = True ,
359352 ** step_kwargs )
360353
361- # This shouldn't be here
362- # idx, = np.where(np.isnan(tzx))
363- # if len(idx) > 0:
364- # print("Encountered nan tzx values ", idx, ecef_x[idx],
365- # ecef_y[idx], ecef_z[idx], tzx[idx], tzy[idx], tzz[idx],
366- # latitude[idx], longitude[idx], altitude[idx])
367-
368354 if centered_diff :
369355 # Negative step
370356 ecef_xz2 , ecef_yz2 , ecef_zz2 = (ecef_x - step_size * tzx ,
@@ -435,15 +421,6 @@ def calculate_mag_drift_unit_vectors_ecef(latitude, longitude, altitude,
435421 & (loop_num + 1 >= min_loops ):
436422 # Reached terminating conditions
437423 repeat_flag = False
438- # print('Using position ', location_info(ecef_xz, ecef_yz, ecef_zz,
439- # datetimes,
440- # return_geodetic=True,
441- # ecef_input=True,
442- # **step_kwargs)[3:], loop_num)
443- # idx, = np.where(np.isnan(ecef_xz))
444- # if len(idx) > 0:
445- # print("Encountered nan apex locations ", ecef_xz[idx],
446- # ecef_yz[idx], ecef_zz[idx])
447424 else :
448425 # Store info into calculation vectors to refine next loop
449426 tzx , tzy , tzz = tzx2 , tzy2 , tzz2
@@ -691,7 +668,6 @@ def step_along_mag_unit_vector(x, y, z, date, direction, num_steps=1,
691668
692669 if direction == 'meridional' :
693670 centered_diff = True
694- # print("False")
695671 else :
696672 centered_diff = False
697673
0 commit comments