C $Header: /u/gcmpack/MITgcm/pkg/dic/DIC_VARS.h,v 1.10 2014/11/04 17:15:52 jmc Exp $ C $Name: $ #include "DIC_OPTIONS.h" C *==========================================================* C | DIC_VARS.h C | o Carbon Variables C *==========================================================* COMMON /CARBON_NEEDS/ & AtmospCO2, AtmosP, pH, pCO2, FluxCO2, & wind, & light, & riv_nutr, riv_dis, & FIce, Silica, Kwexch_Pre _RL AtmospCO2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) _RL AtmosP(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) _RL pH(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) _RL pCO2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) _RL FluxCO2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) _RL wind(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) _RL light(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) _RL riv_nutr(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) _RL riv_dis(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) _RL FIce(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) _RL Silica(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) _RL Kwexch_Pre(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) COMMON /CARBON_CHEM/ & ak0,ak1,ak2,akw,akb,aks,akf, & ak1p,ak2p,ak3p,aksi, fugf, & ff,ft,st,bt, Ksp_TP_Calc _RL ak0(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) _RL ak1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) _RL ak2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) _RL akw(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) _RL akb(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) _RL aks(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) _RL akf(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) _RL ak1p(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) _RL ak2p(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) _RL ak3p(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) _RL aksi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) _RL ff(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) C Fugacity Factor added by Val Bennington Nov. 2010 _RL fugf(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) _RL ft(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) _RL st(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) _RL bt(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) _RL Ksp_TP_Calc(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) COMMON /OXYGEN_CHEM/ & oA0,oA1,oA2,oA3,oA4,oA5, & oB0,oB1,oB2,oB3, & oC0 _RL oA0,oA1,oA2,oA3,oA4,oA5 _RL oB0,oB1,oB2,oB3 _RL oC0 C permil : is conversion factor for mol/m3 to mol/kg C assumes uniform (surface) density C Pa2Atm : for conversion of atmospheric pressure C when coming from atmospheric model COMMON /GLOBAL_SURF_MEAN/ & gsm_alk,gsm_s,gsm_t,gsm_dic, & gsm_c14,permil,Pa2Atm _RL gsm_alk _RL gsm_s _RL gsm_t _RL gsm_DIC _RL gsm_C14 _RL permil _RL Pa2Atm C schmidt number coefficients COMMON /DIC_SCHMIDT_NO/ & sca1, sca2, sca3, sca4, & sox1, sox2, sox3, sox4 _RL sca1 _RL sca2 _RL sca3 _RL sca4 _RL sox1 _RL sox2 _RL sox3 _RL sox4 C-- COMMON /DIC_FILENAMES/ C DIC_windFile :: file name of wind speeds C light_file C riv_nutr_file C riv_dis_file C DIC_atmospFile :: file name of atmospheric pressure C DIC_iceFile :: file name of seaice fraction C DIC_ironFile :: file name of aeolian iron flux C DIC_silicaFile :: file name of surface silica C DIC_parFile :: file name of photosynthetically available radiation (PAR) C DIC_chlaFile :: file name of chlorophyll climatology C DIC_forcingPeriod :: periodic forcing parameter specific for dic (seconds) C DIC_forcingCycle :: periodic forcing parameter specific for dic (seconds) C DIC_forcingPeriod2 C DIC_forcingCycle2 C dic_pCO2 :: Atmospheric pCO2 to be rad in data.dic C dic_int* :: place holder to read in a integer number, set at run time COMMON /DIC_FILENAMES/ & DIC_windFile, DIC_atmospFile, DIC_iceFile, & light_file, riv_nutr_file, riv_dis_file, & DIC_ironFile, DIC_silicaFile, DIC_parFile, & DIC_chlaFile, & DIC_forcingPeriod, DIC_forcingCycle, & DIC_forcingPeriod2, DIC_forcingCycle2, & dic_pCO2, dic_int1, dic_int2, dic_int3, dic_int4 CHARACTER*(MAX_LEN_FNAM) DIC_windFile CHARACTER*(MAX_LEN_FNAM) light_file CHARACTER*(MAX_LEN_FNAM) riv_nutr_file CHARACTER*(MAX_LEN_FNAM) riv_dis_file CHARACTER*(MAX_LEN_FNAM) DIC_atmospFile CHARACTER*(MAX_LEN_FNAM) DIC_iceFile CHARACTER*(MAX_LEN_FNAM) DIC_ironFile CHARACTER*(MAX_LEN_FNAM) DIC_silicaFile CHARACTER*(MAX_LEN_FNAM) DIC_parFile CHARACTER*(MAX_LEN_FNAM) DIC_chlaFile _RL DIC_forcingPeriod _RL DIC_forcingCycle _RL DIC_forcingPeriod2 _RL DIC_forcingCycle2 _RL dic_pCO2 INTEGER dic_int1 INTEGER dic_int2 INTEGER dic_int3 INTEGER dic_int4 #ifdef DIC_BIOTIC C *==========================================================* C | o Biological Carbon Variables C *==========================================================* COMMON /BIOTIC_NEEDS/ & BIOave, CARave, SURave, SUROave, pCO2ave, pHave, & fluxCO2ave, omegaCave, pfluxave, epfluxave, cfluxave, & DIC_timeAve, & par, alpha, rain_ratio, InputFe, omegaC, CHL, & Kpo4, DOPfraction, zcrit, KRemin, & KDOPremin,zca,R_op,R_cp,R_NP, R_FeP, & O2crit, alpfe, KScav, ligand_stab, ligand_tot, KFE, & freefemax, fesedflux_pcm, FeIntSec, & parfrac, k0, kchl, lit0, & alphaUniform, rainRatioUniform, & alphamax, alphamin, & calpha, crain_ratio, cInputFe, calpfe, feload, cfeload, & nlev, QSW_underice INTEGER nlev C For averages _RL BIOave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) _RL CARave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) _RL SURave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) _RL SUROave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) _RL pCO2ave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) _RL pHave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) _RL fluxCO2ave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) _RL OmegaCave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) _RL pfluxave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) _RL epfluxave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) _RL cfluxave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) _RL DIC_timeAve(nSx,nSy) C values for biogeochemistry C CHL :: chlorophyll climatology [mg/m3] C fesedflux_pcm :: ratio of sediment iron to sinking organic matter C FeIntSec :: Sediment Fe flux, intersect value in: C Fe_flux = fesedflux_pcm*pflux + FeIntSec _RL par(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) _RL alpha(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) _RL rain_ratio(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) _RL InputFe(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) _RL omegaC(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) _RL CHL(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) _RL Kpo4 _RL DOPfraction _RL zcrit _RL KRemin _RL KDOPremin _RL zca _RL R_op _RL R_cp _RL R_NP _RL R_FeP _RL O2crit _RL alpfe _RL fesedflux_pcm _RL FeIntSec _RL KScav _RL ligand_stab _RL ligand_tot _RL KFe _RL freefemax C values for light limited bio activity C k0 :: Light attentuation coefficient for water [1/m] C kchl :: Light attentuation coefficient fct of chlorophyll [m2/mg] _RL k0, kchl, parfrac, lit0 _RL alphaUniform _RL rainRatioUniform _RL alphamax, alphamin _RL calpha _RL crain_ratio _RL cInputFe(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) _RL calpfe _RL cfeload _RL feload(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) LOGICAL QSW_underice #endif /* DIC_BIOTIC */ CEH3 ;;; Local Variables: *** CEH3 ;;; mode:fortran *** CEH3 ;;; End: ***