Commit 24e430d1 authored by Rainer Weinberger's avatar Rainer Weinberger

added example cosmo_box_star_formation_3d and TREECOOL files (Puchwein et al....

added example cosmo_box_star_formation_3d and TREECOOL files (Puchwein et al. 2019MNRAS.485...47 and Faucher-Giguere et al. 2009ApJ...703.1416) in directory ./data
parent 73af73d4
This diff is collapsed.
0.0 3.76244e-14 2.0821e-14 1.12165e-16 2.4774e-25 2.21352e-25 5.0084e-27
0.005 3.83213e-14 2.12572e-14 1.14812e-16 2.52581e-25 2.26035e-25 5.12205e-27
0.01 3.90303e-14 2.17025e-14 1.17562e-16 2.57518e-25 2.30815e-25 5.23959e-27
0.015 4.0129e-14 2.23918e-14 1.21829e-16 2.65152e-25 2.38216e-25 5.42173e-27
0.02 4.16007e-14 2.33198e-14 1.27598e-16 2.75388e-25 2.48203e-25 5.66713e-27
0.025 4.31222e-14 2.42879e-14 1.33706e-16 2.86025e-25 2.58638e-25 5.92543e-27
0.03 4.47007e-14 2.52975e-14 1.40176e-16 2.97072e-25 2.69539e-25 6.19738e-27
0.035 4.63453e-14 2.63568e-14 1.47032e-16 3.086e-25 2.81009e-25 6.48397e-27
0.04 4.80462e-14 2.74619e-14 1.54289e-16 3.20577e-25 2.92999e-25 6.78559e-27
0.045 4.98113e-14 2.86143e-14 1.61972e-16 3.33015e-25 3.05526e-25 7.10308e-27
0.05 5.1649e-14 2.98232e-14 1.70114e-16 3.45991e-25 3.18705e-25 7.43765e-27
0.055 5.35503e-14 3.10841e-14 1.78734e-16 3.59472e-25 3.32485e-25 7.78983e-27
0.06 5.55239e-14 3.23989e-14 1.87857e-16 3.73468e-25 3.46884e-25 8.16047e-27
0.065 5.75769e-14 3.37776e-14 1.97519e-16 3.88063e-25 3.6203e-25 8.5509e-27
0.07 5.97016e-14 3.52155e-14 2.07747e-16 4.03223e-25 3.77865e-25 8.96189e-27
0.075 6.19076e-14 3.67145e-14 2.18571e-16 4.18957e-25 3.94413e-25 9.39441e-27
0.08 6.42009e-14 3.82859e-14 2.30028e-16 4.35363e-25 4.11817e-25 9.84984e-27
0.085 6.65751e-14 3.99246e-14 2.42151e-16 4.52402e-25 4.3001e-25 1.03291e-26
0.09 6.9041e-14 4.16331e-14 2.54978e-16 4.70085e-25 4.4902e-25 1.08334e-26
0.095 7.16025e-14 4.34246e-14 2.68549e-16 4.88515e-25 4.69012e-25 1.13642e-26
0.1 7.42548e-14 4.52916e-14 2.82897e-16 5.07649e-25 4.89903e-25 1.19225e-26
0.105 7.701e-14 4.72369e-14 2.98071e-16 5.27498e-25 5.11723e-25 1.25097e-26
0.11 7.98695e-14 4.9276e-14 3.14118e-16 5.48177e-25 5.34672e-25 1.31276e-26
0.115 8.28307e-14 5.13987e-14 3.31072e-16 5.69631e-25 5.58637e-25 1.37771e-26
0.12 8.59068e-14 5.36075e-14 3.48986e-16 5.91873e-25 5.83653e-25 1.44598e-26
0.125 8.90972e-14 5.59217e-14 3.67918e-16 6.15044e-25 6.09961e-25 1.51778e-26
0.13 9.24012e-14 5.83284e-14 3.87904e-16 6.39069e-25 6.37406e-25 1.5932e-26
0.135 9.58331e-14 6.08297e-14 4.09e-16 6.6396e-25 6.66022e-25 1.67241e-26
0.14 9.94039e-14 6.34491e-14 4.31274e-16 6.89882e-25 6.96109e-25 1.75564e-26
0.145 1.03087e-13 6.617e-14 4.54764e-16 7.16738e-25 7.27457e-25 1.843e-26
0.15 1.06911e-13 6.89941e-14 4.79529e-16 7.44539e-25 7.60104e-25 1.93466e-26
0.155 1.10885e-13 7.19495e-14 5.05647e-16 7.73474e-25 7.94415e-25 2.03088e-26
0.16 1.14984e-13 7.50163e-14 5.33157e-16 8.03432e-25 8.30139e-25 2.13177e-26
0.165 1.19226e-13 7.81946e-14 5.62124e-16 8.34412e-25 8.67296e-25 2.2375e-26
0.17 1.23641e-13 8.15181e-14 5.9263e-16 8.66638e-25 9.06333e-25 2.34836e-26
0.175 1.28194e-13 8.49624e-14 6.24716e-16 8.99975e-25 9.46935e-25 2.46445e-26
0.18 1.32903e-13 8.85268e-14 6.58451e-16 9.34417e-25 9.89109e-25 2.58595e-26
0.185 1.37801e-13 9.22523e-14 6.93928e-16 9.70214e-25 1.03339e-24 2.71319e-26
0.19 1.42848e-13 9.61096e-14 7.31182e-16 1.00721e-24 1.07938e-24 2.84623e-26
0.195 1.48065e-13 1.00098e-13 7.70282e-16 1.04539e-24 1.12709e-24 2.98526e-26
0.2 1.5348e-13 1.04261e-13 8.11311e-16 1.08502e-24 1.1771e-24 3.13056e-26
0.205 1.59057e-13 1.08567e-13 8.54314e-16 1.12592e-24 1.22898e-24 3.28223e-26
0.21 1.64817e-13 1.13014e-13 8.99354e-16 1.16809e-24 1.28273e-24 3.44044e-26
0.215 1.70784e-13 1.17645e-13 9.46514e-16 1.21178e-24 1.33896e-24 3.60544e-26
0.22 1.76923e-13 1.22429e-13 9.95835e-16 1.25682e-24 1.39721e-24 3.77733e-26
0.225 1.83257e-13 1.27358e-13 1.04737e-15 1.30318e-24 1.45746e-24 3.95624e-26
0.23 1.89805e-13 1.32482e-13 1.1012e-15 1.35113e-24 1.52037e-24 4.14241e-26
0.235 1.96533e-13 1.37766e-13 1.15735e-15 1.4005e-24 1.58543e-24 4.33589e-26
0.24 2.03467e-13 1.432e-13 1.21586e-15 1.45124e-24 1.65258e-24 4.53676e-26
0.245 2.10616e-13 1.48839e-13 1.27681e-15 1.50361e-24 1.72257e-24 4.74524e-26
0.25 2.17951e-13 1.54642e-13 1.3402e-15 1.55741e-24 1.79481e-24 4.96129e-26
0.255 2.25498e-13 1.60595e-13 1.40606e-15 1.61262e-24 1.86918e-24 5.18497e-26
0.26 2.33258e-13 1.66759e-13 1.47444e-15 1.66946e-24 1.94657e-24 5.4164e-26
0.265 2.41203e-13 1.73085e-13 1.54532e-15 1.72774e-24 2.02624e-24 5.65548e-26
0.27 2.49362e-13 1.79553e-13 1.61871e-15 1.78739e-24 2.10802e-24 5.90216e-26
0.275 2.57724e-13 1.86235e-13 1.69462e-15 1.84866e-24 2.19296e-24 6.15648e-26
0.28 2.66267e-13 1.93073e-13 1.77301e-15 1.91132e-24 2.28012e-24 6.41825e-26
0.285 2.75018e-13 2.00039e-13 1.85382e-15 1.97526e-24 2.36924e-24 6.68724e-26
0.29 2.83999e-13 2.07216e-13 1.93707e-15 2.04084e-24 2.46155e-24 6.96343e-26
0.295 2.93102e-13 2.14538e-13 2.02262e-15 2.10762e-24 2.55596e-24 7.24643e-26
0.3 3.02401e-13 2.21968e-13 2.11041e-15 2.17554e-24 2.6521e-24 7.5359e-26
0.305 3.11903e-13 2.29599e-13 2.20036e-15 2.24497e-24 2.75142e-24 7.83164e-26
0.31 3.21502e-13 2.37362e-13 2.29232e-15 2.31542e-24 2.85271e-24 8.13307e-26
0.315 3.3124e-13 2.45206e-13 2.38612e-15 2.38675e-24 2.95544e-24 8.43965e-26
0.32 3.41165e-13 2.53235e-13 2.48163e-15 2.45942e-24 3.06128e-24 8.75097e-26
0.325 3.5116e-13 2.61371e-13 2.57862e-15 2.53289e-24 3.16882e-24 9.06625e-26
0.33 3.61251e-13 2.69554e-13 2.67684e-15 2.60693e-24 3.27734e-24 9.38466e-26
0.335 3.71488e-13 2.77895e-13 2.77606e-15 2.68202e-24 3.38876e-24 9.70556e-26
0.34 3.81738e-13 2.86311e-13 2.87595e-15 2.75755e-24 3.50146e-24 1.00278e-25
0.345 3.92027e-13 2.94728e-13 2.97615e-15 2.83325e-24 3.61454e-24 1.03504e-25
0.35 4.02401e-13 3.03264e-13 3.07631e-15 2.90959e-24 3.73013e-24 1.06722e-25
0.355 4.12718e-13 3.11833e-13 3.17597e-15 2.98591e-24 3.84645e-24 1.09919e-25
0.36 4.22998e-13 3.20347e-13 3.27461e-15 3.06186e-24 3.96235e-24 1.13079e-25
0.365 4.33286e-13 3.28932e-13 3.37175e-15 3.13794e-24 4.08026e-24 1.16188e-25
0.37 4.43424e-13 3.37495e-13 3.46672e-15 3.21339e-24 4.19816e-24 1.19227e-25
0.375 4.53438e-13 3.45939e-13 3.55884e-15 3.28784e-24 4.31476e-24 1.22175e-25
0.38 4.63352e-13 3.5439e-13 3.64744e-15 3.36172e-24 4.43261e-24 1.25015e-25
0.385 4.73011e-13 3.62758e-13 3.73165e-15 3.43427e-24 4.54961e-24 1.27721e-25
0.39 4.82441e-13 3.7093e-13 3.81065e-15 3.50505e-24 4.66422e-24 1.3027e-25
0.395 4.91671e-13 3.79047e-13 3.88427e-15 3.57461e-24 4.77926e-24 1.32656e-25
0.4 5.00591e-13 3.87039e-13 3.95257e-15 3.64244e-24 4.8928e-24 1.34881e-25
0.405 5.09267e-13 3.948e-13 4.01573e-15 3.70831e-24 5.00329e-24 1.36945e-25
0.41 5.17733e-13 4.02487e-13 4.07414e-15 3.7729e-24 5.11399e-24 1.3886e-25
0.415 5.25902e-13 4.10041e-13 4.1282e-15 3.83579e-24 5.22289e-24 1.40634e-25
0.42 5.33842e-13 4.17334e-13 4.17807e-15 3.89669e-24 5.32802e-24 1.42269e-25
0.425 5.41633e-13 4.2453e-13 4.22358e-15 3.95639e-24 5.43307e-24 1.4376e-25
0.43 5.48998e-13 4.31557e-13 4.2644e-15 4.01381e-24 5.53562e-24 1.45096e-25
0.435 5.56077e-13 4.38264e-13 4.30025e-15 4.06878e-24 5.63337e-24 1.46267e-25
0.44 5.62935e-13 4.44821e-13 4.33092e-15 4.12207e-24 5.73032e-24 1.47268e-25
0.445 5.69306e-13 4.51181e-13 4.35615e-15 4.17272e-24 5.82448e-24 1.4809e-25
0.45 5.75296e-13 4.57167e-13 4.37567e-15 4.22035e-24 5.91289e-24 1.48724e-25
0.455 5.80996e-13 4.6296e-13 4.38935e-15 4.26588e-24 6.00002e-24 1.49166e-25
0.46 5.86202e-13 4.68531e-13 4.39696e-15 4.30856e-24 6.08407e-24 1.49408e-25
0.465 5.90973e-13 4.73684e-13 4.3983e-15 4.34784e-24 6.16146e-24 1.49444e-25
0.47 5.95411e-13 4.78604e-13 4.3933e-15 4.38467e-24 6.23711e-24 1.49273e-25
0.475 5.99308e-13 4.83282e-13 4.3818e-15 4.41838e-24 6.3094e-24 1.48888e-25
0.48 6.02729e-13 4.87503e-13 4.36368e-15 4.44836e-24 6.37419e-24 1.48285e-25
0.485 6.05783e-13 4.9146e-13 4.33895e-15 4.47566e-24 6.43689e-24 1.47466e-25
0.49 6.08265e-13 4.95162e-13 4.30754e-15 4.49965e-24 6.49611e-24 1.46428e-25
0.495 6.10246e-13 4.98382e-13 4.26941e-15 4.5197e-24 6.54716e-24 1.45169e-25
0.5 6.11843e-13 5.01314e-13 4.22468e-15 4.53692e-24 6.59584e-24 1.43694e-25
0.505 6.12855e-13 5.03994e-13 4.17338e-15 4.55078e-24 6.64119e-24 1.42004e-25
0.51 6.13357e-13 5.06176e-13 4.11558e-15 4.5606e-24 6.67784e-24 1.401e-25
0.515 6.13477e-13 5.08061e-13 4.05149e-15 4.56758e-24 6.71202e-24 1.37989e-25
0.52 6.13017e-13 5.09709e-13 3.98125e-15 4.57128e-24 6.74328e-24 1.35676e-25
0.525 6.1206e-13 5.10859e-13 3.90504e-15 4.57099e-24 6.76553e-24 1.33166e-25
0.53 6.10742e-13 5.11711e-13 3.82314e-15 4.56795e-24 6.78531e-24 1.30469e-25
0.535 6.0887e-13 5.12361e-13 3.73584e-15 4.56188e-24 6.80292e-24 1.27593e-25
0.54 6.0653e-13 5.12524e-13 3.64337e-15 4.55199e-24 6.81138e-24 1.24547e-25
0.545 6.03867e-13 5.12396e-13 3.54613e-15 4.53954e-24 6.81745e-24 1.21341e-25
0.55 6.00695e-13 5.12115e-13 3.44448e-15 4.52444e-24 6.82238e-24 1.17989e-25
0.555 5.97115e-13 5.11372e-13 3.33879e-15 4.50585e-24 6.81826e-24 1.14501e-25
0.56 5.93252e-13 5.10365e-13 3.2295e-15 4.485e-24 6.81201e-24 1.10893e-25
0.565 5.88951e-13 5.09259e-13 3.11706e-15 4.46201e-24 6.80577e-24 1.07177e-25
0.57 5.84315e-13 5.0772e-13 3.00188e-15 4.43592e-24 6.79043e-24 1.03366e-25
0.575 5.79553e-13 5.05921e-13 2.88443e-15 4.40814e-24 6.7724e-24 9.94779e-26
0.58 5.74327e-13 5.04102e-13 2.76522e-15 4.37853e-24 6.75586e-24 9.55269e-26
0.585 5.6885e-13 5.01885e-13 2.64469e-15 4.3463e-24 6.73036e-24 9.15271e-26
0.59 5.63301e-13 4.99419e-13 2.5233e-15 4.31269e-24 6.70185e-24 8.7495e-26
0.595 5.57385e-13 4.97035e-13 2.40159e-15 4.27802e-24 6.67718e-24 8.34469e-26
0.6 5.5127e-13 4.94306e-13 2.27998e-15 4.24119e-24 6.64424e-24 7.93969e-26
0.605 5.4511e-13 4.91369e-13 2.15895e-15 4.20332e-24 6.60877e-24 7.53611e-26
0.61 5.3875e-13 4.88596e-13 2.03896e-15 4.16532e-24 6.57919e-24 7.13549e-26
0.615 5.32273e-13 4.85527e-13 1.92042e-15 4.12568e-24 6.54173e-24 6.73916e-26
0.62 5.25828e-13 4.82286e-13 1.80373e-15 4.08546e-24 6.50194e-24 6.34853e-26
0.625 5.19261e-13 4.79277e-13 1.68929e-15 4.04572e-24 6.46977e-24 5.96499e-26
0.63 5.12652e-13 4.76017e-13 1.57745e-15 4.00481e-24 6.43004e-24 5.58963e-26
0.635 5.06142e-13 4.72606e-13 1.46852e-15 3.96369e-24 6.38786e-24 5.2236e-26
0.64 4.99581e-13 4.69496e-13 1.36282e-15 3.92362e-24 6.35514e-24 4.86801e-26
0.645 4.93041e-13 4.6617e-13 1.26063e-15 3.88277e-24 6.31507e-24 4.52365e-26
0.65 4.86655e-13 4.62701e-13 1.16216e-15 3.84196e-24 6.27207e-24 4.19132e-26
0.655 4.80274e-13 4.59597e-13 1.06765e-15 3.8027e-24 6.24047e-24 3.87178e-26
0.66 4.73961e-13 4.56306e-13 9.77272e-16 3.76296e-24 6.20163e-24 3.56551e-26
0.665 4.67838e-13 4.52855e-13 8.91151e-16 3.72335e-24 6.15882e-24 3.27294e-26
0.67 4.61761e-13 4.4983e-13 8.09431e-16 3.68571e-24 6.12947e-24 2.9945e-26
0.675 4.55782e-13 4.46633e-13 7.32188e-16 3.64775e-24 6.09287e-24 2.73037e-26
0.68 4.5001e-13 4.43238e-13 6.59421e-16 3.60985e-24 6.0506e-24 2.48065e-26
0.685 4.44306e-13 4.40324e-13 5.91164e-16 3.57423e-24 6.02402e-24 2.24545e-26
0.69 4.38712e-13 4.37247e-13 5.27374e-16 3.53837e-24 5.9902e-24 2.02466e-26
0.695 4.33323e-13 4.33906e-13 4.67931e-16 3.50229e-24 5.94822e-24 1.81801e-26
0.7 4.28006e-13 4.31098e-13 4.12768e-16 3.46874e-24 5.92445e-24 1.62531e-26
0.705 4.22791e-13 4.28129e-13 3.61758e-16 3.43493e-24 5.89353e-24 1.44613e-26
0.71 4.17762e-13 4.24803e-13 3.14704e-16 3.40042e-24 5.85107e-24 1.27991e-26
0.715 4.12791e-13 4.22061e-13 2.71484e-16 3.36861e-24 5.82969e-24 1.12618e-26
0.72 4.07901e-13 4.19156e-13 2.31925e-16 3.33644e-24 5.80152e-24 9.8429e-27
0.725 4.03235e-13 4.15595e-13 1.95774e-16 3.30245e-24 5.75127e-24 8.53438e-27
0.73 3.98514e-13 4.12787e-13 1.62846e-16 3.27146e-24 5.72877e-24 7.32975e-27
0.735 3.93836e-13 4.09782e-13 1.32991e-16 3.23981e-24 5.69886e-24 6.22278e-27
0.74 3.89305e-13 4.05942e-13 1.06079e-16 3.20528e-24 5.64006e-24 5.2085e-27
0.745 3.84682e-13 4.02955e-13 8.22042e-17 3.17402e-24 5.61438e-24 4.28865e-27
0.75 3.80033e-13 3.9977e-13 6.15944e-17 3.14184e-24 5.58225e-24 3.46929e-27
0.755 3.7538e-13 3.95765e-13 4.44715e-17 3.10623e-24 5.52171e-24 2.75794e-27
0.76 3.70629e-13 3.92516e-13 3.0887e-17 3.07347e-24 5.49367e-24 2.15897e-27
0.765 3.65741e-13 3.89054e-13 2.06961e-17 3.03929e-24 5.46039e-24 1.67078e-27
0.77 3.60711e-13 3.84702e-13 1.3496e-17 3.00081e-24 5.39707e-24 1.28373e-27
0.775 3.55487e-13 3.80985e-13 8.69481e-18 2.96431e-24 5.36356e-24 9.84237e-28
0.78 3.50126e-13 3.77102e-13 5.63078e-18 2.92659e-24 5.32674e-24 7.55904e-28
0.785 3.44716e-13 3.72393e-13 3.70294e-18 2.88519e-24 5.26032e-24 5.81517e-28
0.79 3.3925e-13 3.68308e-13 2.48255e-18 2.84629e-24 5.22131e-24 4.47613e-28
0.795 3.33815e-13 3.64184e-13 1.70973e-18 2.80737e-24 5.18159e-24 3.45843e-28
0.8 3.28484e-13 3.59324e-13 1.21531e-18 2.76578e-24 5.11338e-24 2.69491e-28
0.805 3.23176e-13 3.55056e-13 8.84193e-19 2.72686e-24 5.07029e-24 2.11318e-28
0.81 3.17893e-13 3.50761e-13 6.53144e-19 2.68798e-24 5.02741e-24 1.6628e-28
0.815 3.12648e-13 3.45703e-13 4.87505e-19 2.64603e-24 4.9558e-24 1.31098e-28
0.82 3.07359e-13 3.4115e-13 3.66698e-19 2.60615e-24 4.9071e-24 1.03491e-28
0.825 3.0204e-13 3.36556e-13 2.77156e-19 2.56602e-24 4.85904e-24 8.17011e-29
0.83 2.96708e-13 3.31186e-13 2.10272e-19 2.52258e-24 4.78224e-24 6.44917e-29
0.835 2.91291e-13 3.26269e-13 1.60006e-19 2.48081e-24 4.72701e-24 5.08899e-29
0.84 2.8581e-13 3.21311e-13 1.22011e-19 2.43868e-24 4.67297e-24 4.01304e-29
0.845 2.80289e-13 3.15582e-13 9.31725e-20 2.39313e-24 4.5904e-24 3.16163e-29
0.85 2.74664e-13 3.1028e-13 7.12932e-20 2.34908e-24 4.52881e-24 2.4899e-29
0.855 2.68965e-13 3.04954e-13 5.46017e-20 2.30469e-24 4.46921e-24 1.95887e-29
0.86 2.63257e-13 2.98712e-13 4.18682e-20 2.25635e-24 4.37461e-24 1.53997e-29
0.865 2.57409e-13 2.92929e-13 3.21503e-20 2.20953e-24 4.30192e-24 1.21002e-29
0.87 2.51476e-13 2.87174e-13 2.47075e-20 2.16256e-24 4.23305e-24 9.49923e-30
0.875 2.45535e-13 2.80426e-13 1.90054e-20 2.11129e-24 4.12512e-24 7.45182e-30
0.88 2.39468e-13 2.74242e-13 1.46331e-20 2.06204e-24 4.04292e-24 5.84138e-30
0.885 2.33338e-13 2.68097e-13 1.12719e-20 2.01276e-24 3.96438e-24 4.57425e-30
0.89 2.27174e-13 2.61268e-13 8.68398e-21 1.96042e-24 3.85884e-24 3.57772e-30
0.895 2.2096e-13 2.54865e-13 6.69136e-21 1.90978e-24 3.77282e-24 2.79488e-30
0.9 2.14712e-13 2.485e-13 5.15238e-21 1.8592e-24 3.68957e-24 2.17935e-30
0.905 2.08469e-13 2.41648e-13 3.96403e-21 1.80656e-24 3.58664e-24 1.6963e-30
0.91 2.02191e-13 2.35119e-13 3.04862e-21 1.75522e-24 3.49803e-24 1.31831e-30
0.915 1.95911e-13 2.28635e-13 2.33985e-21 1.70408e-24 3.41158e-24 1.02167e-30
0.92 1.89667e-13 2.21755e-13 1.79111e-21 1.65138e-24 3.30828e-24 7.89231e-31
0.925 1.83417e-13 2.15178e-13 1.36828e-21 1.60001e-24 3.21758e-24 6.0799e-31
0.93 1.77197e-13 2.08681e-13 1.04232e-21 1.5491e-24 3.12959e-24 4.668e-31
0.935 1.71041e-13 2.01855e-13 7.91355e-22 1.49701e-24 3.02662e-24 3.57068e-31
0.94 1.6491e-13 1.95334e-13 5.99161e-22 1.44637e-24 2.93544e-24 2.72229e-31
0.945 1.58837e-13 1.88923e-13 4.52039e-22 1.39643e-24 2.84745e-24 2.06735e-31
0.95 1.52855e-13 1.82243e-13 3.39658e-22 1.34567e-24 2.74615e-24 1.56322e-31
0.955 1.46927e-13 1.75871e-13 2.54299e-22 1.29645e-24 2.656e-24 1.17727e-31
0.96 1.41081e-13 1.69638e-13 1.89565e-22 1.24815e-24 2.56947e-24 8.8252e-32
0.965 1.35349e-13 1.63189e-13 1.40628e-22 1.19934e-24 2.47116e-24 6.58262e-32
0.97 1.29695e-13 1.57052e-13 1.03871e-22 1.15218e-24 2.38346e-24 4.88671e-32
0.975 1.24146e-13 1.51075e-13 7.63165e-23 1.1061e-24 2.29973e-24 3.60782e-32
0.98 1.18727e-13 1.44932e-13 5.57465e-23 1.05977e-24 2.20567e-24 2.64796e-32
0.985 1.13406e-13 1.39099e-13 4.05127e-23 1.01517e-24 2.12169e-24 1.93293e-32
0.99 1.07857e-13 1.33017e-13 2.91731e-23 9.68663e-25 2.03551e-24 1.39791e-32
0.995 1.00976e-13 1.24981e-13 2.05667e-23 9.08836e-25 1.91232e-24 9.89734e-33
1.0 9.28557e-14 1.15548e-13 1.41855e-23 8.38466e-25 1.77249e-24 6.85374e-33
1.005 8.37879e-14 1.04913e-13 9.55468e-24 7.59438e-25 1.61556e-24 4.63391e-33
1.01 7.40917e-14 9.31059e-14 6.26749e-24 6.73026e-25 1.43274e-24 3.05136e-33
1.015 6.40348e-14 8.0892e-14 4.00081e-24 5.83548e-25 1.24677e-24 1.95486e-33
1.02 5.39463e-14 6.85984e-14 2.47791e-24 4.93603e-25 1.06112e-24 1.21498e-33
1.025 4.41478e-14 5.63508e-14 1.48254e-24 4.04893e-25 8.70519e-25 7.29463e-34
1.03 3.48885e-14 4.47794e-14 8.5303e-25 3.21073e-25 6.92448e-25 4.211e-34
1.035 2.64282e-14 3.41487e-14 4.6867e-25 2.44228e-25 5.29443e-25 2.32071e-34
1.04 1.89767e-14 2.46605e-14 2.42781e-25 1.75991e-25 3.82589e-25 1.20595e-34
1.045 1.26822e-14 1.65843e-14 1.16647e-25 1.18077e-25 2.57593e-25 5.81084e-35
1.05 7.64449e-15 1.00704e-14 5.03609e-26 7.15023e-26 1.56831e-25 2.51545e-35
1.055 3.90686e-15 5.17956e-15 1.8344e-26 3.66893e-26 8.07264e-26 9.18775e-36
1.06 1.45165e-15 1.93818e-15 4.84395e-27 1.36933e-26 3.0249e-26 2.43198e-36
1.065 2.10205e-16 2.82933e-16 4.96681e-28 1.99298e-27 4.42766e-27 2.49922e-37
1.07 0 0 0 0 0 0
......@@ -17,7 +17,7 @@ FileBase ics % Base-filename of output files
OutputDir ../ % Directory for output
GlassFile dummy_glass.dat % File with unperturbed glass or
% Cartesian grid
TileFac 4 % Number of times the glass file is
TileFac 2 % Number of times the glass file is
% tiled in each dimension (must be
% an integer)
#!/bin/bash # this line only there to enable syntax highlighting in this file
# Config options for cosmo_box_star_formation_3d #
#------------------------------------------------ mesh
REGULARIZE_MESH_CM_DRIFT # Mesh regularization; Move mesh generating point towards center of mass to make cells rounder.
REGULARIZE_MESH_CM_DRIFT_USE_SOUNDSPEED # Limit mesh regularization speed by local sound speed
REGULARIZE_MESH_FACE_ANGLE # Use maximum face angle as roundness criterion in mesh regularization
ENFORCE_JEANS_STABILITY_OF_CELLS # this imposes an adaptive floor for the temperature
#------------------------------------------------ Time integration options
TREE_BASED_TIMESTEPS # non-local timestep criterion (take 'signal speed' into account)
#------------------------------------------------ Gravity treatment
SELFGRAVITY # gravitational intraction between simulation particles/cells
HIERARCHICAL_GRAVITY # use hierarchical splitting of the time integration of the gravity
CELL_CENTER_GRAVITY # uses geometric centers to calculate gravity of cells, only possible with HIERARCHICAL_GRAVITY
ALLOW_DIRECT_SUMMATION # Performed direct summation instead of tree-based gravity if number of active particles < DIRECT_SUMMATION_THRESHOLD (= 3000 unless specified differently here)
DIRECT_SUMMATION_THRESHOLD=500 # Overrides maximum number of active particles for which direct summation is performed instead of tree based calculation
#------------------------------------------------ Gravity softening
NSOFTTYPES=2 # Number of different softening values to which particle types can be mapped.
MULTIPLE_NODE_SOFTENING # If a tree node is to be used which is softened, this is done with the softenings of its different mass components
INDIVIDUAL_GRAVITY_SOFTENING=32 # bitmask with particle types where the softenig type should be chosen with that of parttype 1 as a reference type
ADAPTIVE_HYDRO_SOFTENING # Adaptive softening of gas cells depending on their size
#------------------------------------------------ TreePM Options
PMGRID=256 # Enables particle mesh; number of cells used for grid in each dimension
RCUT=5.0 # This can be used to override the maximum radius in which the short-range tree-force is evaluated (in case the TreePM algorithm is used). The default value is 4.5, given in mesh-cells.
#------------------------------------------------ Single/Double Precision
DOUBLEPRECISION=1 # Mode of double precision: not defined: single; 1: full double precision 2: mixed, 3: mixed, fewer single precisions; unless short of memory, use 1.
DOUBLEPRECISION_FFTW # FFTW calculation in double precision
OUTPUT_COORDINATES_IN_DOUBLEPRECISION # will always output coordinates in double precision
NGB_TREE_DOUBLEPRECISION # if this is enabled, double precision is used for the neighbor node extension
#------------------------------------------------ On the fly FOF groupfinder
FOF # enable FoF output
FOF_PRIMARY_LINK_TYPES=2 # 2^type for the primary dark matter type
FOF_SECONDARY_LINK_TYPES=1+16+32 # 2^type for the types linked to nearest primaries
#------------------------------------------------ Subfind
SUBFIND # enables substructure finder
SAVE_HSML_IN_SNAPSHOT # stores hsml, density, and velocity dispersion values in the snapshot files
SUBFIND_CALC_MORE # calculates also the velocity dispersion in the local density estimate (this is automatically enabled by several other options, e.g. SAVE_HSML_IN_SNAPSHOT)
SUBFIND_EXTENDED_PROPERTIES # adds calculation of further quantities related to angular momentum in different components
#------------------------------------------------ Things for special behaviour
PROCESS_TIMES_OF_OUTPUTLIST # goes through times of output list prior to starting the simulaiton to ensure that outputs are written as close to the desired time as possible (as opposed to at next possible time if this flag is not active)
#------------------------------------------------ Output/Input options
REDUCE_FLUSH # only flush output to log-files in predefined intervals
OUTPUT_CPU_CSV # output of a cpu.csv file on top of cpu.txt
HAVE_HDF5 # needed when HDF5 I/O support is desired (recommended)
#------------------------------------------------ Testing and Debugging options
DEBUG # enables core-dumps
GENERATE_GAS_IN_ICS # Generates gas from dark matter only ICs.
#------------------------------------------------ cooling and star formation
COOLING # Simple primordial cooling
USE_SFR # Star formation model, turning dense gas into collisionless partices
This diff is collapsed.
""" @package ./examples/cosmo_box_star_formation_3d/
Code that checks results of structure formation simulation
created by Rainer Weinberger, last modified 28.02.2019
""" load libraries """
import sys ## system calls
import numpy as np ## load numpy
import h5py ## load h5py; needed to read snapshots
import matplotlib.pyplot as plt ## plot stuff
simulation_directory = str(sys.argv[1])
print("examples/cosmo_box_star_formation_3d/ checking simulation output in directory " + simulation_directory)
""" script settings """
createReferenceSolution = False
createFigures = False
CompareToReferenceRun = True
FloatType = np.float64 # double precision: np.float64, for single use np.float32
Boxsize = 7.5 ##Mpc/h
CellsPerDimension = 32 ## 3d sim
NumberOfCells = CellsPerDimension * CellsPerDimension * CellsPerDimension
HubbleParam = 0.6774
UnitMass = 1.0e10/HubbleParam
Volume = (Boxsize/HubbleParam) * (Boxsize/HubbleParam) * (Boxsize/HubbleParam)
""" read in star formation rate output snapshot """
directory = simulation_directory+"/output/"
filename = "sfr.txt"
data = np.loadtxt(directory+filename)
if createReferenceSolution:
np.savetxt("./examples/cosmo_box_star_formation_3d/sfr_ref_L8n32_reduced.txt", data[::50,:])
scalefactor = data[:,0]
redshift_plus_one = 1.0 / data[:,0]
sfrDensity = data[:,3] / Volume
cum_mass_stars = data[:,5] * UnitMass / Volume
""" compare to reference output """
data_ref = np.loadtxt("./examples/cosmo_box_star_formation_3d/sfr_ref_L8n32_reduced.txt")
scalefactor_ref = data_ref[:,0]
redshift_plus_one_ref = 1.0 / data_ref[:,0]
cum_mass_stars_ref = data_ref[:,5]* UnitMass / Volume
scalefactor_to_probe = np.array([0.25,0.3333,0.5,0.66667,1])
delta_mass_stars = np.zeros(scalefactor_to_probe.shape)
avg_mass_stars = np.zeros(scalefactor_to_probe.shape)
if createFigures:
## figure -- dark matter and stellar positions
fig, ax = plt.subplots( 2, 2, figsize=np.array([6.9,6.9]), sharex=True, sharey=True )
fig.subplots_adjust(left = 0.09, bottom = 0.09,right = 0.98, top = 0.98, hspace=0.0, wspace=0.0)
for i, a in enumerate(scalefactor_to_probe):
delta_mass_stars[i] = np.interp(a, scalefactor, cum_mass_stars)
delta_mass_stars[i] -= np.interp(a, scalefactor_ref, cum_mass_stars_ref)
avg_mass_stars[i] = 0.5 * np.interp(a, scalefactor, cum_mass_stars)
avg_mass_stars[i] += 0.5 * np.interp(a, scalefactor_ref, cum_mass_stars_ref)
filename = "snap_%03d.hdf5" % (i)
data = h5py.File(directory+filename, "r")
print("could not open "+directory+filename)
pos = np.array(data["PartType1"]["Coordinates"], dtype = FloatType) / HubbleParam / 1000.
posstars = np.array(data["PartType4"]["Coordinates"], dtype = FloatType) / HubbleParam / 1000.
if(pos.shape[0] > 32**3):
i_select = np.random.uniform(low=0.0, high=pos.shape[0], size=32**3).astype(
i_select = np.arange(pos.shape[0])
## Try to show the different gas phases and their spatial distribution
z = 1./ a - 1
if i == 0:
ax_col =
ax_row =
print('col %d, row %d'%(ax_col, ax_row))
ax[ax_row][ax_col].text(5.5,10.0,'redshift %.1f'%z)
ax[ax_row][ax_col].scatter(pos[i_select, 0], pos[i_select, 1], marker='.', s=0.01, alpha=0.5, rasterized=True)
ax[ax_row][ax_col].scatter(posstars[:, 0], posstars[:, 1], marker='*',c='r', s=0.01, alpha=1.0, rasterized=True)
fig.savefig(simulation_directory+'/dark_matter_and_stars_evolution.pdf', dpi=300)
""" figure of stellar mass density """
ax = plt.axes([0.15,0.15,0.8,0.8])
ax.plot(redshift_plus_one, cum_mass_stars, label="new")
ax.plot(redshift_plus_one_ref, cum_mass_stars_ref, label="ref")
ax.set_xticklabels(["0", "1", "2", "4", "8"])
ax.set_ylabel(r"Stellar mass denstiy [M$_\odot$ Mpc$^{-3}$]")
if CompareToReferenceRun:
""" check if deviations within tolerance """
abs_tolerance = 5e7 ## depends on box size...
if np.max( np.abs(delta_mass_stars ) ) > abs_tolerance:
print("Error: stellar mass deviaties from reference more than %g (Msun/Mpc^3): "% abs_tolerance )
print(1./scalefactor_to_probe - 1)
print("relative deviation:")
print(delta_mass_stars / avg_mass_stars)
print("absolute deviation (1e6 Msun/Mpc^3):")
""" if everything is ok """
""" @package examples/cosmo_box_star_formation_3d/
Code creates the output list; ICs need to be present already
created by Rainer Weinberger, last modified 28.02.2019
""" load libraries """
import sys # system calls
import numpy as np # scientific computing package
import h5py # hdf5 format
import os # operating system interface
from subprocess import call # execute bash commands
## create new ics with 'music' or 'ngenic' or just 'copy' existing
## ones to the run directory
""" input """
simulation_directory = str(sys.argv[1])
print("examples/cosmo_box_star_formation_3d/ " + simulation_directory)
""" initial conditions: either copy or create with code """
if ic_creation == 'copy':
call(['cp', './examples/cosmo_box_star_formation_3d/L8n32/ics', simulation_directory+'/ics'])
elif ic_creation == 'music':
status = call(['hg', 'clone', '', simulation_directory+'/music'])
if status != 0:
print('CREATE: ERROR: hg clone failed!')
cwd = os.getcwd()
status = call(['make'])
if status != 0:
print('CREATE: ERROR: make failed!')
status = call(['./MUSIC',cwd+'/examples/cosmo_box_star_formation_3d/param_music.txt'])
if status != 0:
print('CREATE: ERROR: execution failed!')
elif ic_creation == 'ngenic':
status = call(['git', 'clone', '', simulation_directory+'/ngenic'])
if status != 0:
print('CREATE: ERROR: git clone failed!')
cwd = os.getcwd()
status = call(['make'])
if status != 0:
print('CREATE: ERROR: make failed!')
status = call(['mpiexec','-np','1','./N-GenIC',cwd+'/examples/cosmo_box_star_formation_3d/param_ngenic.txt'])
if status != 0:
print('CREATE: ERROR: execution failed!')
print("CREATE: ERROR: no valid option for ic creation! choose 'copy', 'music' or 'ngenic'")
""" set output times """
outputTimes = np.array([0.2,0.25,0.33,0.5,0.66,1], dtype=np.float64)
ones = np.ones(outputTimes.shape,
""" write output list file """
data = np.array([outputTimes, ones]).T
np.savetxt(simulation_directory+"/output_list.txt",data, fmt="%g %1.f" )
""" normal exit """
%% parameters for cosmo_box_star_formation_3d
%---- Relevant files
InitCondFile run/examples/cosmo_box_star_formation_3d/ics
OutputDir run/examples/cosmo_box_star_formation_3d/output
SnapshotFileBase snap
OutputListFilename run/examples/cosmo_box_star_formation_3d/output_list.txt
%---- File formats
ICFormat 1
SnapFormat 3
%---- CPU-time limits
TimeLimitCPU 90000
CpuTimeBetRestartFile 12000
FlushCpuTimeDiff 120
ResubmitOn 0
ResubmitCommand my-scriptfile
%----- Memory alloction
MaxMemSize 2500
%---- Caracteristics of run
TimeBegin 0.0078125 % Begin of the simulation; z=127
TimeMax 1.0 % End of the simulation z=0
%---- Basic code options that set the type of simulation
ComovingIntegrationOn 1
PeriodicBoundariesOn 1
CoolingOn 1
StarformationOn 1
%---- Cosmological parameters (Planck cosmology)
Omega0 0.3089
OmegaLambda 0.6911
OmegaBaryon 0.0486
HubbleParam 0.6774
BoxSize 7500.0
%---- Output frequency and output parameters
OutputListOn 1
TimeBetSnapshot 0.0
TimeOfFirstSnapshot 0.0
TimeBetStatistics 0.01
NumFilesPerSnapshot 1
NumFilesWrittenInParallel 1
%---- Accuracy of time integration
TypeOfTimestepCriterion 0
ErrTolIntAccuracy 0.012
CourantFac 0.3
MaxSizeTimestep 0.005
MinSizeTimestep 1.0e-9
%---- Treatment of empty space and temperature limits
InitGasTemp 244.8095
MinGasTemp 5.0
MinimumDensityOnStartUp 1.0e-20
LimitUBelowThisDensity 0.0
LimitUBelowCertainDensityToThisValue 0.0
MinEgySpec 0.0
%---- Tree algorithm, force accuracy, domain update frequency
TypeOfOpeningCriterion 1
ErrTolTheta 0.7
ErrTolForceAcc 0.0025
MultipleDomains 8
TopNodeFactor 2.5
ActivePartFracForNewDomainDecomp 0.01
%---- Initial density estimate
DesNumNgb 64
MaxNumNgbDeviation 4
%---- System of units
UnitLength_in_cm 3.085678e21 % 1.0 kpc
UnitMass_in_g 1.989e43 % 1.0e10 solar masses
UnitVelocity_in_cm_per_s 1e5 % 1 km/sec
GravityConstantInternal 0
%---- Gravitational softening lengths
SofteningComovingType0 6.0
SofteningComovingType1 6.0
SofteningMaxPhysType0 3.0
SofteningMaxPhysType1 3.0
GasSoftFactor 2.5
SofteningTypeOfPartType0 0
SofteningTypeOfPartType1 1
SofteningTypeOfPartType2 1
SofteningTypeOfPartType3 1
SofteningTypeOfPartType4 1
SofteningTypeOfPartType5 1
MinimumComovingHydroSoftening 0.25
AdaptiveHydroSofteningSpacing 1.2
%----- Mesh regularization options
CellShapingSpeed 0.5
CellMaxAngleFactor 2.25
ReferenceGasPartMass 0
TargetGasMassFactor 1