anchor.diff:    return
   end subroutine calc_baro_v
   
  -! -----------------------------------------------------------------------------
  -
  +! ---------------------------------------------------------------------------
   !> Calculate the free surface anomaly using the velocities
   !! timestepped with the tendencies excluding the free surface
   !! pressure gradient.
  -subroutine calc_eta_star(ub, vb, eta, etastar, freesurfFac, nx, ny, dx, dy, dt)
  +
  +subroutine calc_eta_star(ub, vb, eta, etastar, &
  +    freesurfFac, nx, ny, dx, dy, dt)
     implicit none
   
     double precision, intent(in)  :: ub(nx+1, ny)
  
positive.code: end program MIM
  
  
  ! ------------------------- Beginning of the subroutines ----------------------
  !
  ! -----------------------------------------------------------------------------
  
  !> Evaluate the Bornoulli Potential for n-layer physics.
  !! B is evaluated at the tracer point, for each grid box.
  
negative.code:   return
  end subroutine read_input_fileH_2D
  
  ! ---------------------------------------------------------------------------
  
  subroutine read_input_fileU(name, array, default, nx, ny, layers)
    implicit none
  