median: replace _ by numbers

src/common.fypp

@@ -169,17 +169,17 @@ ${prefix + joinstr.join([line.strip() for line in txt.split("\n")]) + suffix}$
 #!
 #! E.g.,
 #!  select_subvector('j', 5, 2)
-#!   -> (j_, :, j___, j____, j_____)
+#!   -> (j1, :, j3, j4, j5)
 #!
 #:def select_subvector(varname, rank, idim)
   #:assert rank > 0
   #:call join_lines(joinstr=", ", prefix="(", suffix=")")
     #:for i in range(1, idim)
-      ${varname}$${ "_" * (i) }$
+      ${varname}$${i}$
     #:endfor
     :
     #:for i in range(idim + 1, rank + 1)
-      ${varname}$${ "_" * (i) }$
+      ${varname}$${i}$
     #:endfor
   #:endcall
 #:enddef
@@ -197,17 +197,17 @@ ${prefix + joinstr.join([line.strip() for line in txt.split("\n")]) + suffix}$
 #!
 #! E.g.,
 #!  reduce_subvector('j', 5, 2)
-#!   -> (j_, j___, j____, j_____)
+#!   -> (j1, j3, j4, j5)
 #!
 #:def reduce_subvector(varname, rank, idim)
   #:assert rank > 0
   #:if rank > 1
     #:call join_lines(joinstr=", ", prefix="(", suffix=")")
       #:for i in range(1, idim)
-        ${varname}$${ "_" * (i) }$
+        ${varname}$${i}$
       #:endfor
       #:for i in range(idim + 1, rank + 1)
-        ${varname}$${ "_" * (i) }$
+        ${varname}$${i}$
       #:endfor
     #:endcall
   #:endif

src/stdlib_stats_median.fypp

@@ -64,7 +64,7 @@ contains
         integer :: c, n
         #:if rank > 1
         #:for fj in range(1, rank+1)
-        integer :: j${"_" * fj}$
+        integer :: j${fj}$
         #:endfor
         #:endif
         ${t1}$, allocatable :: x_tmp(:)
@@ -84,7 +84,7 @@ contains
           case(${fi}$)
             ! Loop over every dimension of the array except "dim"
             #:for fj in list(range(1, fi)) + list(range(fi+1, rank+1))
-              do j${"_" * fj}$ = 1, size(x, ${fj}$)
+              do j${fj}$ = 1, size(x, ${fj}$)
             #:endfor
                 x_tmp(:) = x${select_subvector('j', rank, fi)}$
                 call sort(x_tmp)
@@ -162,7 +162,7 @@ contains
         integer(kind = int64) :: c, n
         #:if rank > 1
         #:for fj in range(1, rank+1)
-        integer :: j${"_" * fj}$
+        integer :: j${fj}$
         #:endfor
         #:endif
         ${t1}$, allocatable :: x_tmp(:)
@@ -176,7 +176,7 @@ contains
           case(${fi}$)
             ! Loop over every dimension of the array except "dim"
             #:for fj in list(range(1, fi)) + list(range(fi+1, rank+1))
-              do j${"_" * fj}$ = 1, size(x, ${fj}$)
+              do j${fj}$ = 1, size(x, ${fj}$)
             #:endfor
                 x_tmp = pack(x${select_subvector('j', rank, fi)}$, &
                               mask${select_subvector('j', rank, fi)}$)