1 files changed, 41 insertions, 0 deletions
diff --git a/src/macro/MOMYADM_guts.h b/src/macro/MOMYADM_guts.h
new file mode 100644
index 0000000..00e246d
--- /dev/null
+++ b/src/macro/MOMYADM_guts.h
@@ -0,0 +1,41 @@
+/*@@
+  @header   MOMYADM_guts.h
+  @date     Aug 98
+  @author   Gabrielle Allen
+  @desc
+  Macro to calculate the spacetime part of the 
+  y-Momentum Constraint. That is:
+
+       Del_j K_y^j - Del_y trK 
+
+  @enddesc
+@@*/
+
+#ifndef MOMYADM_GUTS
+#define MOMYADM_GUTS
+
+#include "../../packages/CactusEinstein/Einstein/src/macro/UPPERMET_guts.h"
+#include "../../packages/CactusEinstein/Einstein/src/macro/CDK_guts.h"
+
+#ifdef FCODE 
+
+      MOMYADM_MOMYADM = (CDXCDK_CDXCDKXY - CDYCDK_CDYCDKXX)*UPPERMET_UXX 
+     &  + (CDXCDK_CDXCDKYY - CDYCDK_CDYCDKXY)*UPPERMET_UXY + 
+     &  (CDXCDK_CDXCDKYZ - 2*CDYCDK_CDYCDKXZ + CDZCDK_CDZCDKXY)*UPPERMET_UXZ + 
+     &  (-CDYCDK_CDYCDKYZ + CDZCDK_CDZCDKYY)*UPPERMET_UYZ + 
+     &  (-CDYCDK_CDYCDKZZ + CDZCDK_CDZCDKYZ)*UPPERMET_UZZ
+
+#endif
+
+#ifdef CCODE
+
+      MOMYADM_MOMYADM = (CDXCDK_CDXCDKXY - CDYCDK_CDYCDKXX)*UPPERMET_UXX 
+       + (CDXCDK_CDXCDKYY - CDYCDK_CDYCDKXY)*UPPERMET_UXY + 
+       (CDXCDK_CDXCDKYZ - 2*CDYCDK_CDYCDKXZ + CDZCDK_CDZCDKXY)*UPPERMET_UXZ + 
+       (-CDYCDK_CDYCDKYZ + CDZCDK_CDZCDKYY)*UPPERMET_UYZ + 
+       (-CDYCDK_CDYCDKZZ + CDZCDK_CDZCDKYZ)*UPPERMET_UZZ;
+
+#endif
+
+#endif
+