! source file: /raid23/csomes/UVic/2.9/c_n_isotopes/last_glacial_experiments/dlgfe45p15/updates/sed.h ! sed include file ! ipmax = maximum number of sediment points ! nzmax = maximum number of mixed layers ! ibmax = maximum number of buried layers integer ipmax, nzmax, ibmax, imap, jmap, map_sed, ipsed, kmax integer kmin parameter (ipmax=imt*jmt, nzmax=8, ibmax=20) ! imap = 1d maping index for 2d map (i index) ! jmap = 1d maping index for 2d map (j index) ! map_sed = 2d map of sediment indices ! ipsed = number of sediment points (must be >= to ipmax) ! kmax = number of mixed layers (must be less than nzmax) ! kmin = minimum kmt level for sediments common /sed_i/ imap(ipmax), jmap(ipmax), map_sed(imt,jmt) common /sed_i/ ipsed, kmax, kmin ! mixed layer fields ! carb = forms of carbon, 1=co2,2=hco3,3=co3 (mol l-1) ! dcpls = diffusion coefficient above, 1=co2,2=hco3,3=co3 ! dcmin = diffusion coefficient below, 1=co2,2=hco3,3=co3 ! pore = porosity ! form = pore**3. ! o2 = bottom water oxygen (mol l-1) ! orggg = organic carbon fraction ! orgml = organic carbon mass (g cm-2) ! calgg = calcite fraction ! calml = calcite mass (g CaCO3 cm-2) ! dopls = diffusion coefficient for oxygen above ! domin = diffusion coefficient for oxygen below ! dbpls = diffusion coefficient for organic carbon above ! dbmin = diffusion coefficient for organic carbon below ! buried layer fields ! buried_mass = buried mass (g CaCO3 cm-2) ! buried_calfrac = buried calcite fraction ! depth_age = age of buried level (years) ! sediment surface or integrated fields ! zrct = maximum respiration depth (cm) ! water_z_p = ocean depth (m) ! k1 = first dissociation constant for carbonic acid ! k2 = second dissociation constant for carbonic acid ! k3 = dissociation constant for boric acid ! csat = carbonate saturation at sediment surface (mol l-1) ! rc = respiration coefficient ! sed_ml_mass = mixed layer mass (g cm-2) ! ttrorg = dissolution rate organic carbon (mol cm-2 dtsed-1) ! ttrcal = dissolution rate of calcite (mol cm-2 dtsed-1) ! c_advect = burial flux [mol cm-2 dtsed-1] ! zsed = mixed layer grid depth (cm) ! delz = mixed layer grid layer thickness (cm) ! rain_org_p = rain rate of organic carbon (mol cm-2 dtsed-1) ! rain_cal_p = rain rate of calcite (mol cm-2 dtsed-1) ! co3_p = carbonate at sediment surface (mol l-1) ! scalars ! dissc = calcite dissolution constant ! dissn = calcite dissolution constant ! weath = weathering flux (kg s-1) ! weathflx = weathering flux (umol s-1) ! sed_year = year for sediment profile ! sedsa = surface area of potential deep sediments (cm2) ! carblith = change in lithosphere carbon (umol) ! set to zero if time is initialized real carb, dcpls, dcmin, pore, form, o2, orggg, orgml, calgg real calml, dopls, domin, dbpls, dbmin, buried_mass real buried_calfrac, depth_age, zrct, water_z_p, k1, k2, k3, csat real rc, sed_ml_mass, ttrorg, ttrcal, c_advect, zsed, delz real rain_org_p, rain_cal_p, co3_p, dissc, dissn, weath, weathflx real dicwflx, sed_year, sedsa, carblith common /sed_r/ carb(nzmax,3,ipmax), dcpls(nzmax,3,ipmax) common /sed_r/ dcmin(nzmax,3,ipmax) common /sed_r/ pore(nzmax,ipmax), form(nzmax,ipmax) common /sed_r/ o2(nzmax,ipmax), orggg(nzmax,ipmax) common /sed_r/ orgml(nzmax,ipmax), calgg(nzmax,ipmax) common /sed_r/ calml(nzmax,ipmax), dopls(nzmax,ipmax) common /sed_r/ domin(nzmax,ipmax), dbpls(nzmax,ipmax) common /sed_r/ dbmin(nzmax,ipmax), buried_mass(ibmax,ipmax) common /sed_r/ buried_calfrac(ibmax,ipmax), depth_age(ibmax) common /sed_r/ zrct(ipmax), water_z_p(ipmax), k1(ipmax), k2(ipmax) common /sed_r/ k3(ipmax), csat(ipmax), rc(ipmax), ttrorg(ipmax) common /sed_r/ ttrcal(ipmax), sed_ml_mass(ipmax), c_advect(ipmax) common /sed_r/ zsed(nzmax), delz(nzmax), rain_org_p(ipmax) common /sed_r/ rain_cal_p(ipmax), co3_p(ipmax) common /sed_r/ dissc, dissn, weath, weathflx, dicwflx, sed_year common /sed_r/ sedsa, carblith ! ntatss = time step counter for time averaging ! ta_ttrcal = time average dissolution rate of calcite ! ta_rain_cal = time average rain rate of calcite ! ta_cal = time average mixed layer calcite fraction ! ta_calmass = time average mixed layer calcite mass ! ta_calmass_bur = time average buried calcite mass ! ta_co3 = time average carbonate at sediment surface ! ta_co3sat = time average carbonate saturation at surface ! ta_rainr = time average rain ratio integer ntatss common /sed_i/ ntatss real ta_ttrcal, ta_rain_cal, ta_cal, ta_calmass, ta_calmass_bur real ta_co3, ta_co3sat, ta_rainr common /sed_r/ ta_ttrcal(imt,jmt), ta_rain_cal(imt,jmt) common /sed_r/ ta_cal(imt,jmt), ta_calmass(imt,jmt) common /sed_r/ ta_calmass_bur(imt,jmt), ta_co3(imt,jmt) common /sed_r/ ta_co3sat(imt,jmt), ta_rainr(imt,jmt) ! ntatis = number of time averaged time step integrals ! tai_ttrcal = time step integral of ttrcal ! tai_rain_cal = time step integral of rain_cal ! tai_cal = time step integral of cal ! tai_weathflx = time step integral of weathflx ! tai_calmass = time step integral of mixed layer calcite mass ! tai_calmass_bur = time step integral of buried calcite mass ! tai_co3 = time step integral of co3 ! tai_co3sat = time step integral of co3sat ! tai_rainr = time step integral of rain ratio ! tai_carblith = time step integral of carblith ! tai_cfo2s = average total flux ocean to sediments ! tai_cfl2o = average total flux land to ocean (weathering) integer ntatis common /sed_i/ ntatis real tai_ttrcal, tai_rain_cal, tai_cal, tai_weathflx, tai_calmass real tai_calmass_bur, tai_co3, tai_co3sat, tai_rainr, tai_carblith real tai_cfo2s, tai_cfl2o common /sed_r/ tai_ttrcal, tai_rain_cal, tai_cal, tai_weathflx common /sed_r/ tai_calmass, tai_calmass_bur, tai_co3, tai_co3sat common /sed_r/ tai_rainr, tai_carblith, tai_cfo2s, tai_cfl2o