New initialization method is added. It can be off by the flag start_by_grad. Tested with temperate and boreal sites, 100yrs run.