From 32049a6ef493f7a4409e8e543f1cc8433ab421da Mon Sep 17 00:00:00 2001
From: Ben Dudson <dudson2@llnl.gov>
Date: Fri, 19 May 2023 20:31:48 -0700
Subject: [PATCH] Output metric tensor in outputVars

Updates the Jacobian (J) that is otherwise not saved after modifying in init()
---
 sd1d.cxx | 3 +++
 1 file changed, 3 insertions(+)

diff --git a/sd1d.cxx b/sd1d.cxx
index eb0fd70..58d84ab 100644
--- a/sd1d.cxx
+++ b/sd1d.cxx
@@ -1707,6 +1707,9 @@ class SD1D : public PhysicsModel {
 
     state["SD1D_REVISION"].force(sd1d::version::revision);
 
+    // Ensure that metrics are updated
+    mesh->getCoordinates()->outputVars(state);
+
     // Write normalisation constants
 
     BoutReal Pnorm = SI::qe * Tnorm * Nnorm; // Pressure normalisation