anchor.diff:          umen = u00/en
           totvol = 0.0D0
           totarea = 0.0D0
  +
  +
   
           ! Integrate to tabulate h and psi as functions of geographical
           ! latitude
  
positive.code:         momeq_exact%f=0
  
  
        case(21, 22, 23, 52, 53) !  ! Galewsky et al test case - From J. Thuburns code
  
          if(testcase==23)then
            u00 = 200.0
            lat0 = pi/7.0
            lat1 = pi/2.0 - lat0
          else
            if(testcase==52 .or. testcase==53)then ! Jet in Southern Hemisphere
               u00 = 80.0
               lat0 = -5.d0*deg2rad
               lat1 = -45.d0*deg2rad
  
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.
  