Sources

The optional sources parameter is a listdictionary of external sources of carbon, alkalinity, and energy. If an element of sources is a numericfloat, then it is interpreted as a constant rate. If an element of sources is a function, then it is interpreted as the time-integrated forcing as a function of number of elapsed years, i.e.,

in units equal to the units given for the rates below times yr. For example, sources=list(Cas=function(t){if(t<=0){0}else{1}})sources={"Cas":lambda t:0 if t<=0 else 1} means that 1 GtC will be added to the atmosphere and surface ocean between the first and second time steps. By default, sources is NULLNone, which sets all of the external sources to zero. The values that can be set within the sources listdictionary are:

• rad: External top-of-atmosphere radiative forcing in units of W/m². If sources$radsources["rad"] is, instead, a function, then it is interpreted as the time-integrated radiative forcing in units of W/m² x yr and as a function of time in yr.
• Cas, Cdeep, Cveg1, Cveg2, …, CvegN, Csoil1, Csoil2, …, CsoilM: External source of carbon to the atmosphere plus surface ocean, deep ocean, vegetation pools 1 through N (where N is the number specified in picontrol, which is 3 by default), and soil pools 1 through M (where M is the number specified in picontrol, which is 2 by default), respectively, in GtC/yr. If, e.g., sources$Cdeepsources["Cdeep"] is, instead, a function, then it is interpreted as the time-integrated source of carbon to the deep ocean in units of GtC and as a function of time in yr.
• Asurf, Adeep: External source of alkalinity to the surface ocean and deep ocean, respectively, in mol/yr. If, e.g., sources$Adeepsources["Adeep"] is, instead, a function, then it is interpreted as the time-integrated source of alkalinity to the deep ocean in units of mol and as a function of time in yr.