Commit bf0d946c authored by NUNEZ Arturo's avatar NUNEZ Arturo

Automatic commit vendredi 13 juillet 2018, 16:30:01 (UTC+0200)

parent cbca894c
This diff is collapsed.
This diff is collapsed.
......@@ -1267,7 +1267,7 @@
},
{
"cell_type": "code",
"execution_count": 18,
"execution_count": 6,
"metadata": {
"collapsed": false,
"hide_input": false
......@@ -1279,7 +1279,11 @@
"text": [
"HALO B\n",
"loading Dark matter..\n",
"density [9693.97 9871.5625 9805.939 ]\n"
"x min = 0.00324946572073, max = 19881.4238281\n",
"y min = 0.000428158818977, max = 19881.4335938\n",
"z min = 0.00304498802871, max = 19881.4335938\n",
"density [9693.97 9871.5625 9805.939 ]\n",
"163.4765625\n"
]
}
],
......@@ -1289,15 +1293,47 @@
"#path = \"/media/arturo/ARTUROTECA/OUTPUTS/HaloB/output_00417\"\n",
"print simname\n",
"myDMO = wkbl.Galaxy_Hound(pathsim,getcen=False)\n",
"print \"x min = {0}, max = {1}\".format(myDMO.dm.pos3d[:,0].min(),myDMO.dm.pos3d[:,0].max())\n",
"print \"y min = {0}, max = {1}\".format(myDMO.dm.pos3d[:,1].min(),myDMO.dm.pos3d[:,1].max())\n",
"print \"z min = {0}, max = {1}\".format(myDMO.dm.pos3d[:,2].min(),myDMO.dm.pos3d[:,2].max())\n",
"ok,rho,_= CF.getDensity(np.array(myDMO.dm.pos3d.reshape(len(myDMO.dm.pos3d)*3),dtype=np.float32), myDMO.dm.mass)\n",
"centro_rho = myDMO.dm.pos3d[np.where(rho == rho.max())][0]\n",
"print \"density\",centro_rho\n",
"myDMO.center_shift(centro_rho)\n",
"myDMO.r_virial(600,n=25)\n",
"myDMO.r200\n",
"print myDMO.r200\n",
"myDMO.redefine(25)"
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"4.96e+11\n"
]
}
],
"source": [
"\n",
"print \"{0:.2e}\".format(np.sum(myDMO.dm.mass[myDMO.dm.r<myDMO.r200]))"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": []
},
{
"cell_type": "code",
"execution_count": 19,
This diff is collapsed.
This source diff could not be displayed because it is too large. You can view the blob instead.
This diff is collapsed.
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment