ipython solution control. Create .geo from template, convert to PAFEC, parallel solve, import.

fn = 'v9_j55'
if solve:
    createAnalysisFile(dataRoot + fn + '_template.geo', dataRoot + fn + '.geo', 
                   prepend = 'dia_t = 1.067E-3;\n')
    output = !gmsh -2 {dataRoot + fn + '.geo'}  # this works!
    args = setArgs()
    args['FREQUENCIES.FOR.ANALYSIS'] = ffa(3, 15000, 50000, 192)
    args['restrain'] = [{'Heater':1}]
    convert.gmsh2pafec(dataRoot + fn, case = 66, verb = 0, args = args)
    solveJob(fn)
if loadData:
    res['v9_j55: heaterClamp'] = getResult(dataRoot + fn, Vnode = 68)