# create the parametric study
aStudy=ParStudy (par=('CU1','CU2','CUR2'), directory=ON)

# define the parameter
aStudy.define(DISCRETE, par='CU1', domain=(-10.,10.))
aStudy.define(DISCRETE, par='CU2', domain=(0.,2.5,5.))
aStudy.define(DISCRETE, par='CUR2', domain=(0.,0.131,0.262))


# sample the parameter
aStudy.sample(INTERVAL, par=('CU1','CU2','CUR2'), interval=1)

# combine the samples to give the designs
aStudy.combine(MESH)

# generate analysis data
aStudy.generate(template='tpl_CRF1vsCU1') 

# execute all analysis jobs sequentially
aStudy.execute(ALL,execOptions=' double')

#for QA purposes
aStudy.output(file=ODB)

aStudy.gather(results='CU1_ANTIALIASING',step=2,frameValue=LAST,request=HISTORY,variable='CU1_ANTIALIASING',element=1,int=-1);
aStudy.gather(results='CRF1_ANTIALIASING',step=2,frameValue=LAST,request=HISTORY,variable='CRF1_ANTIALIASING',element=1,int=-1);
aStudy.report(FILE,file='tpl_rb_CRF1vsCU1.psr',results=('CU1_ANTIALIASING','CRF1_ANTIALIASING'))