Commit 95af9357 authored by m.tavelli's avatar m.tavelli

TOV star implemented in 3D

parent a7851d6f
This diff is collapsed.
......@@ -40,6 +40,40 @@
TYPE(tEquations) :: EQN
!Variables for NSTOV module
#ifdef SPHERICAL
INTEGER, PARAMETER :: NSTOV_nODE = 3
#else
INTEGER, PARAMETER :: NSTOV_nODE = 4
#endif
INTEGER, PARAMETER :: NSTOV_nODE_p = 3, NSTOV_ATMO = 0 ! 1 => Landau atmo
REAL(8), PARAMETER :: NSTOV_rho_c = 1.28e-3
REAL(8), PARAMETER :: NSTOV_kappa = 100
REAL(8), PARAMETER :: p_floor = 1.0e-16, rho_floor = 1.0e-10
REAL(8) :: NSTOV_rho_atmo != 1e-10
REAL(8), PARAMETER :: NSTOV_p_atmo = 1e-15 , NSTOV_t_atm=1.0
!
REAL(8) :: Mbh = 1.0, aom = 0.0
REAL(8), PARAMETER :: P_eps = 1e-4
integer, parameter :: MYRANK=0
!
TYPE tNSTOVVar
INTEGER :: Computed
REAL(8) :: Int
!INTEGER, PARAMETER :: nODE = 3
REAL(8) :: Mass, radius,r1,dr1,r2,dr2,rW,drW,rUni,C,p_R ,rho_R,lapse_C
INTEGER :: n_B,iradius,ir1,ir2,irW
REAL(8), DIMENSION (:,:), ALLOCATABLE :: q,dq
REAL(8), DIMENSION (:), ALLOCATABLE :: r,dr
REAL(8), DIMENSION (:), ALLOCATABLE :: qloc
END TYPE tNSTOVVar
TYPE(tNSTOVVar) :: NSTOVVar,NSTOVVar_bar,NSTOVVar_barNew
!
REAL(8), PARAMETER :: CoordTol =1e-11
! 3-point Gaussian quadrature
REAL, PARAMETER :: sGP3(3) = (/ 0.5-sqrt(15.)/10.,0.5,0.5+sqrt(15.)/10. /)
REAL, PARAMETER :: wGP3(3) = (/ 5./18., 8./18., 5./18. /)
......
This diff is collapsed.
......@@ -4036,11 +4036,11 @@ RECURSIVE SUBROUTINE pderefinecriteria(refine_flag, max_luh,min_luh,x)
refine_flag=0
end if
if(nDim .eq. 3) then
if(sqrt(x(1)**2+x(2)**2+x(3)**2)<15.0) then
refine_flag=2
end if
end if
!if(nDim .eq. 3) then
!if(sqrt(x(1)**2+x(2)**2+x(3)**2)<15.0) then
! refine_flag=2
!end if
!end if
#endif
END SUBROUTINE pderefinecriteria
......
{
"project_name": "FOCCZ4",
"compiler_flags": "-DCCZ4EINSTEIN -DGLMROT -DRNSTOV",
"paths": {
"peano_kernel_path": "./Peano",
"exahype_path": "./ExaHyPE",
"output_directory": "./ApplicationExamples/FOCCZ4/FOCCZ4",
"log_file": "whatever.log"
},
"architecture": "noarch",
"computational_domain": {
"dimension": 3,
"end_time": 101.0,
"offset": [
-5.0,
-5.0,
-5.0
],
"width": [
10.0,
10.0,
10.0
]
},
"shared_memory": {
"cores": 10,
"properties_file": "sharedmemory.properties",
"autotuning_strategy": "dummy",
"background_job_consumers": 9
},
"distributed_memory": {
"timeout": 6000,
"load_balancing_type": "static",
"buffer_size": 6400,
"load_balancing_strategy": "hotspot",
"node_pool_strategy": "fair",
"ranks_per_node": 10
},
"optimisation": {
"fuse_algorithmic_steps": "all",
"fuse_algorithmic_steps_rerun_factor": 0.99,
"fuse_algorithmic_steps_diffusion_factor": 0.99,
"spawn_predictor_as_background_thread": true,
"spawn_update_as_background_thread": true,
"spawn_amr_background_threads": true,
"disable_vertex_exchange_in_time_steps": true,
"time_step_batch_factor": 0.0,
"disable_metadata_exchange_in_batched_time_steps": false,
"double_compression": 0.0,
"spawn_double_compression_as_background_thread": true
},
"solvers": [
{
"type": "Limiting-ADER-DG",
"name": "FOCCZ4Solver",
"order": 3,
"maximum_mesh_size": 2.0,
"maximum_mesh_depth": 0,
"time_stepping": "global",
"aderdg_kernel": {
"language": "C",
"nonlinear": true,
"terms": [
"flux",
"ncp",
"source"
],
"space_time_predictor": {},
"optimised_terms": [
"fusedsource"
],
"optimised_kernel_debugging": [],
"implementation": "generic"
},
"point_sources": 0,
"limiter": {
"dmp_observables": 3,
"dmp_relaxation_parameter": 1e+3,
"dmp_difference_scaling": 1e+4,
"patch_size": "max",
"implementation": "generic"
},
"fv_kernel": {
"language": "C",
"terms": [
"flux",
"ncp",
"source"
],
"scheme": "musclhancock",
"slope_limiter" : "minmod",
"implementation": "generic"
},
"variables": [
{
"name": "MYQ",
"multiplicity": 96
}
],
"parameters": {
"reference": "CCZ4TOV"
},
"plotters": [
{
"type": "user::defined",
"name": "TecplotWriter",
"time": 0.0,
"repeat": 1.0,
"output": "./output/tecplot",
"variables": 96
}
]
}
]
}
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment