POLEAVG.f

      SUBROUTINE poleavg(IM,JM,JSTA,JEND,SMALL,COSL,SPVAL,VAR)
! This program averages scalor fields at pole points
! 13-05-06 Shrinivas Moorthi : Added protection from divide by zero for icount
      implicit none

!     INCLUDE 'mpif.h'
      INTEGER,intent(in) :: im,jm,jsta,jend
      REAL,intent(in)    :: small,spval
      REAL,intent(in)    :: cosl(im,jsta:jend)
      REAL,intent(inout) :: var(im,jsta:jend)
      INTEGER i,jj,icount
      REAL    work, tem
!
      jj = 1
      IF(jj>=jsta .and. jj<=jend)then
        IF(cosl(1,jj) < small)then
          work   = 0.
          icount = 0
          DO i=1,im
            IF(var(i,jj) /= spval)THEN
              work   = var(i,jj) + work
              icount = icount    + 1
            END IF  
          END DO
          if (icount > 0) then
            tem = work/icount
            DO i=1,im
              var(i,jj) = tem
            END DO
          endif
        END IF      
      END IF
      jj = jm
      IF(jj>=jsta .and. jj<=jend)then
        IF(cosl(1,jj) < small)then
          work   = 0.
          icount = 0
          DO i=1,im
            IF(var(i,jj) /= spval)THEN
              work   = var(i,jj) + work
              icount = icount    + 1
            END IF
          END DO
          if (icount > 0) then
            tem = work/icount
            DO i=1,im
              var(i,jj) = tem
            END DO
          endif
        END IF
      END IF
      RETURN
      END