Skyplot1001: jjh2000(jdat,par) - File Exchange

No BSD License

Highlights from
Skyplot1001

ErrMsg(ErrorString) function ErrMsg(ErrorString)
pickrd function pickrd
plotsky.m plotsky.m
A=pmatequ(T1,T2) function A=pmatequ(T1,T2)
AG=gms2alt(AH,AM,AS) function AG=gms2alt(AH,AM,AS)
LMST=lmst(mjdat,lambda) function LMST=lmst(mjdat,lambda)
MPN=pnmatrix(T0,T) function MPN=pnmatrix(T0,T)
N=moonaddn(N,COEFFN,P,Q,R,S,C...
PQR=gaussvec(lan,inc,aop) function PQR=gaussvec(lan,inc,aop)
[A,D,W,sense]=orient(PNr,Syst...
[ASTR,AH,AM,AS]=alt2gms(AG)
[DLAM,DS,GAM1C,SINPI]=addsol(... function [DLAM,DS,GAM1C,SINPI]=addsol(DLAM,DS,GAM1C,SINPI,COEFFL,COEFFS,COEFFG,COEFFP,P,Q,R,S,CO,SI)
[HorAlign]=horalign(x,y,stol) function [HorAlign]=horalign(x,y,stol)
[R,D,ELONG,PHI,K]=illum(x,y,z...
[R_EQU,FL]=shape(PNr)
[WD,Phi,k,stadium,mag]=physmo... function [WD,Phi,k,stadium,mag]=physmoon(T)
[X,Y]=moonterm(P,Q,R,S,CO,SI)
[ax,ay,az,LONG,LAT,DEC]=rotat...
[h,az]=equ2hor(dec,tau,phi)
[r,theta,phi]=cart2pol(x,y,z)
[rec,dec,r]=moonequ(T)
[rec,dec]=planequ(T)
[vx,vy,vz]=aberrat(T)
[x,y,vx,vy]=ellip(M,a,ecc)
[x,y,z,vx,vy,vz]=posvel(l,b,r...
[x,y,z]=nutequ(T,x,y,z)
[x,y,z]=orbecl(xx,yy,PQR)
[x,y,z]=pol2cart(r,theta,phi)
[x,y,z]=position(PNr,T)
[xP,yP,zP,xS,yS,zS,x,y,z,delt...
[xp,yp,zp]=precart(A,x,y,z)
[xt,yt,zt]=ecl2equ(T,x,y,z)
addtheo(C1,S1,C2,S2)
eccan=eccanom(man,ecc) function eccan=eccanom(man,ecc)
frac(x)
jjh2000(jdat,par) function jjh = jjh2000(jdat,par)
mag=bright(PNr,R,DELTA,PHI,DE... function mag=bright(PNr,R,DELTA,PHI,DEC,DLONG)
mjuldat(jahr,monat,tag,stunde... function mjdat = mjuldat(jahr,monat,tag,stunde)
perterm(T,C5,S5,C,S,IX,I,IT,D...
phi=atn2(y,x)
physsub(PNr,T)
poswi=posang(x,y,z,dx,dy,dz)
pplotsky(arg) function pplotsky(arg)
sinus=sinfrac(phi) function sinus=sinfrac(phi)
View all files

from Skyplot1001 by Markus Penzkofer
Virtual Planetarium showing stars, the planets, the sun and the moon.

jjh2000(jdat,par)

function jjh = jjh2000(jdat,par)

% function jjh = jjh2000(jdat,par)
%
% Julian Century in respect of J2000
% from Julian Date (par = 0) or
% from Modified Julian Date (par = 1)
% algorithms by Montenbruck/Pfleger: "Astronomy with the Personal Computer"
%
% 05.05.04 M.Penzkofer

if (par == 0)
	jjh = (jdat-2451545)/36525;
else
   jjh = (jdat-51544.5)/36525;
end

Highlights from Skyplot1001

Highlights from
Skyplot1001