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

# define the parameter
aStudy.define(DISCRETE, par='CU1', domain=(-10.,-5.,0.,5.,10.))
aStudy.define(DISCRETE, par='CU2', domain=(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_rb_CRF2vsCU2') 

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

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

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