sit=sit*permil

     &        sit/(1.0 + x/ksilocal) -

     &        sit/ksilocal/(1.0+x/ksilocal)**2. +

     I                       t,s,diclocal,pt,sit,ta,

     &          sit/(1.0 + x/ksilocal) -

        sit=sit/permil

     I                       t,s,diclocal,pt,sit,ta,

        _RL  t, s, pt, sit, ta

        sit=sit*permil

        siooh3g = sit*ksilocal / (ksilocal + hguess)

        sit=sit/permil

     I                       t,s,diclocal,pt,sit,ta,

        _RL  t, s, pt, sit, ta

        sit=sit*permil

        siooh3g = sit*ksilocal / (ksilocal + hguess)

        _RL  t, s, pt, sit, ta

        sit=sit/permil