“where 语句”和“merge”内在函数功能强大,可以对数组中的选定位置进行操作,但它们不会将项目移动到数组的前面。使用带有显式索引的老式代码(可以打包成一个函数),例如:
k=1
do i=1, n
if (ymiss (i) == 1) then
y(k) = y(i)
k = k + 1
end if
end do
您可以使用“pack”内在函数通过数组内在函数来完成您想要的事情。将 ymiss 转换为逻辑数组: 0 --> .false., 1 --> .true.. 然后使用类似代码(在没有第二个索引的情况下测试):
y(1:ydim(t), t) = pack (y (:,t), ymiss (:,t))
编辑以添加示例代码,显示 Fortran 内在函数“where”、“count”和“pack”的使用。单独的“where”不能解决问题,但是“pack”可以。在此示例中,我使用“< -90”作为 NaN。OP 不需要步骤“y (ydim+1:LEN) = -99.0”,不需要使用这些元素。
program test1
integer, parameter :: LEN = 6
real, dimension (1:LEN) :: y = [3.0, 1.0, -99.0, 6.0, 2.0, -99.0 ]
real, dimension (1:LEN) :: y2
logical, dimension (1:LEN) :: ymiss
integer :: ydim
y2 = y
write (*, '(/ "The input array:" / 6(F6.1) )' ) y
where (y < -90.0)
ymiss = .false.
elsewhere
ymiss = .true.
end where
ydim = count (ymiss)
where (ymiss) y2 = y
write (*, '(/ "Masking with where does not rearrange:" / 6(F6.1) )' ) y2
y (1:ydim) = pack (y, ymiss)
y (ydim+1:LEN) = -99.0
write (*, '(/ "After using pack, and ""erasing"" the end:" / 6(F6.1) )' ) y
stop
end program test1
输出是:
输入数组:3.0 1.0 -99.0 6.0 2.0 -99.0
掩蔽 where 不重新排列:3.0 1.0 -99.0 6.0 2.0 -99.0
使用包后,“擦除”结束:3.0 1.0 6.0 2.0 -99.0 -99.0