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")); }
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); }
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; }