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