8 USE mod_gauss_points_2d_p1
9 USE mod_gauss_points_2d_p2
12 LOGICAL,
OPTIONAL,
INTENT(IN) :: edge_stab
14 IF (
PRESENT(edge_stab))
THEN 15 mesh%edge_stab=edge_stab
17 mesh%edge_stab=.false.
20 nws =
SIZE(mesh%jjs,1)
21 IF (nw==3 .AND. nws==2)
THEN 22 CALL gauss_points_2d_p1(mesh)
23 ELSE IF (nw==6 .AND. nws==3)
THEN 26 CALL gauss_points_2d_p2(mesh)
27 ELSE IF (nw==0 .AND. nws==0)
THEN 33 ELSE IF (nw==4 .AND. nws==3)
THEN 34 WRITE(*,*)
' Finite element not yet programmed ' 37 ELSE IF (nw==10 .AND. nws==6)
THEN 38 WRITE(*,*)
' Finite element not yet programmed ' 42 WRITE(*,*)
' Finite element not yet programmed '
subroutine, public gen_gauss(mesh, edge_stab)