#define GLOBALORIGIN
#include "view.h"
#undef GLOBALORIGIN

#include "../xgrafix/xgrafix.h"

FILE *InputDeck;
void start();

main(argc, argv)
int argc;
char *argv[];
{
  start(argc, argv);
  InitWindows(argc, argv);

  while(True)
  {
    if(theStepFlag || theRunFlag) 
    {
      t += dt;
      Refresh_Screen();
    }
    Event_Manager();  
  }
}   /* End of MAIN */

/***********************************************************/

void start(iargc, pargv)
int iargc;
char *pargv[];
{
  char aachar[80];
  int i, j, isp;
  
  nsp= 2;
  /* Allocate particle arrays */
  
  x  = (float **)malloc(nsp*sizeof(float *));
  vx = (float **)malloc(nsp*sizeof(float *));
  vy = (float **)malloc(nsp*sizeof(float *));
  vz = (float **)malloc(nsp*sizeof(float *));
  num= (float **)malloc(nsp*sizeof(float *));
  
  for(i=0; i<nsp; i++)
  {
    x[i]  = (float *)malloc(MAXLEN*sizeof(float));
    vx[i] = (float *)malloc(MAXLEN*sizeof(float));
    vy[i] = (float *)malloc(MAXLEN*sizeof(float));
    vz[i] = (float *)malloc(MAXLEN*sizeof(float));
    num[i]= (float *)malloc(MAXLEN*sizeof(float));
  }
  
  if (nsp && (!x[nsp-1] || !vx[nsp-1] || !vy[nsp-1] || !vz[nsp-1] || !num[nsp-1]))
  {
    puts("START: x, v malloc failed");
    exit(1);
  }
  
  /*********************************************************/
  /* Load in all species' particles, properly distributed. */
  /* Start from a dump file if given otherwise load the    */
  /* system with the initial distribution                  */

  if (iargc>1)
  {
    strcpy(aachar, pargv[1]);
    if(access(aachar, F_OK)) strcat(aachar, ".dmp");
  }
  Restore(aachar);

  puts(aachar);
  printf("nsp = %d, np1= %d, np2= %d\n\n", nsp, np[0], np[1]);

  for(i=0; i<nsp; i++)
    for(j=0; j<np[i]; j++) num[i][j] = (float)j;
}   /* End START */

/***************************************************************/