private XXLGroupProcess process; // for accessing useful config data
#endregion Fields
#region Constructors
/// <summary>
///
/// </summary>
/// <param name="_process"></param>
/// <param name="filename"></param>
/// <param name="_snapnum"></param>
public XXLGroupFile(XXLGroupProcess _process, string _pathname, int _ifile, short _snapnum, bool _hasVeldisp)
{
this.process = _process;
this.snapnum = _snapnum;
this.ifile = _ifile;
string filename = _process.GetSubTab(_pathname, _snapnum, ifile);
this.hasVelDisp = _hasVeldisp;
float redshift = process.globals.redshifts[snapnum];
bool hasVelDisp = process.hasVelDisp;
// no need for records/bytesRead, not in fortran format
using (BinaryReader reader = new BinaryReader(new FileStream(filename, FileMode.Open)))
{
numGroups = reader.ReadInt32();
totalGroups = reader.ReadInt64();
numIDs = reader.ReadInt32();
totalIDs = reader.ReadInt64();
numFiles = reader.ReadInt32();
numSubgroups = reader.ReadInt32();
totalSubgroups = reader.ReadInt64();
groups = new GroupInfo_MXXL[numGroups];
for (int i = 0; i < numGroups; i++)
{
groups[i] = new GroupInfo_MXXL();
groups[i].snapnum = snapnum;
groups[i].redshift = redshift;
}
subhalos = new SubhaloInfo_MXXL[numSubgroups];
for (int i = 0; i < numSubgroups; i++)
{
subhalos[i] = new SubhaloInfo_MXXL();
subhalos[i].snapnum = snapnum;
subhalos[i].redshift = redshift;
}
for (int i = 0; i < numGroups; i++)
groups[i].Length = reader.ReadInt32();
for (int i = 0; i < numGroups; i++)
groups[i].offset = reader.ReadInt32();
for (int i = 0; i < numGroups; i++)
groups[i].GroupNr = reader.ReadInt64();
for (int i = 0; i < numGroups; i++)
{
groups[i].cmx = reader.ReadSingle();
groups[i].cmy = reader.ReadSingle();
groups[i].cmz = reader.ReadSingle();
}
for (int i = 0; i < numGroups; i++)
{
groups[i].vx = reader.ReadSingle();
groups[i].vy = reader.ReadSingle();
groups[i].vz = reader.ReadSingle();
}
for (int i = 0; i < numGroups; i++)
{
groups[i].x = reader.ReadSingle();
groups[i].y = reader.ReadSingle();
groups[i].z = reader.ReadSingle();
groups[i].ix = process.GetZone(groups[i].x);
groups[i].iy = process.GetZone(groups[i].y);
groups[i].iz = process.GetZone(groups[i].z);
groups[i].phkey = process.GetPHKey(groups[i].x, groups[i].y, groups[i].z);
}
for (int i = 0; i < numGroups; i++)
groups[i].M_Mean200 = reader.ReadSingle();
for (int i = 0; i < numGroups; i++)
groups[i].M_Crit200 = reader.ReadSingle();
for (int i = 0; i < numGroups; i++)
groups[i].M_TopHat200 = reader.ReadSingle();
for (int i = 0; i < numGroups; i++)
groups[i].VelDisp = reader.ReadSingle();
if (hasVelDisp)
{
for (int i = 0; i < numGroups; i++)
groups[i].VelDisp_Mean200 = reader.ReadSingle();
for (int i = 0; i < numGroups; i++)
groups[i].VelDisp_Crit200 = reader.ReadSingle();
for (int i = 0; i < numGroups; i++)
groups[i].VelDisp_TopHat200 = reader.ReadSingle();
}
for (int i = 0; i < numGroups; i++)
groups[i].numSubs = reader.ReadInt32();
for (int i = 0; i < numGroups; i++)
groups[i].firstSubIndex = reader.ReadInt32();
// -------------------------------------
// now we're at the subhalo table start
// -------------------------------------
for (int i = 0; i < numSubgroups; i++)
subhalos[i].length = reader.ReadInt32();
for (int i = 0; i < numSubgroups; i++)
subhalos[i].offset = reader.ReadInt32();
for (int i = 0; i < numSubgroups; i++)
subhalos[i].GrNr = reader.ReadInt64();
for (int i = 0; i < numSubgroups; i++)
subhalos[i].SubNr = reader.ReadInt64();
for (int i = 0; i < numSubgroups; i++)
{
subhalos[i].x = reader.ReadSingle();
subhalos[i].y = reader.ReadSingle();
subhalos[i].z = reader.ReadSingle();
subhalos[i].ix = process.GetZone(subhalos[i].x);
subhalos[i].iy = process.GetZone(subhalos[i].y);
subhalos[i].iz = process.GetZone(subhalos[i].z);
subhalos[i].phkey = process.GetPHKey(subhalos[i].x, subhalos[i].y, subhalos[i].z);
}
for (int i = 0; i < numSubgroups; i++)
{
subhalos[i].vx = reader.ReadSingle();
subhalos[i].vy = reader.ReadSingle();
subhalos[i].vz = reader.ReadSingle();
}
for (int i = 0; i < numSubgroups; i++)
{
subhalos[i].cmx = reader.ReadSingle();
subhalos[i].cmy = reader.ReadSingle();
subhalos[i].cmz = reader.ReadSingle();
}
for (int i = 0; i < numSubgroups; i++)
{
subhalos[i].sx = reader.ReadSingle();
subhalos[i].sy = reader.ReadSingle();
subhalos[i].sz = reader.ReadSingle();
}
for (int i = 0; i < numSubgroups; i++)
subhalos[i].veldisp = reader.ReadSingle();
for (int i = 0; i < numSubgroups; i++)
subhalos[i].velMax = reader.ReadSingle();
for (int i = 0; i < numSubgroups; i++)
subhalos[i].velMaxRad = reader.ReadSingle();
for (int i = 0; i < numSubgroups; i++)
subhalos[i].hmradius = reader.ReadSingle();
for (int i = 0; i < numSubgroups; i++)
for(int k = 0; k < 6; k++)
subhalos[i].shape[k] = reader.ReadSingle();
for (int i = 0; i < numSubgroups; i++)
subhalos[i].bindingEnergy = reader.ReadSingle();
for (int i = 0; i < numSubgroups; i++)
subhalos[i].potentialEnergy = reader.ReadSingle();
for (int i = 0; i < numSubgroups; i++)
for (int k = 0; k < 9; k++)
subhalos[i].profile[k] = reader.ReadSingle();
}
}
private GroupProcess RefreshGroupsProcess()
{
GroupProcess process;
if ("milliMXXL".Equals(lbSimulation.SelectedItem.ToString())
|| "centiMXXL".Equals(lbSimulation.SelectedItem.ToString())
|| "Millennium-XXL".Equals(lbSimulation.SelectedItem.ToString()))// chkIsMXXL.Checked)
process = new XXLGroupProcess(globalParameters);
else if("Millennium-II".Equals(lbSimulation.SelectedItem.ToString()))
process = new GroupProcess(globalParameters);
else
{
DebugOut.PrintLine("Currently no support for millimil or Millennium FOF groups and subhalos.");
return null;
}
process.inPath = txtFOFSubhaloSourceDir.Text;
process.outPath = txtFOFSubhaloTargetDir.Text;
process.writeFOFs = chkWriteFOF.Checked;
process.writeSubHalos = chkWriteSubHalos.Checked;
process.writeSubhaloIDs = chkIDs.Checked;
process.hasVelDisp = chkDisp.Checked;
process.fofTable = txtFOFTable.Text;
process.subhaloTable = txtSubhaloTable.Text;
process.firstSnap = Int16.Parse(txtFirstGroupsSnapshot.Text);
process.lastSnap = Int16.Parse(txtLastGroupsSnapshot.Text);
string sim = lbSimulation.SelectedItem.ToString();
return process;
}
