release: bump version to 0.29.2

frontend/public/earth/js/celestial.js

@@ -2,6 +2,7 @@ import * as THREE from "three";
 import * as Astronomy from "astronomy-engine";
 
 import { CELESTIAL_CONFIG, EARTH_CONFIG } from "./constants.js";
+import { latLonToVector3 } from "./utils.js";
 
 const textureLoader = new THREE.TextureLoader();
 const defaultSunDirection = new THREE.Vector3(1, 0.2, 0.4).normalize();
@@ -32,6 +33,59 @@ let runtimeFollowConfig = {
 };
 const scratchEuler = new THREE.Euler(0, 0, 0, "YXZ");
 
+function normalizeDegrees180(value) {
+  let normalized = value;
+  while (normalized <= -180) normalized += 360;
+  while (normalized > 180) normalized -= 360;
+  return normalized;
+}
+
+function computeGreenwichMeanSiderealDegrees(date) {
+  const jd = date.getTime() / 86400000 + 2440587.5;
+  const t = (jd - 2451545.0) / 36525.0;
+  const gmst =
+    280.46061837 +
+    360.98564736629 * (jd - 2451545.0) +
+    0.000387933 * t * t -
+    (t * t * t) / 38710000;
+  return THREE.MathUtils.euclideanModulo(gmst, 360);
+}
+
+function computeSubsolarLocalDirection(date) {
+  const vector = Astronomy.GeoVector(Astronomy.Body.Sun, date, false);
+  const radius = Math.sqrt(
+    vector.x * vector.x +
+    vector.y * vector.y +
+    vector.z * vector.z,
+  );
+  if (!Number.isFinite(radius) || radius === 0) {
+    return defaultSunDirection.clone();
+  }
+
+  const rightAscensionDeg = THREE.MathUtils.radToDeg(
+    Math.atan2(vector.y, vector.x),
+  );
+  const declinationDeg = THREE.MathUtils.radToDeg(
+    Math.asin(THREE.MathUtils.clamp(vector.z / radius, -1, 1)),
+  );
+  const gmstDeg = computeGreenwichMeanSiderealDegrees(date);
+  const subsolarLonDeg = normalizeDegrees180(rightAscensionDeg - gmstDeg);
+
+  return latLonToVector3(declinationDeg, subsolarLonDeg, 1).normalize();
+}
+
+function getPhysicalSunDirection(date = new Date()) {
+  const localSunDirection = computeSubsolarLocalDirection(date);
+  if (!linkedEarth) {
+    return localSunDirection;
+  }
+
+  return localSunDirection
+    .clone()
+    .applyQuaternion(linkedEarth.quaternion)
+    .normalize();
+}
+
 function getCelestialEuler() {
   const { x, y, z } = runtimeOrientationEuler;
   return new THREE.Euler(x, y, z, "YXZ");
@@ -279,13 +333,15 @@ function updateSpritePositions() {
 }
 
 function updateLighting() {
-  const calibratedSunDirection = applyCelestialOrientation(sunDirection);
+  const physicalSunDirection = getPhysicalSunDirection(
+    new Date(lastUpdatedAt || Date.now()),
+  );
 
   if (linkedSunLight) {
     linkedSunLight.color.setHex(CELESTIAL_CONFIG.sunLightColor);
     linkedSunLight.intensity = CELESTIAL_CONFIG.sunLightIntensity;
     linkedSunLight.position
-      .copy(calibratedSunDirection)
+      .copy(physicalSunDirection)
       .multiplyScalar(CELESTIAL_CONFIG.sunLightDistance);
   }
 
@@ -293,7 +349,7 @@ function updateLighting() {
     linkedBackLight.color.setHex(CELESTIAL_CONFIG.backLightColor);
     linkedBackLight.intensity = CELESTIAL_CONFIG.backLightIntensity;
     linkedBackLight.position
-      .copy(calibratedSunDirection)
+      .copy(physicalSunDirection)
       .multiplyScalar(-CELESTIAL_CONFIG.sunLightDistance * 0.7);
   }
 }
@@ -413,7 +469,7 @@ export function updateCelestialLayer(date = new Date(), camera = null) {
 }
 
 export function getSunDirection() {
-  return applyCelestialOrientation(sunDirection);
+  return getPhysicalSunDirection(new Date(lastUpdatedAt || Date.now()));
 }
 
 export function getMoonDirection() {