From 2b40052ac2986c1c6d49e2aff8a7fe6d6c5a617c Mon Sep 17 00:00:00 2001
From: Daniel Shapero <shapero.daniel@gmail.com>
Date: Fri, 19 Jan 2024 20:19:27 -0800
Subject: [PATCH] Add plots

---
 .github/workflows/test.yml | 11 +++++++++--
 make_plot.py               | 13 +++++++++++++
 2 files changed, 22 insertions(+), 2 deletions(-)
 create mode 100644 make_plot.py

diff --git a/.github/workflows/test.yml b/.github/workflows/test.yml
index 09c9c17..92ffe48 100644
--- a/.github/workflows/test.yml
+++ b/.github/workflows/test.yml
@@ -24,8 +24,15 @@ jobs:
           pip install git+https://github.com/icepack/icepack.git
       - name: Check out git repository
         uses: actions/checkout@v3
-      - name: run tests
+      - name: Run tests
         env:
           EARTHDATA_USERNAME: ${{ secrets.EARTHDATA_USERNAME }}
           EARTHDATA_PASSWORD: ${{ secrets.EARTHDATA_PASSWORD }}
-        run: python initialize.py
+        run: |
+          python initialize.py
+          python make_plot.py
+      - name: Upload result
+        uses: actions/upload-artifact@v3
+        with:
+          name: kangerdlugssuaq
+          path: kangerdlugssuaq.png
diff --git a/make_plot.py b/make_plot.py
new file mode 100644
index 0000000..819ea40
--- /dev/null
+++ b/make_plot.py
@@ -0,0 +1,13 @@
+import maptlotlib.pyplot as plt
+import firedrake
+
+filename = "kangerdlugssuaq-initial.h5"
+with firedrake.CheckpointFile(filename, "r") as chk:
+    mesh = chk.load_mesh()
+    q = chk.load_function(mesh, name="log_friction")
+
+fig, axes = plt.subplots()
+axes.set_aspect("equal")
+colors = firedrake.tripcolor(q, vmin=-8, vmax=+8, cmap="RdBu", axes=axes)
+fig.colorbar(colors)
+fig.savefig("kangerdlugssuaq.png", bbox_inches="tight")