pro test0630,rse,ve,rse366,ve366

pv=[[-0.236014486,0.880327681,0.381576410,-16994363,-3847358,-1667831],[-0.401394416,0.828404431,0.359062331,-15991627,-6507444,-2821516],[-0.554224354,0.750954032,0.325489267,-14501990,-8940977,-3875321],[-0.702472054,0.639247556,0.277064278,-12365596,-11307895,-4902750],[-0.814628866,0.516959543,0.224045151,-10009763,-13080140,-5670388],[-0.901788006,0.379050686,0.164264377,-7384226,-14433388,-6256747],[-0.956904262,0.245032539,0.106162899,-4837685,-15287397,-6628125],[-0.990534827,0.089294225,0.038643397,-1877496,-15775976,-6838853],[-0.994449629,-0.068900404,-0.029932074,1090519,-15788339,-6844415],[-0.964701280,-0.240285823,-0.104236354,4302781,-15278833,-6624329],[-0.907662275,-0.388679470,-0.168567079,7067051,-14326740,-6210290],[-0.824182331,-0.525521892,-0.227888579,9587034,-12980922,-5627869],[-0.716825624,-0.647092157,-0.280597434,11832835,-11272666,-4887144],[-0.588778745,-0.749860968,-0.325145531,13707376,-9236147,-4003500],[-0.444024285,-0.831061459,-0.360348047,15177133,-6969805,-3022151],[-0.444024285,-0.831061459,-0.360348047,15177133,-6969805,-3022151],[-0.104293055,-0.922596976,-0.400027773,16823787,-1662210,-720444],[0.064601937,-0.926152832,-0.401572976,16880495,952535,412427],[0.231757404,-0.903565151,-0.391778799,16466916,3557489,1543092],[0.392314074,-0.855391652,-0.370889833,15570479,6050506,2622858],[0.541784919,-0.783131751,-0.339567675,14256570,8371454,3629079],[0.688479885,-0.678065564,-0.294012444,12330733,10677315,4629618],[0.801383993,-0.562169009,-0.243767754,10199069,12446340,5395299],[0.891530565,-0.430250676,-0.186579731,7784468,13880995,6018289],[0.961143877,-0.270848004,-0.117466657,4829017,15013299,6508913],[0.995296391,-0.117524377,-0.050999689,1986086,15578922,6753454],[1.000615416,0.039316120,0.016997746,-934064,15716001,6814205],[0.976544788,0.195204462,0.084584731,-3870941,15380065,6667547],[0.923579290,0.345364742,0.149679772,-6696278,14579637,6320615],[0.843161997,0.485405544,0.210397368,-9356828,13356250,5791099],[0.725410796,0.622637164,0.269892416,-11986315,11513090,4990793],[0.595183828,0.727777696,0.315472602,-13995633,9464435,4103693],[0.446776108,0.810938349,0.351532271,-15614476,7118294,3086288],[0.284515970,0.869300180,0.376831343,-16750367,4521892,1959933],[0.113520427,0.900944658,0.390552329,-17364848,1793410,778382],[-0.061089105,0.904918175,0.392280244,-17465924,-1006624,-436525]]
help,pv
sz=SIZE(pv)
x=FLTARR(sz(2)) & y=x & z=x & rse=x
vx=x & vy=vx & vz=vx & ve=x & vs=x 
for i=0,sz(2)-1 do begin
    x(i)=pv(0,i)
    y(i)=pv(1,i)
    z(i)=pv(2,i)
    vx(i)=pv(3,i)
    vy(i)=pv(4,i)
    vz(i)=pv(5,i)
endfor
;print,x
;print,'y:',y
;print,'z:',z
;print,vx
vx=vx*1.0E-9
vy=vy*1.0E-9
vz=vz*1.0E-9
snpc0=rse & snpc1=rse & snpc2=rse
sdp=25.38
for i=0,sz(2)-1 do begin
    rse(i)=SQRT(x(i)*x(i)+y(i)*y(i)+z(i)*z(i))
    ve(i)=SQRT(vx(i)*vx(i)+vy(i)*vy(i)+vz(i)*vz(i))
    vs(i)=rse(i)*2*!pi/sdp                                   ;in AU/day
    ;ve=2*!pi*rse/365.
    snpc0(i)=sdp+(ve(i)*sdp)/vs(i)
    snpc1(i)=sdp+(ve(i)*27.2753)/vs(i)
    snpc2(i)=sdp/(1-ve(i)/vs(i))
endfor
print,'snpc2:',snpc2
print,'snpc1:',snpc1
psout,031,'Almanac36.ps'
!p.multi=[0,1,5]
plot,x,xrange=[0,sz(2)],xsty=1,yrange=[-1.0,1.0],line=1,ystyle=1,ytitle='x,y,z',title='Astronomic Alamanac 2008'
oplot,y,line=2
oplot,z,line=3
plot,rse,xrange=[0,sz(2)],xsty=1,yrange=[0.96,1.04],ystyle=1,ytitle='Rse in AU'
oplot,[0,sz(2)],[1,1]
plot,vs,xrange=[0,sz(2)],xsty=1,yrange=[0.23,0.26],ytitle="Sun's sidereal rotation speed in AU/day"
plot,ve,xrange=[0,sz(2)],xsty=1,yrange=[0.016,0.018],ystyle=1,ytitle="Earth's Ve in AU/day"
;plot,snp0,xrange=[0,sz(2)],yrange=[27.0,27.6],ystyle=1,ytitle='Synodic period'
plot,snpc0,xrange=[0,sz(2)],xsty=1,yrange=[27.15,27.36],ystyle=1,ytitle='Synodic period in days',line=1
;oplot,snpc0,line=1
oplot,snpc1,line=2
oplot,snpc2,line=3
oplot,[0,sz(2)],[27.2753,27.2753]
oplot,[0,sz(2)],[27.2000,27.2000]
oplot,[0,sz(2)],[27.3400,27.3400]
!p.multi=0
psdone
print,'Almanac36.ps'
rse366=fltarr(366)
ve366=rse366
rse366(0:10)=rse(0)
ve366(0:10)=ve(0)
rse366(11:20)=rse(1)
ve366(11:20)=ve(1)
rse366(21:30)=rse(2)
ve366(21:30)=ve(2)
rse366(31:40)=rse(3)
ve366(31:40)=ve(3)
rse366(41:50)=rse(4)
ve366(41:50)=ve(4)
rse366(51:61)=rse(5)
ve366(51:61)=ve(5)
rse366(62:71)=rse(6)
ve366(62:71)=ve(6)
rse366(72:81)=rse(7)
ve366(72:81)=ve(7)
rse366(82:91)=rse(8)
ve366(82:91)=ve(8)
rse366(92:101)=rse(9)
ve366(92:101)=ve(9)
rse366(102:112)=rse(10)
ve366(102:112)=ve(10)
rse366(113:122)=rse(11)
ve366(113:122)=ve(11)
rse366(123:131)=rse(12)
ve366(123:131)=ve(12)
rse366(132:141)=rse(13)
ve366(132:141)=ve(13)
rse366(142:151)=rse(14)
ve366(142:151)=ve(14)
rse366(152:162)=rse(15)
ve366(152:162)=ve(15)
rse366(163:172)=rse(16)
ve366(163:172)=ve(16)
rse366(173:182)=rse(17)
ve366(173:182)=ve(17)
rse366(183:192)=rse(18)
ve366(183:192)=ve(18)
rse366(193:202)=rse(19)
ve366(193:202)=ve(19)
rse366(203:213)=rse(20)
ve366(203:213)=ve(20)
rse366(214:223)=rse(21)
ve366(214:223)=ve(21)
rse366(224:233)=rse(22)
ve366(224:233)=ve(22)
rse366(234:243)=rse(23)
ve366(234:243)=ve(23)
rse366(244:253)=rse(24)
ve366(244:253)=ve(24)
rse366(254:263)=rse(25)
ve366(54:263)=ve(25)
rse366(264:274)=rse(26)
ve366(264:274)=ve(26)
rse366(275:284)=rse(27)
ve366(275:284)=ve(27)
rse366(285:294)=rse(28)
ve366(285:294)=ve(28)
rse366(295:304)=rse(29)
ve366(295:304)=ve(29)
rse366(305:314)=rse(30)
ve366(305:314)=ve(30)
rse366(315:324)=rse(31)
ve366(315:324)=ve(31)
rse366(325:334)=rse(32)
ve366(325:334)=ve(32)
rse366(335:344)=rse(33)
ve366(335:344)=ve(33)
rse366(345:354)=rse(34)
ve366(345:354)=ve(34)
rse366(355:365)=rse(35)
ve366(355:365)=ve(35)

end