anchor.diff:    call wrap_fields_3D(v, nx, ny, layers)
     call wrap_fields_3D(h, nx, ny, layers)
   
  +  end if
     ! Now the model is ready to start.
     ! - We have h, u, v at the zeroth time step, and the tendencies at
     !   two older time steps.
  
positive.code:       spongeHTimeScale, spongeUTimeScale, spongeVTimeScale, &
        spongeH, spongeU, spongeV, &
        nx, ny, layers, RedGrav, DumpWind, &
        MPI_COMM_WORLD, num_procs, ilower, iupper, &
        hypre_grid)
    print *, 'Execution ended normally'
    stop 0
  end program aronnax
  
  
negative.code:           dx, dy, wetmask, hfacW, hfacS, dt, &
            maxits, eps, rjac, freesurfFac, thickness_error, &
            g_vec, nx, ny, layers, n, &
            MPI_COMM_WORLD, myid, num_procs, ilower, iupper, &
            hypre_grid, hypre_b, ierr)
            ! , &
            ! MPI_COMM_WORLD, num_procs, ilower, iupper, &
  