Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision 		Previous revision
examples_aswns_nonbaro [2022/08/16 14:30]
theoastro [A first Test] 		examples_aswns_nonbaro [2022/08/16 20:52] (current)
theoastro [Computing non-barotropic stars]
Line 1: Line 1:
 ===== Computing non-barotropic stars ===== ===== Computing non-barotropic stars =====
- +==== Theoretical Background and Model Parameters ====
-This test follows the computation focusses on the computation of a non-barotropic configuration +
- +
-==== Model Parameters ==== +
- +
-==== Theoretical Background ====+
  
 For a careful theoretical understanding, one has to investigate the exact model parameters that are employed. The equation of state follows the form  (Eq. 45 of [[https://arxiv.org/pdf/1908.11258.pdf|Camelio et al., 2019]]):  For a careful theoretical understanding, one has to investigate the exact model parameters that are employed. The equation of state follows the form  (Eq. 45 of [[https://arxiv.org/pdf/1908.11258.pdf|Camelio et al., 2019]]): 
Line 23: Line 18:
 {{ASWNS_EOS_4.png|EOS4}} {{ASWNS_EOS_4.png|EOS4}}
  
-==== Changing the code input  ====+==== Code input  ====
  
 The simple example that is given in the git repository looks at a cold, rigidly rotating neutron star.  The simple example that is given in the git repository looks at a cold, rigidly rotating neutron star. 
Line 80: Line 75:
 ==== Plotting your results ==== ==== Plotting your results ====
  
 +ASWNS also provides a simple script to transform the binary output (that is by default stored in `star.out') to ascii file format. If you want to do this, please start a python environment and run
 +  
 +  from binfile_reader import * 
 +  bin2dat(binstar('star.out'),'star.dat')  
 + 
 +This would create an output file with the name star.dat in which you can see the columns theta, r, rest mass density, pressure, internal energy density, fluid velocity, conformal factor, lapse, ZAMO energy, entropy per baryon, and the temperature. 
 +
 +The most straight forward way to make a 2d plot from the data would be to use gnuplot or any other low-level routine and type: 
 +
 +  gnuplot 
 +  gnuplot> set xrange [0:10]
 +  gnuplot> set yrange [0:10]
 +  gnuplot> splot 'star.dat' u ($2)*cos($1):($2)*sin($1):($3) w l
 +
 +{{Simple2d_plot.png}}
 +
 +
 +Of course, it would also be possible to produce simple 1d plots and/or use python for plotting, where an example is given below, where we show the density along one radial direction (by using just the first 1000 entries) for approximately theta=0: 
 +
 +  import numpy as np
 +  import matplotlib.pyplot as plt
 +  
 +  th, r, rho = np.loadtxt('star.dat', usecols=[0,1,2], unpack=True, comments='#')
 +  fig = plt.figure()
 +  plt.plot(r[0:1000], rho[0:1000])
 +  plt.xlabel('radius')
 +  plt.ylabel('density')
 +  plt.show();
 +
 +
 +{{ASWNS_nonbaro_1d.png}}
 + 
Last modified: le 2022/08/16 14:30

Page Tools