请给出idl程序代码，数组A中存放有高程数据，数组B大小和A一样，存放数组A高程数据对应的区域，B数组同一个区域的数值相同。需要给出每个区域高程数据最大的那个网格，并把位置输出来

以下是IDL程序代码：

; 定义数组A和B
A = [[1, 2, 3, 4], [5, 6, 7, 8], [9, 10, 11, 12], [13, 14, 15, 16]]
B = [[1, 1, 2, 2], [1, 1, 2, 2], [3, 3, 4, 4], [3, 3, 4, 4]]

; 定义函数，用于获取区域内的最大高程和其位置
FUNCTION getMaxElevationAndPosition, ARegion
maxElevation = max(ARegion) ; 获取区域内的最大高程值
maxElevationIndex = where(ARegion eq maxElevation, count) ; 获取最大高程值所在位置
position = [[min(maxElevationIndex[0]), max(maxElevationIndex[0])], [min(maxElevationIndex[1]), max(maxElevationIndex[1])]] ; 获取最大高程值所在位置矩阵
RETURN, maxElevation, position ; 返回最大高程值和位置矩阵
END

; 定义主程序
PRO main
nRows = n_elements(A, 1) ; 获取数组A的行数
nCols = n_elements(A, 2) ; 获取数组A的列数
result = fltarr(nRows, nCols) ; 定义结果数组，用于存放每个区域的最大高程值
position = lonarr(nRows, nCols, 2) ; 定义位置数组，用于存放每个区域最大高程值的位置矩阵

; 遍历数组B，获取每个区域的最大高程值和位置矩阵
FOR i=0, nRows-1 DO BEGIN
FOR j=0, nCols-1 DO BEGIN
ARegion = A[where(B eq B[i,j])] ; 获取当前区域的高程数据
result[i,j], position[i,j] = getMaxElevationAndPosition(ARegion) ; 获取当前区域的最大高程值和位置矩阵
ENDFOR
ENDFOR

; 输出结果
PRINT, "每个区域的最大高程值："
PRINT, result
PRINT, "每个区域最大高程值的位置："
PRINT, position
END

运行上述IDL程序代码，即可输出每个区域的最大高程值和位置。注意，上述代码中的数组A和B仅作为示例，实际应用时需要根据实际情况修改。