This is a _temporary_ repository to be able to start to work on the

code right now. I have put in the public version of Whisky to start from.

Everybody with commit rights should get commit messages (and the other
way around). It should not be a problem to add people to that list, just
ask. I don't want to get into political problems because someone feels
excluded, but I also don't want to give everyone access per se.

Frank



git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/GRHydro_InitData/trunk@3 ac85fae7-cede-4708-beff-ae01c7fa1c26

1 files changed, 84 insertions, 0 deletions

diff --git a/interface.ccl b/interface.ccl
new file mode 100644
index 0000000..48a5456
--- /dev/null
+++ b/interface.ccl
@@ -0,0 +1,84 @@
+# Interface definition for thorn Whisky_Init_Data
+# $Header$
+
+implements: whisky_init_data
+inherits: whisky grid
+
+USES INCLUDE: EOS_Base.inc
+USES INCLUDE: SpaceMask.h
+USES INCLUDE: EOS_Base.h
+
+
+SUBROUTINE SpatialDet(CCTK_REAL IN gxx, CCTK_REAL IN gxy, \
+                      CCTK_REAL IN gxz, CCTK_REAL IN gyy, \
+		      CCTK_REAL IN gyz, CCTK_REAL IN gzz, \
+		      CCTK_REAL OUT det)
+
+
+SUBROUTINE Prim2ConPoly(CCTK_INT IN handle, \
+                        CCTK_REAL IN gxx, CCTK_REAL IN gxy, CCTK_REAL IN gxz, \
+			CCTK_REAL IN gyy, CCTK_REAL IN gyz, CCTK_REAL IN gzz, \
+			CCTK_REAL IN det, CCTK_REAL OUT dens, \
+			CCTK_REAL OUT sx, CCTK_REAL OUT sy, \
+			CCTK_REAL OUT sz, CCTK_REAL OUT tau, \
+			CCTK_REAL IN rho, CCTK_REAL IN velx, \
+			CCTK_REAL IN vely, \
+			CCTK_REAL IN velz, CCTK_REAL OUT epsilon, \
+			CCTK_REAL OUT press, CCTK_REAL OUT w_lorentz)
+
+
+SUBROUTINE Prim2ConGen(CCTK_INT IN handle, \
+                          CCTK_REAL IN gxx, CCTK_REAL IN gxy, \
+                          CCTK_REAL IN gxz, CCTK_REAL IN gyy, \
+                          CCTK_REAL IN gyz, CCTK_REAL IN gzz, \
+                          CCTK_REAL IN det, CCTK_REAL OUT dens, \
+                          CCTK_REAL OUT sx, CCTK_REAL OUT sy, \
+                          CCTK_REAL OUT sz, CCTK_REAL OUT tau, \
+                          CCTK_REAL IN rho, CCTK_REAL IN velx, \
+                          CCTK_REAL IN vely, \
+                          CCTK_REAL IN velz, CCTK_REAL IN epsilon, \
+                          CCTK_REAL OUT press, CCTK_REAL OUT w_lorentz)
+
+SUBROUTINE Con2PrimPoly(CCTK_INT IN handle, CCTK_REAL OUT dens, \
+                        CCTK_REAL OUT sx, CCTK_REAL OUT sy, \
+			CCTK_REAL OUT sz, CCTK_REAL OUT tau, \
+			CCTK_REAL OUT rho, CCTK_REAL OUT velx, \
+			CCTK_REAL OUT vely, CCTK_REAL OUT velz, \
+			CCTK_REAL OUT epsilon, CCTK_REAL OUT press, \
+			CCTK_REAL OUT w_lorentz, CCTK_REAL IN uxx, \
+			CCTK_REAL IN uxy, CCTK_REAL IN uxz, CCTK_REAL IN uyy, \
+			CCTK_REAL IN uyz, CCTK_REAL IN uzz, CCTK_REAL IN det, \
+			CCTK_REAL IN x, CCTK_REAL IN y, CCTK_REAL IN z, \
+			CCTK_REAL IN r, CCTK_REAL IN rho_min, \
+			CCTK_INT IN whisky_reflevel, CCTK_REAL IN whisky_C2P_failed)
+
+
+USES FUNCTION SpatialDet
+USES FUNCTION Prim2ConPoly
+USES FUNCTION Prim2ConGen
+USES FUNCTION Con2PrimPoly
+
+protected:
+
+CCTK_REAL simple_wave_grid_functions TYPE=GF TIMELEVELS=1 tags='checkpoint="no"'
+{
+  simple_tmp
+  c_s
+} "1D arrays for the simple-wave routine"
+
+CCTK_REAL simple_wave_scalars TYPE=scalar
+{
+     simple_rho_0
+     simple_eps_0
+} "values at v=0"
+
+CCTK_REAL simple_wave_output TYPE=GF TIMELEVELS=1 tags='checkpoint="no"'
+{
+  simple_rho
+  simple_eps
+#  simple_entropy
+} "output variables for the simple-wave routine"
+
+private:
+
+int whisky_init_data_reflevel type = SCALAR tags='checkpoint="no"' "Refinement level Whisky is working on right now"