!-----------------------------------------------------------------------------
Real Function Ave(X, N, Miss)
Integer :: N, i, icount
Real :: X(N), Mean, Miss
Mean=0.; icount=0
Do i=1, N
   if (X(i)/=Miss) then
      icount=icount+1; Mean=Mean+X(i)
   endif
Enddo
if (icount>0) then
   Mean=Mean/real(icount)
else
   Mean=Miss
endif
Ave=Mean
End Function Ave
!-----------------------------------------------------------------------------