search

xavierdechamps / Shallow_water_FV

Solve the 2D shallow water equations with a finite volume method based on a Q-scheme
10 stars 1 forks source link

Rectangles as cell element type #5

Closed xavierdechamps closed 2 years ago

xavierdechamps commented 2 years ago

Up to now only triangular elements were readable and usable by the program. It would be nice to be able to have the possibility to combine rectangles and triangles inside the same mesh. The main work would lie inside subroutines gmsh_operations.f90 and get_normal_to_cell.f90.

xavierdechamps commented 2 years ago

The bulk of the work was indeed inside subroutines gmsh_operations.f90 and get_normal_to_cell.f90 but the main programs were also affected by the change in the parameters of the subroutines. Now quadrangular / triangular elements are possible in the mesh. Nothing special to be done by the user with the exception of the mesh construction. Changes effective in commit 62d7227096ca46c5ff680b5f80f78da427249107