/*@@
  @header   STRESSENERGY_guts.h
  @date     Nov 98
  @author   Gabrielle Allen
  @desc
  Macro to calculate the hydro quantities 
  scheme
  @enddesc
@@*/

#ifndef STRESSENERGY_GUTS
#define STRESSENERGY_GUTS

/* These geometrical quantities are supplied for CalcTmunu */

#include "DETG_guts.h"
#include "UPPERMET_guts.h"

      uxx = UPPERMET_UXX; uyy = UPPERMET_UYY; uzz = UPPERMET_UZZ
      uxy = UPPERMET_UXY; uxz = UPPERMET_UXZ; uyz = UPPERMET_UYZ

      detg = DETG_DETCG

/*   Zero the stress-energy tensor */

      STRESSENERGY_TTT = 0.0D0; STRESSENERGY_TTX = 0.0D0
      STRESSENERGY_TTY = 0.0D0; STRESSENERGY_TTZ = 0.0D0
      STRESSENERGY_TXX = 0.0D0; STRESSENERGY_TXY = 0.0D0 
      STRESSENERGY_TXZ = 0.0D0; STRESSENERGY_TYY = 0.0D0
      STRESSENERGY_TYZ = 0.0D0; STRESSENERGY_TZZ = 0.0D0

      
/*   Calculate the stress-energy tensor */
#include "CalcTmunu.inc"


/*  Calculate the trace of the stress-energy tensor */
#include "TRT_guts.h"

#endif