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!

src/mo_fortran_tools.F90

@@ -2261,6 +2261,22 @@ CONTAINS
       base_shape(2) = in_shape(2)
       CALL insert_dimension_r_dp_3_2_s(ptr_out, ptr_in(1, 1), &
                                        base_shape, new_dim_rank)
+      IF (in_stride(1) > 1 .OR. in_stride(2) > in_shape(1) &
+          .OR. base_shape(1) /= in_shape(1)) THEN
+        out_stride(1) = in_stride(1)
+        out_stride(2) = 1
+        out_shape(1:out_rank - 1) = in_shape
+        DO i = out_rank, new_dim_rank + 1, -1
+          out_shape(i) = out_shape(i - 1)
+          out_stride(i) = out_stride(i - 1)
+        END DO
+        out_stride(new_dim_rank) = 1
+        out_shape(new_dim_rank) = 1
+        out_shape = (out_shape - 1)*out_stride + 1
+        ptr_out => ptr_out(:out_shape(1):out_stride(1), &
+             &             :out_shape(2):out_stride(2), &
+             &             :out_shape(3):out_stride(3))
+      END IF
     ELSE
       out_shape(1:out_rank - 1) = SHAPE(ptr_in)
       DO i = out_rank, new_dim_rank + 1, -1
@@ -2322,6 +2338,22 @@ CONTAINS
       base_shape(2) = in_shape(2)
       CALL insert_dimension_r_sp_3_2_s(ptr_out, ptr_in(1, 1), &
                                        base_shape, new_dim_rank)
+      IF (in_stride(1) > 1 .OR. in_stride(2) > in_shape(1) &
+          .OR. base_shape(1) /= in_shape(1)) THEN
+        out_stride(1) = in_stride(1)
+        out_stride(2) = 1
+        out_shape(1:out_rank - 1) = in_shape
+        DO i = out_rank, new_dim_rank + 1, -1
+          out_shape(i) = out_shape(i - 1)
+          out_stride(i) = out_stride(i - 1)
+        END DO
+        out_stride(new_dim_rank) = 1
+        out_shape(new_dim_rank) = 1
+        out_shape = (out_shape - 1)*out_stride + 1
+        ptr_out => ptr_out(:out_shape(1):out_stride(1), &
+             &             :out_shape(2):out_stride(2), &
+             &             :out_shape(3):out_stride(3))
+      END IF
     ELSE
       out_shape(1:out_rank - 1) = SHAPE(ptr_in)
       DO i = out_rank, new_dim_rank + 1, -1
@@ -2383,6 +2415,22 @@ CONTAINS
       base_shape(2) = in_shape(2)
       CALL insert_dimension_i4_3_2_s(ptr_out, ptr_in(1, 1), &
                                      base_shape, new_dim_rank)
+      IF (in_stride(1) > 1 .OR. in_stride(2) > in_shape(1) &
+          .OR. base_shape(1) /= in_shape(1)) THEN
+        out_stride(1) = in_stride(1)
+        out_stride(2) = 1
+        out_shape(1:out_rank - 1) = in_shape
+        DO i = out_rank, new_dim_rank + 1, -1
+          out_shape(i) = out_shape(i - 1)
+          out_stride(i) = out_stride(i - 1)
+        END DO
+        out_stride(new_dim_rank) = 1
+        out_shape(new_dim_rank) = 1
+        out_shape = (out_shape - 1)*out_stride + 1
+        ptr_out => ptr_out(:out_shape(1):out_stride(1), &
+             &             :out_shape(2):out_stride(2), &
+             &             :out_shape(3):out_stride(3))
+      END IF
     ELSE
       out_shape(1:out_rank - 1) = SHAPE(ptr_in)
       DO i = out_rank, new_dim_rank + 1, -1
@@ -2442,6 +2490,22 @@ CONTAINS
       base_shape(2) = in_shape(2)
       CALL insert_dimension_l_3_2_s(ptr_out, ptr_in(1, 1), &
                                     base_shape, new_dim_rank)
+      IF (in_stride(1) > 1 .OR. in_stride(2) > in_shape(1) &
+          .OR. base_shape(1) /= in_shape(1)) THEN
+        out_stride(1) = in_stride(1)
+        out_stride(2) = 1
+        out_shape(1:out_rank - 1) = in_shape
+        DO i = out_rank, new_dim_rank + 1, -1
+          out_shape(i) = out_shape(i - 1)
+          out_stride(i) = out_stride(i - 1)
+        END DO
+        out_stride(new_dim_rank) = 1
+        out_shape(new_dim_rank) = 1
+        out_shape = (out_shape - 1)*out_stride + 1
+        ptr_out => ptr_out(:out_shape(1):out_stride(1), &
+             &             :out_shape(2):out_stride(2), &
+             &             :out_shape(3):out_stride(3))
+      END IF
     ELSE
       out_shape(1:out_rank - 1) = SHAPE(ptr_in)
       DO i = out_rank, new_dim_rank + 1, -1