anchor.diff:          enddo
           deallocate(hgg)!, psigg)
   
  +
  +        ! Flow over Andes mountain - adds topography
  +        if(testcase == 54)then
  +          allocate(alt_table(nlat_alt*nlon_alt,3))
  +          call getunit(iunit2)
  +          call andean_mountain_data(alt_table, iunit2)
  +          do i=1, mesh%nv
  +            lon=mesh%v(i)%lon
  +            lat=mesh%v(i)%lat
  +            bt%f(i) = smooth_andean_mountain([lat,lon], alt_table)
  +            !print*, i, lon*rad2deg, lat*rad2deg, bt%f(i)
  +            ! Correct h to allow orography
  +            h%f(i)=h%f(i)-bt%f(i)
  +          end do
  +          deallocate(alt_table)
  +        end if
  +
           !Set velocity field
           do l=1,mesh%ne
             utmp=0._r8
  
positive.code:     gravi, &
      sec2day, &
      day2sec, &
      rotatn
  
    !Global variables and operators for shallow water model
    use swm_data !Everything
  
negative.code:     div_cell_Cgrid, &
      grad_edge_Cgrid
  
    !Use topography and bilinear interpolation routines
    use refinter, only: &
     andean_mountain_data, &
     smooth_andean_mountain
  
    !use eispack, only: &
    !  rg
  