diff --git a/experiments/ClimaCore/heat-diffusion/run.jl b/experiments/ClimaCore/heat-diffusion/run.jl
index ce7d1d3cf8..c99985dc0e 100644
--- a/experiments/ClimaCore/heat-diffusion/run.jl
+++ b/experiments/ClimaCore/heat-diffusion/run.jl
@@ -376,5 +376,4 @@ Plots.png(
 #src     end
 #src end
 #src
-#src dirname = "heat"
-#src linkfig("output/$(dirname)/heat_end.png", "Heat End Simulation")
+#src linkfig(ARTIFACTS_DIR * "/heat_end.png", "Heat End Simulation")
diff --git a/output/README.md b/output/README.md
deleted file mode 100644
index b4eb40088b..0000000000
--- a/output/README.md
+++ /dev/null
@@ -1 +0,0 @@
-Run generated output directory.
diff --git a/plotting/README.md b/plotting/README.md
deleted file mode 100644
index 56a1cc6029..0000000000
--- a/plotting/README.md
+++ /dev/null
@@ -1 +0,0 @@
-Quick plotting utilities
diff --git a/plotting/plot-field.py b/plotting/plot-field.py
deleted file mode 100644
index d6d6a62dd5..0000000000
--- a/plotting/plot-field.py
+++ /dev/null
@@ -1,56 +0,0 @@
-#!/usr/bin/env python
-# coding: utf-8
-
-# To set up pyvista use Miniconda and pip e.g.
-#
-#  wget https://repo.anaconda.com/miniconda/Miniconda3-latest-MacOSX-x86_64.sh
-#  chmod +x Miniconda3-latest-MacOSX-x86_64.sh
-#  ./Miniconda3-latest-MacOSX-x86_64.sh -b -p miniconda3
-# source miniconda3/bin/activate
-# conda create --name vtk_env python=3.8
-# conda activate vtk_env
-# conda install nodejs matplotlib
-# pip install jupyter pyvista ipyvtk_simple
-# 
-# can be be run as a script e.g.
-#
-# python plot-field.py
-#
-# or in notebook - 
-#
-# jupyter notebook
-
-import vtk
-import pyvista as pv
-from matplotlib import cm
-import numpy as np
-
-fcode='θ'
-# fcode='state_debug'
-
-mesh22 = pv.read('Coupler_UnitTest_atmosphere_long_num0004.pvtu')
-
-xc=mesh22.get_array('xc')
-yc=mesh22.get_array('yc')
-zc=mesh22.get_array('zc')
-
-
-R=np.sqrt(np.square(xc)+np.square(yc)+np.square(zc))
-
-# Mask to surface and "top" half of sphere.
-tol=1e-1 
-sval=( np.abs(R[0]-R)>tol ) | ( zc>0 )
-
-l=np.array(list(range(len( R )))); l[ np.abs(R[0]-R)<tol ];len(l);len( R[np.abs(R[0]-R)<tol] )
-ptnan=l[sval]
-ptplot=l[~sval]
-
-mesh22[fcode][ptnan]=np.nan
-
-msurf=mesh22.extract_points(ptplot,adjacent_cells=True,include_cells=True)
-
-plotter = pv.Plotter()
-# plotter.add_mesh(msurf,style='points',show_edges=True,scalars='θ',culling=True,cmap='jet',lighting=True,nan_opacity=0,point_size=10)
-plotter.add_mesh(msurf,style='surface',show_edges=False,scalars=fcode,culling=True,cmap='jet',lighting=True,nan_opacity=0)
-# plotter.add_mesh(msurf,style='points',scalars='θ',culling=True,cmap='jet',lighting=True,nan_opacity=0)
-plotter.show()