def plotFinal(solPlot,aSolPlot,tn,test,tpars,problems,
              mlMesh,mlScalarTransport,testOut):
    """
    plot out solution and mesh at last step in a couple of formats
    """
    if tpars['DG'] == False:
        solPlot.hardcopy(testOut+'_sol.eps', eps=1,enhanced=1,color=1)
        aSolPlot.hardcopy(testOut+'_asol.eps', eps=1,enhanced=1,color=1)
    #end if
    mlMesh.meshList[-1].writeMeshEnsight(test,test)
    mlScalarTransport.modelList[-1].u.name='u'
    #mlScalarTransport.modelList[-1].writeBoundaryTermsEnsight(test)
    mlScalarTransport.modelList[-1].u.writeFunctionEnsight(test,append=False)

    return solPlot,aSolPlot