0 HST::[jordan@quadrans] /data/lrp/orbit_model_updates 04-10 11:30> linear_model_projector.py modern_altitude_profile.rpt Table of differences between IJ-Solar05-Ref and CASM model projected to given dates: Projecting model 'apr16' to 19.090 = 14573.5 linear 14.7 - solar 18.8 = -4.2 Projecting model 'apr16' to 19.182 = 14665.5 linear 126.3 - solar 130.5 = -4.2 Projecting model 'apr16' to 19.275 = 14758.5 linear 231.3 - solar 235.5 = -4.2 Projecting model 'apr16' to 20.091 = 14939.5 linear 114.4 - solar 118.5 = -4.1 Projecting model 'apr16' to 20.153 = 15001.5 linear 64.4 - solar 68.4 = -4.0 Projecting model 'apr16' to 20.275 = 15123.5 linear 337.6 - solar 341.5 = -3.9 Projecting model 'apr16' to 21.091 = 15305.5 linear 214.1 - solar 217.5 = -3.4 Table of differences between IJ-Solar05-Ref and CASM model projected to given dates: Projecting model 'oct16' to 19.090 = 14573.5 linear 15.8 - solar 18.8 = -3.0 Projecting model 'oct16' to 19.182 = 14665.5 linear 127.2 - solar 130.5 = -3.3 Projecting model 'oct16' to 19.275 = 14758.5 linear 232.0 - solar 235.5 = -3.5 Projecting model 'oct16' to 20.091 = 14939.5 linear 114.6 - solar 118.5 = -3.9 Projecting model 'oct16' to 20.153 = 15001.5 linear 64.4 - solar 68.4 = -4.0 Projecting model 'oct16' to 20.275 = 15123.5 linear 337.3 - solar 341.5 = -4.1 Projecting model 'oct16' to 21.091 = 15305.5 linear 213.3 - solar 217.5 = -4.2 Table of differences between IJ-Solar05-Ref and CASM model projected to given dates: Projecting model 'nov16' to 19.090 = 14573.5 linear 15.8 - solar 18.8 = -3.1 Projecting model 'nov16' to 19.182 = 14665.5 linear 127.3 - solar 130.5 = -3.2 Projecting model 'nov16' to 19.275 = 14758.5 linear 232.2 - solar 235.5 = -3.3 Projecting model 'nov16' to 20.091 = 14939.5 linear 115.0 - solar 118.5 = -3.4 Projecting model 'nov16' to 20.153 = 15001.5 linear 65.0 - solar 68.4 = -3.4 Projecting model 'nov16' to 20.275 = 15123.5 linear 338.1 - solar 341.5 = -3.4 Projecting model 'nov16' to 21.091 = 15305.5 linear 214.3 - solar 217.5 = -3.2 Table of differences between IJ-Solar05-Ref and CASM model projected to given dates: Projecting model 'apr17' to 19.090 = 14573.5 linear 17.7 - solar 18.8 = -1.2 Projecting model 'apr17' to 19.182 = 14665.5 linear 129.3 - solar 130.5 = -1.2 Projecting model 'apr17' to 19.275 = 14758.5 linear 234.3 - solar 235.5 = -1.2 Projecting model 'apr17' to 20.091 = 14939.5 linear 117.4 - solar 118.5 = -1.1 Projecting model 'apr17' to 20.153 = 15001.5 linear 67.4 - solar 68.4 = -1.0 Projecting model 'apr17' to 20.275 = 15123.5 linear 340.6 - solar 341.5 = -0.8 Projecting model 'apr17' to 21.091 = 15305.5 linear 217.1 - solar 217.5 = -0.4 Table of differences between IJ-Solar05-Ref and CASM model projected to given dates: Projecting model 'apr18' to 19.090 = 14573.5 linear 18.1 - solar 18.8 = -0.7 Projecting model 'apr18' to 19.182 = 14665.5 linear 129.8 - solar 130.5 = -0.7 Projecting model 'apr18' to 19.275 = 14758.5 linear 234.7 - solar 235.5 = -0.8 Projecting model 'apr18' to 20.091 = 14939.5 linear 117.8 - solar 118.5 = -0.7 Projecting model 'apr18' to 20.153 = 15001.5 linear 67.8 - solar 68.4 = -0.6 Projecting model 'apr18' to 20.275 = 15123.5 linear 341.0 - solar 341.5 = -0.4 Projecting model 'apr18' to 21.091 = 15305.5 linear 217.5 - solar 217.5 = -0.0 Table of differences between IJ-Solar05-Ref and CASM model projected to given dates: Projecting model 'ppsep18' to 19.090 = 14573.5 linear 18.5 - solar 18.8 = -0.4 Projecting model 'ppsep18' to 19.182 = 14665.5 linear 130.1 - solar 130.5 = -0.4 Projecting model 'ppsep18' to 19.275 = 14758.5 linear 235.1 - solar 235.5 = -0.4 Projecting model 'ppsep18' to 20.091 = 14939.5 linear 118.2 - solar 118.5 = -0.3 Projecting model 'ppsep18' to 20.153 = 15001.5 linear 68.2 - solar 68.4 = -0.2 Projecting model 'ppsep18' to 20.275 = 15123.5 linear 341.4 - solar 341.5 = -0.0 Projecting model 'ppsep18' to 21.091 = 15305.5 linear 217.9 - solar 217.5 = 0.4 Table of differences between IJ-Solar05-Ref and CASM model projected to given dates: Projecting model 'mar19' to 19.090 = 14573.5 linear 19.8 - solar 18.8 = 0.9 Projecting model 'mar19' to 19.182 = 14665.5 linear 131.2 - solar 130.5 = 0.7 Projecting model 'mar19' to 19.275 = 14758.5 linear 235.9 - solar 235.5 = 0.4 Projecting model 'mar19' to 20.091 = 14939.5 linear 118.6 - solar 118.5 = 0.1 Projecting model 'mar19' to 20.153 = 15001.5 linear 68.4 - solar 68.4 = -0.0 Projecting model 'mar19' to 20.275 = 15123.5 linear 341.3 - solar 341.5 = -0.1 Projecting model 'mar19' to 21.091 = 15305.5 linear 217.3 - solar 217.5 = -0.2 Table of differences between IJ-Solar05-Ref and CASM model projected to given dates: Projecting model 'LeastSquares-IJ-SOLAR05-Ref' to 19.090 = 14573.5 linear 19.0 - solar 18.8 = 0.1 Projecting model 'LeastSquares-IJ-SOLAR05-Ref' to 19.182 = 14665.5 linear 130.5 - solar 130.5 = 0.0 Projecting model 'LeastSquares-IJ-SOLAR05-Ref' to 19.275 = 14758.5 linear 235.4 - solar 235.5 = -0.1 Projecting model 'LeastSquares-IJ-SOLAR05-Ref' to 20.091 = 14939.5 linear 118.4 - solar 118.5 = -0.1 Projecting model 'LeastSquares-IJ-SOLAR05-Ref' to 20.153 = 15001.5 linear 68.3 - solar 68.4 = -0.1 Projecting model 'LeastSquares-IJ-SOLAR05-Ref' to 20.275 = 15123.5 linear 341.4 - solar 341.5 = -0.0 Projecting model 'LeastSquares-IJ-SOLAR05-Ref' to 21.091 = 15305.5 linear 217.7 - solar 217.5 = 0.2 Details of the least-squares fit: numpy.polyfit(x, y, 1) = array([-6.61371391e+00, 9.60439428e+04]) = mod(360, 96043.9428257198) = asc_node 283.9428257198015 deg Ref time: 0 HJD Precession Rate: -6.613713913628986 deg/day Ascending Node: 283.9428257198015 deg where x = numpy.array(dates) y = numpy.array(ascn) To get an ascending node from this model at a 'reasonable' reference HJD, take the asc_node, precession_rate, and your intended reference HJD and plug it into... (asc_node + (reference_hjd * precession_rate ) ) % 360 e.g.: 19.090: ( 283.9428 + ( 14573.50 * -6.613713914 ) ) % 360 = 18.983 e.g.: 19.182: ( 283.9428 + ( 14665.50 * -6.613713914 ) ) % 360 = 130.521 e.g.: 19.275: ( 283.9428 + ( 14758.50 * -6.613713914 ) ) % 360 = 235.446 e.g.: 20.091: ( 283.9428 + ( 14939.50 * -6.613713914 ) ) % 360 = 118.364 e.g.: 20.153: ( 283.9428 + ( 15001.50 * -6.613713914 ) ) % 360 = 68.314 e.g.: 20.275: ( 283.9428 + ( 15123.50 * -6.613713914 ) ) % 360 = 341.440 e.g.: 21.091: ( 283.9428 + ( 15305.50 * -6.613713914 ) ) % 360 = 217.744 [14522, 14702, 14882, 15062, 15242]