Simulate the Solar System FOR REAL

{"name":"Simulate the Solar System [Brownsville]","description":"Data-driven approach","photo":"https://res.cloudinary.com/dzv7hjwfn/image/upload/rsoiws15wjy50gbdsiiv","view_count":1,"fork_count":0,"expires_at":null,"created_at":"2024-10-09T17:27:47.649884Z","updated_at":"2024-10-09T17:28:14.325719Z","user_id":60,"status":"PRIVATE","parent_id":3535,"id":3536,"user":{"photo":null,"email":"jacob@jippity.pro","recess_id":"34481dd2-2bdf-14f9-99da-131de59dc90f","created_at":"2024-02-26T18:31:54.473216Z","reset_token":null,"stripe_customer_id":"cus_RPbovGjLueUufv","seats":100,"username":"2718281828459045","id":60,"password_hash":"$2b$12$ZpelZmCi5wEA/VuXSdOJXeylQ6vOUks0sfwQxIKdjwBAeamaWIMOW","updated_at":"2024-02-26T18:31:54.473218Z","reset_token_expiry":null,"subscription_status":"active"},"code":{"script.js":"/*\n========================================\nTHE JAVASCRIPT SOLAR SYSTEM SIMULATION\n========================================\n\nPlease visit the URL to learn more:\nhttp://www.planetaryorbits.com/tutorial-javascript-orbit-simulation.html\n----------------------------------------\n*/\n\nlet sunImg, mercuryImg, venusImg, earthImg, marsImg, jupiterImg, saturnImg, uranusImg, neptuneImg;\n\nfunction preload() {\n sunImg = loadImage(\"@asset('sun')\");\n mercuryImg = loadImage(\"@asset('mercury')\");\n venusImg = loadImage(\"@asset('venus')\");\n earthImg = loadImage(\"@asset('earth')\");\n marsImg = loadImage(\"@asset('mars')\");\n jupiterImg = loadImage(\"@asset('jupiter')\");\n saturnImg = loadImage(\"@asset('saturn')\");\n uranusImg = loadImage(\"@asset('uranus')\");\n neptuneImg = loadImage(\"@asset('neptune')\");\n}\n\nfunction setup() {\n createCanvas(windowWidth, windowHeight);\n frameRate(60);\n T6.contextbackground = createGraphics(windowWidth, windowHeight);\n T6.contextforeground = createGraphics(windowWidth, windowHeight);\n init_T6();\n}\n\nfunction draw() {\n run_T6();\n image(T6.contextbackground, 0, 0);\n image(T6.contextforeground, 0, 0);\n}\n\nvar T6 = {\n contextbackground: null,\n contextforeground: null,\n width: window.innerWidth,\n height: window.innerHeight,\n fps: 60,\n now: null,\n then: Date.now(),\n interval: null,\n delta: null,\n nowTPS: null,\n thenTPS: Date.now(),\n avgTPSCount: 0,\n TPSCount: 0,\n deltaTPS: 0,\n nowFPS: null,\n thenFPS: Date.now(),\n avgFPSCount: 0,\n FPSCount: 0,\n deltaFPS: 0,\n julianCenturyInJulianDays: 36525,\n julianEpochJ2000: 2451545.0,\n julianDate: null,\n current: null,\n newDate: null,\n DAY: null,\n MONTH: null,\n YEAR: null,\n planetElements: [\n [0.38709927, 0.20563593, 7.00497902, 252.25032350, 77.45779628, 48.33076593],\n [0.72333566, 0.00677672, 3.39467605, 181.97909950, 131.60246718, 76.67984255],\n [1.00000261, 0.01671123, -0.00001531, 100.46457166, 102.93768193, 0.0],\n [1.52371034, 0.09339410, 1.84969142, -4.55343205, -23.94362959, 49.55953891],\n [5.20288700, 0.04838624, 1.30439695, 34.39644051, 14.72847983, 100.47390909],\n [9.53667594, 0.05386179, 2.48599187, 49.95424423, 92.59887831, 113.66242448],\n [19.18916464, 0.04725744, 0.77263783, 313.23810451, 170.95427630, 74.01692503],\n [30.06992276, 0.00859048, 1.77004347, -55.12002969, 44.96476227, 131.78422574]\n ],\n planetRates: [\n [0.00000037, 0.00001906, -0.00594749, 149472.67411175, 0.16047689, -0.1253408],\n [0.00000390, -0.00004107, -0.00078890, 58517.81538729, 0.00268329, -0.27769418],\n [0.00000562, -0.00004392, -0.01294668, 35999.37244981, 0.32327364, 0.0],\n [0.00001847, 0.00007882, -0.00813131, 19140.30268499, 0.44441088, -0.29257343],\n [-0.00011607, -0.00013253, -0.00183714, 3034.74612775, 0.21252668, 0.20469106],\n [-0.00125060, -0.00050991, 0.00193609, 1222.49362201, -0.41897216, -0.28867794],\n [-0.00196176, -0.00004397, -0.00242939, 428.48202785, 0.40805281, 0.04240589],\n [0.00026291, 0.00005105, 0.00035372, 218.45945325, -0.32241464, -0.00508664]\n ],\n orbitalElements: null,\n xMercury: null,\n yMercury: null,\n xVenus: null,\n yVenus: null,\n xEarth: null,\n yEarth: null,\n xMars: null,\n yMars: null,\n xJupiter: null,\n yJupiter: null,\n xSaturn: null,\n ySaturn: null,\n xUranus: null,\n yUranus: null,\n xNeptune: null,\n yNeptune: null,\n scale: 150,\n jupiterScaleDivider: 2.5,\n saturnScaleDivider: 3.5,\n uranusScaleDivider: 6.2,\n neptuneScaleDivider: 8.7\n}\n\nfunction init_T6() {\n T6.contextbackground = createGraphics(windowWidth, windowHeight);\n T6.contextforeground = createGraphics(windowWidth, windowHeight);\n\n T6.interval = 1000 / T6.fps;\n\n T6.current = new Date(Date.UTC(2000, 0, 2, 0, 0, 0));\n T6.DAY = T6.current.getDate();\n T6.MONTH = T6.current.getMonth() + 1;\n T6.YEAR = T6.current.getFullYear();\n T6.julianDate = getJulianDate_T6(T6.YEAR, T6.MONTH, T6.DAY);\n T6.T = (T6.julianDate - T6.julianEpochJ2000) / T6.julianCenturyInJulianDays;\n\n renderBackground_T6();\n run_T6();\n}\n\nfunction run_T6() {\n T6.now = Date.now();\n T6.delta = T6.now - T6.then;\n\n if (T6.delta > T6.interval) {\n T6.then = T6.now - (T6.delta % T6.interval);\n\n T6.FPSCount++;\n T6.nowFPS = Date.now();\n T6.deltaFPS = T6.nowFPS - T6.thenFPS;\n if (T6.deltaFPS > 1000) {\n T6.avgFPSCount = T6.FPSCount;\n T6.FPSCount = 0;\n T6.thenFPS = T6.nowFPS - (T6.deltaFPS % 1000);\n }\n\n updateDate_T6(1);\n\n T6.orbitalElements = plotPlanet_T6(T6.T, 0); T6.xMercury = T6.orbitalElements[0]; T6.yMercury = T6.orbitalElements[1];\n T6.orbitalElements = plotPlanet_T6(T6.T, 1); T6.xVenus = T6.orbitalElements[0]; T6.yVenus = T6.orbitalElements[1];\n T6.orbitalElements = plotPlanet_T6(T6.T, 2); T6.xEarth = T6.orbitalElements[0]; T6.yEarth = T6.orbitalElements[1];\n T6.orbitalElements = plotPlanet_T6(T6.T, 3); T6.xMars = T6.orbitalElements[0]; T6.yMars = T6.orbitalElements[1];\n T6.orbitalElements = plotPlanet_T6(T6.T, 4); T6.xJupiter = T6.orbitalElements[0]; T6.yJupiter = T6.orbitalElements[1];\n T6.orbitalElements = plotPlanet_T6(T6.T, 5); T6.xSaturn = T6.orbitalElements[0]; T6.ySaturn = T6.orbitalElements[1];\n T6.orbitalElements = plotPlanet_T6(T6.T, 6); T6.xUranus = T6.orbitalElements[0]; T6.yUranus = T6.orbitalElements[1];\n T6.orbitalElements = plotPlanet_T6(T6.T, 7); T6.xNeptune = T6.orbitalElements[0]; T6.yNeptune = T6.orbitalElements[1];\n }\n\n T6.TPSCount++;\n T6.nowTPS = Date.now();\n T6.deltaTPS = T6.nowTPS - T6.thenTPS;\n if (T6.deltaTPS > 1000) {\n T6.avgTPSCount = T6.TPSCount;\n T6.TPSCount = 0;\n T6.thenTPS = T6.nowTPS - (T6.deltaTPS % 1000);\n }\n\n renderForeground_T6();\n}\n\nfunction renderBackground_T6() {\n T6.contextbackground.background(9, 9, 9);\n\n if (sunImg) {\n T6.contextbackground.image(sunImg, T6.width / 2 - 25, T6.height / 2 - 25, 50, 50);\n } else {\n console.error(\"Sun image not loaded\");\n }\n}\n\nfunction renderForeground_T6() {\n T6.contextforeground.clear();\n\n if (!mercuryImg || !venusImg || !earthImg || !marsImg || !jupiterImg || !saturnImg || !uranusImg || !neptuneImg) {\n console.error(\"One or more planet images not loaded\");\n return;\n }\n\n // Helper function to draw images with preserved aspect ratio\n function drawImageWithAspectRatio(img, x, y, desiredWidth, desiredHeight) {\n const aspectRatio = img.width / img.height;\n let width = desiredWidth;\n let height = desiredHeight;\n\n // Adjust width and height to preserve aspect ratio\n if (width / height > aspectRatio) {\n width = height * aspectRatio;\n } else {\n height = width / aspectRatio;\n }\n\n T6.contextforeground.image(img, x - width / 2, y - height / 2, width, height);\n }\n\n drawImageWithAspectRatio(mercuryImg, (T6.width / 2) + (T6.xMercury * T6.scale), (T6.height / 2) - (T6.yMercury * T6.scale), 10, 10);\n drawImageWithAspectRatio(venusImg, (T6.width / 2) + (T6.xVenus * T6.scale), (T6.height / 2) - (T6.yVenus * T6.scale), 20, 20);\n drawImageWithAspectRatio(earthImg, (T6.width / 2) + (T6.xEarth * T6.scale), (T6.height / 2) - (T6.yEarth * T6.scale), 20, 20);\n drawImageWithAspectRatio(marsImg, (T6.width / 2) + (T6.xMars * T6.scale), (T6.height / 2) - (T6.yMars * T6.scale), 14, 14);\n drawImageWithAspectRatio(jupiterImg, (T6.width / 2) + (T6.xJupiter * T6.scale / T6.jupiterScaleDivider), (T6.height / 2) - (T6.yJupiter * T6.scale / T6.jupiterScaleDivider), 40, 40);\n drawImageWithAspectRatio(saturnImg, (T6.width / 2) + (T6.xSaturn * T6.scale / T6.saturnScaleDivider), (T6.height / 2) - (T6.ySaturn * T6.scale / T6.saturnScaleDivider), 36, 36);\n drawImageWithAspectRatio(uranusImg, (T6.width / 2) + (T6.xUranus * T6.scale / T6.uranusScaleDivider), (T6.height / 2) - (T6.yUranus * T6.scale / T6.uranusScaleDivider), 28, 28);\n drawImageWithAspectRatio(neptuneImg, (T6.width / 2) + (T6.xNeptune * T6.scale / T6.neptuneScaleDivider), (T6.height / 2) - (T6.yNeptune * T6.scale / T6.neptuneScaleDivider), 28, 28);\n\n T6.contextforeground.fill(105, 105, 105);\n T6.contextforeground.textSize(20);\n T6.contextforeground.textAlign(LEFT);\n T6.contextforeground.text(\"FPS = \" + T6.avgFPSCount, 10, 30);\n T6.contextforeground.textAlign(RIGHT);\n T6.contextforeground.text(\"TPS = \" + T6.avgTPSCount, T6.width - 10, 30);\n T6.contextforeground.textAlign(LEFT);\n T6.contextforeground.text(\"GREGORIAN DATE = \" + T6.YEAR + \".\" + T6.MONTH + \".\" + T6.DAY, 10, T6.height - 10);\n}\n\nfunction getJulianDate_T6(Year, Month, Day) {\n var inputDate = new Date(Year, Month, Math.floor(Day));\n var switchDate = new Date(\"1582\", \"10\", \"15\");\n\n var isGregorianDate = inputDate >= switchDate;\n\n if (Year < 0) {\n Year++;\n }\n\n if (Month == 1 || Month == 2) {\n Year = Year - 1;\n Month = Month + 12;\n }\n\n var A = Math.floor(Year / 100);\n var B = 2 - A + Math.floor(A / 4);\n\n if (!isGregorianDate) { B = 0; }\n\n return ((Math.floor(365.25 * Year)) + (Math.floor(30.6001 * (Month + 1))) + Day + 1720994.5 + B);\n}\n\nfunction updateDate_T6(increment) {\n T6.newDate = new Date(T6.current.getFullYear(), T6.current.getMonth(), T6.current.getDate() + increment);\n T6.current = T6.newDate;\n\n newdd = T6.newDate.getDate();\n newmm = T6.newDate.getMonth() + 1;\n newyyyy = T6.newDate.getFullYear();\n\n T6.YEAR = newyyyy;\n T6.MONTH = newmm;\n T6.DAY = newdd;\n\n T6.julianDate = getJulianDate_T6(T6.YEAR, T6.MONTH, T6.DAY);\n T6.T = (T6.julianDate - T6.julianEpochJ2000) / T6.julianCenturyInJulianDays;\n}\n\nfunction plotPlanet_T6(TGen, planetNumber) {\n aGen = T6.planetElements[planetNumber][0] + (T6.planetRates[planetNumber][0] * TGen);\n eGen = T6.planetElements[planetNumber][1] + (T6.planetRates[planetNumber][1] * TGen);\n iGen = T6.planetElements[planetNumber][2] + (T6.planetRates[planetNumber][2] * TGen);\n iGen = iGen % 360;\n WGen = T6.planetElements[planetNumber][5] + (T6.planetRates[planetNumber][5] * TGen);\n WGen = WGen % 360;\n wGen = T6.planetElements[planetNumber][4] + (T6.planetRates[planetNumber][4] * TGen);\n wGen = wGen % 360;\n if (wGen < 0) { wGen = 360 + wGen; }\n LGen = T6.planetElements[planetNumber][3] + (T6.planetRates[planetNumber][3] * TGen);\n LGen = LGen % 360;\n if (LGen < 0) { LGen = 360 + LGen; }\n\n MGen = LGen - (wGen);\n if (MGen < 0) { MGen = 360 + MGen; }\n\n EGen = EccAnom_T6(eGen, MGen, 6);\n\n trueAnomalyArgGen = (Math.sqrt((1 + eGen) / (1 - eGen))) * (Math.tan(toRadians_T6(EGen) / 2));\n\n K = Math.PI / 180.0;\n if (trueAnomalyArgGen < 0) {\n nGen = 2 * (Math.atan(trueAnomalyArgGen) / K + 180);\n }\n else {\n nGen = 2 * (Math.atan(trueAnomalyArgGen) / K)\n }\n\n rGen = aGen * (1 - (eGen * (Math.cos(toRadians_T6(EGen)))));\n\n xGen = rGen * (Math.cos(toRadians_T6(WGen)) * Math.cos(toRadians_T6(nGen + wGen - WGen)) - Math.sin(toRadians_T6(WGen)) * Math.sin(toRadians_T6(nGen + wGen - WGen)) * Math.cos(toRadians_T6(iGen)));\n yGen = rGen * (Math.sin(toRadians_T6(WGen)) * Math.cos(toRadians_T6(nGen + wGen - WGen)) + Math.cos(toRadians_T6(WGen)) * Math.sin(toRadians_T6(nGen + wGen - WGen)) * Math.cos(toRadians_T6(iGen)));\n zGen = rGen * (Math.sin(toRadians_T6(nGen + wGen - WGen)) * Math.sin(toRadians_T6(iGen)));\n\n return [xGen, yGen];\n}\n\n//TAKEN FROM: http://www.jgiesen.de/kepler/kepler.html\n//CREDIT: Juergen Giesen\n//Used to solve for E\nfunction EccAnom_T6(ec, m, dp) {\n // arguments:\n // ec=eccentricity, m=mean anomaly,\n // dp=number of decimal places\n\n var pi = Math.PI, K = pi / 180.0;\n var maxIter = 30, i = 0;\n var delta = Math.pow(10, -dp);\n var E, F;\n\n m = m / 360.0;\n m = 2.0 * pi * (m - Math.floor(m));\n\n if (ec < 0.8) E = m; else E = pi;\n\n F = E - ec * Math.sin(m) - m;\n\n while ((Math.abs(F) > delta) && (i < maxIter)) {\n E = E - F / (1.0 - ec * Math.cos(E));\n F = E - ec * Math.sin(E) - m;\n i = i + 1;\n }\n\n E = E / K;\n\n return Math.round(E * Math.pow(10, dp)) / Math.pow(10, dp);\n}\n\nfunction toRadians_T6(deg) {\n return deg * (Math.PI / 180);\n}\n\nfunction round_T6(value, decimals) {\n return Number(Math.round(value + 'e' + decimals) + 'e-' + decimals);\n}","style.css":"body {\n color:#c32de05e\n}","index.html":"\n\n\n\n \n \n Simulate the Solar System FOR REAL\n \n\n\n\n \n

\n\n \n

\n \n

\n\n

\n \n \n\n\n"},"chat_history":[],"assets":[{"name":"sun","cloudinary_id":"m2fnyrzanyad18imd8nr","resource_type":"image","id":24046,"project_id":3536,"created_at":"2024-10-09T17:27:47.682495Z","deleted":false},{"name":"neptune","cloudinary_id":"zun3dpkcjgdgd4jwqwc1","resource_type":"image","id":24045,"project_id":3536,"created_at":"2024-10-09T17:27:47.681166Z","deleted":false},{"name":"uranus","cloudinary_id":"zkpxbh9u1z2gqqd50rku","resource_type":"image","id":24044,"project_id":3536,"created_at":"2024-10-09T17:27:47.679781Z","deleted":false},{"name":"saturn","cloudinary_id":"cp7elsfrletlgcrpc2mw","resource_type":"image","id":24043,"project_id":3536,"created_at":"2024-10-09T17:27:47.678253Z","deleted":false},{"name":"jupiter","cloudinary_id":"knztmajsviml6ydjrjbx","resource_type":"image","id":24042,"project_id":3536,"created_at":"2024-10-09T17:27:47.676720Z","deleted":false},{"name":"mars","cloudinary_id":"dyzkeopx9xzgqdzwhhun","resource_type":"image","id":24041,"project_id":3536,"created_at":"2024-10-09T17:27:47.675027Z","deleted":false},{"name":"earth","cloudinary_id":"rsoiws15wjy50gbdsiiv","resource_type":"image","id":24040,"project_id":3536,"created_at":"2024-10-09T17:27:47.673061Z","deleted":false},{"name":"venus","cloudinary_id":"i6exysrb54inn0dbzz5e","resource_type":"image","id":24039,"project_id":3536,"created_at":"2024-10-09T17:27:47.670975Z","deleted":false},{"name":"mercury","cloudinary_id":"v9monl3gqcd7jdhzewxo","resource_type":"image","id":24038,"project_id":3536,"created_at":"2024-10-09T17:27:47.657957Z","deleted":false}],"messages_sent":0,"comments":[],"steps":[]}