; how many simulations
NSIMU     = 2l

; how many time steps per simulation
NTIMES    = 2l

; how many halos (just a guess....this should be the absolute MAX value!)
MAXHALOS  = 40000l

; how many bins in profile  (just a guess....this should be the absolute MAX value!)
MAXBINS   = 50l

; which files
DATApath     = '/Users/aknebe/DATA/DEBUG/LGR2Mpc_wind_497/'

halofileAHF  = [$
		; z=0
	        [  DATApath+'test-v229.LGR2Mpc_wind_217.z0.954.AHF_halos',$
              DATApath+'test-v229.LGR2Mpc_wind_218.z0.948.AHF_halos' ],$
		; z=1
	        [  DATApath+'test-v229.LGR2Mpc_wind_217.z0.954.AHF_halos',$
              DATApath+'test-v229.LGR2Mpc_wind_218.z0.948.AHF_halos' ]$
	]	

proffileAHF  = [$
		; z=0
	        [  DATApath+'test-v229.LGR2Mpc_wind_217.z0.954.AHF_profiles',$
              DATApath+'test-v229.LGR2Mpc_wind_218.z0.948.AHF_profiles' ],$
		; z=1
	        [  DATApath+'test-v229.LGR2Mpc_wind_217.z0.954.AHF_profiles',$
              DATApath+'test-v229.LGR2Mpc_wind_218.z0.948.AHF_profiles' ]$
	]	

simutag  = [ '217','218' ]
ztag     = ["z0.0","z1.0"]
zred     = [0.0,0.0]
line     = [0,1]
col      = [0,70]
colTab   = 2

pmass    = 1678784.125

; positions in *_halos file
NHFILE   = 60
NHEXTRA  = 4
NHCOLUMNS= NHFILE+NHEXTRA
Hnpart   = 0
Hnvpart  = 1
HXc      = 2
HYc      = 3
HZc      = 4
HVXc     = 5
HVYc     = 6
HVZc     = 7
HMvir    = 8
HRvir    = 9
HVmax    = 10
HRmax    = 11
HsigV    = 12
Hlambda  = 13
HLx      = 14
HLy      = 15
HLz      = 16
Ha       = 17
HEax     = 18
HEay     = 19
HEaz     = 20
Hb       = 21
HEbx     = 22
HEby     = 23
HEbz     = 24
Hc       = 25
HEcx     = 26
HEcy     = 27
HEcz     = 28
Hovdens  = 29
HRedge   = 30
Hnbins   = 31
HEkin    = 32
HEpot    = 33
Hmbp_off = 34
Hcom_off = 35
Hr2      = 36
HlambdaE = 37
; NHEXTRA:
Hc5      = (NHFILE-1)+1
HVslope  = (NHFILE-1)+2
HDslope  = (NHFILE-1)+3
HT       = (NHFILE-1)+4

; positions in *_profiles file
NPCOLUMNS= 24 
Prad     = 0
Pnpart   = 1
Pnvpart  = 2
Povdens  = 3
Pdens    = 4
Pvcirc   = 5
Psigv    = 6
PLx      = 7
PLy      = 8
PLz      = 9
Pa       = 10
PEax     = 11
PEay     = 12
PEaz     = 13
Pb       = 14
PEbx     = 15
PEby     = 16
PEbz     = 17
Pc       = 18
PEcx     = 19
PEcy     = 20
PEcz     = 21
PEkin    = 22
PEpot    = 23

; some useful constants/conversion factors
kpc2cm        = 3.086e21
kpc2km        = 3.086e16
km2cm         = 1e5
cm2km         = 1e-5
cm2kpc        = 1/kpc2cm
Gkm3Msun1s2   = 13.271801e10
GMpckm2Msun1s2= 4.3006485E-9
Gcm3Msun1s2   = Gkm3Msun1s2 * km2cm^3
a0            = 1.2e-8
Hubble        = 100.
Grav          = GMpckm2Msun1s2
