Added functions for vertical flux, molar density and delta 13 c. Compute for each 12c and 13c the modell values