don't use rank() in argument to avoid oneapi 2025.x ICE

Author: scivision

src/attr_read.inc

@@ -1,6 +1,7 @@
 integer :: ier, attr_class
+integer, parameter :: rA = rank(A)
 integer(HID_T) :: attr_id, space_id
-integer(HSIZE_T) :: attr_dims(rank(A))
+integer(HSIZE_T) :: attr_dims(rA)
 logical :: is_scalar
 
 attr_dims = shape(A, HSIZE_T)
@@ -11,8 +12,8 @@ call estop(ier, "attr_read:H5Aopen_by_name", self%filename, obj_name, attr_name)
 call H5Aget_space_f(attr_id, space_id, ier)
 call estop(ier, "attr_read:H5Aget_space", self%filename, obj_name, attr_name)
 
-if(rank(A) == 0) then
-  call hdf_rank_check(self, obj_name // ":" // attr_name, space_id, rank(A), is_scalar)
+if(rA == 0) then
+  call hdf_rank_check(self, obj_name // ":" // attr_name, space_id, rA, is_scalar)
 else
   call hdf_shape_check(self, obj_name // ":" // attr_name, space_id, attr_dims)
 endif

src/read_scalar.f90

@@ -22,7 +22,7 @@
 call H5Dget_space_f(dset_id, file_space_id, ier)
 call estop(ier, "read_scalar:H5Dget_space", self%filename, dname)
 
-call hdf_rank_check(self, dname, file_space_id, rank(A), is_scalar)
+call hdf_rank_check(self, dname, file_space_id, 0, is_scalar)
 
 if (is_scalar) then
   call hdf_get_slice(dims, dset_id, file_space_id, mem_space_id, [1], [1])

test/test_minimal.f90

@@ -24,7 +24,7 @@ program test_minimal
 if (i/=0) error stop 'minimal: could not create file'
 print *, 'minimal: created '//filename
 
-call h5ltmake_dataset_f(lid, "A", rank(p), shape(p, kind=HSIZE_T), h5kind_to_type(kind(p),H5_INTEGER_KIND), p, i)
+call h5ltmake_dataset_f(lid, "A", 0, shape(p, kind=HSIZE_T), h5kind_to_type(kind(p),H5_INTEGER_KIND), p, i)
 if (i/=0) error stop 'minimal: could not create dataset A'
 print *, 'minimal: created variable'