<s> import os <EOL> import sys <EOL> import math <EOL> import numpy as np <EOL> from random import * <EOL> from numpy import log10 <EOL> from optparse import OptionParser , OptionGroup <EOL> ID , D_GC , PVMAX , VMAX , RMAX , MTIDAL , RTIDAL , X , Y , Z , Vx , Vy , Vz , M300pc , M600pc = range ( <NUM_LIT:15> ) <EOL> Rad_ext , J_nfw , J_approx = range ( <NUM_LIT:8> , <NUM_LIT:11> ) <EOL> mm , rr , tt , pp , dd , jj = range ( <NUM_LIT:6> ) <EOL> x_earth = <NUM_LIT> <EOL> J_D = <NUM_LIT> <EOL> R_vir_VL = <NUM_LIT> <EOL> R_vir_AQ = <NUM_LIT> <EOL> km_per_kpc = <NUM_LIT> <EOL> omega_matter = <NUM_LIT> <EOL> omega_lambda = <NUM_LIT> <EOL> omega_baryon = <NUM_LIT> <EOL> omega_rad = <NUM_LIT> <EOL> hubble0 = <NUM_LIT> <EOL> ns_init = <NUM_LIT> <EOL> sigma_8 = <NUM_LIT> <EOL> w_de = - <NUM_LIT:1.0> <EOL> ID , X , Y , Z , VX , VY , VZ = range ( <NUM_LIT:7> ) <EOL> H_ID = <NUM_LIT:7> <EOL> N200a , M200a , R200a , s200a , N500a , M500a , R500a , s500a = range ( <NUM_LIT:7> , <NUM_LIT:15> ) <EOL> N200c , M200c , R200c , s200c , N500c , M500c , R500c , s500c = range ( <NUM_LIT:15> , <NUM_LIT> ) <EOL> Np = <NUM_LIT> ** <NUM_LIT:3> <EOL> Nh = <NUM_LIT> <EOL> scalefactor = <NUM_LIT:1.> <EOL> a = scalefactor <EOL> box = <NUM_LIT> <EOL> Ndm = <NUM_LIT> <EOL> h0 = <NUM_LIT> <EOL> om = <NUM_LIT> <EOL> Mpc2cm = <NUM_LIT> <EOL> H0_cgs = <NUM_LIT> <EOL> x_unit = box / Ndm <EOL> time_unit = <NUM_LIT:2> / ( <NUM_LIT:3> * H0_cgs * h0 * ( om ** <NUM_LIT> ) ) * a ** <NUM_LIT:2> <EOL> vel_unit = ( x_unit / time_unit ) / ( <NUM_LIT> ) <EOL> if __name__ == '<STR_LIT:__main__>' : <EOL> print '<STR_LIT>' <EOL> print ID , X , Y , Z , VX , VY , VZ , N200a , M200a , R200a , s200a , N500a , M500a , R500a , s500a , N200c , M200c , R200c , s200c , N500c , M500c , R500c , s500c <EOL> print '<STR_LIT>' , Np <EOL> print '<STR_LIT>' , Nh </s>
