anchor.diff:          umen = u00/en
           totvol = 0.0D0
           totarea = 0.0D0
  +
  +
   
           ! Integrate to tabulate h and psi as functions of geographical
           ! latitude
  
positive.code:         enddo
          deallocate(hgg)!, psigg)
  
          !Set velocity field
          do l=1,mesh%ne
            utmp=0._r8
  
negative.code:     character (len=256):: tctmp, atime
  
      !File units
      integer (i4):: iunit, iunit2
      logical::  ifile
  
      ! Andes smooth topography - testcases 54 and 55
      real(r8), allocatable :: alt_table(:,:)
  
      maxvel=1.
      maxh=1.
  