Thanks so much for your advice. I will try to analyze the eel3d code. However, it currently displays that my computer does not have GNU library. What puzzles me is that GNU is installed on my computer. Next, I will try to compile and re-install the GNU compiler.
Besides, I encountered a problem when I changed the grid size and calculation domain. When I increased the calculation domain and changed the grid size. To make sure that there is more or less one Lagrangian point per Eulerian cell on the finest mesh level in my changed domain configuration, I adjusted the size of N (actual number of grid cells on coarsest grid level)to meet the above requirements. For example, I chose the grid with 40*16 grid cells in the coarsest-level patch, and I changed the value of N to 16(In this case, I chose the four-level Cartesian grid). But there was an error in .log file:
// Lower/Upper side of integration box is not aligned with sides on coarsest level in dimension 0/1. Modifying coordinate to nearest box side.//
No matter what changes I made to the coordinates, this error always appears. What caused it?
Best,
Dewu