!!---------------------------------------------------------------------
!!
!!                         COMMON passivetrc.pisces.h
!!                       ******************************
!!
!!  purpose :
!!  ---------
!!     INCLUDE COMMON FILE for PISCES biological model
!!
!!  modifications :
!!   -------------
!!      original    : 00-02 (O. Aumont)
!!
!!
#if defined key_trc_pisces
!!
!!----------------------------------------------------------------------
!!
!! Variable for chemistry of the CO2 cycle
!!
!! ---------------------------------------------------------------------
!!
      REAL(wp), DIMENSION(jpi,jpj,jpk) :: &
         akb3, ak13, ak23, aksp, akw3, &
         akp13, akp23, akp33, aksi3, aks3, akf3, &
         hi, borat

      REAL :: &
         atcco2, atcox

!!
!!----------------------------------------------------------------------
!!
!!  Variable for chemistry of the CO2 cycle
!!
!! ---------------------------------------------------------------------
!!
      REAL(wp), DIMENSION(10)  :: &
         devk1, devk2, devk3, devk4, devk5     
   
      REAL(wp) :: & 
         akcc1, akcc2, akcc3, akcc4,akcc5, akcc6, akcc7, akcc8, akcc9 , &
         bor1, bor2, c00, c01, c02, c03, c04, c05, c10, c11, &
         c12, c13, c20, c21, c22, c23, cb0, cb1, cb2, cb3, &
         cb4, cb5, cb6, cb7, cb8, cb9, cb10, cb11, c14, &
         cw3, cw4, cw5, cw6, cw0, cw1, cw2, ox0, ox1, ox2, ox3, ox4,ox5, &
         salchl, rgas, oxyco, ca0, ca1, ca2, ca3, ca4, ca5, ca6, &
         cp10, cp11, cp12, cp13, cp14, cp15, cp16, cp20, cp21, &
         cp22, cp23, cp24, cp25, cp26, cp30, cp31, cp32, cp33, &
         cp34, cp35, cs10, cs11, cs12, cs13, cs14, cs15, cs16, &
         cs17, cs18, cs19, cs20, cs21, &
         st1, st2, ft1, ft2, ks0, ks1, ks2, ks3, ks4, ks5, &
         ks6, ks7, ks8, ks9, ks10, ks11, ks12, kf0, kf1, &
         kf2, kf3, kf4

      REAL(wp), DIMENSION(jpi,jpj,3) :: &
         chemc

!!
!!----------------------------------------------------------------------
!!
!!  Variable for chemistry of Fe and SIO3
!!
!! ---------------------------------------------------------------------
!!
      REAL(wp), DIMENSION(jpi,jpj,jpk) :: &
         sio3eq, fekeq

!!
!!----------------------------------------------------------------------
!!
!!  Time variables
!!
!! ---------------------------------------------------------------------
!!

      INTEGER :: &
         iabsyr, nrdttrc, ndayflxtr

      REAL(wp) :: &
         tspyr, absyr, xtvit , &
         rfact, rfactr, rfact2, rfact2r


!!
!!----------------------------------------------------------------------
!!
!!  Gas exchange
!!
!! ---------------------------------------------------------------------
!!
      REAL(wp), DIMENSION(jpi,jpj) :: & 
         strn
!!---------------------------------------
!!
!!  Biological parameters 
!!
!! --------------------------------------
!! 
      INTEGER :: &
         jkopt

      REAL(wp) :: &
         caco3r, kdca, nca, part, rno3, o2ut, po4r,  &
         sco2, dispo0, conc0,sumdepsi,rivalkinput,sedfeinput, &
         calcon, rivpo4input,nitdepinput,oxymin,spocri, &
         nitrif,rdenit,o2nit,concnnh4,concdnh4, &
         pislope,excret,wsbio,resrat,mprat,wchl,wchld, &
         mzrat,grazrat,xprefc,xprefp,unass,xkgraz,xkmort, &
         xksi1,xksi2,sicmax,xremik,xremip,xkdoc1, &
         xkdoc2,grosip,resrat2,excret2,mprat2,mzrat2,xprefz, &
         xkgraz2,grazrat2,xlam1,conc1,conc2,conc3, &
         unass2,xprefpoc,epsher,epsher2,pislope2,mpratm, &
         sigma1, sigma2, zprefc, zprefp, zprefd, wchl2, ferat3, &
         xsirem, fecnm, fecdm, chlcnm, chlcdm



!!
!!---------------------------------------------
!!
!!  Biological fluxes for light
!!
!!---------------------------------------------
!!
      REAL(wp), DIMENSION(jpi,jpj,jpk) :: &
         etot, etot3, emoy

      REAL(wp), DIMENSION(jpi,jpj) :: &
         zmeu

      REAL(wp), DIMENSION(3,61) :: & 
         xkrgb


!!
!!----------------------------------------------------------
!!
!!  Biological fluxes for primary production
!!
!!----------------------------------------------------------
!!
      REAL(wp), DIMENSION(jpi,jpj,jpk) :: &
         prmax, tgfunc, tgfunc2, &
         prcaca, prorca, prorca2, prorca3, &
         prorca4, prorca5, prorca6, prorca7, &
         pronew, pronew2, proreg, proreg2, &
         xnanono3, xdiatno3, xnanonh4, xdiatnh4, &
         xlimphy, xlimdia, concdfe, concnfe, znegtr

      REAL(wp), DIMENSION(jpi,jpj) :: &      
         xksimax, xksi



!!
!!------------------------------------------
!!
!! Sinks for phytoplankton
!!
!!------------------------------------------
!!
      REAL(wp), DIMENSION(jpi,jpj,jpk) :: &      
         tortp, tortnf, tortnch, &
         respp, respp2, respnch, respdch, &
         tortp2, tortdf, tortdch, tortds, &
         respds, respdf, respnf


!!
!!------------------------------------
!!
!!  SMS for zooplankton
!!
!!-------------------------------------
!!
      REAL(wp), DIMENSION(jpi,jpj,jpk) :: &   
         respz, tortz, grazp, grazpf, &
         grazpch, grazm, grazmf, grazsd, &
         grazsf, grazss, grazsch, grarem, &
         grafer,respz2,tortz2,grazd, grazz,grazn, &
         grazpoc,graznf, graznch, grazs, grazf, &
         grazdch, grazpof, grarem2, grafer2, grapoc2, &
         grapoc, grazffe, grazfff


!!
!!---------------------------------------------
!!
!!  SMS for the organic matter
!!
!!---------------------------------------------
!!
      REAL(wp), DIMENSION(jpi,jpj,jpk) :: &   
         xagg, xaggfe, zdiss, xaggdoc, xaggdfe, xbactfer, &
         xscave, olimi, orem, orem2, ofer, ofer2, &
         osil, xaggdoc2, nitrfac, xlimbac, & 
         wsbio4, wsbio3, wscal

      REAL(wp) :: &
         wsbio2


 
!!
!!---------------------------------------------------------
!!
!! External sources of nutrients in ocean
!!
!!---------------------------------------------------------
!!
      REAL(wp), DIMENSION(jpi,jpj,jpk) :: &         
         onitr, denitr, ironsed

      REAL(wp), DIMENSION(jpi,jpj) :: & 
         dust, cotdep, nitdep, rivinp

      REAL(wp), DIMENSION(jpi,jpj,12) ::  &  
         dustmo

      REAL(wp) :: &
         areacot

      LOGICAL :: &
         bdustfer, briver, bndepo, bsedinput



#if defined key_trc_kriest
!!
!!---------------------------------------------------------
!!
!!  Kriest parameter for aggregation
!!
!!---------------------------------------------------------
!!
      REAL(wp) :: &
         xkr_eta, xkr_zeta, xkr_sfact, xkr_stick,  &
         xkr_mass_min, xkr_mass_max, xkr_massp, xkr_frac, &
         xkr_dnano, xkr_ddiat, xkr_dmeso, xkr_daggr, &
         xkr_nnano, xkr_ndiat, xkr_nmeso, xkr_naggr, &
         xkr_wsbio_min, xkr_wsbio_max
      
      REAL(wp), DIMENSION(jpk) :: &
         xnumm
#endif

#endif