diff --git a/src/forcad_nurbs_curve.f90 b/src/forcad_nurbs_curve.f90
index 34744fcf1..6a99ca824 100644
--- a/src/forcad_nurbs_curve.f90
+++ b/src/forcad_nurbs_curve.f90
@@ -1927,8 +1927,8 @@ pure subroutine ansatz(this, ie, ig, Tgc, dTgc_dXg, dL)
         type(nurbs_curve) :: th, th_e
         real(rk), allocatable :: dTtth_dXksi(:), Ttth(:), dTgc_dXt(:), dXg_dXt(:)
         real(rk) :: Xt, dXt_dXksi
-        real(rk), allocatable :: dXg_dXksi(:) !> Jacobian matrix
-        real(rk) :: det_dXg_dXksi !> Determinant of the Jacobian matrix
+        real(rk), allocatable :: dXg_dXksi(:) !! Jacobian matrix
+        real(rk) :: det_dXg_dXksi !! Determinant of the Jacobian matrix
 
         call gauss_leg([0.0_rk, 1.0_rk], this%degree, Xksi, Wksi)
 
diff --git a/src/forcad_nurbs_surface.f90 b/src/forcad_nurbs_surface.f90
index 4865dd9a4..b58223e70 100644
--- a/src/forcad_nurbs_surface.f90
+++ b/src/forcad_nurbs_surface.f90
@@ -2550,8 +2550,8 @@ pure subroutine ansatz(this, ie, ig, Tgc, dTgc_dXg, dA)
         integer, allocatable :: elem_th(:,:), elem_c(:,:), elem_ce(:)
         type(nurbs_surface) :: th, th_e
         real(rk), allocatable :: dTtth_dXksi(:,:), Ttth(:), dTgc_dXt(:,:), Xt(:), dXt_dXksi(:,:), dXg_dXt(:,:)
-        real(rk), allocatable :: dXg_dXksi(:,:) !> Jacobian matrix
-        real(rk) :: det_dXg_dXksi !> Determinant of the Jacobian matrix
+        real(rk), allocatable :: dXg_dXksi(:,:) !! Jacobian matrix
+        real(rk) :: det_dXg_dXksi !! Determinant of the Jacobian matrix
 
         call gauss_leg([0.0_rk, 1.0_rk], [0.0_rk, 1.0_rk], this%degree, Xksi, Wksi)
 
diff --git a/src/forcad_nurbs_volume.f90 b/src/forcad_nurbs_volume.f90
index 97eeae690..c2d2303f8 100644
--- a/src/forcad_nurbs_volume.f90
+++ b/src/forcad_nurbs_volume.f90
@@ -3280,8 +3280,8 @@ pure subroutine ansatz(this, ie, ig, Tgc, dTgc_dXg, dV)
         integer, allocatable :: elem_th(:,:), elem_c(:,:), elem_ce(:)
         type(nurbs_volume) :: th, th_e
         real(rk), allocatable :: dTtth_dXksi(:,:), Ttth(:), dTgc_dXt(:,:), Xt(:), dXt_dXksi(:,:), dXg_dXt(:,:)
-        real(rk), allocatable :: dXg_dXksi(:,:) !> Jacobian matrix
-        real(rk) :: det_dXg_dXksi !> Determinant of the Jacobian matrix
+        real(rk), allocatable :: dXg_dXksi(:,:) !! Jacobian matrix
+        real(rk) :: det_dXg_dXksi !! Determinant of the Jacobian matrix
 
         call gauss_leg([0.0_rk, 1.0_rk], [0.0_rk, 1.0_rk], [0.0_rk, 1.0_rk], this%degree, Xksi, Wksi)