/*
 * full_tpc_paper.java
 */

import com.comsol.model.*;
import com.comsol.model.util.*;

/** Model exported on Jul 26 2014, 22:55 by COMSOL 4.3.2.189. */
public class full_tpc_paper {

  public static void main(String[] args) {
    run();
  }

  public static Model run() {
    Model model = ModelUtil.create("Model");

    model.modelPath("/data/fopiroot/fopiroot_dev/tpc_field_COMSOL");

    model.name("full_tpc_paper2.mph");

    model.param().set("Els", "-341.0 [V/cm]");
    model.param().set("Edrift", "-312.9 [V/cm]");
    model.param().set("dLsG1", "4.7 [mm]");
    model.param().set("Vls", "-160.27 [V]");
    model.param().set("Vd", "-22622.67 [V]");
    model.param().set("Vg", "0.0 [V]");
    model.param().set("shortsi", "0");
    model.param().set("shortso", "0");
    model.param()
         .set("Vstepi", "(Vd-Vls)/(nofStrips-shortsi-0.5*(shortsi==0))");
    model.param()
         .set("Vstepo", "(Vd-Vls)/(nofStrips-shortso-0.5*(shortso==0))");
    model.param().set("Bforce", "0.0 [N]");
    model.param().set("Emod", "36.0 [MPa]");
    model.param().set("AreaMoment", "83333.33333333333 [mm^4]");
    model.param().set("nofStrips", "482");
    model.param().set("Cz", "0.25[mm]+0.001[mm]+dLsG1+(nofStrips*1.5[mm])");
    model.param().set("Clength", "105[mm]");
    model.param().set("CmaxR", "155[mm]");
    model.param()
         .set("CminZ", "(Bforce*Clength^3)/(6*Emod*AreaMoment)*( (( 50[mm]-CmaxR)/Clength)^3-3*(( 50[mm]-CmaxR)/Clength)^2)+Cz");
    model.param().set("shi", "-(Cz-CminZ)");
    model.param().set("sho", "0.0");
    model.param().set("Anode_rot", "0.0");
    model.param().set("Ez_nom", "(Vd-Vg)/(Cz)");

    model.modelNode().create("mod1");

    model.func().create("an1", "Analytic");
    model.func("an1")
         .set("plotargs", new String[][]{{"x", "50[mm]", "155[mm]"}});
    model.func("an1").set("argunit", "m");
    model.func("an1")
         .set("expr", "(Bforce*Clength^3)/(6*Emod*AreaMoment)*( (( x -CmaxR)/Clength)^3-3*(( x -CmaxR)/Clength)^2)+Cz");

    model.geom().create("geom1", 2);
    model.geom("geom1").axisymmetric(true);
    model.geom("geom1").lengthUnit("mm");
    model.geom("geom1").run();


    model.variable().create("var1");
    model.variable("var1").set("Ez", "(Vd-Vg)/(Cz-1[mm])");
    model.variable().create("var2");
    model.variable("var2").model("mod1");
    model.variable("var2").set("Ez2", "(Vd-Vg)/(Cz-1[mm])");
    model.variable("var2").set("Ez3", "Ez_nom");

    model.view().create("view4", "geom1");
    model.view().create("view5", "geom1");
    model.view().create("view6", "geom1");

    model.material().create("mat1");
    model.material("mat1").propertyGroup()
         .create("Enu", "Young's modulus and Poisson's ratio");
    model.material("mat1").propertyGroup()
         .create("linzRes", "Linearized resistivity");
    model.material().create("mat2");
    model.material("mat2").propertyGroup("def").func()
         .create("eta", "Piecewise");
    model.material("mat2").propertyGroup("def").func()
         .create("Cp", "Piecewise");
    model.material("mat2").propertyGroup("def").func()
         .create("rho", "Analytic");
    model.material("mat2").propertyGroup("def").func()
         .create("k", "Piecewise");
    model.material("mat2").propertyGroup("def").func()
         .create("cs", "Analytic");
    model.material("mat2").selection().set(new int[]{1});
    model.material().create("mat3");
    model.material("mat3").propertyGroup("def").func()
         .create("k", "Piecewise");
    model.material("mat3").propertyGroup("def").func()
         .create("C", "Piecewise");
    model.material("mat3").propertyGroup("def").func()
         .create("rho", "Piecewise");
    model.material("mat3").selection().named("sel_Insu");

    model.physics().create("es", "Electrostatics", "geom1");

    model.mesh().create("mesh1", "geom1");

    model.result().table().create("tbl1", "Table");
    model.result().table().create("tbl2", "Table");
    model.result().table().create("tbl3", "Table");

    model.variable("var1").name("Variables 1a");
    model.variable("var2").name("Variables 2a");

    model.view("view1").axis().set("xmin", "-3304.452880859375");
    model.view("view1").axis().set("ymin", "-1463.3902587890625");
    model.view("view1").axis().set("xmax", "3509.452880859375");
    model.view("view1").axis().set("ymax", "2191.34130859375");
    model.view("view4").name("close to strips inner");
    model.view("view4").set("locked", true);
    model.view("view4").axis().set("xmin", "44");
    model.view("view4").axis().set("ymin", "435.0");
    model.view("view4").axis().set("manualgrid", "on");
    model.view("view4").axis().set("yspacing", "0.5");
    model.view("view4").axis().set("xmax", "57");
    model.view("view4").axis().set("ymax", "439");
    model.view("view5").name("close to cathode corner");
    model.view("view5").set("locked", true);
    model.view("view5").axis().set("xmin", "30");
    model.view("view5").axis().set("ymin", "708");
    model.view("view5").axis().set("manualgrid", "on");
    model.view("view5").axis().set("xmax", "90");
    model.view("view5").axis().set("ymax", "735");
    model.view("view6").name("close to last strip");
    model.view("view6").set("locked", true);
    model.view("view6").axis().set("xmin", "30");
    model.view("view6").axis().set("ymin", "-10");
    model.view("view6").axis().set("manualgrid", "on");
    model.view("view6").axis().set("xmax", "90");
    model.view("view6").axis().set("ymax", "17");

    model.material("mat1").name("Copper");
    model.material("mat1").propertyGroup("def")
         .set("relpermeability", new String[]{"1", "0", "0", "0", "1", "0", "0", "0", "1"});
    model.material("mat1").propertyGroup("def")
         .set("electricconductivity", new String[]{"5.998e7[S/m]", "0", "0", "0", "5.998e7[S/m]", "0", "0", "0", "5.998e7[S/m]"});
    model.material("mat1").propertyGroup("def")
         .set("thermalexpansioncoefficient", new String[]{"17e-6[1/K]", "0", "0", "0", "17e-6[1/K]", "0", "0", "0", "17e-6[1/K]"});
    model.material("mat1").propertyGroup("def")
         .set("heatcapacity", "385[J/(kg*K)]");
    model.material("mat1").propertyGroup("def")
         .set("relpermittivity", new String[]{"1", "0", "0", "0", "1", "0", "0", "0", "1"});
    model.material("mat1").propertyGroup("def")
         .set("density", "8700[kg/m^3]");
    model.material("mat1").propertyGroup("def")
         .set("thermalconductivity", new String[]{"400[W/(m*K)]", "0", "0", "0", "400[W/(m*K)]", "0", "0", "0", "400[W/(m*K)]"});
    model.material("mat1").propertyGroup("Enu")
         .set("youngsmodulus", "110e9[Pa]");
    model.material("mat1").propertyGroup("Enu").set("poissonsratio", "0.35");
    model.material("mat1").propertyGroup("linzRes").set("rho0", "");
    model.material("mat1").propertyGroup("linzRes").set("alpha", "");
    model.material("mat1").propertyGroup("linzRes").set("Tref", "");
    model.material("mat1").propertyGroup("linzRes")
         .set("rho0", "1.72e-8[ohm*m]");
    model.material("mat1").propertyGroup("linzRes")
         .set("alpha", "0.0039[1/K]");
    model.material("mat1").propertyGroup("linzRes").set("Tref", "298[K]");
    model.material("mat1").propertyGroup("linzRes").addInput("temperature");
    model.material("mat2").name("Air");
    model.material("mat2").propertyGroup("def").func("eta").set("arg", "T");
    model.material("mat2").propertyGroup("def").func("eta")
         .set("pieces", new String[][]{{"200.0", "1600.0", "-8.38278E-7+8.35717342E-8*T^1-7.69429583E-11*T^2+4.6437266E-14*T^3-1.06585607E-17*T^4"}});
    model.material("mat2").propertyGroup("def").func("Cp").set("arg", "T");
    model.material("mat2").propertyGroup("def").func("Cp")
         .set("pieces", new String[][]{{"200.0", "1600.0", "1047.63657-0.372589265*T^1+9.45304214E-4*T^2-6.02409443E-7*T^3+1.2858961E-10*T^4"}});
    model.material("mat2").propertyGroup("def").func("rho")
         .set("argders", new String[][]{{"pA", "d(pA*0.02897/8.314/T,pA)"}, {"T", "d(pA*0.02897/8.314/T,T)"}});
    model.material("mat2").propertyGroup("def").func("rho")
         .set("plotargs", new String[][]{{"pA", "0", "1"}, {"T", "0", "1"}});
    model.material("mat2").propertyGroup("def").func("rho")
         .set("dermethod", "manual");
    model.material("mat2").propertyGroup("def").func("rho")
         .set("args", new String[]{"pA", "T"});
    model.material("mat2").propertyGroup("def").func("rho")
         .set("expr", "pA*0.02897/8.314/T");
    model.material("mat2").propertyGroup("def").func("k").set("arg", "T");
    model.material("mat2").propertyGroup("def").func("k")
         .set("pieces", new String[][]{{"200.0", "1600.0", "-0.00227583562+1.15480022E-4*T^1-7.90252856E-8*T^2+4.11702505E-11*T^3-7.43864331E-15*T^4"}});
    model.material("mat2").propertyGroup("def").func("cs")
         .set("argders", new String[][]{{"T", "d(sqrt(1.4*287*T),T)"}});
    model.material("mat2").propertyGroup("def").func("cs")
         .set("plotargs", new String[][]{{"T", "0", "1"}});
    model.material("mat2").propertyGroup("def").func("cs")
         .set("dermethod", "manual");
    model.material("mat2").propertyGroup("def").func("cs")
         .set("args", new String[]{"T"});
    model.material("mat2").propertyGroup("def").func("cs")
         .set("expr", "sqrt(1.4*287*T)");
    model.material("mat2").propertyGroup("def")
         .set("relpermeability", new String[]{"1", "0", "0", "0", "1", "0", "0", "0", "1"});
    model.material("mat2").propertyGroup("def")
         .set("relpermittivity", new String[]{"1", "0", "0", "0", "1", "0", "0", "0", "1"});
    model.material("mat2").propertyGroup("def")
         .set("dynamicviscosity", "eta(T[1/K])[Pa*s]");
    model.material("mat2").propertyGroup("def")
         .set("ratioofspecificheat", "1.4");
    model.material("mat2").propertyGroup("def")
         .set("electricconductivity", new String[]{"0[S/m]", "0", "0", "0", "0[S/m]", "0", "0", "0", "0[S/m]"});
    model.material("mat2").propertyGroup("def")
         .set("heatcapacity", "Cp(T[1/K])[J/(kg*K)]");
    model.material("mat2").propertyGroup("def")
         .set("density", "rho(pA[1/Pa],T[1/K])[kg/m^3]");
    model.material("mat2").propertyGroup("def")
         .set("thermalconductivity", new String[]{"k(T[1/K])[W/(m*K)]", "0", "0", "0", "k(T[1/K])[W/(m*K)]", "0", "0", "0", "k(T[1/K])[W/(m*K)]"});
    model.material("mat2").propertyGroup("def")
         .set("soundspeed", "cs(T[1/K])[m/s]");
    model.material("mat2").propertyGroup("def").addInput("temperature");
    model.material("mat2").propertyGroup("def").addInput("pressure");
    model.material("mat3").name("Polyimide tape (Kapton HN) [solid]");
    model.material("mat3").propertyGroup("def").func("k").set("arg", "T");
    model.material("mat3").propertyGroup("def").func("k")
         .set("pieces", new String[][]{{"5.0", "140.0", "-0.001372384+0.005601653*T^1+2.082966E-6*T^2-5.05445E-9*T^3"}, {"140.0", "300.0", "-0.007707532+0.005769136*T^1+5.622796E-7*T^2-4.329984E-10*T^3"}});
    model.material("mat3").propertyGroup("def").func("C").set("arg", "T");
    model.material("mat3").propertyGroup("def").func("C")
         .set("pieces", new String[][]{{"4.0", "30.0", "2.809666-1.394349*T^1+0.2106639*T^2+0.004752016*T^3-3.279998E-4*T^4+4.282249E-6*T^5"}, {"30.0", "300.0", "-86.86946+7.816917*T^1-0.03664788*T^2+9.9128E-5*T^3-1.08441E-7*T^4"}});
    model.material("mat3").propertyGroup("def").func("rho").set("arg", "T");
    model.material("mat3").propertyGroup("def").func("rho")
         .set("pieces", new String[][]{{"293.0", "303.0", "1420.0"}});
    model.material("mat3").propertyGroup("def")
         .set("thermalconductivity", new String[]{"k(T[1/K])[W/(m*K)]", "0", "0", "0", "k(T[1/K])[W/(m*K)]", "0", "0", "0", "k(T[1/K])[W/(m*K)]"});
    model.material("mat3").propertyGroup("def")
         .set("heatcapacity", "C(T[1/K])[J/(kg*K)]");
    model.material("mat3").propertyGroup("def")
         .set("density", "rho(T[1/K])[kg/m^3]");
    model.material("mat3").propertyGroup("def")
         .set("relpermittivity", new String[]{"3.83", "0", "0", "0", "3.83", "0", "0", "0", "3.83"});
    model.material("mat3").propertyGroup("def").addInput("temperature");

    model.mesh("mesh1").run();

    model.result().table("tbl1").comments("Surface Average 1 (V)");
    model.result().table("tbl2").comments("Point Evaluation 1 (es.Ez)");
    model.result().table("tbl3").comments("Line Maximum 1 (es.Ez)");

    model.study().create("std1");
    model.study("std1").feature().create("param", "Parametric");
    model.study("std1").feature().create("stat", "Stationary");
    model.study().create("std2");
    model.study("std2").feature().create("param", "Parametric");
    model.study("std2").feature().create("stat", "Stationary");

    model.sol().create("sol1");
    model.sol("sol1").study("std1");
    model.sol("sol1").attach("std1");
    model.sol("sol1").feature().create("st1", "StudyStep");
    model.sol("sol1").feature().create("v1", "Variables");
    model.sol("sol1").feature().create("s1", "Stationary");
    model.sol("sol1").feature("s1").feature().create("p1", "Parametric");
    model.sol("sol1").feature("s1").feature().create("fc1", "FullyCoupled");
    model.sol("sol1").feature("s1").feature().remove("fcDef");

    model.study("std1").feature("stat").set("initstudyhide", "on");
    model.study("std1").feature("stat").set("initsolhide", "on");
    model.study("std1").feature("stat").set("notstudyhide", "on");
    model.study("std1").feature("stat").set("notsolhide", "on");
    model.study("std2").feature("stat").set("initstudyhide", "on");
    model.study("std2").feature("stat").set("initsolhide", "on");
    model.study("std2").feature("stat").set("notstudyhide", "on");
    model.study("std2").feature("stat").set("notsolhide", "on");

    model.sol().create("sol2");
    model.sol("sol2").study("std2");
    model.sol("sol2").attach("std2");
    model.sol("sol2").feature().create("st1", "StudyStep");
    model.sol("sol2").feature().create("v1", "Variables");
    model.sol("sol2").feature().create("s1", "Stationary");
    model.sol("sol2").feature("s1").feature().create("p1", "Parametric");
    model.sol("sol2").feature("s1").feature().create("fc1", "FullyCoupled");
    model.sol("sol2").feature("s1").feature().remove("fcDef");

    model.study("std1").feature("stat").set("initstudyhide", "on");
    model.study("std1").feature("stat").set("initsolhide", "on");
    model.study("std1").feature("stat").set("notstudyhide", "on");
    model.study("std1").feature("stat").set("notsolhide", "on");
    model.study("std2").feature("stat").set("initstudyhide", "on");
    model.study("std2").feature("stat").set("initsolhide", "on");
    model.study("std2").feature("stat").set("notstudyhide", "on");
    model.study("std2").feature("stat").set("notsolhide", "on");

    model.result().dataset().create("cpt1", "CutPoint2D");
    model.result().dataset().create("cln1", "CutLine2D");
    model.result().dataset().create("cln2", "CutLine2D");
    model.result().dataset().create("cln4", "CutLine2D");
    model.result().dataset().create("cln5", "CutLine2D");
    model.result().dataset().create("cln6", "CutLine2D");
    model.result().dataset().create("cln7", "CutLine2D");
    model.result().dataset().create("cpt2", "CutPoint2D");
    model.result().dataset("cln5").set("data", "dset2");
    model.result().dataset("cln7").set("data", "dset2");
    model.result().numerical().create("av1", "AvSurface");
    model.result().numerical().create("pev1", "EvalPoint");
    model.result().numerical().create("max1", "MaxLine");
    model.result().numerical("av1").selection().set(new int[]{1});
    model.result().numerical("av1").set("probetag", "none");
    model.result().numerical("pev1").set("probetag", "none");
    model.result().numerical("max1").set("probetag", "none");
    model.result().create("pg_pot", "PlotGroup2D");
    model.result().create("pg_er", "PlotGroup2D");
    model.result().create("pg_ez", "PlotGroup2D");
    model.result().create("pg_e", "PlotGroup2D");
    model.result().create("pg4", "PlotGroup1D");
    model.result().create("pg8", "PlotGroup1D");
    model.result().create("pg11", "PlotGroup1D");
    model.result().create("pg9", "PlotGroup1D");
    model.result().create("pg10", "PlotGroup1D");
    model.result().create("pg12", "PlotGroup1D");
    model.result().create("pg_er1", "PlotGroup2D");
    model.result().create("pg_ez1", "PlotGroup2D");
    model.result().create("pg_e1", "PlotGroup2D");
    model.result().create("pg13", "PlotGroup1D");
    model.result("pg_pot").feature().create("surf1", "Surface");
    model.result("pg_pot").feature().create("con1", "Contour");
    model.result("pg_er").feature().create("surf1", "Surface");
    model.result("pg_ez").feature().create("surf1", "Surface");
    model.result("pg_e").feature().create("surf1", "Surface");
    model.result("pg4").feature().create("lngr1", "LineGraph");
    model.result("pg4").feature().create("lngr2", "LineGraph");
    model.result("pg8").feature().create("lngr2", "LineGraph");
    model.result("pg11").feature().create("lngr1", "LineGraph");
    model.result("pg9").feature().create("lngr1", "LineGraph");
    model.result("pg10").feature().create("lngr1", "LineGraph");
    model.result("pg12").feature().create("lngr1", "LineGraph");
    model.result("pg_er1").feature().create("surf1", "Surface");
    model.result("pg_ez1").feature().create("surf1", "Surface");
    model.result("pg_e1").feature().create("surf1", "Surface");
    model.result("pg13").feature().create("lngr1", "LineGraph");
    model.result("pg13").feature().create("lngr2", "LineGraph");
    model.result("pg13").feature().create("lngr3", "LineGraph");
























    model.study("std1").name("Vls sweep");
    model.study("std1").feature("param").set("probesel", "none");
    model.study("std1").feature("param").set("pname", new String[]{"Vls"});
    model.study("std1").feature("param")
         .set("plistarr", new String[]{"range(-150,-5,-170)"});
    model.study("std1").feature("param").set("plot", "on");
    model.study("std1").feature("param").set("sweeptype", "filled");
    model.study("std2").name("Real Pot with Shorts");
    model.study("std2").feature("param")
         .set("pname", new String[]{"shortsi", "shortso", "Vg", "Vls", "Vd"});
    model.study("std2").feature("param")
         .set("plistarr", new String[]{"4", "3", "-3458.6", "-3528.6", "-25940"});

    model.sol("sol1").attach("std1");
    model.sol("sol1").feature("st1").name("Compile Equations: Stationary");
    model.sol("sol1").feature("st1").set("studystep", "stat");
    model.sol("sol1").feature("v1").set("control", "stat");
    model.sol("sol1").feature("s1").set("control", "stat");
    model.sol("sol1").feature("s1").feature("p1").set("control", "param");
    model.sol("sol1").feature("s1").feature("p1")
         .set("plistarr", new String[]{"range(-150,-5,-170)"});
    model.sol("sol1").feature("s1").feature("p1").set("sweeptype", "filled");
    model.sol("sol1").feature("s1").feature("p1").set("plot", "on");
    model.sol("sol1").feature("s1").feature("p1").set("probesel", "none");
    model.sol("sol1").feature("s1").feature("p1")
         .set("pname", new String[]{"Vls"});
    model.sol("sol1").runAll();
    model.sol("sol2").attach("std2");
    model.sol("sol2").feature("st1").name("Compile Equations: Stationary");
    model.sol("sol2").feature("st1").set("study", "std2");
    model.sol("sol2").feature("st1").set("studystep", "stat");
    model.sol("sol2").feature("v1").set("control", "stat");
    model.sol("sol2").feature("s1").set("control", "stat");
    model.sol("sol2").feature("s1").feature("p1").set("control", "param");
    model.sol("sol2").feature("s1").feature("p1")
         .set("plistarr", new String[]{"4", "3", "-3458.6", "-3528.6", "-25940"});
    model.sol("sol2").feature("s1").feature("p1")
         .set("pname", new String[]{"shortsi", "shortso", "Vg", "Vls", "Vd"});
    model.sol("sol2").runAll();

    model.result().dataset("cpt1").set("method", "grid");
    model.result().dataset("cpt1").set("gridx", "range(50,0.5,155)");
    model.result().dataset("cpt1")
         .set("gridy", "range(0,0.5,dLsG1+nofStrips*1.5)");
    model.result().dataset("cln1").name("left");
    model.result().dataset("cln1")
         .set("genpoints", new String[][]{{"51", "0"}, {"51", "1000"}});
    model.result().dataset("cln1").set("genparaactive", "on");
    model.result().dataset("cln1").set("genparadist", "range(0,-0.05,-50)");
    model.result().dataset("cln2").name("right");
    model.result().dataset("cln2")
         .set("genpoints", new String[][]{{"155+12.5e-3", "0"}, {"155+12.5e-3", "1000"}});
    model.result().dataset("cln4").name("Last Strip");
    model.result().dataset("cln4")
         .set("genpoints", new String[][]{{"49", "dLsG1+0.5"}, {"156", "dLsG1+0.5"}});
    model.result().dataset("cln5").name("Last Strip 1");
    model.result().dataset("cln5").set("spacevars", new String[]{"cln4x"});
    model.result().dataset("cln5")
         .set("genpoints", new String[][]{{"49", "dLsG1+0.5"}, {"156", "dLsG1+0.5"}});
    model.result().dataset("cln6").name("middle");
    model.result().dataset("cln6")
         .set("genpoints", new String[][]{{"49", "350"}, {"156", "350"}});
    model.result().dataset("cln7").name("middle 1");
    model.result().dataset("cln7").set("spacevars", new String[]{"cln6x"});
    model.result().dataset("cln7")
         .set("genpoints", new String[][]{{"49", "350"}, {"156", "350"}});
    model.result().dataset("cpt2").set("pointx", "105");
    model.result().dataset("cpt2").set("pointy", "350");
    model.result().numerical("av1").set("unit", "V/cm");
    model.result().numerical("av1").set("table", "tbl1");
    model.result().numerical("av1")
         .set("descr", "Electric field, z component");
    model.result().numerical("av1").set("expr", "es.Ez");
    model.result().numerical("pev1").set("data", "cpt2");
    model.result().numerical("pev1").set("expr", "es.Ez");
    model.result().numerical("pev1")
         .set("descr", "Electric field, z component");
    model.result().numerical("pev1").set("unit", "V/cm");
    model.result().numerical("pev1").set("table", "tbl2");
    model.result().numerical("max1").set("data", "cln1");
    model.result().numerical("max1")
         .set("descr", "Electric field, z component");
    model.result().numerical("max1").set("unit", "V/m");
    model.result().numerical("max1").set("table", "tbl3");
    model.result().numerical("max1").set("expr", "es.Ez");
    model.result().numerical("max1").set("solnum", new String[]{"4"});
    model.result().numerical("max1").set("innerinput", "manual");
    model.result().numerical("av1").setResult();
    model.result().numerical("pev1").setResult();
    model.result().numerical("max1").setResult();
    model.result("pg_pot").name("Potential");
    model.result("pg_pot").set("solnum", "4");
    model.result("pg_pot").feature("con1").set("number", "80");
    model.result("pg_pot").feature("con1").set("color", "black");
    model.result("pg_pot").feature("con1").set("coloring", "uniform");
    model.result("pg_pot").feature("con1").set("colorlegend", false);
    model.result("pg_er").name("Radial Field");
    model.result("pg_er").set("xlabelactive", true);
    model.result("pg_er").set("window", "window2");
    model.result("pg_er").set("ylabelactive", true);
    model.result("pg_er").set("titletype", "custom");
    model.result("pg_er").set("solutionintitle", false);
    model.result("pg_er").set("typeintitle", false);
    model.result("pg_er").set("ylabel", "Z (mm)");
    model.result("pg_er").set("windowtitleactive", true);
    model.result("pg_er").set("xlabel", "Radius (mm)");
    model.result("pg_er").set("solnum", "4");
    model.result("pg_er").set("windowtitle", "Radial Field");
    model.result("pg_er").feature("surf1")
         .set("descr", "Electric field, r component");
    model.result("pg_er").feature("surf1").set("rangecolormax", "20");
    model.result("pg_er").feature("surf1").set("expr", "es.Er");
    model.result("pg_er").feature("surf1").set("rangecoloractive", "on");
    model.result("pg_er").feature("surf1").set("resolution", "fine");
    model.result("pg_er").feature("surf1").set("rangecolormin", "-20");
    model.result("pg_er").feature("surf1").set("unit", "V/cm");
    model.result("pg_ez").name("Drift Field");
    model.result("pg_ez").set("xlabelactive", true);
    model.result("pg_ez").set("window", "window3");
    model.result("pg_ez").set("ylabelactive", true);
    model.result("pg_ez").set("titletype", "custom");
    model.result("pg_ez").set("solutionintitle", false);
    model.result("pg_ez").set("typeintitle", false);
    model.result("pg_ez").set("ylabel", "Z (mm)");
    model.result("pg_ez").set("windowtitleactive", true);
    model.result("pg_ez").set("xlabel", "Radius (mm)");
    model.result("pg_ez").set("solnum", "4");
    model.result("pg_ez").set("windowtitle", "Drift Field");
    model.result("pg_ez").feature("surf1")
         .set("descr", "Electric field, z component");
    model.result("pg_ez").feature("surf1").set("rangecolormax", "311.01");
    model.result("pg_ez").feature("surf1").set("expr", "es.Ez");
    model.result("pg_ez").feature("surf1").set("rangecoloractive", "on");
    model.result("pg_ez").feature("surf1").set("resolution", "fine");
    model.result("pg_ez").feature("surf1").set("rangecolormin", "311");
    model.result("pg_ez").feature("surf1").set("unit", "V/cm");
    model.result("pg_e").name("Abs Field");
    model.result("pg_e").set("xlabelactive", true);
    model.result("pg_e").set("window", "window4");
    model.result("pg_e").set("ylabelactive", true);
    model.result("pg_e").set("titletype", "custom");
    model.result("pg_e").set("solutionintitle", false);
    model.result("pg_e").set("typeintitle", false);
    model.result("pg_e").set("ylabel", "Z (mm)");
    model.result("pg_e").set("windowtitleactive", true);
    model.result("pg_e").set("xlabel", "Radius (mm)");
    model.result("pg_e").set("solnum", "4");
    model.result("pg_e").set("windowtitle", "Absolute Field");
    model.result("pg_e").feature("surf1").set("descr", "Electric field norm");
    model.result("pg_e").feature("surf1").set("rangecolormax", "400");
    model.result("pg_e").feature("surf1").set("expr", "es.normE");
    model.result("pg_e").feature("surf1").set("rangecoloractive", "on");
    model.result("pg_e").feature("surf1").set("resolution", "fine");
    model.result("pg_e").feature("surf1").set("rangecolormin", "200");
    model.result("pg_e").feature("surf1").set("unit", "V/cm");
    model.result("pg4").name("Potential along strips");
    model.result("pg4").set("ylabel", "Electric potential (V)");
    model.result("pg4").set("xlabel", "Arc length");
    model.result("pg4").set("ylabelactive", false);
    model.result("pg4").set("xlabelactive", false);
    model.result("pg4").feature("lngr1").set("data", "cln1");
    model.result("pg4").feature("lngr1").set("solrepresentation", "solnum");
    model.result("pg4").feature("lngr1").set("legend", true);
    model.result("pg4").feature("lngr2").set("data", "cln2");
    model.result("pg4").feature("lngr2").set("solrepresentation", "solnum");
    model.result("pg4").feature("lngr2").set("legend", true);
    model.result("pg8").name("Ez along z");
    model.result("pg8").set("solnum", new String[]{"4"});
    model.result("pg8").set("ylabel", "dE<SUB>z</SUB>/E<SUB>z</SUB>");
    model.result("pg8").set("titletype", "custom");
    model.result("pg8").set("windowtitleactive", true);
    model.result("pg8").set("ylabelactive", true);
    model.result("pg8").set("typeintitle", false);
    model.result("pg8").set("xlabel", "Z (cm)");
    model.result("pg8").set("windowtitle", "Ez along Z");
    model.result("pg8").set("window", "window8");
    model.result("pg8").set("xlabelactive", true);
    model.result("pg8").set("innerinput", "manual");
    model.result("pg8").feature("lngr2").name("Ez");
    model.result("pg8").feature("lngr2").set("data", "cln1");
    model.result("pg8").feature("lngr2").set("titletype", "custom");
    model.result("pg8").feature("lngr2").set("solrepresentation", "solnum");
    model.result("pg8").feature("lngr2").set("unit", "1");
    model.result("pg8").feature("lngr2").set("typeintitle", false);
    model.result("pg8").feature("lngr2")
         .set("expr", "abs(  (es.Ez- (Vg-Vd)/(Cz) ) / ( (Vg-Vd)/(Cz) )  )");
    model.result("pg8").feature("lngr2")
         .set("descr", "Relative deviation from nominal drift field");
    model.result("pg8").feature("lngr2").set("unitintitle", false);
    model.result("pg8").feature("lngr2").set("solnum", new String[]{"4"});
    model.result("pg8").feature("lngr2").set("innerinput", "manual");
    model.result("pg8").feature("lngr2").set("descractive", true);
    model.result("pg8").feature("lngr2").set("legend", true);
    model.result("pg11").name("Er along z");
    model.result("pg11").feature("lngr1").name("Er");
    model.result("pg11").feature("lngr1").set("data", "cln1");
    model.result("pg11").feature("lngr1").set("solrepresentation", "solnum");
    model.result("pg11").feature("lngr1").set("unit", "V/m");
    model.result("pg11").feature("lngr1").set("expr", "es.Er");
    model.result("pg11").feature("lngr1")
         .set("descr", "Electric field, r component");
    model.result("pg11").feature("lngr1").set("solnum", new String[]{"4"});
    model.result("pg11").feature("lngr1").set("innerinput", "manual");
    model.result("pg11").feature("lngr1").set("legend", true);
    model.result("pg9").name("Er along last strip");
    model.result("pg9").set("data", "cln4");
    model.result("pg9").set("ylabel", "Electric field, r component (V/cm)");
    model.result("pg9").set("windowtitleactive", true);
    model.result("pg9").set("xlabel", "Arc length");
    model.result("pg9").set("windowtitle", "Er along ls");
    model.result("pg9").set("window", "window7");
    model.result("pg9").set("ylabelactive", false);
    model.result("pg9").set("xlabelactive", false);
    model.result("pg9").feature("lngr1").set("unit", "V/cm");
    model.result("pg9").feature("lngr1").set("expr", "es.Er");
    model.result("pg9").feature("lngr1")
         .set("descr", "Electric field, r component");
    model.result("pg9").feature("lngr1").set("resolution", "fine");
    model.result("pg9").feature("lngr1").set("legend", true);
    model.result("pg10").name("Ez along last strip");
    model.result("pg10").set("ylabel", "Electric field, z component (V/cm)");
    model.result("pg10").set("windowtitleactive", true);
    model.result("pg10").set("xlabel", "Arc length");
    model.result("pg10").set("windowtitle", "Ez along ls");
    model.result("pg10").set("window", "window8");
    model.result("pg10").set("ylabelactive", false);
    model.result("pg10").set("xlabelactive", false);
    model.result("pg10").feature("lngr1").set("data", "cln4");
    model.result("pg10").feature("lngr1").set("solrepresentation", "solnum");
    model.result("pg10").feature("lngr1").set("unit", "V/cm");
    model.result("pg10").feature("lngr1").set("expr", "es.Ez");
    model.result("pg10").feature("lngr1")
         .set("descr", "Electric field, z component");
    model.result("pg10").feature("lngr1").set("legend", true);
    model.result("pg12").name("Er alonng middle");
    model.result("pg12").set("data", "cln6");
    model.result("pg12").set("ylabel", "Electric field, r component (V/cm)");
    model.result("pg12").set("xlabel", "Arc length");
    model.result("pg12").set("ylog", true);
    model.result("pg12").set("ylabelactive", false);
    model.result("pg12").set("xlabelactive", false);
    model.result("pg12").feature("lngr1").set("unit", "V/cm");
    model.result("pg12").feature("lngr1").set("expr", "es.Er");
    model.result("pg12").feature("lngr1")
         .set("descr", "Electric field, r component");
    model.result("pg12").feature("lngr1").set("resolution", "extrafine");
    model.result("pg_er1").name("Radial Field 1");
    model.result("pg_er1").set("data", "dset2");
    model.result("pg_er1").set("xlabelactive", true);
    model.result("pg_er1").set("window", "window2");
    model.result("pg_er1").set("ylabelactive", true);
    model.result("pg_er1").set("titletype", "custom");
    model.result("pg_er1").set("solrepresentation", "solutioninfo");
    model.result("pg_er1").set("solutionintitle", false);
    model.result("pg_er1").set("typeintitle", false);
    model.result("pg_er1").set("ylabel", "Z (mm)");
    model.result("pg_er1").set("windowtitleactive", true);
    model.result("pg_er1").set("looplevel", new String[]{"1"});
    model.result("pg_er1").set("xlabel", "Radius (mm)");
    model.result("pg_er1").set("windowtitle", "Radial Field");
    model.result("pg_er1").feature("surf1")
         .set("descr", "Electric field, r component");
    model.result("pg_er1").feature("surf1").set("rangecolormax", "20");
    model.result("pg_er1").feature("surf1").set("expr", "es.Er");
    model.result("pg_er1").feature("surf1").set("rangecoloractive", "on");
    model.result("pg_er1").feature("surf1").set("resolution", "fine");
    model.result("pg_er1").feature("surf1").set("rangecolormin", "-20");
    model.result("pg_er1").feature("surf1").set("unit", "V/cm");
    model.result("pg_ez1").name("Drift Field 1");
    model.result("pg_ez1").set("data", "dset2");
    model.result("pg_ez1").set("xlabelactive", true);
    model.result("pg_ez1").set("window", "window3");
    model.result("pg_ez1").set("ylabelactive", true);
    model.result("pg_ez1").set("titletype", "custom");
    model.result("pg_ez1").set("solrepresentation", "solutioninfo");
    model.result("pg_ez1").set("solutionintitle", false);
    model.result("pg_ez1").set("typeintitle", false);
    model.result("pg_ez1").set("ylabel", "Z (mm)");
    model.result("pg_ez1").set("windowtitleactive", true);
    model.result("pg_ez1").set("looplevel", new String[]{"1"});
    model.result("pg_ez1").set("xlabel", "Radius (mm)");
    model.result("pg_ez1").set("windowtitle", "Drift Field");
    model.result("pg_ez1").feature("surf1")
         .set("descr", "Electric field, z component");
    model.result("pg_ez1").feature("surf1").set("rangecolormax", "400");
    model.result("pg_ez1").feature("surf1").set("expr", "es.Ez");
    model.result("pg_ez1").feature("surf1").set("rangecoloractive", "on");
    model.result("pg_ez1").feature("surf1").set("resolution", "fine");
    model.result("pg_ez1").feature("surf1").set("rangecolormin", "200");
    model.result("pg_ez1").feature("surf1").set("unit", "V/cm");
    model.result("pg_e1").name("Abs Field 1");
    model.result("pg_e1").set("data", "dset2");
    model.result("pg_e1").set("xlabelactive", true);
    model.result("pg_e1").set("window", "window4");
    model.result("pg_e1").set("ylabelactive", true);
    model.result("pg_e1").set("titletype", "custom");
    model.result("pg_e1").set("solrepresentation", "solutioninfo");
    model.result("pg_e1").set("solutionintitle", false);
    model.result("pg_e1").set("typeintitle", false);
    model.result("pg_e1").set("ylabel", "Z (mm)");
    model.result("pg_e1").set("windowtitleactive", true);
    model.result("pg_e1").set("looplevel", new String[]{"1"});
    model.result("pg_e1").set("xlabel", "Radius (mm)");
    model.result("pg_e1").set("windowtitle", "Absolute Field");
    model.result("pg_e1").feature("surf1")
         .set("descr", "Electric field norm");
    model.result("pg_e1").feature("surf1").set("rangecolormax", "400");
    model.result("pg_e1").feature("surf1").set("expr", "es.normE");
    model.result("pg_e1").feature("surf1").set("rangecoloractive", "on");
    model.result("pg_e1").feature("surf1").set("resolution", "fine");
    model.result("pg_e1").feature("surf1").set("rangecolormin", "200");
    model.result("pg_e1").feature("surf1").set("unit", "V/cm");
    model.result("pg13").name("Ez diff along middle");
    model.result("pg13").set("data", "cln6");
    model.result("pg13").set("solnum", new String[]{"3", "4", "5", "6"});
    model.result("pg13")
         .set("ylabel", "abs( (es.Ez - (Vg-Vd)/Cz  )   /  ( ( Vg-Vd)/Cz )  ) (1)");
    model.result("pg13").set("xlabel", "Reversed arc length");
    model.result("pg13").set("innerinput", "manual");
    model.result("pg13").set("ylog", true);
    model.result("pg13").set("ylabelactive", false);
    model.result("pg13").set("xlabelactive", false);
    model.result("pg13").feature("lngr1").set("unit", "1");
    model.result("pg13").feature("lngr1").set("xdata", "reversedarc");
    model.result("pg13").feature("lngr1")
         .set("expr", "abs( (es.Ez - (Vg-Vd)/Cz  )   /  ( ( Vg-Vd)/Cz )  )");
    model.result("pg13").feature("lngr1")
         .set("descr", "abs( (es.Ez - (Vg-Vd)/Cz  )   /  ( ( Vg-Vd)/Cz )  )");
    model.result("pg13").feature("lngr1").set("legend", true);
    model.result("pg13").feature("lngr2").active(false);
    model.result("pg13").feature("lngr2").set("unit", "V/m");
    model.result("pg13").feature("lngr2").set("expr", "(Vg-Vd)/(Cz)");
    model.result("pg13").feature("lngr2").set("descr", "(Vg-Vd)/(Cz)");
    model.result("pg13").feature("lngr3").active(false);
    model.result("pg13").feature("lngr3").set("unit", "cV/m");
    model.result("pg13").feature("lngr3").set("expr", "es.Ez");
    model.result("pg13").feature("lngr3")
         .set("descr", "Electric field, z component");
    model.result("pg13").feature("lngr3").set("legend", true);

    model.name("full_tpc_paper2.mph");

    return model;
  }

}