type c_diff = Real(final quantity = "diffusion coeffient", final unit = "m2/s", displayUnit = "m2/s");