Add the TA calculation and call the CO2DYN program here. Need to add a switch for pCO2 calculation method in the nml file

[Peisheng-Huang]

libaed2/src/aed2_carbon.F90

Line 299 in 0cdfaad