moment_dev: removed some explicit statements

Author: jvdp1

src/stdlib_experimental_stats_moment.fypp

@@ -27,7 +27,7 @@ contains
           return
         end if
 
-        n = real(size(x, kind = int64), ${k1}$)
+        n = size(x, kind = int64)
         mean = sum(x) / n
 
         #:if t1[0] == 'r'
@@ -57,7 +57,7 @@ contains
           return
         end if
 
-        n = real(size(x, kind = int64), dp)
+        n = size(x, kind = int64)
         mean = sum(real(x, dp)) / n
 
         res = sum((real(x, dp) - mean)**order) / n
@@ -86,11 +86,11 @@ contains
           return
         end if
 
-        res = 0._${k1}$
+        res = 0
         select case(dim)
           #:for fi in range(1, rank+1)
           case(${fi}$)
-            n = real(size(x, dim), ${k1}$)
+            n = size(x, dim)
             mean = sum(x, dim) / n
             do i = 1, size(x, dim)
               #:if t1[0] == 'r'
@@ -129,11 +129,11 @@ contains
           return
         end if
 
-        res = 0._dp
+        res = 0
         select case(dim)
           #:for fi in range(1, rank+1)
           case(${fi}$)
-            n = real(size(x, dim), dp)
+            n = size(x, dim)
             mean = sum(real(x, dp), dim) / n
             do i = 1, size(x, dim)
               res = res + (real(x${select_subarray(rank, [(fi, 'i')])}$, dp) - mean)**order
@@ -161,7 +161,7 @@ contains
         real(${k1}$) :: n
         ${t1}$ :: mean
 
-        n = real(count(mask, kind = int64), ${k1}$)
+        n = count(mask, kind = int64)
         mean = sum(x, mask) / n
 
         #:if t1[0] == 'r'
@@ -186,7 +186,7 @@ contains
 
         real(dp) :: n, mean
 
-        n = real(count(mask, kind = int64), dp)
+        n = count(mask, kind = int64)
         mean = sum(real(x, dp), mask) / n
 
         res = sum((real(x, dp) - mean)**order, mask) / n
@@ -210,11 +210,11 @@ contains
         real(${k1}$) :: n${reduced_shape('x', rank, 'dim')}$
         ${t1}$ :: mean${reduced_shape('x', rank, 'dim')}$
 
-        res = 0._${k1}$
+        res = 0
         select case(dim)
           #:for fi in range(1, rank+1)
           case(${fi}$)
-            n = real(count(mask, dim), ${k1}$)
+            n = count(mask, dim)
             mean = sum(x, dim, mask) / n
             do i = 1, size(x, dim)
               #:if t1[0] == 'r'
@@ -250,11 +250,11 @@ contains
         real(dp) :: n${reduced_shape('x', rank, 'dim')}$
         real(dp) :: mean${reduced_shape('x', rank, 'dim')}$
 
-        res = 0._dp
+        res = 0
         select case(dim)
           #:for fi in range(1, rank+1)
           case(${fi}$)
-            n = real(count(mask, dim), dp)
+            n = count(mask, dim)
             mean = sum(real(x, dp), dim, mask) / n
             do i = 1, size(x, dim)
               res = res + merge((real(x${select_subarray(rank, [(fi, 'i')])}$, dp) - mean)**order,&