/*@@
   @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