SUBROUTINE SHC2GRID( LMAX, SHC, NLON, NLAT, GRID, P )

      REAL    SHC((1+LMAX)*(1+LMAX))     ! INPUT  spherical harmonics

      CALL SHCROTATE(LMAX,SHC,ANGLE)

         CALL SHC2ZONE( LMAX, SHC, XLAT1, XLAT2, HS, HN, P )

      CALL SHCROTATE(LMAX,SHC,-ANGLE)

      SUBROUTINE SHC4GRID( LMAX, SHC, NLON, NLAT, GRID, P )

      REAL    SHC((1+LMAX)*(1+LMAX))     ! OUTPUT spherical harmonics

         CALL SHC4ZONE(LMAX,SHC,XLAT1,XLAT2,HS,HN,P)

      CALL SHCROTATE(LMAX,SHC,-ANGLE)

      SUBROUTINE SHCROTATE(LMAX,SHC,ANGLE)

      REAL     SHC((1+LMAX)*(1+LMAX)) ! spherical harmonic coeffs.

               C = SHC(K+M)

               S = SHC(K-M)

               SHC(K+M) = COSA*C + SINA*S

               SHC(K-M) =-SINA*C + COSA*S

      SUBROUTINE SHC2ZONE(LMAX,SHC,XLAT1,XLAT2,HS,HN,P)

      REAL    SHC((1+LMAX)*(1+LMAX)) ! INPUT  spherical harmonic coeffs.

            DA   = Q*SHC(J+M)

            DB   = Q*SHC(J-M)

            DA   = Q*SHC(J+M)

            DB   = Q*SHC(J-M)

            DA   = Q*SHC(J+M)

            DB   = Q*SHC(J-M)

      SUBROUTINE SHC4ZONE(LMAX,SHC,XLAT1,XLAT2,HS,HN,P)

      REAL    SHC((1+LMAX)*(1+LMAX))    ! spherical harmonic coeffs.

            SHC(K+M) = SHC(K+M) + P(I1) * CE

            SHC(K-M) = SHC(K-M) + P(I1) * SE

            SHC(K+M) = SHC(K+M) + P(I2) * CO

            SHC(K-M) = SHC(K-M) + P(I2) * SO

            SHC(K+M) = SHC(K+M) + P(I) * CE

            SHC(K-M) = SHC(K-M) + P(I) * SE