Variant of copy for 2d dpsp

What is the new feature

New copy variants for 2d,4d,5d fields from double-precision to single-precision.

Mandatory steps before review

• Gitlab CI passes (Hint: use make format for linting)
• New feature is covered by additional unit tests
• Mark the merge request as ready by removing Draft:

Mandatory steps before merge

• Test coverage does not decrease
• Reviewed by a maintainer
• Incorporate review suggestions
• Prior to merging, please remove any boilerplate from the MR description, retaining only the Please describe your feature in a couple of words and describe important implementation details of the feature section to maintain
• Remember to edit the commit message and select the proper changelog category (feature/bugfix/other)

You are not supposed to merge this request by yourself, the maintainers of fortan-support take care of this action!

test/fortran/test_fortran_tools.f90

@@ -442,6 +442,34 @@ CONTAINS
     CALL ASSERT_EQUAL(assert_real_spdp_3d_array(dest, src), .TRUE.)
   END SUBROUTINE
 
+  SUBROUTINE Test_copy_4d_dpsp
+    REAL(dp) :: src(5, 5, 5, 5) = 1.0
+    REAL(sp) :: dest(5, 5, 5, 5)
+
+    CALL TAG_TEST("Test_copy_4d_dpsp_ones")
+    CALL copy(src, dest, .FALSE.)
+    CALL ASSERT_EQUAL(assert_real_spdp_4d_array(dest, src), .TRUE.)
+
+    CALL RANDOM_NUMBER(src)
+    CALL TAG_TEST("Test_copy_4d_dpsp_random")
+    CALL copy(src, dest, .FALSE.)
+    CALL ASSERT_EQUAL(assert_real_spdp_4d_array(dest, src), .TRUE.)
+  END SUBROUTINE
+
+  SUBROUTINE Test_copy_5d_dpsp
+    REAL(dp) :: src(5, 5, 5, 5, 5) = 1.0
+    REAL(sp) :: dest(5, 5, 5, 5, 5)
+
+    CALL TAG_TEST("Test_copy_5d_dpsp_ones")
+    CALL copy(src, dest, .FALSE.)
+    CALL ASSERT_EQUAL(assert_real_spdp_5d_array(dest, src), .TRUE.)
+
+    CALL RANDOM_NUMBER(src)
+    CALL TAG_TEST("Test_copy_5d_dpsp_random")
+    CALL copy(src, dest, .FALSE.)
+    CALL ASSERT_EQUAL(assert_real_spdp_5d_array(dest, src), .TRUE.)
+  END SUBROUTINE
+
   SUBROUTINE Test_copy_2d_spdp
     REAL(sp) :: src(10, 10) = 1.0
     REAL(dp) :: dest(10, 10)
@@ -1909,7 +1937,7 @@ CONTAINS
       DO j = 1, SIZE(array1, 2)
         DO k = 1, SIZE(array1, 3)
           DO l = 1, SIZE(array1, 4)
-            IF (array1(i, j, k, l) /= array2(i, j, k, l)) THEN
+            IF (array1(i, j, k, l) /= REAL(array2(i, j, k, l), KIND=sp)) THEN
               assert_real_spdp_4d_array = .FALSE.
               EXIT
             END IF
@@ -1930,7 +1958,7 @@ CONTAINS
         DO k = 1, SIZE(array1, 3)
           DO l = 1, SIZE(array1, 4)
             DO m = 1, SIZE(array1, 5)
-              IF (array1(i, j, k, l, m) /= array2(i, j, k, l, m)) THEN
+              IF (array1(i, j, k, l, m) /= REAL(array2(i, j, k, l, m), KIND=sp)) THEN
                 assert_real_spdp_5d_array = .FALSE.
                 EXIT
               END IF