calc_c_closure function

Calculate canopy closure date