Ejemplo n.º 1
0
nemo_main()
{
    int seed;

    alpha = getdparam("alpha");
    rcut = getdparam("rcut");
    mdisk = getdparam("mdisk");
    Qtoomre = getdparam("Qtoomre");
    gammas = getdparam("gamma");
    ndisk = getiparam("nbody");
    z0 = getdparam("z0");             /* was called epsi in mkbaredisk */
    cmode = getiparam("mode");
    zmode = getiparam("zmode");
    seed = init_xrandom(getparam("seed"));
    Qtab = getbparam("tab");
    inittables();
    makedisk();
    writesnap(getparam("out"), getparam("headline"));
}
Ejemplo n.º 2
0
int main(int argc, string argv[])
{
  float tmp1, tmp2;

  initparam(argv, defv);
  if (sscanf(getparam("mdisk"), "%f,%f", &tmp1, &tmp2) == 2) {
    mdisk1 = tmp1;
    mdisk2 = tmp2;
    if (sscanf(getparam("alpha"), "%f,%f", &tmp1, &tmp2) != 2)
      error("%s: must specify two alpha values\n", getargv0());
    alpha1 = tmp1;
    alpha2 = tmp2;
    if (sscanf(getparam("epsilon"), "%f,%f", &tmp1, &tmp2) != 2)
      error("%s: must specify two epsilon values\n", getargv0());
    epsilon1 = tmp1;
    epsilon2 = tmp2;
  } else {
    mdisk1 = getdparam("mdisk");
    mdisk2 = 0.0;
    alpha1 = getdparam("alpha");
    alpha2 = 12.0;			 /* nonzero value stops div by zero */
    epsilon1 = getdparam("epsilon");
    epsilon2 = -1.0;
  }
  zdisk = getdparam("zdisk");
  mu = getdparam("mu");
  r_mu = getdparam("r_mu");
  eta = getdparam("eta");
  rcut = getdparam("rcut");
  readgsp();
  setprof();
  layout_body(bodyfields, Precision, NDIM);
  ndisk = getiparam("ndisk");
  disk = (bodyptr) allocate(ndisk * SizeofBody);
  init_random(getiparam("seed"));
  if (getiparam("nlist") > 0)
    listdisk(getiparam("nlist"));
  makedisk();
  writemodel();
  return (0);
}
Ejemplo n.º 3
0
int main(int argc, string argv[])
{
  double alpha, rcut;
  int model;

  initparam(argv, defv);
  alpha = getdparam("alpha");
  rcut = getdparam("rcut");
  model = (alpha <= 0.0 ? -1 : getiparam("model"));
  ndisk = getiparam("ndisk");
  init_random(getiparam("seed"));
  readgsp();
  setprof(model, alpha, rcut);
  layout_body(bodyfields, Precision, NDIM);
  disk = (bodyptr) allocate(ndisk * SizeofBody);
  makedisk(getbparam("randspin"));
  if (! getbparam("randspin"))			// if spins not random
    bodyfields[2] = NULL;			// don't write AuxVec field
  writemodel();
  fflush(NULL);
  return 0;
}