C
!     CALL RTHS(RHS,IMX0,JMX0,JMAXMH,C0)

      SUBROUTINE RTHS(RHS,IMX0,JMX0,JMAXMH,C0)
      implicit none
C     ****
C     ****     INSERT RHS IN STENCIL C0(10)
C     ****
C     ****     INPUT:
C     ****       RHS((IMX0+1)/2,(JMX0+1)/2)   RIGHT HAND SIDE OF PDE
C     ****
C     ****       CS(JMX0) COSINE ARRAY OF JMX0 LATITUDE VALUES
C     ****
C     ****       IMX0, JMX0   DIMENSION OF FINEST GRID IS IMX0 BY JMX0
C     ****
C     ****     OUTPUT:
C     ****       C0(IMX0,JMX0,10)
C     ****         IN C0(I,J,10) FOR I=1,IMX0 AND J=1,JMX0
C     ****
! Args:
      integer,intent(in) :: imx0,jmx0,jmaxmh
      real,intent(in) :: RHS(IMX0,1)
      real,intent(out) :: C0(IMX0,JMX0,10)
!
! Local:
      integer :: j,jj,i
C     ****
C     ****     CALCULATE NUMERICAL QUADRATURES OF WEIGHTED RHS AND
C     ****       WEIGHTS.
C     ****
      DO 4 J = 1,JMX0
        JJ = JMAXMH-JMX0+J
        DO 4 I = 1,IMX0
          C0(I,J,10) = RHS(I,JJ)
    4 CONTINUE
      RETURN
      END
C