aboutsummaryrefslogtreecommitdiff
path: root/src/Initial.F77
diff options
context:
space:
mode:
Diffstat (limited to 'src/Initial.F77')
-rw-r--r--src/Initial.F7751
1 files changed, 51 insertions, 0 deletions
diff --git a/src/Initial.F77 b/src/Initial.F77
new file mode 100644
index 0000000..dca2e8b
--- /dev/null
+++ b/src/Initial.F77
@@ -0,0 +1,51 @@
+ /*@@
+ @file WaveToy.F
+ @date
+ @author Tom Goodale
+ @desc
+ Solve the 3D Wave Equation
+ @enddesc
+ @@*/
+
+#include "cctk.h"
+#include "cctk_parameters.h"
+#include "cctk_arguments.h"
+
+
+
+ /*@@
+ @routine WaveToy_Initial
+ @date
+ @author Tom Goodale
+ @desc
+ Set up initial properties for Wave Toy evolver
+ @enddesc
+ @calls
+ @calledby
+ @history
+
+ @endhistory
+
+@@*/
+
+ subroutine WaveToy_Initial(CCTK_FARGUMENTS)
+
+ implicit none
+
+ DECLARE_CCTK_FARGUMENTS
+ DECLARE_CCTK_PARAMETERS
+
+ CCTK_REAL :: min_delta,dx,dy,dz
+
+ dx = cctk_delta_space(1)
+ dy = cctk_delta_space(2)
+ dz = cctk_delta_space(3)
+
+c Calculate timestep
+c ------------------
+ min_delta = min(dx,dy,dz)
+ cctk_delta_time = dtfac*min_delta
+
+ end subroutine wavetoy_initial
+
+