diff options
author | cott <cott@8e189c6b-2ab8-4400-aa02-70a9cfce18b9> | 2010-06-29 14:59:02 +0000 |
---|---|---|
committer | cott <cott@8e189c6b-2ab8-4400-aa02-70a9cfce18b9> | 2010-06-29 14:59:02 +0000 |
commit | 8b44bae2014fefc6baca909e7c363c877fb5e359 (patch) | |
tree | 4a8a7437badbaa4559717db00bcb277b91937940 |
* EOS_Omni thorn:
Initial implementation providing a polytrope and a gamma-law EOS for testing.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEOS/EOS_Omni/EOS_Omni@1 8e189c6b-2ab8-4400-aa02-70a9cfce18b9
-rw-r--r-- | README | 24 | ||||
-rw-r--r-- | doc/documentation.tex | 144 | ||||
-rw-r--r-- | interface.ccl | 89 | ||||
-rw-r--r-- | param.ccl | 26 | ||||
-rw-r--r-- | schedule.ccl | 8 | ||||
-rw-r--r-- | src/EOS_Omni_Module.F90 | 22 | ||||
-rw-r--r-- | src/EOS_Omni_SingleVarCalls.F90 | 333 | ||||
-rw-r--r-- | src/EOS_Omni_Startup.F90 | 18 | ||||
-rw-r--r-- | src/make.code.defn | 7 |
9 files changed, 671 insertions, 0 deletions
@@ -0,0 +1,24 @@ +Cactus Code Thorn EOS_Omni +Author(s) : Christian D. Ott <cott@tapir.caltech.edu> +Maintainer(s): Christian D. Ott <cott@tapir.caltech.edu> +Licence : LGPL +-------------------------------------------------------------------------- + +This Thorn provides a one-fits-all EOS handler. +All EOS are actually implemented right inside this single +thorn without crazy interfaces and lots of function calls. + +eoskey controls the type of EOS being used: + +eoskey: + +1 --- polytropic EOS +2 --- gamma-law EOS +3 --- hybrid EOS +4 --- finite-T microphysical NSE EOS + + +keyerr settings: + +-1 -- called EOS routine with keytemp=1, but + routine does not allow such a call diff --git a/doc/documentation.tex b/doc/documentation.tex new file mode 100644 index 0000000..3912bd2 --- /dev/null +++ b/doc/documentation.tex @@ -0,0 +1,144 @@ +% *======================================================================* +% Cactus Thorn template for ThornGuide documentation +% Author: Ian Kelley +% Date: Sun Jun 02, 2002 +% $Header$ +% +% Thorn documentation in the latex file doc/documentation.tex +% will be included in ThornGuides built with the Cactus make system. +% The scripts employed by the make system automatically include +% pages about variables, parameters and scheduling parsed from the +% relevant thorn CCL files. +% +% This template contains guidelines which help to assure that your +% documentation will be correctly added to ThornGuides. More +% information is available in the Cactus UsersGuide. +% +% Guidelines: +% - Do not change anything before the line +% % START CACTUS THORNGUIDE", +% except for filling in the title, author, date, etc. fields. +% - Each of these fields should only be on ONE line. +% - Author names should be separated with a \\ or a comma. +% - You can define your own macros, but they must appear after +% the START CACTUS THORNGUIDE line, and must not redefine standard +% latex commands. +% - To avoid name clashes with other thorns, 'labels', 'citations', +% 'references', and 'image' names should conform to the following +% convention: +% ARRANGEMENT_THORN_LABEL +% For example, an image wave.eps in the arrangement CactusWave and +% thorn WaveToyC should be renamed to CactusWave_WaveToyC_wave.eps +% - Graphics should only be included using the graphicx package. +% More specifically, with the "\includegraphics" command. Do +% not specify any graphic file extensions in your .tex file. This +% will allow us to create a PDF version of the ThornGuide +% via pdflatex. +% - References should be included with the latex "\bibitem" command. +% - Use \begin{abstract}...\end{abstract} instead of \abstract{...} +% - Do not use \appendix, instead include any appendices you need as +% standard sections. +% - For the benefit of our Perl scripts, and for future extensions, +% please use simple latex. +% +% *======================================================================* +% +% Example of including a graphic image: +% \begin{figure}[ht] +% \begin{center} +% \includegraphics[width=6cm]{MyArrangement_MyThorn_MyFigure} +% \end{center} +% \caption{Illustration of this and that} +% \label{MyArrangement_MyThorn_MyLabel} +% \end{figure} +% +% Example of using a label: +% \label{MyArrangement_MyThorn_MyLabel} +% +% Example of a citation: +% \cite{MyArrangement_MyThorn_Author99} +% +% Example of including a reference +% \bibitem{MyArrangement_MyThorn_Author99} +% {J. Author, {\em The Title of the Book, Journal, or periodical}, 1 (1999), +% 1--16. {\tt http://www.nowhere.com/}} +% +% *======================================================================* + +% If you are using CVS use this line to give version information +% $Header$ + +\documentclass{article} + +% Use the Cactus ThornGuide style file +% (Automatically used from Cactus distribution, if you have a +% thorn without the Cactus Flesh download this from the Cactus +% homepage at www.cactuscode.org) +\usepackage{../../../../doc/latex/cactus} + +\begin{document} + +% The author of the documentation +\author{Christian D. Ott \textless cott@tapir.caltech.edu\textgreater} + +% The title of the document (not necessarily the name of the Thorn) +\title{EOS\_Omni} + +% the date your document was last changed, if your document is in CVS, +% please use: +% \date{$ $Date: 2004-01-07 12:12:39 -0800 (Wed, 07 Jan 2004) $ $} +\date{June 27 2010} + +\maketitle + +% Do not delete next line +% START CACTUS THORNGUIDE + +% Add all definitions used in this documentation here +% \def\mydef etc + +% Add an abstract for this thorn's documentation +\begin{abstract} + +\end{abstract} + +% The following sections are suggestive only. +% Remove them or add your own. + +\section{Introduction} + +\section{Physical System} + +\section{Numerical Implementation} + +\section{Using This Thorn} + +\subsection{Obtaining This Thorn} + +\subsection{Basic Usage} + +\subsection{Special Behaviour} + +\subsection{Interaction With Other Thorns} + +\subsection{Examples} + +\subsection{Support and Feedback} + +\section{History} + +\subsection{Thorn Source Code} + +\subsection{Thorn Documentation} + +\subsection{Acknowledgements} + + +\begin{thebibliography}{9} + +\end{thebibliography} + +% Do not delete next line +% END CACTUS THORNGUIDE + +\end{document} diff --git a/interface.ccl b/interface.ccl new file mode 100644 index 0000000..3eb5752 --- /dev/null +++ b/interface.ccl @@ -0,0 +1,89 @@ +# Interface definition for thorn EOS_Omni + +implements: EOS_Omni + +inherits: + +void FUNCTION EOS_Omni_press(CCTK_INT IN eoskey, \ + CCTK_INT IN havetemp, \ + CCTK_INT IN npoints, \ + CCTK_REAL IN ARRAY rho, \ + CCTK_REAL INOUT ARRAY eps, \ + CCTK_REAL INOUT ARRAY temp, \ + CCTK_REAL IN ARRAY ye, \ + CCTK_REAL OUT ARRAY press, \ + CCTK_INT OUT ARRAY keyerr, \ + CCTK_INT OUT anyerr) + +PROVIDES FUNCTION EOS_Omni_press WITH EOS_Omni_EOS_Press LANGUAGE Fortran + + +void FUNCTION EOS_Omni_DPressByDEps(CCTK_INT IN eoskey, \ + CCTK_INT IN havetemp, \ + CCTK_INT IN npoints, \ + CCTK_REAL IN ARRAY rho, \ + CCTK_REAL INOUT ARRAY eps, \ + CCTK_REAL INOUT ARRAY temp, \ + CCTK_REAL IN ARRAY ye, \ + CCTK_REAL OUT ARRAY DPressByDEps, \ + CCTK_INT OUT ARRAY keyerr, \ + CCTK_INT OUT anyerr) + +PROVIDES FUNCTION EOS_Omni_DPressByDEps WITH EOS_Omni_EOS_DPressByDEps LANGUAGE Fortran + + +void FUNCTION EOS_Omni_DPressByDRho(CCTK_INT IN eoskey, \ + CCTK_INT IN havetemp, \ + CCTK_INT IN npoints, \ + CCTK_REAL IN ARRAY rho, \ + CCTK_REAL INOUT ARRAY eps, \ + CCTK_REAL INOUT ARRAY temp, \ + CCTK_REAL IN ARRAY ye, \ + CCTK_REAL OUT ARRAY DPressByDRho, \ + CCTK_INT OUT ARRAY keyerr, \ + CCTK_INT OUT anyerr) + +PROVIDES FUNCTION EOS_Omni_DPressByDRho WITH EOS_Omni_EOS_DPressByDRho LANGUAGE Fortran + + +void FUNCTION EOS_Omni_cs2(CCTK_INT IN eoskey, \ + CCTK_INT IN havetemp, \ + CCTK_INT IN npoints, \ + CCTK_REAL IN ARRAY rho, \ + CCTK_REAL INOUT ARRAY eps, \ + CCTK_REAL INOUT ARRAY temp, \ + CCTK_REAL IN ARRAY ye, \ + CCTK_REAL OUT ARRAY cs2, \ + CCTK_INT OUT ARRAY keyerr, \ + CCTK_INT OUT anyerr) + +PROVIDES FUNCTION EOS_Omni_cs2 WITH EOS_Omni_EOS_cs2 LANGUAGE Fortran + + +void FUNCTION EOS_Omni_EpsFromPress(CCTK_INT IN eoskey, \ + CCTK_INT IN havetemp, \ + CCTK_INT IN npoints, \ + CCTK_REAL IN ARRAY rho, \ + CCTK_REAL INOUT ARRAY eps, \ + CCTK_REAL INOUT ARRAY temp, \ + CCTK_REAL IN ARRAY ye, \ + CCTK_REAL IN ARRAY press, \ + CCTK_REAL OUT ARRAY xeps, \ + CCTK_INT OUT ARRAY keyerr, \ + CCTK_INT OUT anyerr) + +PROVIDES FUNCTION EOS_Omni_EpsFromPress WITH EOS_Omni_EOS_eps_from_press LANGUAGE Fortran + + +void FUNCTION EOS_Omni_RestMassDensityFromEpsPress(CCTK_INT IN eoskey, \ + CCTK_INT IN havetemp, \ + CCTK_INT IN npoints, \ + CCTK_REAL OUT ARRAY rho, \ + CCTK_REAL IN ARRAY eps, \ + CCTK_REAL INOUT ARRAY temp, \ + CCTK_REAL IN ARRAY ye, \ + CCTK_REAL IN ARRAY press, \ + CCTK_INT OUT ARRAY keyerr, \ + CCTK_INT OUT anyerr) + +PROVIDES FUNCTION EOS_Omni_RestMassDensityFromEpsPress WITH EOS_Omni_EOS_RestMassDensityFromEpsPress LANGUAGE Fortran diff --git a/param.ccl b/param.ccl new file mode 100644 index 0000000..6b8acf7 --- /dev/null +++ b/param.ccl @@ -0,0 +1,26 @@ +# Parameter definitions for thorn EOS_Omni + + +restricted: + +# poly EOS +CCTK_REAL poly_gamma "Adiabatic Index for poly EOS" +{ + : :: "" +} 2.0 + +REAL poly_k "Polytropic constant in c=G=Msun=1" +{ + : :: "" +} 100.0 + +# gamma-law EOS +CCTK_REAL gl_gamma "Adiabatic Index for gamma-law EOS" +{ + : :: "" +} 2.0 + +REAL gl_k "Polytropic constant in c=G=Msun=1 for gamma-law EOS" +{ + : :: "" +} 100.0 diff --git a/schedule.ccl b/schedule.ccl new file mode 100644 index 0000000..a9cbd79 --- /dev/null +++ b/schedule.ccl @@ -0,0 +1,8 @@ +# Schedule definitions for thorn EOS_Omni + + +schedule EOS_Omni_Startup AT WRAGH +{ + LANG: Fortran + OPTIONS: global +} "Set up conversion factors and other fun stuff" diff --git a/src/EOS_Omni_Module.F90 b/src/EOS_Omni_Module.F90 new file mode 100644 index 0000000..702d841 --- /dev/null +++ b/src/EOS_Omni_Module.F90 @@ -0,0 +1,22 @@ +module EOS_Omni_Module + + implicit none + + real*8,parameter :: rho_gf = 1.61930347d-18 + real*8,parameter :: press_gf = 1.80171810d-39 + real*8,parameter :: eps_gf = 1.11265006d-21 + real*8,parameter :: time_gf = 2.03001708d+05 + real*8,parameter :: mass_gf = 5.02765209d-34 + real*8,parameter :: length_gf = 6.77140812d-06 + + real*8,parameter :: inv_rho_gf = 6.17549470205236d17 + real*8,parameter :: inv_press_gf = 5.55025783445257d38 + real*8,parameter :: inv_eps_gf = 8.98755175549085d20 + real*8,parameter :: inv_time_gf = 4.92606692747629d-6 + real*8,parameter :: inv_mass_gf = 1.98899999860571d33 + real*8,parameter :: inv_length_gf = 147679.77092481d0 + + real*8 :: poly_k_cgs = 0.0d0 + real*8 :: gl_k_cgs = 0.0d0 + +end module EOS_Omni_Module diff --git a/src/EOS_Omni_SingleVarCalls.F90 b/src/EOS_Omni_SingleVarCalls.F90 new file mode 100644 index 0000000..aa0f79d --- /dev/null +++ b/src/EOS_Omni_SingleVarCalls.F90 @@ -0,0 +1,333 @@ +#include "cctk.h" +#include "cctk_Parameters.h" +#include "cctk_Arguments.h" +#include "cctk_Functions.h" + +! eoskey: +! 1 --- polytropic EOS +! 2 --- gamma-law EOS +! 3 --- hybrid EOS +! 4 --- finite-T microphysical NSE EOS + +subroutine EOS_Omni_EOS_Press(eoskey,keytemp,npoints,& + rho,eps,temp,ye,press,keyerr,anyerr) + + use EOS_Omni_Module + implicit none + DECLARE_CCTK_PARAMETERS + + CCTK_INT, intent(in) :: eoskey,keytemp,npoints + CCTK_INT, intent(out) :: keyerr(npoints) + CCTK_INT, intent(out) :: anyerr + CCTK_REAL, intent(in) :: rho(npoints),ye(npoints) + CCTK_REAL, intent(inout) :: eps(npoints), temp(npoints) + CCTK_REAL, intent(out) :: press(npoints) + + ! local vars + integer :: i + character(256) :: warnstring + + anyerr = 0 + keyerr(:) = 0 + + select case (eoskey) + case (1) + ! polytropic EOS + if(keytemp.eq.1) then + do i=1,npoints + eps(i) = press_gf * poly_k_cgs * & + (rho(i)*inv_rho_gf)**(poly_gamma) / & + (poly_gamma - 1.0d0) / rho(i) + enddo + endif + do i=1,npoints + press(i) = press_gf * poly_k_cgs * & + (rho(i)*inv_rho_gf)**poly_gamma + enddo + case (2) + ! gamma-law EOS + if(keytemp.eq.1) then + do i=1,npoints + eps(i) = press_gf * gl_k_cgs * & + (rho(i)*inv_rho_gf)**(gl_gamma) / & + (gl_gamma - 1.0d0) / rho(i) + enddo + endif + do i=1,npoints + press(i) = (gl_gamma - 1.0d0) * rho(i) * eps(i) + enddo + + case DEFAULT + write(warnstring,*) "eoskey ",eoskey," not implemented!" + call CCTK_WARN(0,warnstring) + end select + +end subroutine EOS_Omni_EOS_Press + + + +subroutine EOS_Omni_EOS_DPressByDRho(eoskey,keytemp,npoints,& + rho,eps,temp,ye,dpdrhoe,keyerr,anyerr) + + use EOS_Omni_Module + implicit none + DECLARE_CCTK_PARAMETERS + + CCTK_INT, intent(in) :: eoskey,keytemp,npoints + CCTK_INT, intent(out) :: keyerr(npoints) + CCTK_INT, intent(out) :: anyerr + CCTK_REAL, intent(in) :: rho(npoints),ye(npoints) + CCTK_REAL, intent(inout) :: eps(npoints), temp(npoints) + CCTK_REAL, intent(out) :: dpdrhoe(npoints) + + ! local vars + integer :: i + character(256) :: warnstring + + anyerr = 0 + keyerr(:) = 0 + + select case (eoskey) + case (1) + ! polytropic EOS + if(keytemp.eq.1) then + do i=1,npoints + eps(i) = press_gf * poly_k_cgs * & + (rho(i)*inv_rho_gf)**(poly_gamma) / & + (poly_gamma - 1.0d0) / rho(i) + enddo + endif + do i=1,npoints + dpdrhoe(i) = press_gf * poly_k_cgs * & + poly_gamma * inv_rho_gf * & + (rho(i)*inv_rho_gf) ** (poly_gamma - 1.d0) + enddo + case (2) + ! gamma-law EOS + if(keytemp.eq.1) then + do i=1,npoints + eps(i) = press_gf * gl_k_cgs * & + (rho(i)*inv_rho_gf)**(gl_gamma) / & + (gl_gamma - 1.0d0) / rho(i) + enddo + endif + do i=1,npoints + dpdrhoe(i) = (gl_gamma-1.0d0) * & + eps(i) + enddo + + case DEFAULT + write(warnstring,*) "eoskey ",eoskey," not implemented!" + call CCTK_WARN(0,warnstring) + end select + +end subroutine EOS_Omni_EOS_DPressByDRho + +subroutine EOS_Omni_EOS_DPressByDEps(eoskey,keytemp,npoints,& + rho,eps,temp,ye,dpdepsrho,keyerr,anyerr) + + use EOS_Omni_Module + implicit none + DECLARE_CCTK_PARAMETERS + + CCTK_INT, intent(in) :: eoskey,keytemp,npoints + CCTK_INT, intent(out) :: keyerr(npoints) + CCTK_INT, intent(out) :: anyerr + CCTK_REAL, intent(in) :: rho(npoints),ye(npoints) + CCTK_REAL, intent(inout) :: eps(npoints), temp(npoints) + CCTK_REAL, intent(out) :: dpdepsrho(npoints) + + ! local vars + integer :: i + character(256) :: warnstring + + anyerr = 0 + keyerr(:) = 0 + + select case (eoskey) + case (1) + ! polytropic EOS + if(keytemp.eq.1) then + do i=1,npoints + eps(i) = press_gf * poly_k_cgs * & + (rho(i)*inv_rho_gf)**(poly_gamma) / & + (poly_gamma - 1.0d0) / rho(i) + enddo + endif + do i=1,npoints + dpdepsrho(i) = 0.0d0 + enddo + case (2) + ! gamma-law EOS + if(keytemp.eq.1) then + do i=1,npoints + eps(i) = press_gf * gl_k_cgs * & + (rho(i)*inv_rho_gf)**(gl_gamma) / & + (gl_gamma - 1.0d0) / rho(i) + enddo + endif + do i=1,npoints + dpdepsrho(i) = (gl_gamma - 1.0d0) * & + rho(i) + enddo + + case DEFAULT + write(warnstring,*) "eoskey ",eoskey," not implemented!" + call CCTK_WARN(0,warnstring) + end select + +end subroutine EOS_Omni_EOS_DPressByDEps + +subroutine EOS_Omni_EOS_cs2(eoskey,keytemp,npoints,& + rho,eps,temp,ye,cs2,keyerr,anyerr) + + use EOS_Omni_Module + implicit none + DECLARE_CCTK_PARAMETERS + + CCTK_INT, intent(in) :: eoskey,keytemp,npoints + CCTK_INT, intent(out) :: keyerr(npoints) + CCTK_INT, intent(out) :: anyerr + CCTK_REAL, intent(in) :: rho(npoints),ye(npoints) + CCTK_REAL, intent(inout) :: eps(npoints), temp(npoints) + CCTK_REAL, intent(out) :: cs2(npoints) + + ! local vars + integer :: i + character(256) :: warnstring + real*8 :: xpress,xdpdrhoe,xdpderho + + anyerr = 0 + keyerr(:) = 0 + + select case (eoskey) + case (1) + ! polytropic EOS + if(keytemp.eq.1) then + do i=1,npoints + eps(i) = press_gf * poly_k_cgs * & + (rho(i)*inv_rho_gf)**(poly_gamma) / & + (poly_gamma - 1.0d0) / rho(i) + enddo + endif + do i=1,npoints + xpress = press_gf*poly_k_cgs * & + (rho(i)*inv_rho_gf)**(poly_gamma) + cs2(i) = poly_gamma * xpress / rho(i) / & + (1 + eps(i) + xpress/rho(i)) + enddo + case (2) + ! gamma-law EOS + if(keytemp.eq.1) then + do i=1,npoints + eps(i) = press_gf * gl_k_cgs * & + (rho(i)*inv_rho_gf)**(gl_gamma) / & + (gl_gamma - 1.0d0) / rho(i) + enddo + endif + do i=1,npoints + xpress = (gl_gamma-1.0d0)*rho(i)*eps(i) + xdpdrhoe = (gl_gamma-1.0d0)*eps(i) + xdpderho = (gl_gamma-1.0d0)*rho(i) + cs2(i) = (xdpdrhoe + xpress * xdpderho / (rho(i)**2)) / & + (1.0d0 + eps(i) + xpress/rho(i)) + enddo + + case DEFAULT + write(warnstring,*) "eoskey ",eoskey," not implemented!" + call CCTK_WARN(0,warnstring) + end select + +end subroutine EOS_Omni_EOS_cs2 + + +subroutine EOS_Omni_EOS_eps_from_press(eoskey,keytemp,npoints,& + rho,eps,temp,ye,press,xeps,keyerr,anyerr) + + use EOS_Omni_Module + implicit none + DECLARE_CCTK_PARAMETERS + + CCTK_INT, intent(in) :: eoskey,keytemp,npoints + CCTK_INT, intent(out) :: keyerr(npoints) + CCTK_INT, intent(out) :: anyerr + CCTK_REAL, intent(in) :: rho(npoints),ye(npoints),press(npoints) + CCTK_REAL, intent(inout) :: eps(npoints), temp(npoints) + CCTK_REAL, intent(out) :: xeps(npoints) + + ! local vars + integer :: i + character(256) :: warnstring + + if(keytemp.eq.1) then + anyerr = 1 + keyerr(:) = -1 + else + anyerr = 0 + keyerr(:) = 0 + endif + + select case (eoskey) + case (1) + ! polytropic EOS + do i=1,npoints + xeps(i) = press(i) / (poly_gamma - 1.0d0) / rho(i) + eps(i) = xeps(i) + enddo + case (2) + ! gamma-law EOS + do i=1,npoints + xeps(i) = press(i) / (gl_gamma - 1.0d0) / rho(i) + eps(i) = xeps(i) + enddo + case DEFAULT + write(warnstring,*) "eoskey ",eoskey," not implemented!" + call CCTK_WARN(0,warnstring) + end select + +end subroutine EOS_Omni_EOS_eps_from_press + + +subroutine EOS_Omni_EOS_RestMassDensityFromEpsPress(eoskey,keytemp,npoints,& + rho,eps,temp,ye,press,keyerr,anyerr) + + use EOS_Omni_Module + implicit none + DECLARE_CCTK_PARAMETERS + + CCTK_INT, intent(in) :: eoskey,keytemp,npoints + CCTK_INT, intent(out) :: keyerr(npoints) + CCTK_INT, intent(out) :: anyerr + CCTK_REAL, intent(in) :: ye(npoints),press(npoints),eps(npoints) + CCTK_REAL, intent(inout) :: rho(npoints),temp(npoints) + + + ! local vars + integer :: i + character(256) :: warnstring + + if(keytemp.eq.1) then + anyerr = 1 + keyerr(:) = -1 + else + anyerr = 0 + keyerr(:) = 0 + endif + + select case (eoskey) + case (1) + ! polytropic EOS + do i=1,npoints + rho(i) = press(i) / ((poly_gamma - 1.0d0)*eps(i)) + enddo + case (2) + ! gamma-law EOS + do i=1,npoints + rho(i) = press(i) / ((gl_gamma - 1.0d0)*eps(i)) + enddo + case DEFAULT + write(warnstring,*) "eoskey ",eoskey," not implemented!" + call CCTK_WARN(0,warnstring) + end select + +end subroutine EOS_Omni_EOS_RestMassDensityFromEpsPress diff --git a/src/EOS_Omni_Startup.F90 b/src/EOS_Omni_Startup.F90 new file mode 100644 index 0000000..7a9a2f7 --- /dev/null +++ b/src/EOS_Omni_Startup.F90 @@ -0,0 +1,18 @@ +#include "cctk.h" +#include "cctk_Parameters.h" +#include "cctk_Arguments.h" +#include "cctk_Functions.h" + +subroutine EOS_Omni_Startup(CCTK_ARGUMENTS) + + use EOS_Omni_Module + implicit none + + DECLARE_CCTK_PARAMETERS + DECLARE_CCTK_ARGUMENTS + + poly_k_cgs = poly_k * rho_gf**poly_gamma / press_gf + + gl_k_cgs = gl_k * rho_gf**gl_gamma / press_gf + +end subroutine EOS_Omni_Startup diff --git a/src/make.code.defn b/src/make.code.defn new file mode 100644 index 0000000..835441d --- /dev/null +++ b/src/make.code.defn @@ -0,0 +1,7 @@ +# Main make.code.defn file for thorn EOS_Omni + +# Source files in this directory +SRCS = EOS_Omni_Module.F90 EOS_Omni_Startup.F90 EOS_Omni_SingleVarCalls.F90 + +# Subdirectories containing source files +SUBDIRS = |