SUBROUTINE GRD2MSC(A,IM,JM,IGRD,B,MXCHPS,NCHP)

      _RL A(IM,JM), B(MXCHPS)

         B(I) = A(IGRD(I),1)

      SUBROUTINE MSC2GRD(IGRD,CHFR,B,MXCHPS,NCHP,FRACG,A,IM,JM)

      _RL A(IM,JM), B(MXCHPS), CHFR(MXCHPS), FRACG(IM,JM)

         VT1(IGRD(I),1) = VT1(IGRD(I),1) + B(I)*CHFR(I)

        B(I,J) = A(I+OFFSET-1,J)                                                

        B(I,J) = A(I+OFFSET-1,J)                                                

        B(Lena+I,J) = A(Lena+OFFSET-1,J)                                            

      SUBROUTINE STRIPINT(A,B,IA,IB,L,K)                                           

      INTEGER A(IA,L), B(IB,L)                                                

      SUBROUTINE VTRI0 ( A,B,C,F,Y,K,irun)

      _RL A(irun,K),B(irun,K),C(irun,K),Y(irun,K+1)

        A(I,L) = 1. / ( A(I,L) - B(I,LM1) * C(I,L) )

       Y(I,L) = (F(I,L) + B(I,L) * Y(I,L+1)) * A(I,L)

      SUBROUTINE VTRI1 ( A,B,Y,K,irun)

      _RL A(irun,K),B(irun,K),Y(irun,K+1)

        Y(I,L) = B(I,L) * Y(I,L+1) * A(I,L)

      SUBROUTINE VTRI2 ( A,B,C,F,Y,K,irun)

      _RL A(irun,K),B(irun,K),C(irun,K),F(irun,K)

        Y(I,L) = (F(I,L) + B(I,L) * Y(I,L+1)) *  A(I,L)

           B(I,J)=QS1*6.11 _d +00*EXP(17.2694 _d +00

           A2(I,J)=-D1*UG(I,J)*TICE(I,J,bi,bj)-D1W*UG(I,J)*B(I,J)

            B(I,J)=QS1*(C1*TICE(I,J,bi,bj)**4+C2*TICE(I,J,bi,bj)**3

            A2(I,J)=-D1*UG(I,J)*TICE(I,J,bi,bj)-D1I*UG(I,J)*B(I,J)

            B(I,J)=XKS/(HSNOW(I,J,bi,bj)/HICE(I,J)+XKS/XKI)/HICE(I,J)

            A3(I,J)=4.0 _d +00*D3*(TICE(I,J,bi,bj)**3)+B(I,J)+D1*UG(I,J)

            B(I,J)=B(I,J)*(TB-TICE(I,J,bi,bj))

     &                     +(A1(I,J)+A2(I,J)+B(I,J))/A3(I,J) 

      _RL B          (1-OLx:sNx+OLx, 1-OLy:sNy+OLy)