anchor.diff:          momeq_exact%f=0
   
   
  -      case(21, 22, 23, 52, 53) !  ! Galewsky et al test case - From J. Thuburns code
  +      case(21, 22, 23, 52, 53, 54) !  ! 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
  +          if(testcase==52 .or. testcase==53 .or. testcase==54)then ! Jet in Southern Hemisphere
                u00 = 80.0
                lat0 = -5.d0*deg2rad
                lat1 = -45.d0*deg2rad
  
positive.code:     div_cell_Cgrid, &
      grad_edge_Cgrid
  
  
    !use eispack, only: &
    !  rg
  
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.
  