###################################################################
# Mesh convergence study for plate pentration verification        #
# problem in ABAQUS/Explicit                                      # 
# Study parameters:                                               # 
#    contype: Contact type                                        #
###################################################################

# Create the study
cont = ParStudy(par='contype')

# Define the parameter
cont.define(DISCRETE, par='contype',
	    domain=('CONTACT PAIR, PENALTY','CONTACT PAIR'))

# Sample the parameter
cont.sample(INTERVAL, par='contype', interval=1)

# Combine the samples into design points
cont.combine(MESH, name='con')

# Generate analysis job data
cont.generate(template='pp_con_study')

# Execute all analysis jobs sequentially
cont.execute(ALL)

# Study output at end of step 1
cont.output(step=1)

# Gather the results for the velocity of the projectile
cont.gather(results='vel', variable='V', node=99999)

# Report the gathered results
cont.report(PRINT, results=('vel.2'), par=('contype'))

# read same set of results from ODB file

# Study output at end of step 1
cont.output(step=1, file=ODB)

# Gather the results for the velocity of the projectile
cont.gather(results='o vel', variable='V2', node=99999, request=HISTORY)

# Report the gathered results
cont.report(PRINT, results=('o vel'), par=('contype'))

# FOR QA PURPOSES
qa = open('pp_con_study.psr','w')
keys = cont.resultData.keys()
keys.sort()
for key in keys:
    for item in cont.resultData[key]:
        qa.write( '%6.0f, ' % item )
    qa.write( '\n' )
qa.close()