Dear Wolfgang Thanks your reply and I'm sorry for this late answer, apparently I'm having problems with receiving emails from the mailing list.
I wasn't sure about which adaptive refinement indicator to use, so I followed the one suggested in step 9 as they mention it could be more reliably when the function doesn't have second derivatives (although I think is not my case) but mostly because it doesn't require boundary conditions details. Now, I'm facing two new problems. One is regarding the mesh input from gmsh. Apparently deal.ii is not happy with two volumes next to each other. I think there is a function for deleting duplicated nodes in deal.ii that I think it might help and I'm in the process of learning how to use it. On the other hand after refinement there's a problem accessing the point values for cells that happens to have hanging nodes. Apparently the problem comes from the function GridTools::find_active_cell_around_point. I had opened two threads in the mailing list with these questions. Once again I appreciate your help. Best Regards Javier _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
