anchor.diff:     IMPLICIT NONE
   
   !  Input & output variables:
  -   REAL(kind=r8), INTENT(in) :: temp
  +   REAL(kind=r8), INTENT(in) :: Tc
      REAL(kind=r8) :: sco2
   
  -   sco2 = 1920.4 - 135.6*temp  + 5.2122*temp**2 - 0.10939*temp**3  + 0.00093777*temp**4
  +   sco2 = 1920.4 - 135.6*Tc  + 5.2122*Tc**2 - 0.10939*Tc**3  + 0.00093777*Tc**4
   
      RETURN
   END FUNCTION sco2
  
positive.code:   !     ================
    !     kw660   = gas transfer velocity (piston velocity) for CO2 [m/s] 
    !               without T-dependant Schmidt number correction
    !               but accounting for sea ice fraction (See OCMIP2 design & HOWTO documents for details)
    !     xco2    = atmospheric mole fraction of CO2 [ppm]
    !     Patm    = atmospheric pressure at surface [atm]
    !     dz1     = depth of the top vertical layer of the model [m]
  
negative.code: END FUNCTION scco2
  
  !>    Compute Schmidt number for O2 in seawater from temperature
  FUNCTION sco2(Tc)
  
  !  Compute Schmidt number of O2 in seawater w/ formulation from Wanninkhof (Limnol. Oceanogr.: Methods 12, 2014, 351–362)
  !  Input is temperature in deg C.
  