&                DUM2,          DUM3,            AA,      BB,

      AA = 0.5 - (0.633 + 0.330*VGCHIL(ITYP(ChNo)))*VGCHIL(ITYP(ChNo))

      BB = 0.877 * ( ONE - 2.*AA )

      CC =  ( AA + BB*SUNANG(ChNo) ) / SUNANG(ChNo)

      DUM2 = (ONE-PDIR(ChNo)) * ( BB*(ONE/3.+PIE/4.) + AA*1.5 )

      _RL AA(irun,NLEV),BB(irun,NLEV),SSDEV(irun,NLEV)

       AA(I,NLEV) = 1. / (1. + CLH * BB(I,NLEV) )

       BB(I,NLEV) = BB(I,NLEV) * AA(I,NLEV) * plk(I,nlev)

       SBAR(I,NLEV) = AA(I,NLEV) * (SH(I,NLEV) - SHSAT(I,NLEV))

       SSDEV(I,NLEV)=CT(I)*(AA(I,NLEV)*DSH(I,NLEV)

       BETAW1(I,NLEV) = BETAW(I,NLEV) + AA(I,NLEV) * FCC(I,NLEV)

        AA(I,L) = 1. / (1. + CLH * BB(I,L) )

        BB(I,L) = BB(I,L) * AA(I,L)

        SBAR(I,L) = AA(I,L) * (SH(I,L) - SHSAT(I,L))

       SSDEV(I,1) = XL(I,1)*(AA(I,1)*DSH(I,1)-

        SSDEV(I,L) = XL(I,L-1)*(AA(I,L)*DSH(I,L-1)-

        SSDEV(I,L) = XL(I,L)*(AA(I,L)*DSH(I,L)-

        BETAW1(I,L) = BETAW(I,L) + AA(I,L) * FCC(I,L) * BETAL(I,L)

      _RL AA(irun,nlev), BB(irun,nlev), CC(irun,nlev+1)

        AA(I,L) = 1. + CC(I,L) + BB(I,L)

         AA(I,L) = AA(I,L) - XX2(I,L)

      CALL VTRI0(AA,BB,CC,XX1,XX1,NLEV,irun)

      CALL VTRI1(AA,BB,XX2,NLEV,irun)

      IF(ITYPE.EQ.3) CALL VTRI2 (AA,BB,CC,XX2,XX2,NLEV,irun)

      _RL AA(IRUN,5),TEMP(IRUN)

        AA(I,K) = AA1(K)

        IF( INT2(I).EQ.1 )  AA(I,K) = AA2(K)

        IF( INT3(I).EQ.1 )  AA(I,K) = AA3(K)

        IF( INT4(I).EQ.1 )  AA(I,K) = AA4(K)

       VDZSEA(I)  =  ( AA(I,4) + AA(I,5) * VUSTAR(I) ) * VUSTAR(I)

       VZSEA(I)  =  AA(I,2) + ( AA(I,3) + VDZSEA(I) ) * VUSTAR(I)

       TEMP(I) = AA(I,1) / VUSTAR(I)

        VDZSEA(I)  =  3. * VDZSEA(I) -(AA(I,4)*VUSTAR(I) - AA(I,3))