vz = vz[cut] M = M[cut] Mstar = Mstar[cut] Mgas = Mgas[cut] host = host[cut] fhires = fhires[cut] ngas = ngas[cut] nstar = nstar[cut] npart = npart[cut] Rvir = Rvir[cut] Vsig = Vsig[cut] Zs = Zs[cut] Ns = [] for theZ in Zs: N = SF.N_for_z(theZ) Ns.append(N) count = 0 jump_use = [] iSinX = [] iSinY = [] iSinZ = [] ShalX = [] ShalY = [] ShalZ = [] TheLen = len(Zs) DinX = []
print ID print 'all host' print host print 'at first, this halo is the following: ' print ID[0], Zs[0], nDM[0], nstar[0], ngas[0], host[0], fhires[0] print 'at the end, this halo becomes: ' print ID[-1], Zs[-1], nDM[-1], nstar[-1], ngas[-1], host[-1], fhires[-1] count = 0 type = -1 while (count < len(ID)): if (host[count] != -1): print 'this halo stopped being a main halo at ', count, Zs[ count], ' N = ', SF.N_for_z(Zs[count]) override = raw_input('do you agree?') if (override != 'n'): type = 0 print 'terminating' break count += 1 if (host[-1] != -1): type = 0 if (count > 0): count -= 1 good = True