aboutsummaryrefslogtreecommitdiff
path: root/src/nuc_eos/eosmodule.F90
diff options
context:
space:
mode:
Diffstat (limited to 'src/nuc_eos/eosmodule.F90')
-rw-r--r--src/nuc_eos/eosmodule.F90111
1 files changed, 0 insertions, 111 deletions
diff --git a/src/nuc_eos/eosmodule.F90 b/src/nuc_eos/eosmodule.F90
deleted file mode 100644
index 174d817..0000000
--- a/src/nuc_eos/eosmodule.F90
+++ /dev/null
@@ -1,111 +0,0 @@
-#include "cctk.h"
-#include "cctk_Parameters.h"
-#include "cctk_Arguments.h"
-#include "cctk_Functions.h"
-
- module eosmodule
-
- implicit none
-
- integer,save :: nrho,ntemp,nye
-
- integer,save :: warn_from !warn from given reflevel
-
- real*8 :: energy_shift = 0.0d0
-
- real*8 :: precision = 1.0d-9
-
-! min-max values:
- real*8 :: eos_rhomin,eos_rhomax
- real*8 :: eos_yemin,eos_yemax
- real*8 :: eos_tempmin,eos_tempmax
-
- real*8 :: t_max_hack = 240.0d0
-
-! basics
- integer, parameter :: nvars = 19
- real*8,pointer :: alltables(:,:,:,:)
- real*8,allocatable :: epstable(:,:,:)
- ! index variable mapping:
- ! 1 -> logpress
- ! 2 -> logenergy
- ! 3 -> entropy
- ! 4 -> munu
- ! 5 -> cs2
- ! 6 -> dedT
- ! 7 -> dpdrhoe
- ! 8 -> dpderho
- ! 9 -> muhat
- ! 10 -> mu_e
- ! 11 -> mu_p
- ! 12 -> mu_n
- ! 13 -> xa
- ! 14 -> xh
- ! 15 -> xn
- ! 16 -> xp
- ! 17 -> abar
- ! 18 -> zbar
- ! 19 -> gamma
-
- real*8,pointer :: logrho(:)
- real*8,pointer :: logtemp(:)
- real*8,pointer :: ye(:)
-
-! constants
- real*8,save :: mev_to_erg = 1.60217733d-6
- real*8,save :: amu_cgs = 1.66053873d-24
- real*8,save :: amu_mev = 931.49432d0
- real*8,save :: pi = 3.14159265358979d0
- real*8,save :: ggrav = 6.672d-8
- real*8,save :: temp_mev_to_kelvin = 1.1604447522806d10
- real*8,save :: clight = 2.99792458d10
- real*8,save :: kb_erg = 1.380658d-16
- real*8,save :: kb_mev = 8.61738568d-11
-
-
- end module eosmodule
-
-! This function is called from within readtable.c
-! It creates pointers to the arrays given as parameters into arrays
-! of the eos module.
-subroutine setup_eosmodule(nrho_, ntemp_, nye_, alltables_, logrho_, logtemp_, ye_, energy_shift_)
- use eosmodule
- DECLARE_CCTK_PARAMETERS
-
- CCTK_INT :: nrho_, ntemp_, nye_
- CCTK_REAL, TARGET :: alltables_(nrho_, ntemp_, nye_, 19)
- CCTK_REAL, TARGET :: logrho_(nrho_)
- CCTK_REAL, TARGET :: logtemp_(ntemp_)
- CCTK_REAL, TARGET :: ye_(nye_)
- CCTK_REAL :: energy_shift_
-
- nrho = nrho_
- ntemp = ntemp_
- nye = nye_
-
- alltables => alltables_
- logrho => logrho_
- logtemp => logtemp_
- ye => ye_
-
- allocate(epstable(nrho,ntemp,nye))
- epstable(1:nrho,1:ntemp,1:nye) = 10.0d0**alltables(1:nrho,1:ntemp,1:nye,2)
-
- energy_shift = energy_shift_
- if(do_energy_shift.ne.1) then
- energy_shift = 0.0d0
- endif
-
- ! set min-max values:
- eos_rhomin = 10.0d0**logrho(1)
- eos_rhomax = 10.0d0**logrho(nrho)
-
- eos_yemin = ye(1)
- eos_yemax = ye(nye)
-
- eos_tempmin = 10.0d0**logtemp(1)
- eos_tempmax = 10.0d0**logtemp(ntemp)
-
-
-end subroutine setup_eosmodule
-