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 16:15]
theoastro [Plotting your results] 		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 94: Line 89:
   gnuplot> splot 'star.dat' u ($2)*cos($1):($2)*sin($1):($3) w l   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 16:15

Page Tools