Add plots

.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

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")