1 files changed, 40 insertions, 0 deletions
diff --git a/src/macro/DZDK_guts.h b/src/macro/DZDK_guts.h
new file mode 100644
index 0000000..f388e37
--- /dev/null
+++ b/src/macro/DZDK_guts.h
@@ -0,0 +1,40 @@
+/*@@
+  @header   DZDG_guts.h
+  @date     Jul 98
+  @author   Gabrielle Allen
+  @desc
+  Macro to calculate the first derivatives of the 
+  extrinsic with respect to z
+  @enddesc
+@@*/
+
+#ifndef DZDK_GUTS
+#define DZDK_GUTS
+
+#ifdef FCODE 
+
+      DZDK_OO2DZ = 1D0/(2D0*DZDK_DZ)
+    
+      DZDK_DZDKXX = DZDK_OO2DZ*(DZDK_KXX_KP - DZDK_KXX_KM)
+      DZDK_DZDKXY = DZDK_OO2DZ*(DZDK_KXY_KP - DZDK_KXY_KM)
+      DZDK_DZDKXZ = DZDK_OO2DZ*(DZDK_KXZ_KP - DZDK_KXZ_KM)
+      DZDK_DZDKYY = DZDK_OO2DZ*(DZDK_KYY_KP - DZDK_KYY_KM)
+      DZDK_DZDKYZ = DZDK_OO2DZ*(DZDK_KYZ_KP - DZDK_KYZ_KM)
+      DZDK_DZDKZZ = DZDK_OO2DZ*(DZDK_KZZ_KP - DZDK_KZZ_KM)
+
+#endif
+
+#ifdef CCODE
+
+      DZDK_OO2DZ = 1/(2*GH->dz0);
+    
+      DZDK_DZDKXX = DZDK_OO2DZ*(DZDK_KXX_KP - DZDK_KXX_KM);
+      DZDK_DZDKXY = DZDK_OO2DZ*(DZDK_KXY_KP - DZDK_KXY_KM);
+      DZDK_DZDKXZ = DZDK_OO2DZ*(DZDK_KXZ_KP - DZDK_KXZ_KM);
+      DZDK_DZDKYY = DZDK_OO2DZ*(DZDK_KYY_KP - DZDK_KYY_KM);
+      DZDK_DZDKYZ = DZDK_OO2DZ*(DZDK_KYZ_KP - DZDK_KYZ_KM);
+      DZDK_DZDKZZ = DZDK_OO2DZ*(DZDK_KZZ_KP - DZDK_KZZ_KM);
+
+#endif
+
+#endif