! Loop over list of filenames and make plot (changing '.nc' to '.png')
repeat/range=1:`numfiles`/name=in (let filename = filenames[i=`in`] ; \ let filepng = {spawn:"echo ($filen) | sed s:.nc:.png:"} ; \
def sym filen basinmasks/`filename` ; \
use ($filen) ; \
def sym filep `filepng` ; \
set win/new ; \
shade temp[k=1] )