! source file: /Users/oschlies/UVIC/master/testcase/updates/npzd.h
!====================== include file "npzd.h" =========================

!   variables for npzd model

!   nbio    = number of npzd timesteps per ocean timestep
!   trcmin  = minimum tracer for flux calculations
!   alpha   = initial slope P-I curve [(W/m^2)^(-1)/day]
!   kw      = light attenuation due to water [1/m]
!   kc      = light attenuation by phytoplankton [1/(m*mmol m-3)]
!   ki      = light attenuation through sea ice & snow
!   abio    = maximum growth rate parameter [1/day]
!   bbio    = b
!   cbio    = [1/deg_C]
!   k1      = half saturation constant for N uptake [mmol m-3]
!   nup     = specific mortality rate (Phytoplankton) [day-1]
!   gamma1  = assimilation efficiency (zpk)
!   gbio    = maximum grazing rate [da-1y]
!   epsbio  = prey capture rate [(mmol m-3)-2 day-1]
!   nuz     = quadratic mortality (zpk)
!   gamma2  = excretion (zpk)
!   nud0    = remineralization rate [day-1]
!   LFe     = Iron limitation
!   wd      = sinking speed of detritus [m day-1]
!   ztt     = depth to top of grid cell [cm]
!   rkwz    = reciprical of light attenuation times grid depth
!   par     = fraction of photosythetically active radiation
!   dtnpzd  = time step of biology
!   rnpp    = rate of net primary production [nmol cm-3 s-1]
!   rgraz   = rate of grazing [nmol cm-3 s-1]
!   rmorp   = rate of mortality of phytoplankton [nmol cm-3 s-1]
!   rmorz   = rate of mortality of zooplankton [nmol cm-3 s-1]
!   rremi   = rate of remineralization [nmol cm-3 s-1]
!   rexcr   = rate of excretion [nmol cm-3 s-1]
!   rexpo   = rate of export through the bottom [nmol cm-3 s-1]
!   capr    = carbonate to carbon production ratio
!   dcaco3  = remineralisation depth of calcite [cm]
!   rcak    = array used in calculating calcite remineralization
!   rcab    = array used in calculating bottom calcite remineralization
!   nupt0   = specific mortality rate (Phytoplankton) [1/day]
!   wd0     = sinking speed of detritus at surface [m/day]
!   k1p     = half saturation constant for P uptake
!   jdiar   = factor reducing the growth rate of diazotrophs
!   rnpp_D  = npp for diazotraphs [nmol cm-3 s-1]
!   rgraz_D = rgraz for diazotraphs [nmol cm-3 s-1]
!   rmorp_D = rmorp for diazotraphs [nmol cm-3 s-1]
!   rnfix   = rate of nitrogen fixation [nmol cm-3 s-1]
!   rdeni   = rate of denitrification [nmol cm-3 s-1]
!   redctn  = C/N Redfield ratio (includes mol to mmol conversion)
!   redctp  = C/P Redfield ratio (includes mol to mmol conversion)
!   redptn  = P/N Redfield ratio
!   redntp  = N/P Redfield ratio
!   redotn  = O/N Redfield ratio (includes mol to mmol conversion)
!   redotp  = O/P Redfield ratio (includes mol to mmol conversion)
!   rnbio   = reciprical of nbio
!   rdtts   = reciprical of dtts [s-1]
!   dtbio   = npzd time step [s]

      integer ntnpzd, nbio
      parameter (ntnpzd=6)
      common /npzd_i/ nbio(km)

      real trcmin

      parameter (trcmin=5e-12)

      real alpha, kw, kc, ki, abio, bbio, cbio, k1, nup, gamma1, gbio
      real epsbio, nuz, gamma2, nud0, LFe, wd, ztt, rkwz, par, dtnpzd
      real rnpp, rgraz, rmorp, rmorpt, rmorz, rremi, rexcr, rexpo, capr
      real dcaco3, rcak, rcab, nupt0, wd0, k1p, jdiar, rnpp_D, rgraz_D
      real rmorp_D, rnfix, rdeni, redctn, redctp, redptn, redntp, redotn
      real redotp, rnbio, rdtts, dtbio
      real pco2z0

      common /npzd_r/ alpha, kw, kc, ki, abio, bbio, cbio, k1, nup
      common /npzd_r/ gamma1, gbio, epsbio, nuz, gamma2, nud0, LFe
      common /npzd_r/ wd(km), ztt(km), rkwz(km), par, dtnpzd, capr
      common /npzd_r/ dcaco3, rcak(km), rcab(km), nupt0, wd0, k1p
      common /npzd_r/ jdiar, redctn, redctp, redptn, redntp, redotn
      common /npzd_r/ redotp, rnbio(km), rdtts(km), dtbio(km)
      common /npzd_r/ pco2z0(imt,jmt)

      common /npzd_r/ rnpp(imt,kpzd,jsmw:jemw)
      common /npzd_r/ rgraz(imt,kpzd,jsmw:jemw)
      common /npzd_r/ rmorp(imt,kpzd,jsmw:jemw)
      common /npzd_r/ rmorpt(imt,kpzd,jsmw:jemw)
      common /npzd_r/ rmorz(imt,kpzd,jsmw:jemw)
      common /npzd_r/ rexcr(imt,kpzd,jsmw:jemw)
      common /npzd_r/ rremi(imt,km,jsmw:jemw)
      common /npzd_r/ rexpo(imt,km,jsmw:jemw)
      common /npzd_r/ rnpp_D(imt,kpzd,jsmw:jemw)
      common /npzd_r/ rgraz_D(imt,kpzd,jsmw:jemw)
      common /npzd_r/ rmorp_D(imt,kpzd,jsmw:jemw)
      common /npzd_r/ rnfix(imt,kpzd,jsmw:jemw)
      common /npzd_r/ rdeni(imt,km,jsmw:jemw)