mdb.saveAs(pathName = mdbName + '.cae')
func.runJob(modelName)
#Write CPU time to file
func.readMsgFile(modelName, 'results.txt')
#=========== Post ============#
print 'Post processing...'
# #Contour
# func.countourPrint(modelName, defScale, printFormat)
#Energy
func.xyEnergyPlot(modelName)
#R2 at col base
beam.xyColBaseR2(modelName,x,z)
#Displacement at colTop
beam.xyAPMcolPrint(modelName, APMcol)
# #Check largest peeq against criteria
# print '\n' + "Getting data from ODB..."
# elmOverLim = func.getElmOverLim(modelName, var,
# stepName, var_invariant, limit)
# print " done"
# with open('results.txt','a') as f:
# if elmOverLim:
func.readMsgFile(modelName, 'results.txt')
#===================================================#
#===================================================#
# Post #
#===================================================#
#===================================================#
print 'Post processing...'
# #Contour
# func.countourPrint(modelName, defScale, printFormat)
#Animation
#func.animate(modelName, defScale, frameRate= animeFrameRate)
#Energy
func.xyEnergyPlot(modelName)
#R2 at column base
shell.xyR2colBase(modelName, x, z)
shell.xyUcolTop(modelName, column=APMcol)
# #Force and displacement
# shell.xyCenterU2_colBaseR2(modelName,x,z)
print ' done'
print '########### END OF SCRIPT ###########'