c  Mahc_Init_Data: Various Initial Data Routines for Mahc_Evolve
c  Copyright (C) 2001 Mark Miller

/*@@
  @file      test_primvar.F
  @date      March 2000
  @author    Mark Miller
  @desc
     Testing the primative variable solver.
  @enddesc
@@*/

#include "cctk.h"
#include "cctk_Parameters.h"
#include "cctk_Arguments.h"

/*@@
  @routine    test_primvar
  @date       March 2000
  @author     Mark Miller
  @desc
     Testing the primative variable solver.
  @enddesc
  @calls
  @calledby
  @history
  @endhistory
@@*/


      subroutine test_primvar(CCTK_FARGUMENTS)
      implicit  none
      DECLARE_CCTK_FARGUMENTS
      DECLARE_CCTK_PARAMETERS
      INTEGER CCTK_Equals

c     Local Vars
      integer didit,i,j,k,nx,ny,nz
      logical some_hydro_dust,mahc_eos_tables,mahc_eos_polytrope
      logical mahc_eos_ideal
      CCTK_REAL det,uxx,uxy,uxz,uyy,uyz,uzz
      CCTK_REAL dens_send,sx_send,sy_send,sz_send,tau_send
      CCTK_REAL rho_send,velx_send,vely_send,velz_send,eps_send
      CCTK_REAL press_send,w_lorentz_send,x_send,y_send,z_send,r_send

      nx = cctk_lsh(1)
      ny = cctk_lsh(2)
      nz = cctk_lsh(3)
      some_hydro_dust = CCTK_EQUALS(some_hydro,'MahcHydroDust')
      mahc_eos_tables = CCTK_EQUALS(mahc_eos,'tables')
      mahc_eos_polytrope = CCTK_EQUALS(mahc_eos,'polytrope')
      mahc_eos_ideal = CCTK_EQUALS(mahc_eos,'ideal_fluid')

      x_send = 0.0d0
      y_send = 0.0d0
      z_send = 0.0d0
      r_send = 0.0d0

      det = 1.0d0
      uxx = 1.0d0
      uyy = 1.0d0 
      uzz = 1.0d0
      uxy = 0.0d0
      uxz = 0.0d0
      uyz = 0.0d0

      dens_send = 1.29047362d0
      sx_send   = 0.166666658d0
      sy_send   = 0.166666658d0
      sz_send   = 0.166666658d0
      tau_send  = 0.484123939d0

      rho_send = 1.0d0
      velx_send = 0.0d0
      vely_send = 0.0d0
      velz_send = 0.0d0
      eps_send = 1.0d-6
      press_send = 6.666666666666667d-7
      w_lorentz_send = 1.0d0

      write(*,*) 'test_primvar: calling get_prim_var_pt...'
      call get_prim_var_pt(dens_send,sx_send,sy_send,sz_send,
     .   tau_send,rho_send,velx_send,vely_send,velz_send,
     .   eps_send,press_send,w_lorentz_send,mahc_eos_gamma,
     .   uxx,uxy,uxz,uyy,uyz,uzz,det,x_send,y_send,z_send,r_send,
     .   mahc_eps_min,mahc_densmin,
     .      mahc_eos_k,mahc_perc_ptol,mahc_del_ptol,some_hydro_dust,
     .      mahc_eos_tables,mahc_eos_polytrope,mahc_eos_ideal,
     .      mahc_eos_handle,mahc_primvar_solver,mahc_prim_energyadd)

      write(*,*) 'test_primvar: done with get_prim_var_pt.'
      write(*,*) '   primative variables: '
      write(*,*) '   rho   : ',rho_send
      write(*,*) '   velx  : ',velx_send
      write(*,*) '   vely  : ',vely_send
      write(*,*) '   velz  : ',velz_send
      write(*,*) '   eps   : ',eps_send
      write(*,*) '   press : ',press_send
      write(*,*) '   w_lor : ',w_lorentz_send

      STOP




      return
      end