感谢中科院大气所的Huang Yong-jie给出解决方案：

>A = (/ 1, 2, 3, 4, 5, 6, 7, 8, 9/)

>B = (/ 2, 14, 7, 1/)

>C = A(ind(.not.ismissing(get1Dindex(B,A))))

>print(C)

Variable: C
Type: integer
Total Size: 12 bytes
           3 values
Number of Dimensions: 1
Dimensions and sizes:   [3]
Coordinates:
(0)     1
(1)     2
(2)     7