mirror/astronomy
1
0
Fork 0
mirror of https://github.com/cosinekitty/astronomy.git synced 2026-05-25 01:09:09 -04:00
Code Issues Packages Projects Releases Wiki Activity

Eliminating JD: fixed nutation()

Browse Source
This commit is contained in:
Don Cross
2019-04-15 09:56:09 -04:00
parent fe70afa483
commit a3f899e234
2 changed files with 8 additions and 8 deletions
Download Patch File Download Diff File Expand all files Collapse all files

8
generate/template/astronomy.js
View File

@@ -767,8 +767,8 @@ function terra(observer, st) {
};
}
function nutation(jd_tdb, direction, pos) {
const tilt = e_tilt(jd_tdb);
function nutation(time, direction, pos) {
const tilt = e_tilt(time.jd_tt);
const oblm = tilt.mobl * DEG2RAD;
const oblt = tilt.tobl * DEG2RAD;
const psi = tilt.dpsi * ASEC2RAD;
@@ -811,7 +811,7 @@ function geo_pos(time, observer) {
const tilt = e_tilt(time.jd_tt);
const gast = gmst + tilt.ee/3600;
const pos1 = terra(observer, gast).pos;
const pos2 = nutation(time.jd_tt, -1, pos1);
const pos2 = nutation(time, -1, pos1);
const pos3 = precession(time.jd_tt, pos2, T0);
return pos3;
}
@@ -994,7 +994,7 @@ Astronomy.SkyPos = function(gc_vector, observer) { // based on NOVAS place()
let j2000_radec = vector2radec(j2000_vector);
let pos7 = precession(T0, j2000_vector, gc_vector.t.jd_tt);
let ofdate_vector = nutation(gc_vector.t.jd_tt, 0, pos7);
let ofdate_vector = nutation(gc_vector.t, 0, pos7);
let ofdate_radec = vector2radec(ofdate_vector);
let sky = {

8
source/js/astronomy.js
View File

@@ -1457,8 +1457,8 @@ function terra(observer, st) {
};
}
function nutation(jd_tdb, direction, pos) {
const tilt = e_tilt(jd_tdb);
function nutation(time, direction, pos) {
const tilt = e_tilt(time.jd_tt);
const oblm = tilt.mobl * DEG2RAD;
const oblt = tilt.tobl * DEG2RAD;
const psi = tilt.dpsi * ASEC2RAD;
@@ -1501,7 +1501,7 @@ function geo_pos(time, observer) {
const tilt = e_tilt(time.jd_tt);
const gast = gmst + tilt.ee/3600;
const pos1 = terra(observer, gast).pos;
const pos2 = nutation(time.jd_tt, -1, pos1);
const pos2 = nutation(time, -1, pos1);
const pos3 = precession(time.jd_tt, pos2, T0);
return pos3;
}
@@ -1684,7 +1684,7 @@ Astronomy.SkyPos = function(gc_vector, observer) { // based on NOVAS place()
let j2000_radec = vector2radec(j2000_vector);
let pos7 = precession(T0, j2000_vector, gc_vector.t.jd_tt);
let ofdate_vector = nutation(gc_vector.t.jd_tt, 0, pos7);
let ofdate_vector = nutation(gc_vector.t, 0, pos7);
let ofdate_radec = vector2radec(ofdate_vector);
let sky = {