anchor.diff:  program aronnax
     implicit none
   
  -#ifdef useExtSolver
     include 'mpif.h'
   
     integer, parameter :: HYPRE_PARCSR = 5555
  -#endif
   
     integer, parameter :: layerwise_input_length = 10000
     ! Resolution
  
positive.code: 
    namelist /MODEL/ hmean, depthFile, H0, RedGrav
  
  #ifdef useExtSolver
    namelist /PRESSURE_SOLVER/ nProcX, nProcY
  #endif
  
    namelist /SPONGE/ spongeHTimeScaleFile, spongeUTimeScaleFile, &
        spongeVTimeScaleFile, spongeHfile, spongeUfile, spongeVfile
  
negative.code:   integer, dimension(:,:), allocatable :: ilower, iupper
    integer, dimension(:,:), allocatable :: jlower, jupper
    integer*8 :: hypre_grid
    integer :: i, j
    integer   :: offsets(2,5)
  
  
    ! Set default values here
  