Permalink
Browse files

add in sample parameter files

  • Loading branch information...
Evan O'Connor
Evan O'Connor committed Nov 19, 2014
1 parent 507a92f commit 0fd503e56f903e56a285b07dcd5c69995424fd1b
View
@@ -0,0 +1,6 @@
all:
$(MAKE) -C src
clean:
$(MAKE) -C src clean
View
@@ -0,0 +1,25 @@
While this verison of GR1D is open source, its copyright is held by
Evan O’Connor and Christian Ott. In the absence of suitable open
scientific software licenses, we release this version of GR1D to the
community under the Creative Commons attribution-noncommercial-share
alike license:
http://creativecommons.org/licenses/by-nc-sa/3.0/us
Essentially, you may use GR1D, but must make reference to our work,
must not use GR1D for commercial purposes, and any code including or
using our routines or part of them may be made publically available,
and if so, only under the same license.
When using this version of GR1D in published work please make
reference to the papers describing it
Evan O’Connor, title..., arXiv:1411.????, 2014
Evan O’Connor and Christian D. Ott, A New Spherically-Symmetric
General Relativistic Hydrodynamics Code for Stellar Collapse to
Neutron Stars and Black Holes, Class. Quant. Grav., 27 114103, 2010.
View
@@ -0,0 +1,13 @@
F90=gfortran
F90FLAGS=-O3 -g -fopenmp
MODINC="-I ./"
HDF5DIR=
HDF5INCS=-I${HDF5DIR}/include
HDF5LIBS=-L${HDF5DIR}/lib -lhdf5 -lhdf5_fortran -lhdf5 -lz
HAVE_NUC_EOS=1
HAVE_RESTART=1
HAVE_LEAK_ROS=0
HAVE_LAPACK=1
@@ -0,0 +1,74 @@
##### Job parameters #######
jobname = "Abdikamalov Test"
GR = 0 # 1 for GR, 0 for Newtonian, enforced by M1_testcase.F90
outdir = "Data_TC2" # output directory
initial_data = "M1test" # "Sedov", "Shocktube", "OSC", "Collapse", "M1test"
gravity_active = 0 # do we want gravity?
ntmax = 1000000 # maximum timestep
tend = 0.001d0 # maximum time
####### Grid parameters #######
geometry = 2 # 1: planar, 2: spherical
gridtype = "blah" # set by M1_testcase.F90
radial_zones = 100 # number of radial zones, enforced by M1_testcase.F90
ghosts1 = 4 # number of ghost cells
####### Hydro Parameters #######
do_hydro = 0 # 1 for hydro on
iorder_hydro = 2 # RK order
cffac = 0.5d0 # CFL factor
reconstruction_method = "tvd" # "tvd", "pc","ppm"
ppm_origin_TVD = 5
ppm_reconstruct_temp = 0 ###### added by csul
tvd_limiter = "MC" # MC, minmod,
flux_type = "HLLE" # HLLE, HLLC
####### EOS parameters #######
eoskey = 4 # hybrid: 1
# poly: 2
# hot nuclear: 3
# ideal: 4
####### Output parameters #######
ntinfo = 100 # stdout
dtout = 1.0d-4 # time between outputs
dtout_scalar = 1.0d-5 # time between scalar out
ntout = -1 # output every timesteps
ntout_scalar = -1 # scalar output
####### Restart parameters #######
ntout_restart = -1 # restart output frequency in timesteps
dtout_restart = 0.005d0 # restart output frequency in time
do_restart = 0
restart_file_name = "blah"
########## M1 settings ##########
do_M1 = 1 #1 for M1 transport scheme
v_order = 0 #enforced by M1_testcases.F90
extraction_radii = -1.0d0 #set by M1_testcases.F90
number_species = 3 #number of neutrino species, limited
#options will work, must be comensurate with table
number_groups = 18 #number of energy groups, limited
#options will work, must be comensurate with table
opacity_table = "/Volumes/vault/EOS/NuLib_LS220_18.h5"
number_eas = 3 #number of opacity variables to read in
M1closure = 'ME' # Closure choice
testcase = 2 #if you want a test case set this to the
#test case number, along with M1test as "initial_data"
include_epannihil_kernels = 0 #for full thermal treatment - forced by M1_testcase.F90
include_Ielectron_exp = 0 #Inelastic scattering on electrons, explicit? - forced by M1_testcase.F90
include_Ielectron_imp = 0 #Inelastic scattering on electrons, implicit? - forced by M1_testcase.F90
include_energycoupling_exp = 0 #energy coupling, explicit? - forced by M1_testcase.F90
include_energycoupling_imp = 0 #energy coupling, implicit? - forced by M1_testcase.F90
M1_control = 0 #1 means you must specify the settings for each phase
# in the parameter file, collapse, bounce and postbounce.
####### Neutrino parameters #######
fake_neutrinos = 0 # 1 for ANY neutrino scheme
####### Atmosphere parameters #######
atmo_rho_rel_min = 0.0d0
atmo_rho_abs_min = 1.5d3
atmo_fac = 0.1d0
@@ -0,0 +1,74 @@
##### Job parameters #######
jobname = "Smit Test"
GR = 0 # 1 for GR, 0 for Newtonian, enforced by M1_testcase.F90
outdir = "Data_TC4" # output directory
initial_data = "M1test" # "Sedov", "Shocktube", "OSC", "Collapse", "M1test"
gravity_active = 0 # do we want gravity?
ntmax = 1000000 # maximum timestep
tend = 0.0005d0 # maximum time
####### Grid parameters #######
geometry = 2 # 1: planar, 2: spherical
gridtype = "blah" # set by M1_testcase.F90
radial_zones = 800 # number of radial zones, enforced by M1_testcase.F90
ghosts1 = 4 # number of ghost cells
####### Hydro Parameters #######
do_hydro = 0 # 1 for hydro on
iorder_hydro = 2 # RK order
cffac = 0.5d0 # CFL factor
reconstruction_method = "tvd" # "tvd", "pc","ppm"
ppm_origin_TVD = 5
ppm_reconstruct_temp = 0 ###### added by csul
tvd_limiter = "MC" # MC, minmod,
flux_type = "HLLE" # HLLE, HLLC
####### EOS parameters #######
eoskey = 4 # hybrid: 1
# poly: 2
# hot nuclear: 3
# ideal: 4
####### Output parameters #######
ntinfo = 100 # stdout
dtout = 0.5d-4 # time between outputs
dtout_scalar = 0.5d-5 # time between scalar out
ntout = -1 # output every timesteps
ntout_scalar = -1 # scalar output
####### Restart parameters #######
ntout_restart = -1 # restart output frequency in timesteps
dtout_restart = 0.005d0 # restart output frequency in time
do_restart = 0
restart_file_name = "blah"
########## M1 settings ##########
do_M1 = 1 #1 for M1 transport scheme
v_order = 0 #enforced by M1_testcases.F90
extraction_radii = -1.0d0 #set by M1_testcases.F90
number_species = 3 #number of neutrino species, limited
#options will work, must be comensurate with table
number_groups = 18 #number of energy groups, limited
#options will work, must be comensurate with table
opacity_table = "/Volumes/vault/EOS/NuLib_LS220_18.h5"
number_eas = 3 #number of opacity variables to read in
M1closure = 'ME' # Closure choice
testcase = 4 #if you want a test case set this to the
#test case number, along with M1test as "initial_data"
include_epannihil_kernels = 0 #for full thermal treatment - forced by M1_testcase.F90
include_Ielectron_exp = 0 #Inelastic scattering on electrons, explicit? - forced by M1_testcase.F90
include_Ielectron_imp = 0 #Inelastic scattering on electrons, implicit? - forced by M1_testcase.F90
include_energycoupling_exp = 0 #energy coupling, explicit? - forced by M1_testcase.F90
include_energycoupling_imp = 0 #energy coupling, implicit? - forced by M1_testcase.F90
M1_control = 0 #1 means you must specify the settings for each phase
# in the parameter file, collapse, bounce and postbounce.
####### Neutrino parameters #######
fake_neutrinos = 0 # 1 for ANY neutrino scheme
####### Atmosphere parameters #######
atmo_rho_rel_min = 0.0d0
atmo_rho_abs_min = 1.5d3
atmo_fac = 0.1d0
@@ -0,0 +1,74 @@
##### Job parameters #######
jobname = "Mueller Test: velocity, no gravity"
GR = 1 # 1 for GR, 0 for Newtonian, enforced by M1_testcase.F90
outdir = "Data_TC5" # output directory
initial_data = "M1test" # "Sedov", "Shocktube", "OSC", "Collapse", "M1test"
gravity_active = 1 # do we want gravity?
ntmax = 10000000 # maximum timestep
tend = 0.05d0 # maximum time
####### Grid parameters #######
geometry = 2 # 1: planar, 2: spherical
gridtype = "blah" # set by M1_testcase.F90
radial_zones = 300 # number of radial zones, enforced by M1_testcase.F90
ghosts1 = 4 # number of ghost cells
####### Hydro Parameters #######
do_hydro = 0 # 1 for hydro on
iorder_hydro = 2 # RK order
cffac = 0.5d0 # CFL factor
reconstruction_method = "tvd" # "tvd", "pc","ppm"
ppm_origin_TVD = 5
ppm_reconstruct_temp = 0 ###### added by csul
tvd_limiter = "MC" # MC, minmod,
flux_type = "HLLE" # HLLE, HLLC
####### EOS parameters #######
eoskey = 4 # hybrid: 1
# poly: 2
# hot nuclear: 3
# ideal: 4
####### Output parameters #######
ntinfo = 100 # stdout
dtout = 1.0d-3 # time between outputs
dtout_scalar = 1.0d-4 # time between scalar out
ntout = -1 # output every timesteps
ntout_scalar = -1 # scalar output
####### Restart parameters #######
ntout_restart = -1 # restart output frequency in timesteps
dtout_restart = 0.005d0 # restart output frequency in time
do_restart = 0
restart_file_name = "blah"
########## M1 settings ##########
do_M1 = 1 #1 for M1 transport scheme
v_order = -1 #enforced by M1_testcases.F90
extraction_radii = -1.0d0 #set by M1_testcases.F90
number_species = 3 #number of neutrino species, limited
#options will work, must be comensurate with table
number_groups = 18 #number of energy groups, limited
#options will work, must be comensurate with table
opacity_table = "/Volumes/vault/EOS/NuLib_LS220_18.h5"
number_eas = 3 #number of opacity variables to read in
M1closure = 'ME' # Closure choice
testcase = 5 #if you want a test case set this to the
#test case number, along with M1test as "initial_data"
include_epannihil_kernels = 0 #for full thermal treatmetn - forced by M1_testcase.F90
include_Ielectron_exp = 0 #Inelastic scattering on electrons, explicit? - forced by M1_testcase.F90
include_Ielectron_imp = 0 #Inelastic scattering on electrons, implicit? - forced by M1_testcase.F90
include_energycoupling_exp = 0 #energy coupling, explicit? - forced by M1_testcase.F90
include_energycoupling_imp = 0 #energy coupling, implicit? - forced by M1_testcase.F90
M1_control = 0 #1 means you must specify the settings for each phase
# in the parameter file, collapse, bounce and postbounce.
####### Neutrino parameters #######
fake_neutrinos = 0 # 1 for ANY neutrino scheme
####### Atmosphere parameters #######
atmo_rho_rel_min = 0.0d0
atmo_rho_abs_min = 1.5d3
atmo_fac = 0.1d0
@@ -0,0 +1,74 @@
##### Job parameters #######
jobname = "Mueller Test: no velocity, gravity"
GR = 1 # 1 for GR, 0 for Newtonian, enforced by M1_testcase.F90
outdir = "Data_TC6" # output directory
initial_data = "M1test" # "Sedov", "Shocktube", "OSC", "Collapse", "M1test"
gravity_active = 1 # do we want gravity?
ntmax = 10000000 # maximum timestep
tend = 0.05d0 # maximum time
####### Grid parameters #######
geometry = 2 # 1: planar, 2: spherical
gridtype = "blah" # set by M1_testcase.F90
radial_zones = 300 # number of radial zones, enforced by M1_testcase.F90
ghosts1 = 4 # number of ghost cells
####### Hydro Parameters #######
do_hydro = 0 # 1 for hydro on
iorder_hydro = 2 # RK order
cffac = 0.5d0 # CFL factor
reconstruction_method = "tvd" # "tvd", "pc","ppm"
ppm_origin_TVD = 5
ppm_reconstruct_temp = 0 ###### added by csul
tvd_limiter = "MC" # MC, minmod,
flux_type = "HLLE" # HLLE, HLLC
####### EOS parameters #######
eoskey = 4 # hybrid: 1
# poly: 2
# hot nuclear: 3
# ideal: 4
####### Output parameters #######
ntinfo = 100 # stdout
dtout = 1.0d-3 # time between outputs
dtout_scalar = 1.0d-4 # time between scalar out
ntout = -1 # output every timesteps
ntout_scalar = -1 # scalar output
####### Restart parameters #######
ntout_restart = -1 # restart output frequency in timesteps
dtout_restart = 0.005d0 # restart output frequency in time
do_restart = 0
restart_file_name = "blah"
########## M1 settings ##########
do_M1 = 1 #1 for M1 transport scheme
v_order = -1 #enforced by M1_testcases.F90
extraction_radii = -1.0d0 #set by M1_testcases.F90
number_species = 3 #number of neutrino species, limited
#options will work, must be comensurate with table
number_groups = 18 #number of energy groups, limited
#options will work, must be comensurate with table
opacity_table = "/Volumes/vault/EOS/NuLib_LS220_18.h5"
number_eas = 3 #number of opacity variables to read in
M1closure = 'ME' # Closure choice
testcase = 6 #if you want a test case set this to the
#test case number, along with M1test as "initial_data"
include_epannihil_kernels = 0 #for full thermal treatmetn - forced by M1_testcase.F90
include_Ielectron_exp = 0 #Inelastic scattering on electrons, explicit? - forced by M1_testcase.F90
include_Ielectron_imp = 0 #Inelastic scattering on electrons, implicit? - forced by M1_testcase.F90
include_energycoupling_exp = 0 #energy coupling, explicit? - forced by M1_testcase.F90
include_energycoupling_imp = 0 #energy coupling, implicit? - forced by M1_testcase.F90
M1_control = 0 #1 means you must specify the settings for each phase
# in the parameter file, collapse, bounce and postbounce.
####### Neutrino parameters #######
fake_neutrinos = 0 # 1 for ANY neutrino scheme
####### Atmosphere parameters #######
atmo_rho_rel_min = 0.0d0
atmo_rho_abs_min = 1.5d3
atmo_fac = 0.1d0
Oops, something went wrong.

0 comments on commit 0fd503e

Please sign in to comment.