!!--------------------------------------------------------------------- !! !! COMMON passivetrc.hamocc3.h !! ****************************** !! !! purpose : !! --------- !! INCLUDE COMMON FILE for HAMO!!3 biological model !! !! modifications : !! ------------- !! original : 00-02 (O. Aumont) !! !!--------------------------------------------------------------------- !! opa8, ipsl (11/96) !!--------------------------------------------------------------------- !! #if defined key_trc_hamocc3 !! !!---------------------------------------------------------------------- !! !! COMMON/cchem1/ : Variable for chemistry of the CO2 cycle !! !! --------------------------------------------------------------------- !! REAL akb3(jpi,jpj,jpk), ak13(jpi,jpj,jpk), ak23(jpi,jpj,jpk) REAL aksp(jpi,jpj,jpk), co3(jpi,jpj,jpk), hi(jpi,jpj,jpk) REAL borat(jpi,jpj,jpk) COMMON/cchem1/ akb3,ak13,ak23,aksp,co3,hi,borat !! !!---------------------------------------------------------------------- !! !! COMMON/cchem2/ : Variable for chemistry of the CO2 cycle !! !! --------------------------------------------------------------------- !! !! dmin3 : vertical distribution of newly formed biogenic poc !! disc3 : vertical distribution of newly formed biogenic caco3 !! diss3 : vertical distribution of newly formed biogenic sio4 REAL atcco2, atcco20, atcox,h2co3(jpi,jpj),bobsx(jpi,jpj) REAL discl(jpi,jpj),akw3(jpi,jpj,jpk) REAL dissl(jpi,jpj),dmin(jpi,jpj),dminl(jpi,jpj),fugaci(jpi,jpj) REAL diss3(jpi,jpj,jpk),dmin3(jpi,jpj,jpk),disc3(jpi,jpj,jpk) REAL str,strn(jpi,jpj),fbod(jpi,jpj),sedlam,sedpoc(jpi,jpj) REAL sedcal(jpi,jpj),sedsil(jpi,jpj) REAL sedlostpoc, sedlostsil,sedlostcal COMMON/cchem2/ atcco2, atcco20, atcox, h2co3, bobsx, & &, akw3, dmin, fugaci, str, strn COMMON/cotsedim/sedlam,dmin3,dminl,discl,dissl,diss3,disc3,& & sedpoc,sedcal,sedsil,fbod,& & sedlostpoc, sedlostsil,sedlostcal !! !!---------------------------------------------------------------------- !! !! COMMON/cchem3/ : Variable for chemistry of the CO2 cycle !! !! --------------------------------------------------------------------- !! !! VARIABLE TYPE PURPOSE. !! -------- ---- -------- !! !! *THIRD* REAL 1./3. HALF PRECISION !! *TFREE* REAL APPROXIMATE VALUE OF FREEZING POINT OF !! SEAWATER !! *TBOIL* REAL SETS UPPER LIMIT FOR SEAWATER TEMPERATURE !! *SWEET* REAL SETS LOWER LIMIT FOR SALINITY !! *BRINE* REAL SETS UPPER LIMIT FOR SALINITY !! *TEN* REAL 10. HALF PRECISION !! *THOUSI* REAL 0.001 HALF PRECICION !! *SMICR* REAL 1E-6 HALF PRECISION !! *SALCHL* REAL CONVERSION FACTOR TO CALCULATE CHLORID !! CONCENTRATION !! S(O/OO)=1.80655*CL(O/OO) !! KALLE/DIETRICH , P. 60, OR !! WOOSTER ET AL., 1969. !! *TEMZER* REAL ZERO DEG C EXPRESSED AS ABS. TEMPERATURE !! (DEG KELVIN) !! *PERC* REAL 0.01 HALF PRECISION !! *AKCC1* REAL COEFFICIENT FOR EMPIRICALLY DETERMINED !! APPARENT SOLUBILITY PRODUCT K'SP OF !! CALCITE (INGLE, 1800, EQ. 6) !! (CF. BROECKER ET AL., 1982) !! *AKCC2* REAL COEFFICIENT FOR EMPIRICALLY DETERMINED !! APPARENT SOLUBILITY PRODUCT K'SP OF !! CALCITE (INGLE, 1800, EQ. 6) !! (CF. BROECKER ET AL., 1982) !! *AKCC3* REAL COEFFICIENT FOR EMPIRICALLY DETERMINED !! APPARENT SOLUBILITY PRODUCT K'SP OF !! CALCITE (INGLE, 1800, EQ. 6) !! (CF. BROECKER ET AL., 1982) !! *AKCC4* REAL COEFFICIENT FOR EMPIRICALLY DETERMINED !! APPARENT SOLUBILITY PRODUCT K'SP OF !! CALCITE (INGLE, 1800, EQ. 6) !! (CF. BROECKER ET AL., 1982) !! *ARAFRA* REAL FRACTION OF ARAGONITE IN BIOGENIC CACO3 !! PARTICLES (E.G. 0.3 FOR 30 PERCENT) !! *CALFRA* REAL FRACTION OF CALCITE IN BIOGENIC CACO3 !! PARTICLES (E.G. 0.3 FOR 30 PERCENT) !! *ARACAL* REAL FACTOR TO CONVERT APP. CALCITE SOLUBILITY !! PRODUCT (0 DBAR) TO THE APP. SOLUBILITY !! PRODUCT OF ARAGONITE (BERNER, 1976; !! CF. BROECKER ET AL., 1982) !! *DEVK1* REAL COEFFICIENT FOR SEAWATER PRESSURE CORRECTION !! OF 1. DISSOCIATION CONSTANT OF CARBONIC !! ACID AFTER CULBERSON AND PYTKOWICZ, 1968 !! (CF. BROECKER ET AL., 1982) !! *DEVK2* REAL COEFFICIENT FOR SEAWATER PRESSURE CORRECTION !! OF 2. DISSOCIATION CONSTANT OF CARBONIC !! ACID AFTER CULBERSON AND PYTKOWICZ, 1968 !! (CF. BROECKER ET AL., 1982) !! *DEVKB* REAL COEFFICIENT FOR SEAWATER PRESSURE CORRECTION !! OF 1. DISSOCIATION CONSTANT OF BORIC !! ACID AFTER CULBERSON AND PYTKOWICZ, 1968 !! (CF. BROECKER ET AL., 1982) !! *DEVK1T* REAL COEFFICIENT FOR SEAWATER PRESSURE CORRECTION !! OF FIRST DISSOCIATION CONSTANT OF CARBONIC !! ACID AFTER CULBERSON AND PYTKOWICZ, 1968 !! (CF. BROECKER ET AL., 1982) !! *DEVK2T* REAL COEFFICIENT FOR SEAWATER PRESSURE CORRECTION !! OF SECOND DISSOCIATION CONSTANT OF CARBONIC !! ACID AFTER CULBERSON AND PYTKOWICZ, 1968 !! (CF. BROECKER ET AL., 1982) !! *DEVKBT* REAL COEFFICIENT FOR SEAWATER PRESSURE CORRECTION !! OF DISSOCIATION CONSTANT OF BORIC !! ACID AFTER CULBERSON AND PYTKOWICZ, 1968 !! (CF. BROECKER ET AL., 1982) !! *DEVKS* REAL COEFFICIENT FOR PRESSURE CORRECTION OF !! SOLUBILITY PRODUCT OF CALCITE OR ARAGONITE !! AFTER EDMOND AND GIESKES (1970), P. 1285 !! (REFERENCE TO CULBERSON AND PYTKOWICZ, 1968, !! AS DONE IN BROECKER ET AL., 1982, IS !! NOT CORRECT) !! *DEVKST* REAL COEFFICIENT FOR PRESSURE CORRECTION OF !! SOLUBILITY PRODUCT OF CALCITE OR ARAGONITE !! AFTER EDMOND AND GIESKES (1970), P. 1285, !! IN TERM WITH TEMPERATURE !! (REFERENCE TO CULBERSON AND PYTKOWICZ, 1968, !! AS DONE IN BROECKER ET AL., 1982, IS !! NOT CORRECT) !! *RGAS* REAL UNIVERSAL GAS CONSTANT (BOLTZMANN'S CONSTANT !! TIMES AVOGADRO'S CONSTANT = !! 1.3804E-16*6.023*10E+23=8.3143E+7 ERG/GRD*MOL= !! 83.143E+6 ERG/GRD*MOL=8.3143 J/K*MOL) !! MULTIPLIED WITH 10 (TO A!!OUNT FOR !! CHANGE FROM BAR TO DBAR) !! (CF. EDMOND AND GIESKES, 1970. P. 1285, !! BROECKER ET AL., 1982, P. 79) !! *BOR1* REAL TOTAL BORON CONTENT IN G/KG AT CL=19 O/OO !! (S=35) !! (CF. RILEY AND SKIRROW, VOL. 1, P. 648) !! *BOR2* REAL INVERSE OF ATOMIC WEIGHT OF BORON FOR !! CONVERTING SPECIFIC TOTAL BORATE IN !! CONCENTRATIONS !! *OXYCO* REAL INVERS OF NORMAL MOLAL VOLUME OF AN !! IDEAL GAS [CM**-3] !! *C00* REAL VOLUMETRIC SOLUBILITY CONSTANT A1 FOR !! THE SOLUBILITY OF CO2 IN ML/L FROM AIR !! AT ONE ATMOSPHERE (WEISS, 1974) !! *C01* REAL VOLUMETRIC SOLUBILITY CONSTANT A2 FOR !! THE SOLUBILITY OF CO2 IN ML/L FROM AIR !! AT ONE ATMOSPHERE (WEISS, 1974) !! *C02* REAL VOLUMETRIC SOLUBILITY CONSTANT A3 FOR !! THE SOLUBILITY OF CO2 IN ML/L FROM AIR !! AT ONE ATMOSPHERE (WEISS, 1974) !! *C03* REAL VOLUMETRIC SOLUBILITY CONSTANT B1 FOR !! THE SOLUBILITY OF CO2 IN ML/L FROM AIR !! AT ONE ATMOSPHERE (WEISS, 1974) !! *C04* REAL VOLUMETRIC SOLUBILITY CONSTANT B2 FOR !! THE SOLUBILITY OF CO2 IN ML/L FROM AIR !! AT ONE ATMOSPHERE (WEISS, 1974) !! *C05* REAL VOLUMETRIC SOLUBILITY CONSTANT B3 FOR !! THE SOLUBILITY OF CO2 IN ML/L FROM AIR !! AT ONE ATMOSPHERE (WEISS, 1974) !! *C10* REAL COEFF. FOR 1. H2CO3 DISSOC. CONST. AFTER !! EDMOND AND GIESKES (1970) !! *C11* REAL COEFF. FOR F1. H2CO3 DISSOC. CONST. AFTER !! EDMOND AND GIESKES (1970) !! *C12* REAL COEFF. FOR 1. H2CO3 DISSOC. CONST. AFTER !! EDMOND AND GIESKES (1970) !! *C13* REAL COEFF. FOR 1. H2CO3 DISSOC. CONST. AFTER !! EDMOND AND GIESKES (1970) !! *C20* REAL COEFF. FOR 2. H2CO3 DISSOC. CONST. AFTER !! EDMOND AND GIESKES (1970) !! *C21* REAL COEFF. FOR 2. H2CO3 DISSOC. CONST. AFTER !! EDMOND AND GIESKES (1970) !! *C22* REAL COEFF. FOR 2. H2CO3 DISSOC. CONST. AFTER !! EDMOND AND GIESKES (1970) !! *C23* REAL COEFF. FOR 2. H2CO3 DISSOC. CONST. AFTER !! EDMOND AND GIESKES (1970) !! *CB0* REAL COEFF. FOR 1. H3BO3 DISSOC. CONST. AFTER !! EDMOND AND GIESKES (1970) !! *CB1* REAL COEFF. FOR 1. H3BO3 DISSOC. CONST. AFTER !! EDMOND AND GIESKES (1970) !! *CB2* REAL COEFF. FOR 1. H3BO3 DISSOC. CONST. AFTER !! EDMOND AND GIESKES (1970) !! *CB3* REAL COEFF. FOR 1. H3BO3 DISSOC. CONST. AFTER !! EDMOND AND GIESKES (1970) !! *CW0* REAL COEFF. FOR KW (DICKSON AND RILEY, 1979) !! *CW1* REAL COEFF. FOR KW (DICKSON AND RILEY, 1979) !! *CW2* REAL COEFF. FOR KW (DICKSON AND RILEY, 1979) !! (CORRECTED A!!ORDING TO B. BACASTOW, !! PERS. COMMUN., 1988) !! *OX0* REAL VOLUMETRIC SOLUBILITY CONSTANT A1 FOR !! THE SOLUBILITY OF OXYGEN IN ML/L FROM !! MOIST AIR AT ONE ATMOSPHERE (WEISS, 1970) !! *OX1* REAL VOLUMETRIC SOLUBILITY CONSTANT A2 FOR !! THE SOLUBILITY OF OXYGEN IN ML/L FROM !! MOIST AIR AT ONE ATMOSPHERE (WEISS, 1970) !! *OX2* REAL VOLUMETRIC SOLUBILITY CONSTANT A3 FOR !! THE SOLUBILITY OF OXYGEN IN ML/L FROM !! MOIST AIR AT ONE ATMOSPHERE (WEISS, 1970) !! *OX3* REAL VOLUMETRIC SOLUBILITY CONSTANT A4 FOR !! THE SOLUBILITY OF OXYGEN IN ML/L FROM !! MOIST AIR AT ONE ATMOSPHERE (WEISS, 1970) !! *OX4* REAL VOLUMETRIC SOLUBILITY CONSTANT B1 FOR !! THE SOLUBILITY OF OXYGEN IN ML/L FROM !! MOIST AIR AT ONE ATMOSPHERE (WEISS, 1970) !! *OX5* REAL VOLUMETRIC SOLUBILITY CONSTANT B2 FOR !! THE SOLUBILITY OF OXYGEN IN ML/L FROM !! MOIST AIR AT ONE ATMOSPHERE (WEISS, 1970) !! *OX6* REAL VOLUMETRIC SOLUBILITY CONSTANT B3 FOR !! THE SOLUBILITY OF OXYGEN IN ML/L FROM !! MOIST AIR AT ONE ATMOSPHERE (WEISS, 1970) !! *T* REAL DUMMY VARIABLE, ABSOLUTE SEAWATER TEMP. !! *qtt* REAL DUMMY VARIABLE, ABSOLUTE SEAWATER TEMP., !! DIVIDED BY 100. !! *S* REAL DUMMY VARIABLE, SALINITY !! *CL* REAL CHLORINITY (CL(O/OO)=S(O/OO)/1.80655) !! AFTER WOOSTER ET AL., 1969 !! (C.F. KALLE/DIETRICH , P. 60) !! *CEK0* REAL LN(K0), LOGARITHM OF CO2 SOLUBILITY IN !! SEAWATER IN VOLUMETRIC UNITS (EQ. 12 IN !! WEISS, 1974) !! *CK1* REAL PK1-VALUE (K1= 1. H2CO3 DISSOC. CONST..), !! AFTER EDMOND AND GIESKES (1970) !! *CK2* REAL PK2-VALUE (K2= 2. H2CO3 DISSOC. CONST..), !! AFTER EDMOND AND GIESKES (1970) !! *CKB* REAL PKB-VALUE (KB= 1. H3BO3 DISSOC. CONST..) !! AFTER EDMOND AND GIESKES (1970) !! *CKW* REAL PKW-VALUE (KW=H2O DISSOC. CONST.) AFTER !! DICKSON AND RILEY (1979) !! *OXY* REAL LN(C*), LOGARITHM OF O2 SOLUBILITY IN !! SEAWATER IN VOLUMETRIC UNITS (EQ. 4 IN !! WEISS, 1970) !! *AK1* REAL K1, 1. H2CO3 DISSOC. CONSTANT !! (EDMOND AND GIESKES, 1970) !! *AK2* REAL K2, 2. H2CO3 DISSOC. CONSTANT !! (EDMOND AND GIESKES, 1970) !! *AKB* REAL KB, 1. H3BO3 DISSOC. CONSTANT !! (EDMOND AND GIESKES, 1970) !! *AKW* REAL KW, H2O DISSOC. CONSTANT, LIT ? !! *AK0* REAL EXP(LN(K0))=K0 CO2 SOLUBILITY IN SEAWATER !! IN VOLUMETRIC UNITS (ML/L)(WEISS, 1974, !! CF. EQ. 12) !! *RRR* REAL SIGMA-T IN OCEAN MODEL, DUMMY VARIABLE !! (USED FOR CALCULATION OF TOTAL BORAT !! CONCENTRATION) (SIGMA-T=RHO(S,T,0)/1000.) !! *BOR* REAL TOTAL BORAT CONCENTRATION , DUMMY VAR. !! *TENM7* REAL 10.**-7 HALF PRECISION !! *P* REAL APPROXIMATE PRESSURE AT DEPTH OF U-POINTS !! IN BAR, DUMMY VARIABLE !! *AKSP0* REAL CACO3 SOLUBILITY PRODUCT AT P=0 DBAR !! A!!ORDING TO INGLE (1800), EQ. 6; THE !! CITATION OF CULBERSON AND PYTKOWICZ, 1968, !! IN BROECKER ET AL., 1982, IS PRESUMABLY !! NOT CORRECT) !! *AKSP(jpi,jpj,jpk)* REAL CACO3 SOLUBILITY PRODUCT AT IN SITU PRESSURE !! FOLLOWING THE PROCEDURE DESCRIBED IN EDMOND !! AND GIESKES (1970), P. 1285 !! *CP* REAL TERM IN EXPONENT OF EQUATIONS FOR PRESSURE !! CORRECTION OF DISSOC. CONSTANTS (CARB., !! BOR. ACID) AND CALCITE SOLUB. PRODUCT !! (CF. BROECKER ET AL., 1982, EDMOND AND !! GIESKES, 1970) !! *TC* REAL TEMPERATURE AT OCEAN GRID POINTS (DEG C), !! DUMMY VARIABLE !! *KI* INTEGER COUNTS ITERATIONS FOR NEWTON-RAPHSON METHOD !! FOR INITIATION OF [CO3--] AND [H+] !! *H* REAL [H+], DUMMY VARIABLE !! *R* REAL [CO3--] [MOLE/L], DUMMY VARIABLE !! *ALKA* REAL GIVEN ALKALINITY [EQV/L], DUMMY VARIABLE !! *C* REAL GIVEN [SUM(12C)O2] [MOLE/L], DUMMY VARIABLE !! *A* REAL ALKALINITY [EQV/L] AS FUNCTION OF [CO3--] !! AND [H+], DUMMY VARIABLE !! *DCDS* REAL LOCAL DERIVATIVE !! [SUM(CO2)]([CO3--],H+]) -> [CO3--] !! *DADS* REAL LOCAL DERIVATIVE !! [ALK]([CO3--],H+]) -> [CO3--] !! *DCDH* REAL LOCAL DERIVATIVE !! [SUM(CO2)]([CO3--],H+]) -> [H+] !! *DADH* REAL LOCAL DERIVATIVE !! [ALK]([CO3--],H+]) -> [H+] !! *F* REAL FUNCTION [ALK] GIVEN MINUS [ALK] IN TERMS !! OF [CO3--] AND [H+] !! *!!* REAL [SUM(CO2)] [MOLE/L] AS FUNCTION OF !! [CO3--] AND [H+], DUMMY VARIABLE !! *GG* REAL FUNCTION [SUM(CO2)] GIVEN MINUS [SUM(CO2)] !! IN TERMS OF [CO3--] AND [H+] !! *DETI* REAL DETERMINANT WITH LOCAL DERIVATIVES FOR !! NEWTON-RAPHSON ITERATION !! REAL brine, akcc1, akcc2, akcc3, akcc4 REAL arafra, calfra, aracal, devk1, devk2, devkb REAL devk1t, devk2t, devkbt, devkst, devks REAL bor1, bor2, c00, c01, c02, c03, c04, c05, c10, c11 REAL c12, c13, c20, c21, c22, c23, cb0, cb1, cb2, cb3 REAL c14, c15, c16, c17, c24, c25, c26, c27 REAL cb4, cb5, cb6, cb7, cb8, cb9, cb10, cb11 REAL cw3, cw4, cw5, cw6 REAL cw0, cw1, cw2, ox0, ox1, ox2, ox3, ox4 REAL cek0, ckb, ck1, ck2, ckw, ak1, ak2, ak0 REAL bor, aksp0, eak0, co2, ah, ba REAL tenm7, smicr, thousi, perc, third, ten, tfree, tboil REAL sweet, salchl, temzer, sucall, rgas, oxyco, ox5 REAL vi0, vi1, vi2, vi3, vi4, vi5, oxy, vit, rrr REAL soloxy, tco2, qtt REAL chemc(jpi,jpj,3), cpexp COMMON/cchem3/ brine, akcc1, akcc2, akcc3, akcc4 & &, arafra, calfra, aracal, devk1, devk2, devkb & &, devk1t, devk2t, devkbt, devkst, devks & &, bor1, bor2, c00, c01, c02, c03, c04, c05, c10, c11 & &, c12, c13, c20, c21, c22, c23, cb0, cb1, cb2, cb3 & &, c14, c15, c16, c17, c24, c25, c26, c27 & &, cb4, cb5, cb6, cb7, cb8, cb9, cb10, cb11 & &, cw3, cw4, cw5, cw6 & &, cw0, cw1, cw2, ox0, ox1, ox2, ox3, ox4 & &, cek0, ckb, ck1, ck2, ckw, ak1, ak2, ak0 & &, bor, aksp0, eak0, co2, ah, ba & &, tenm7, smicr, thousi, perc, third, ten, tfree, tboil & &, sweet, salchl, temzer, sucall, rgas, oxyco, ox5 & &, vi0, vi1, vi2, vi3, vi4, vi5, oxy, vit, rrr & &, soloxy, tco2, qtt, chemc, cpexp !! !!---------------------------------------------------------------------- !! !! COMMON/cotcon/ : Time variables !! !! --------------------------------------------------------------------- !! INTEGER nyrbeg, nyrend, nyr, nyrfend, nyrf PARAMETER (nyrbeg=1950,nyrend=1990) PARAMETER (nyr=nyrend+1-nyrbeg) PARAMETER (nyrfend=2300,nyrf=nyrfend+1-nyrbeg) INTEGER ntspmn, iabsyr REAL ttsec, ttyear, tspyr, absyr REAL rfact, rfactr INTEGER ilecvit REAL xtvit COMMON/cottim/ ttsec, ttyear, tspyr, ntspmn, absyr, iabsyr & &, rfact, rfactr, xtvit, ilecvit !! !!---------------------------------------------------------------------- !! !! COMMON/cotgas/ : Gas exchange !! !! --------------------------------------------------------------------- !! INTEGER igaswind,icice, nvit1p,nutrstm REAL gasfac,qcumul(jptra),tdepart(jptra),tcumul(jptra) REAL kgwanin(jpi,jpj), wsmo(jpi,jpj,12), kgwanmo(jpi,jpj,12) REAL cicemo(jpi,jpj,12) REAL patm(jpi,jpj),ppres(jpi,jpj,12) COMMON/cotgas/ gasfac, igaswind, nvit1p, icice & &,kgwanin, wsmo, kgwanmo, cicemo & &,patm,ppres,qcumul,tdepart,tcumul !! !!---------------------------------------------------------------------- !! !! COMMON/cotcon/ : Variable for testing convergence !! !! --------------------------------------------------------------------- !! REAL rconv(jpj), rconvs REAL rsumds, rsumdx COMMON/cotcon/ rconv, rconvs, rsumds, rsumdx !! !!---------------------------------------------------------------------- !! !! COMMON/cotham/ : biological parameters specific to HAMOCC3 !! !! --------------------------------------------------------------------- !! !! calcon = MEAN TOTAL [CA++] IN SEAWATER (MOLES/KG) !! (SEE BROECKER A. PENG, 1982, P. 26) !! ([CA++](MOLES/KG)=1.026E-2*(S/35.) AFTER !! CULKIN(1965), CF. BROECKER ET AL. 1982) !! spocri = variable for sediments !! conc0 = MICHAELIS MENTON CONSTANT OF PRODUCTIVITY (HALF SATU- !! RATION CONSTANT; NUTRIENT CONCENTRATION, WHERE NUTRIENT !! UPTAKE VELOCITY IS HALF ITS POTENTIAL MAXIMUM VALUE; !! E.G. DUGDALE (1967), PARSONS AND TAKAHASHI (1973)) !! rno3 = REDFIELD RATIO MOLES N (+ HPO4-) / MOLES C !! (FOR CHANGE IN ALKALINITY DUE TO PRODUCTION/REMINE- !! RALIZATION OF ORGANIC MATTER) !! N:C=16:122, SEE TAKAHASHI ET AL. (1985), P. 6923 !! !! COMPOSITION OF PLANKTONIC MATERIAL AND RESPIRATION: !! C106H263O110N16P + 138 O2 -> !! -> 106 CO2 + 16 NO3- + HPO4(2-) + 122 H2O + 18 H+ , !! SEE DEGENS ET AL. (1984), P. 152 !! *PO4R* REAL RATIO (MOLES P)/(MOLES C) (REDFIELD RATIO P:C) !! caco3r = ratio c-c(organic)/c-cacao3 in biogenic particles !! ----------------------------------------------------- !! achtung caco3r = 0.3 hat gut funktioniert bis auf zu wenig eqpco2 !! in expofa 0.5 auf 0.2 gesetzt !! sio2r = ratio moles c / moles sio4 in biogenic particles !! cristl = CRISTALLIZATION RATE IN CACO3 SUPERSATURATED WATER !! fecpef = FRACTION OF NEWLY FORMED BIOGENIC POC THAT FALLS !! IMMEDIATELY TO THE BOTTOM LAYER !! tobox = FRACTION OF NEWLY FORMED BIOGENIC CACO3 PARTICLES !! THAT FALLS IMMEDIATELY TO THE BOTTOM LAYER !! dispo0 = inertia concentration of caco3 for spontaneous dissolution [moles/l] !! oxymin = minimum [o2] [moles/l] for bacterial decomposition organic matter !! REAL caco3r, rno3, o2ut, sio2r, tobox ,po4r REAL prlat(jpi,jpj) REAL prcaca(jpi,jpj,jpk), silpro(jpi,jpj,jpk), prorca(jpi,jpj,jpk) REAL sco2, dispo0, conc0, silic0 REAL cristl, calcon, plafr13, pdb, c13at REAL grosip, oxymin, fecpef, bigshl, spocri REAL one, half, zero, two, four, fourth REAL pendec,pendes,fluexp REAL cmask(jpi,jpj),areacot INTEGER jpkb COMMON/cotham/ caco3r, rno3, o2ut, sio2r & &, tobox, po4r, prlat & &, prcaca, silpro, prorca & &, sco2, dispo0, conc0, silic0 & &, cristl, calcon, plafr13, pdb, c13at & &, grosip, oxymin, fecpef, bigshl, spocri & &, one, half, zero, two, four, fourth & &, pendec,pendes,fluexp,jpkb & & ,cmask, areacot #endif