We consider the symmetric exclusion process with jumps given by a symmetric, translation invariant, transition probability p(·). The process is put in contact with stochastic reservoirs whose strength is tuned by a parameter θ ∈ IR. Depending on the value of the parameter θ and the range of the transition probability p(·) we obtain the hydrodynamical behavior of the system. The type of hydrodynamic equation depends on whether the underlying probability p(·) has finite or infinite variance and the type of boundary condition depends on the strength of the stochastic reservoirs, that is, it depends on the value of θ. More precisely, when p(·) has finite variance we obtain either a reaction or reaction-diffusion equation with Dirichlet boundary conditions or the heat equation with different types of boundary conditions (of Dirichlet, Robin or Neumann type). When p(·) has infinite variance we obtain a fractional reaction-diffusion equation given by the regional fractional laplacian with Dirichlet boundary conditions but for a particular strength of the reservoirs.