There is a problem when use "Initial temperature model.Adiabatic.half-space cooling" with number of composition greater than 1(ascii model)

Dear Shaohui,

The adiabatic initial temperature model has a Function subsection, with the function expression set to zero by default. This is needed because the material properties needed to compute the profile (density, thermal expansivity etc.) depend on the composition, so we need to let the model know what values to use for the composition (i.e. what is the composition along the adiabatic reference profile?).

By default, this function expression is set to zero, i.e.

subsection Function
    set Function constants = 
    set Function expression = 0
    set Variable names = x,t
end

Since you have two compositional fields, you also have to provide two components in your function expression rather than just one, i.e.

subsection Function
    set Function constants = 
    set Function expression = 0; 0
    set Variable names = x,t
end

(if 0 is a reasonable assumption for your average 1d profile.)
You have to add this as a subsection into your subsection Adiabatic.

This is what the error message is trying to tell you when it says The number of components (2) (you should provide) is not equal to the number of expressions (1) (which is what is currently in the input file).

Best, 
Juliane