Initial conditions

As specified in initialization.c, the initial conditions are straightforward in that all quantities are initialized to zero except the temperature field, which is set to 1 everywhere. It is treated as a passive tracer when $\gamma=1$ (suntans.dat: gamma = 1).