<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<meta http-equiv="content-type" content="text/html; charset=ISO-8859-1">
</head>
<body bgcolor="#ffffff" text="#000000">
Hi all,<br>
<br>
I need to get radial profiles of my sims - like density as a
function or radius. I know how to plot them using the method
pc.add_profile_sphere but I now need to print out the values in a
txt file. For this reason, I tried to use the class BinnedProfile1D.
Below is my script:<br>
<br>
<address>##################################</address>
<address>file = 'DD0000/CommonEnvelope0021'</address>
<address>pf = load(file)</address>
<address>r_min = pf.h.get_smallest_dx()</address>
<address>r_max = 0.4</address>
<address>sphere = pf.h.sphere([0.5, 0.5, 0.5], r_max)</address>
<address><br>
prof1d = L.BinnedProfile1D(sphere, 32, "Radius", r_min,
r_max,log_space=False, lazy_reader=True)</address>
<address>prof1d.add_fields("Density")</address>
<address><br>
prof1d.write_out('test.txt','%0.6e','left')</address>
<address>
##################################</address>
<br>
It run well but returned a density profile equal to 0:<br>
<br>
<address>##############################</address>
<address>[17:25:54] litai05:$ more test.txt </address>
<address># Radius Density</address>
<address>7.812500e-03 0.000000e+00</address>
<address>2.006836e-02 0.000000e+00</address>
<address>3.232422e-02 0.000000e+00</address>
<address>4.458008e-02 0.000000e+00</address>
<address>5.683594e-02 0.000000e+00</address>
<address>6.909180e-02 0.000000e+00</address>
<address>8.134766e-02 0.000000e+00</address>
<address>9.360352e-02 0.000000e+00</address>
<address>1.058594e-01 0.000000e+00</address>
<address>1.181152e-01 0.000000e+00</address>
<address>1.303711e-01 0.000000e+00</address>
<address>1.426270e-01 0.000000e+00</address>
<address>1.548828e-01 0.000000e+00</address>
<address>1.671387e-01 0.000000e+00</address>
<address>1.793945e-01 0.000000e+00</address>
<address>1.916504e-01 0.000000e+00</address>
<address>2.039063e-01 0.000000e+00</address>
<address>2.161621e-01 0.000000e+00</address>
<address>2.284180e-01 0.000000e+00</address>
<address>2.406738e-01 0.000000e+00</address>
<address>2.529297e-01 0.000000e+00</address>
<address>2.651855e-01 0.000000e+00</address>
<address>2.774414e-01 0.000000e+00</address>
<address>2.896973e-01 0.000000e+00</address>
<address>3.019531e-01 0.000000e+00</address>
<address>3.142090e-01 0.000000e+00</address>
<address>3.264648e-01 0.000000e+00</address>
<address>3.387207e-01 0.000000e+00</address>
<address>3.509766e-01 0.000000e+00</address>
<address>3.632324e-01 0.000000e+00</address>
<address>3.754883e-01 0.000000e+00</address>
<address>3.877441e-01 0.000000e+00</address>
<address>
##############################</address>
<br>
Could someone tell me what I am doing wrong ?<br>
<br>
Thanks,<br>
<br>
JC<br>
</body>
</html>