_RL T, D, W, QX, DQX

      T = MIN(TT,TMAX) - TICE

       qx = aw0+T*(aw1+T*(aw2+T*(aw3+T*(aw4+T*(aw5+T*aw6)))))

      dqx = bw0+T*(bw1+T*(bw2+T*(bw3+T*(bw4+T*(bw5+T*bw6)))))

       qx =     w *(aw0+T*(aw1+T*(aw2+T*(aw3+T*(aw4+T*(aw5+T*aw6))))))

     .    + (1.-w)*(ai0+T*(ai1+T*(ai2+T*(ai3+T*(ai4+T*(ai5+T*ai6))))))

      dqx =     w *(bw0+T*(bw1+T*(bw2+T*(bw3+T*(bw4+T*(bw5+T*bw6))))))

     .    + (1.-w)*(bi0+T*(bi1+T*(bi2+T*(bi3+T*(bi4+T*(bi5+T*bi6))))))

       qx = d0+T*(d1+T*(d2+T*(d3+T*(d4+T*(d5+T*d6)))))

      dqx = e0+T*(e1+T*(e2+T*(e3+T*(e4+T*(e5+T*e6)))))

      _RL  S, T, P, PR

      CALL SW_PTMP  (S,T,PR,P,rv)

      SUBROUTINE SW_ADTG  (S,T,P, rv)

      _RL S,T,P

      rv =      a0 + (a1 + (a2 + a3*T)*T)*T

     &     + (b0 + b1*T)*(S-sref)

     &     + ( (c0 + (c1 + (c2 + c3*T)*T)*T) + (d0 + d1*T)*(S-sref) )*P

     &     + (  e0 + (e1 + e2*T)*T )*P*P

      SUBROUTINE SW_PTMP  (S,T,P,PR, rv)

      _RL S,T,P,PR

      call SW_ADTG(S,T,P, adtg_val)

      th     = T + onehalf*del_th

      SUBROUTINE SW_TEMP( S, T, P, PR, rv)

      SR = ((((9.1697E-10*T+2.0816E-8)*T-9.9348E-7) * S + 

     &     (5.2787E-8*T-6.12293E-6)*T+3.47718E-5) *P + 

     &     (1.91075E-4 * SR + (-1.6078E-6*T-1.0981E-5)*T+2.2838E-3) * 

     &     S + ((-5.77905E-7*T+1.16092E-4)*T+1.43713E-3)*T-0.1194975) 

     &     *P + (((-5.3009E-4*T+1.6483E-2)*T+7.944E-2) * SR + 

     &     ((-6.1670E-5*T+1.09987E-2)*T-0.603459)*T+54.6746) * S + 

     &     (((-5.155288E-5*T+1.360477E-2)*T-2.327105)*T+148.4206)*T -

      _RL DPT,LAT,S,T

      Q = PR*((((-2.1687E-16*T+1.8676E-14)*T-4.6206E-13)*P+

     &     ((2.7759E-12*T-1.1351E-10)*(S-35.0)+

     &     ((-5.4481E-14*T+8.733E-12)*T-6.7795E-10)*T+

     &     1.8741E-8))*P+(-4.2393E-8*T+1.8932E-6)*(S-35.0)+

     &     ((6.6228E-10*T-6.836E-8)*T+8.5258E-6)*T+3.5803E-5)

      T = T + 0.5*Q

      x = PR*((((-2.1687E-16*T+1.8676E-14)*T-4.6206E-13)*P+

     &     ((2.7759E-12*T-1.1351E-10)*(S-35.0)+

     &     ((-5.4481E-14*T+8.733E-12)*T-6.7795E-10)*T+

     &     1.8741E-8))*P+(-4.2393E-8*T+1.8932E-6)*(S-35.0)+

     &     ((6.6228E-10*T-6.836E-8)*T+8.5258E-6)*T+3.5803E-5)

      T = T + 0.29289322*(x-Q)

      x = PR*((((-2.1687E-16*T+1.8676E-14)*T-4.6206E-13)*P+

     &     ((2.7759E-12*T-1.1351E-10)*(S-35.0)+

     &     ((-5.4481E-14*T+8.733E-12)*T-6.7795E-10)*T+

     &     1.8741E-8))*P+(-4.2393E-8*T+1.8932E-6)*(S-35.0)+

     &     ((6.6228E-10*T-6.836E-8)*T+8.5258E-6)*T+3.5803E-5)

      _RL FUNCTION SBO_RHO( DPT, LAT, S, T )

      T = T + 1.707106781*(x-Q)

      x = PR*((((-2.1687E-16*T+1.8676E-14)*T-4.6206E-13)*P+

     &     ((2.7759E-12*T-1.1351E-10)*(S-35.0)+

     &     ((-5.4481E-14*T+8.733E-12)*T-6.7795E-10)*T+

     &     1.8741E-8))*P+(-4.2393E-8*T+1.8932E-6)*(S-35.0)+

     &     ((6.6228E-10*T-6.836E-8)*T+8.5258E-6)*T+3.5803E-5)

      T = T + (x-2.0*Q)/6.0

     &     ((-1.6546E-6*T+1.0227E-4)*T-5.72466E-3) * SR +

     &     (((5.3875E-9*T-8.2467E-7)*T+7.6438E-5)*T-4.0899E-3)*T

     &     +8.24493E-1)*S + ((((6.536332E-9*T-1.120083E-6)

     &     *T+1.001685E-4)*T-9.095290E-3)*T+6.793952E-2)*T-28.263737

      REAL DEG2RAD, T, ARG

      T      = H(I+1)+H(I+1)

      H(I+3) = T*H(I+1) - 1.0

      H(I+4) = T*H(I+2)

      H(L)   = T*H(L-2)-H(L-4)

     +T , 1)

     +T , 1)

      _RL  S, T, P, PR

      SW_temp = SW_PTMP (S,T,PR,P)

      _RL FUNCTION SW_ADTG  (S,T,P)

      _RL S,T,P

      SW_ADTG =      a0 + (a1 + (a2 + a3*T)*T)*T

     &     + (b0 + b1*T)*(S-sref)

     &     + ( (c0 + (c1 + (c2 + c3*T)*T)*T) + (d0 + d1*T)*(S-sref) )*P

     &     + (  e0 + (e1 + e2*T)*T )*P*P

      _RL FUNCTION SW_PTMP  (S,T,P,PR)

      _RL S,T,P,PR

      del_th = del_P*sw_adtg(S,T,P)

      th     = T + onehalf*del_th

      _RL FUNCTION SW_TEMP( S, T, P, PR )

         T(j,bj) = 0.0

          T(j,bj) = T0(no_so) + T2(no_so)*P2(i,j,bj)

          LW(j,bj) = LW0 + LW1*(T(j,bj)-t_mlt)

      _RL T(1-OLy:sNy+OLy,nSy)