1 files changed, 73 insertions, 0 deletions

diff --git a/src/cal/get_tle/check_tle.pro b/src/cal/get_tle/check_tle.pro
new file mode 100644
index 0000000..5388d9d
--- /dev/null
+++ b/src/cal/get_tle/check_tle.pro
@@ -0,0 +1,73 @@
+;
+fname='check_tle.dat'
+get_lun,unit0
+openr,unit0,fname
+day=fltarr(1000)
+apogee=fltarr(1000)
+perigee=fltarr(1000)
+semiax=fltarr(1000)
+raan=fltarr(1000)
+incl=fltarr(1000)
+mean_anom=fltarr(1000)
+arg_perig=fltarr(1000)
+eccen=fltarr(1000)
+d=0.0
+a=0.0
+p=0.0
+sem=0.0
+ra=0.0
+inc=0.0
+ma=0.0
+ap=0.0
+ec=0.0
+j=0
+while (not EOF(unit0)) do begin
+   readf,unit0,d,a,p,sem,ra,inc,ma,ap,ec
+   day[j]=d
+   apogee[j]=a
+   perigee[j]=p
+   semiax[j]=sem
+   raan[j]=ra
+   incl[j]=inc
+   mean_anom[j]=ma
+   arg_perig[j]=ap
+   eccen[j]=ec
+   j=j+1
+endwhile
+day=day(0:j-1)
+apogee=apogee(0:j-1)
+perigee=perigee(0:j-1)
+semiax=semiax(0:j-1)
+raan=raan(0:j-1)
+incl=incl(0:j-1)
+mean_anom=mean_anom(0:j-1)
+arg_perig=arg_perig(0:j-1)
+eccen=eccen(0:j-1)
+close,unit0
+free_lun,unit0
+!x.title='Day of Year 1999'
+!y.title='Arbitrary units'
+!p.title='FUSE Orbital Elements'
+!x.style=1
+!y.range=[0,1000]
+!x.range=[170,day(0)+10]
+!p.charsize=1.7
+!p.linestyle=0
+plot,day,apogee
+oplot,day,perigee
+oplot,day,semiax*100.0-714470.0+700.0
+oplot,day,raan*2.0
+oplot,day,incl*3500.0-87445.0+900.0
+oplot,day,mean_anom*2.0
+oplot,day,arg_perig*2.0
+oplot,day,eccen/2E-6
+xyouts,day(0),apogee(0),"Apogee",charsize=1.0
+xyouts,day(0),perigee(0),"Perigee",charsize=1.0
+xyouts,day(0),semiax(0)*100.0-714470.0+700.0,"Semimajor axis*100",charsize=1.0
+xyouts,day(0),raan(0)*2.0,"RAAN*2",charsize=1.0
+xyouts,day(0),incl(0)*3500.0-87445.0+900.0,"incl*3500",charsize=1.0
+xyouts,day(0),mean_anom(0)*2.0,"mean anomaly*2",charsize=1.0
+xyouts,day(0),arg_perig(0)*2.0,"argument of perigee*2",charsize=1.0
+xyouts,day(0),eccen(0)/2E-6,"eccen/2E-6",charsize=1.0
+
+end