These scripts do a good job on showing inputs to model and on showing scatter plots.  Note basemapfigs.py is a significant modification of pism/examples/std-greenland/basemapfigs.py.  Still need script for showing bwat,bwprel outputs as maps.