add comment + combined fypp loops

Author: jvdp1

src/stdlib_stats_median.fypp

@@ -83,12 +83,10 @@ contains
         select case(dim)
           #:for fi in range(1, rank+1)
           case(${fi}$)
-            #:for fj in range(1, 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}$)
             #:endfor
-              #:for fj in range(fi+1, rank+1)
-                do j${"_" * fj}$ = 1, size(x, ${fj}$)
-              #:endfor
                   x_tmp(:) = x${select_subvector('j', rank, fi)}$
                   call sort(x_tmp)
 
@@ -177,12 +175,10 @@ contains
         select case(dim)
           #:for fi in range(1, rank+1)
           case(${fi}$)
-            #:for fj in range(1, 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}$)
             #:endfor
-              #:for fj in range(fi+1, rank+1)
-                do j${"_" * fj}$ = 1, size(x, ${fj}$)
-              #:endfor
                   x_tmp = pack(x${select_subvector('j', rank, fi)}$, &
                                 mask${select_subvector('j', rank, fi)}$)
                   call sort(x_tmp)