;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;;; PROCEDURE Mainbody                 ;;;
;;; PROGRAMMERS EXTRAORDINAIRE         ;;;
;;; Darrin Casebeer &                  ;;;
;;; Kevin Tubbs                        ;;;
;;; Date: 6/01/98-8/08/98              ;;;
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;
;PRO Mainbody is a componet of Galaxy that is
;used to calculate new positions and print to a
;frame which can be used by Xinteranimate.
;
PRO MAINBODY,X,Y,Z,X1,Y1,Z1,SF2,X2,Y2,Z2,VX,VY,VZ,AX,AY,AZ,M1,TIME1,M2,NS,I

  XC=(M1*X1+M2*X2)/(M1+M2)
  YC=(M1*Y1+M2*Y2)/(M1+M2)
  ZC=(M1*Z1+M2*Z2)/(M1+M2)
  XX1=(X1-XC)
  YY1=(Y1-YC)
  ZZ1=(Z1-ZC)
  XX2=(X2-XC)
  YY2=(Y2-YC)
  ZZ2=(Z2-ZC)
  MinMaxRange=40
  SET_PLOT, 'x'
  SURFACE, dist(3,1200),/NODATA,/SAVE,XRANGE=[-MinMaxRange,MinMaxRange], $
  XSTYLE=1,XTITLE='X',YRANGE=[-MinMaxRange,MinMaxRange], YSTYLE=1,$
  YTITLE='Y',CHARSIZE=.5, ZRANGE=[-MinMaxRange,MinMaxRange],$
  ZSTYLE=1,ZTITLE='Z'
  PLOTS,XX1,YY1,ZZ1,PSYM=2,COL=76,SYMSIZE=.5,/T3D
  PLOTS,XX2,YY2,ZZ2,PSYM=2,COL=93,SYMSIZE=.5,/T3D
  FOR I=1,NS DO BEGIN
    XP=(X(I)-XC)
    YP=(Y(I)-YC)
    ZP=2*(Z(I)-ZC)
    PLOTS,XP,YP,ZP,PSYM=2,COL=76,SYMSIZE=.5,/T3D
  ENDFOR
END