anchor.diff:          Kb(i)  = Kb(i)*total2SWS_0p
   
   !       Already on SEAWATER scale: K1p, K2p, K3p, Kb, Ksi, Kw
  -!       K1p(i) = K1p(i)*total2SWS_0p 
  -!       K2p(i) = K2p(i)*total2SWS_0p
  -!       K3p(i) = K3p(i)*total2SWS_0p
  -!       Ksi(i) = Ksi(i)*total2SWS_0p
  -!       Kw(i)  = Kw(i)*total2SWS_0p
  -
  -!       Keep other contants on same scale:
  -!          - Ks (already on Free scale)
  -!          - Kspc, Kspa (no scale)
  +
  +!       Other contants (keep on another scale):
  +!          - K0         (independent of pH scale, already pressure corrected)
  +!          - Ks         (already on Free scale;   already pressure corrected)
  +!          - Kf         (already on Total scale;  already pressure corrected)
  +!          - Kspc, Kspa (independent of pH scale; pressure-corrected below)
   
   !       Perform actual pressure correction (on seawater scale)
           K1(i)   = K1(i)*EXP(lnkpok0(1))
  
positive.code:   REAL(kind=r8) :: total2free_0p, free2SWS_0p, total2SWS_0p
  ! REAL(kind=r8) :: free2SWS, free2SWS_0p
  
    REAL(kind=r8) :: R
  
    REAL(kind=r8) :: pK1o, ma1, mb1, mc1, pK1
  
negative.code:                      St, Ft, Bt,                              &
                       temp, sal, Patm,                         &
                       depth, lat, N,                           &
                       optT, optP, optB, optK1K2, optKf, optGAS)
  
    !   Purpose:
    !     Compute thermodynamic constants
  