Error Codes

       ! DL_MG_ERR_UNSPECIFIED = 1
       errm("Unspecified error:"), &
       ! Dl_MG_ERR_INIT = 2
       errm("dl_mg_init called before dl_mg_free"), &
       ! DL_MG_ERR_NEUMANM_BC = 3
       errm("Neumann boundary condition is not implemented"), &
       !DL_MG_ERR_UNKNOWN_BC = 4
       errm("dl_mg_init: wrong boundary condition specification &
       &in input. It must be one of: DL_MG_BC_PERIODIC, DL_MG_BC_DIRICHLET"),&
       !DL_MG_ERR_MPI_TOPO = 5
       errm("input communicator does not have a cartesian topology"), &
       !DL_MG_ERR_MPI_COMM_DUP = 6
       errm("duplicate communicator does not have a &
       &cartesian topology"), &
       !DL_MG_ERR_MPI_TOPO_BC = 7
       errm("boundary conditions inconsistent with MPI topology"),&
       ! DL_MG_ERR_GRID_EMPTY = 8
       errm("empty inner grid in set_mg_grids"),&
       ! DL_MG_ERR_PROLONG_MAP = 9
       errm("build_prolong_map: inactive site left behind"), &
       ! DL_MG_ERR_DEBYE = 10
       errm("dl_mg_init_nonlin: Debye length squared is negative"), &
       !DL_MG_ERR_NOINIT_POISSON = 11
       errm("dl_mg_solver_poisson is called before dl_mg_init"), &
       ! DL_MG_ERR_NOINIT_PBE = 12
       errm("dl_mg_solver_pbe called before initialisation subroutines"), &
       ! DL_MG_ERR_DERPOT = 13
       errm("dl_mg_solver_pbe: not an error, functionality added"), &
       ! DL_MG_ERR_MOD_BC = 14
       errm("Dirichlet boundary values were modified by the solver"), &
       ! DL_MG_ERR_RES_RATIO = 15
       errm("Consecutive residuals ratio is too large.&
       & It is possible that the round-off limit was reached.&
       & Check the residual ratio distribution histogram (enabled with DL_MG_RES_HIST)."), &
       ! DL_MG_ERR_NITER = 16
       errm("Computation failed to converge after the prescribed number of iterations"), &
       ! DL_MG_ERR_NEWTON_TYPE = 17
       errm("wrong eq_type in Newton's MG"), &
       ! DL_MG_ERR_NEWTON_DAMP = 18
       errm("compute_damping_parameter: failed to find a damping parameter for&
       & Newton method"), &
       ! DL_MG_ERR_NEWTON_DAMP_DERFUN = 19
       errm("compute_damping_parameter: functional derivative is non-negative"), &
       ! DL_MG_ERR_ASSERTION          = 20    ! JCW
       errm("assertion error"), &
       ! DL_MG_ERR_ALLOCATION         = 21    ! JCW
       errm("error in memory allocation"), &
       ! DL_MG_ERR_DEALLOCATION       = 22    ! JCW
       errm("error in memory deallocation"), &
       ! DL_MG_ERR_IO                 = 23    ! JCW
       errm("error in file I/O"),&
       ! DL_MG_ERR_GRID_GEOM          = 24    ! JCW
       errm("error in setting up grid geometry"), &
       ! DL_MG_ERR_MPI_FAIL           = 25    ! JCW
       errm("MPI routine returned a non-zero error code"), &
       ! DL_MG_ERR_EQTYPE = 26
       errm("Wrong equation type"), &
       ! DL_MG_ERR_FD_ORDER           = 27    ! JCW
       errm("requested FD order not supported"), &
       ! DL_MG_ERR_DEFCO_ITER         = 28    ! JCW
       errm("convergence was not achieved within the maximum number of &
       &defect correction iterations"), &
       ! DL_MG_ERR_DEFCO_UNPHYSICAL   = 29    ! JCW
       errm("unphysical solution obtained in defect correction"), &
       ! DL_MG_ERR_DEFCO_DAMPING      = 30    ! JCW
       errm("error damping for defect correction failed"), &
       ! DL_MG_ERR_DEFCO_DERPOT       = 31    ! JCW
       errm("der_pot should not be passed to the solver when the defect &
       &correction is used, since the potential at which the derivative is &
       &evaluated is generated internally&
       & (within the defect correction loop)"), &
       ! DL_MG_ERR_NOMPI              = 32
       errm("Using serial build code in MPI run?"), &
       ! DL_MG_ERR_CION_ZERO          = 33
       errm("Ion concentration is zero in PBE solver. &
       &Use Poisson or linearised PBE"), &
       ! DL_MG_ERR_CION_NEG           = 34
       errm("Shifted ion concentration is negative.") &