C $Header: /home/ubuntu/mnt/e9_copy/MITgcm_contrib/bling/pkg/BLING_VARS.h,v 1.1 2014/05/23 17:33:42 mmazloff Exp $ C $Name: $ C ========================================================== C Carbon chemistry variables C ========================================================== COMMON /CARBON_NEEDS/ & AtmospCO2, AtmosP, pH, pCO2, FluxCO2, & wind, FIce, Silica _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,Nr,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 FIce(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) _RL Silica(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) C ========================================================== C Carbon and oxygen chemistry parameters C ========================================================== COMMON /CARBON_CHEM/ & ak0,ak1,ak2,akw,akb,aks,akf, & ak1p,ak2p,ak3p,aksi, fugf, & ff,ft,st,bt, & Ksp_TP_Calc,Ksp_TP_Arag _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) _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) _RL Ksp_TP_Arag(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 COMMON /GLOBAL_SURF_MEAN/ & permil,Pa2Atm 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 _RL permil _RL Pa2Atm COMMON /SCHMIDT_NO/ & sca1, sca2, sca3, sca4, & sox1, sox2, sox3, sox4 C Schmidt number coefficients _RL sca1,sca2,sca3,sca4 _RL sox1,sox2,sox3,sox4 C ========================================================== C Bling inputs (specified in data.bling) C ========================================================== COMMON /BLING_INPUTS/ & bling_windFile, bling_atmospFile, bling_iceFile, & bling_ironFile, bling_silicaFile, & bling_forcingPeriod, bling_forcingCycle, & bling_pCO2, & riverconc_dic, riverconc_alk, riverconc_nut, & riverconc_dom, riverconc_o2, riverconc_fe C bling_windFile :: file name of wind speeds C bling_atmospFile :: file name of atmospheric pressure C bling_iceFile :: file name of sea ice fraction C bling_ironFile :: file name of aeolian iron flux C bling_silicaFile :: file name of surface silica C bling_forcingPeriod :: period of forcing for biogeochemistry (seconds) C bling_forcingCycle :: periodic forcing parameter for biogeochemistry C bling_pCO2 :: Atmospheric pCO2 to be read in data.bling C riverconc_* :: River concentration CHARACTER*(MAX_LEN_FNAM) bling_windFile CHARACTER*(MAX_LEN_FNAM) bling_atmospFile CHARACTER*(MAX_LEN_FNAM) bling_iceFile CHARACTER*(MAX_LEN_FNAM) bling_ironFile CHARACTER*(MAX_LEN_FNAM) bling_silicaFile _RL bling_forcingPeriod _RL bling_forcingCycle _RL bling_pCO2 _RL riverconc_dic _RL riverconc_alk _RL riverconc_nut _RL riverconc_dom _RL riverconc_o2 _RL riverconc_fe C ========================================================== C Ecosystem variables and parameters C ========================================================== COMMON /BIOTIC_NEEDS/ & InputFe, & omegaC, & omegaAr, & irr_inst, & irr_mem, & P_sm, & P_lg, & Pstar, & Pc_0, & lambda_0, & resp_frac, & chl_min, & CtoP, & NtoP, & O2toP, & CatoP, & NUTfac, & alpha_max, & alpha_min, & theta_Fe_max_hi, & theta_Fe_max_lo, & gamma_biomass, & gamma_irr_mem, & gamma_DOM, & gamma_POM, & k_Fe, & k_O2, & k_NUT, & k_FetoP, & KFeLeq_max, & KFeLeq_min, & Fe_min, & IFeL, & kFe_org, & kFe_inorg, & FetoP_max, & FetoPsed, & Fe_lim_min, & remin_min, & O2_min, & ligand, & b_const, & kappa_eppley, & kappa_remin, & ca_remin_depth, & phi_DOM, & phi_sm, & phi_lg, & wsink0z, & wsink0, & wsinkacc, & parfrac, & alpfe, & k0, & epsln, & QSW_underice _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 omegaAr(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) _RL irr_inst(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) _RL irr_mem(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) _RL P_sm(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) _RL P_lg(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) _RL Pstar _RL Pc_0 _RL lambda_0 _RL resp_frac _RL chl_min _RL CtoP _RL NtoP _RL O2toP _RL CatoP _RL NUTfac _RL alpha_max _RL alpha_min _RL theta_Fe_max_hi _RL theta_Fe_max_lo _RL gamma_biomass _RL gamma_irr_mem _RL gamma_DOM _RL gamma_POM _RL k_Fe _RL k_O2 _RL k_NUT _RL k_FetoP _RL KFeLeq_max _RL KFeLeq_min _RL Fe_min _RL IFeL _RL kFe_org _RL kFe_inorg _RL FetoP_max _RL FetoPsed _RL Fe_lim_min _RL remin_min _RL O2_min _RL ligand _RL b_const _RL kappa_eppley _RL kappa_remin _RL ca_remin_depth _RL phi_DOM _RL phi_sm _RL phi_lg _RL wsink0z _RL wsink0 _RL wsinkacc _RL parfrac _RL alpfe _RL k0 _RL epsln LOGICAL QSW_underice CEH3 ;;; Local Variables: *** CEH3 ;;; mode:fortran *** CEH3 ;;; End: ***