1 files changed, 92 insertions, 0 deletions
diff --git a/src/macro/DZZDG_declare.h b/src/macro/DZZDG_declare.h
new file mode 100644
index 0000000..9ab559f
--- /dev/null
+++ b/src/macro/DZZDG_declare.h
@@ -0,0 +1,92 @@
+/*@@
+  @header   DZZDG_declare.h
+  @date     Jun 98
+  @author   Gabrielle Allen
+  @desc
+  Declarations for macro to calculate the (first and) second derivatives 
+  of the physical metric with respect to z
+
+  The macro uses @seefile DZDG_guts.h and @seefile DZDG_declare.h
+  @enddesc
+@@*/
+
+#ifndef DZZDG_DECLARE
+#define DZZDG_DECLARE
+
+#ifdef FCODE
+
+/* Output variables */
+#undef  DZZDG_DZZDGXX
+#define DZZDG_DZZDGXX dzzdgxx
+#undef  DZZDG_DZZDGXY
+#define DZZDG_DZZDGXY dzzdgxy
+#undef  DZZDG_DZZDGXZ
+#define DZZDG_DZZDGXZ dzzdgxz
+#undef  DZZDG_DZZDGYY
+#define DZZDG_DZZDGYY dzzdgyy
+#undef  DZZDG_DZZDGYZ
+#define DZZDG_DZZDGYZ dzzdgyz
+#undef  DZZDG_DZZDGZZ
+#define DZZDG_DZZDGZZ dzzdgzz
+
+/* Internal variables */
+#undef  DZZDG_FAC
+#define DZZDG_FAC dzzdg_fac
+#undef  DZZDG_OODZ2
+#define DZZDG_OODZ2 dzzdg_oodz2
+#undef  DZZDG_DZZDPSI_O_PSI
+#define DZZDG_DZZDPSI_O_PSI psizz(i,j,k)
+
+/* Declare internal variables */
+      REAL DZZDG_FAC
+      REAL DZZDG_OODZ2
+
+/* Declare output variables */
+      REAL DZZDG_DZZDGXX
+      REAL DZZDG_DZZDGXY
+      REAL DZZDG_DZZDGXZ
+      REAL DZZDG_DZZDGYY
+      REAL DZZDG_DZZDGYZ
+      REAL DZZDG_DZZDGZZ
+
+#endif
+
+#ifdef CCODE
+
+/* Output variables */
+#undef  DZZDG_DZZDGXX
+#define DZZDG_DZZDGXX deldelg3311
+#undef  DZZDG_DZZDGXY
+#define DZZDG_DZZDGXY deldelg3312
+#undef  DZZDG_DZZDGXZ
+#define DZZDG_DZZDGXZ deldelg3313
+#undef  DZZDG_DZZDGYY
+#define DZZDG_DZZDGYY deldelg3322
+#undef  DZZDG_DZZDGYZ
+#define DZZDG_DZZDGYZ deldelg3323
+#undef  DZZDG_DZZDGZZ
+#define DZZDG_DZZDGZZ deldelg3333
+
+/* Internal variables */
+#undef  DZZDG_FAC
+#define DZZDG_FAC dzzdg_fac
+#undef  DZZDG_OODZ2
+#define DZZDG_OODZ2 dzzdg_oodz2
+#undef  DZZDG_DZZDPSI_O_PSI
+#define DZZDG_DZZDPSI_O_PSI psizz[ijk]
+
+/* Declare internal variables */
+double DZZDG_FAC;
+double DZZDG_OODZ2 = 1/(GH->dz0*GH->dz0);
+
+/* Declare output variables */
+double DZZDG_DZZDGXX;
+double DZZDG_DZZDGXY;
+double DZZDG_DZZDGXZ;
+double DZZDG_DZZDGYY;
+double DZZDG_DZZDGYZ;
+double DZZDG_DZZDGZZ;
+
+#endif
+
+#endif