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
general_aswns [2022/08/16 13:53]
theoastro 		general_aswns [2022/08/16 21:17] (current)
theoastro
Line 4: Line 4:
  
 If you open ASWNS.f90, you will find the following parameters  If you open ASWNS.f90, you will find the following parameters 
 +  * **nr** (integer) --> total number of radial grid points 
 +  * **nth** (integer) --> total number of angular grid points 
 +  * **mls** (integer) -->  maximum number of employed Legendre polynomials for the expansion in theta 
 +  * **rhalf** (double) --> end of the inner radial grid/beginning of the outer grid  
 +  * **rmax** (double) --> end of outer (radial) grid 
 +  * **planar_symmetry** (logical) --> employing symmetry through neglecting odd harmonics 
 +  * **qfactor** (double) -->  damping factor for the metric quantities during iteration
 +  * **abs_tol** (double) --> tolerance when solution is considered to be converged
 +  * **verbose** (logical) –> determines if the output is put on the screen
 +  * **wait_iters** (integer) --> number of iterations before checking for convergence after hydro relaxation 
 +  
 +And the following constants are used: 
 +  * nm (double) = 8.42335d-58 --> neutron mass given in Msun
 +  * rhon (double) =  4.339d-4 --> nuclear saturation density given in Msun^-2
 +  * km_to_msol (double) = 0.677218135d0 --> conversion from km to solar masses 
 +  
 +While it isn't to be expected that you need to change the overall constants, it might be worth changing some of the iteration parameters. For this purpose, let us consider running the example for the computation of the non-barotropic configuration and let's vary the following parameters
  
-integer, parameter system::nr = 6000 +  nr = 4000 
-  total number of radial grid points More... +  nth = 15 
-  +  mls = 10 
-integer, parameter system::nth = 50 +which leads to a setup with mass of 2.151388177700834 Msun, but if we increase the resolution  
-  number of angular mesh points (1 = 1D TOV solution) More... + 
-  +  nr = 15000 
-integer, parameter system::mls = 30 +  nth = 100 
-  number of Legendre polinomials for the expansion in theta (0 = 1D) More... +  mls = 50 
-  + 
-logical, parameter system::verbose = .true. +we get 2.1513872961030418Msun, but have to wait much longer. It might also be of interest for you to change the outer boundary location and see which effect this might have on the overall accuracy of the solution.  
-  determine the amount of logging More... + 
-  +===== Creating the Doxygen Documentation ===== 
-double precision, parameter system::rhalf = 30d0 + 
-  end of the inner radial grid / begin of the outer radial grid More... +Another useful feature of ASWNS is the existence of a doxygen documentation. To create these, please first install doxygen, and then run it via  
-  +   sudo apt-get install doxygen 
-double precision, parameter system::rmax = 1000d0 +    doxygen 
-  outer radius (end of the radial grid) More... +This creates HTML pages on which the functions are defined and also additional definitions of variables are given. 
-  +
-logical, parameter system::planar_symmetry = .true. +
-  if planar_symmetry == .true., only even harmonics are used (for the metric) More... +
-  +
-double precision, parameter system::qfactor = .1d0 +
-  damping factor for the metric quantities in the convergence loop More... +
-  +
-double precision, parameter system::abs_tol = 1d-11 +
-  criterion for the absolute convergence of the configuration More... +
-  +
-integer system::wait_iters = 50 +
-  number of iterations before start checking for convergence after hydro relaxation More... +
-  +
-double precision, parameter system::mn = 8.42335d-58 +
-  neutron mass [Msun] More... +
-  +
-double precision, parameter system::rhon = 4.339d-4 +
-  nuclear saturation density [Msun^-2] More... +
-  +
-double precision, parameter system::km_to_msol = 0.677218135d0 +
-  convertion from km to solar masses More... +
-  +
-double precision, parameter system::pi =DACOS(-1d0) +
-  greek pi More... +
-  +
-double precision, dimension(0:nr+1) system::r +
-  radial grid (1 ghost zone per side) More... +
-  +
-double precision, dimension(0:nth+1) system::th +
-  angular grid (1 ghost zone per side) More...+
Last modified: le 2022/08/16 13:53

Page Tools