###############################################################
#Study of the effect of shell element type and mesh density   #
#Parameters used in study:                                    #
#   eltype: the element type                                  #
#   m_density: the mesh density                               #
#      value of 5 indicates mesh density of 5x5               #
#      value of 10 indicates mesh density of 10x10            #
#      value of 20 indicates mesh density of 20x20            #
###############################################################

# create the study
pinch=ParStudy(par=('eltype','m_density','toggle1','toggle2'),
	       verbose=ON, directory=ON)

# define the parameters
pinch.define(DISCRETE, par='eltype', domain=('S4','S3R','S8R'))
pinch.define(DISCRETE, par='m_density', domain=(5,10,20))

# parameters toggle1 and toggle2 are used to switch between triangular
# and quadrilateral elements in the element connectivity definition
pinch.define(DISCRETE, par='toggle1', domain=('','**'))
pinch.define(DISCRETE, par='toggle2', domain=('','**'))

# sample parameters for design set 1:
#   sample S4 and S8R element types
#   sample 5X5 and 10X10 meshes
#   sample toggle1 as ** (i.e. comment out triangular elements)
#   sample toggle2 as    (i.e. include quadrilateral elements)
pinch.sample(INTERVAL, par='eltype', interval=2)
pinch.sample(VALUES, par='m_density', values=(5,10))
pinch.sample(REFERENCE, par='toggle1', reference=1,
	     interval=1, numSymPairs=0)
pinch.sample(REFERENCE, par='toggle2', reference=0, 
	     interval=1, numSymPairs=0)

# combine samples into 4 designs:
#   5X5 and 10X10 meshes of S4 elements
#   5X5 and 10X10 meshes of S8R elements   
pinch.combine(MESH, name='dSet1')

# sample parameters for design set 2:
#   sample S4 element type
#   sample 20X20 mesh
#   sample toggle1 as ** (i.e. comment out triangular elements)
#   sample toggle2 as    (i.e. include quadrilateral elements)
pinch.sample(REFERENCE, par='eltype', reference=0,
	     interval=1, numSymPairs=0)
pinch.sample(INTERVAL, par='m_density',interval=-3)
pinch.sample(REFERENCE, par='toggle1', reference=1,
	     interval=1, numSymPairs=0)
pinch.sample(REFERENCE, par='toggle2', reference=0,
	     interval=1,numSymPairs=0)

# combine samples into 1 design:
#   20X20 mesh of S4 elements
pinch.combine(TUPLE, name='dSet2')

# sample parameters for design set 3:
#   sample S3R element type
#   sample 5X5, 10X10, and 20X20 meshes
#   sample toggle1 as    (i.e. include triangular elements)
#   sample toggle2 as ** (i.e. comment out quadrilateral elements)
pinch.sample(REFERENCE, par='eltype', reference=1,
	     interval=1, numSymPairs=0)
pinch.sample(NUMBER, par='m_density', number=3)
pinch.sample(REFERENCE, par='toggle1', reference=0,
	     interval=1, numSymPairs=0)
pinch.sample(REFERENCE, par='toggle2', reference=1,
	     interval=1, numSymPairs=0)

# combine samples into 3 designs:
#   5X5, 10X10, and 20X20 meshes of S3R elements
#
pinch.combine(MESH, name='dSet3')

# generate analysis data
pinch.generate(template='pinchcyl_parametric')

# execute runs sequentially
pinch.execute()

# parametric study output at end of step 1
pinch.output(step=1)

# gather displacement results at node 2001
pinch.gather(results='N2001_U', variable='U', node=2001)


# screen report of U2 displacement results at node 2001
pinch.report(PRINT,par=('eltype','m_density'),
	     results=('N2001_U.2'), truncation=OFF)

# read same set of results from ODB file

# parametric study output at end of step 1
pinch.output(step=1, file=ODB)

# gather displacement results at node 2001
pinch.gather(results='o N2001_U', variable='U', node=2001)


# screen report of U2 displacement results at node 2001
pinch.report(PRINT,par=('eltype','m_density'),
	     results=('o N2001_U.2'), truncation=OFF)

# for QA purposes
pinch.report(FILE, par=(),
	     file='pinchcyl_parametric.psr',
	     results=('N2001_U.2'))

pinch.report(FILE, par=(),
	     file='pinchcyl_parametric.psr',
	     results=('o N2001_U.2'))