! source file: /Users/oschlies/UVIC/master/source/mom/neptune.F
       subroutine neptune

!=======================================================================
!     calculate neptune (maximum entropy) velocities

!     the option "neptune" provides a subgridscale parameterization
!     for the interaction of eddies and topography

!     reference:
!       Holloway, G., 1992: representing topographic stress for large
!     scale ocean models, J. Phys. Oceanogr., 22, 1033-1046

!     neptune is calculated as an equilibrium streamfunction given by
!     pnep=-f*snep*snep*hnep and is applied through eddy viscosity

!     hnep = model streamfunction depth
!     snep = spnep + (senep-spnep)*(0.5 + 0.5*cos(2.0*latitude))

!     the neptune length scale snep has a value of senep at the
!     equator and smoothly changes to spnep at the poles

!     based on code by: G. Holloway and M. Eby
!=======================================================================

      return
      end