I am trying to use a "parameter" field variable in interactive mode. What I am trying to do is very similar to the thermo-elaticity example, but for some reason I cannot make it work in interactive mode. 

In the thermo-elasticity example the parameter field variable is defined as follows 

def get_temperature_load(ts, coors, region=None):
    """
    Temperature load depends on the `x` coordinate.
    """
    x = coors[:, 0]
    return (x - x.min())**2 - T0

fields = {
    'displacement': ('real', 3, 'Omega', 1),
    'temperature': ('real', 1, 'Omega', 1),
}

variables = {
    'u' : ('unknown field', 'displacement', 0),
    'v' : ('test field', 'displacement', 'u'),
    'T' : ('parameter field', 'temperature',
           {'setter' : 'get_temperature_load'}),
}

equations = {
    'balance_of_forces' :
    """dw_lin_elastic.2.Omega( solid.D, v, u )
     - dw_biot.2.Omega( solid.alpha, v, T )
     = 0""",
}

In interactive mode this is what I am doing the following:

field_p = Field.from_args('fp', nm.float64, 1, omega,

                               approx_order=order_p)


press_load = Function('press_load',get_pressure_load)

press = FieldVariable('pressure','parameter',field_p,special={'setter':press_load},

                          primary_var_name='(set-to-None)')


And the corresponding term:


t20 = Term.new('dw_laplace(matalpha.k, q, press)',

                   integral, omega, matalpha=matalpha, q=q, press=press)


But this does not work.


Any thoughts?


Thanks