anchor.diff:      ! them consistent.
       call enforce_depth_thickness_consistency(h, eta, depth, &
           freesurfFac, thickness_error, nx, ny, layers)
  -
     end if
   
     !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
  
positive.code:   integer :: num_procs, myid
    integer :: ilower(0:num_procs-1,2), iupper(0:num_procs-1,2)
  
    if (myid .eq. 0) then
    ! Calculate the barotropic velocities
    call calc_baro_u(ub, unew, hnew, eta, freesurfFac, nx, ny, layers)
    call calc_baro_v(vb, vnew, hnew, eta, freesurfFac, nx, ny, layers)
  
negative.code:   integer*8 :: hypre_x
    integer   ::  nentries, nvalues, stencil_indices(5)
    double precision, dimension(:), allocatable :: values
    integer   :: temp(2)
  
    integer :: MPI_COMM_WORLD
    integer :: myid
  