導入：なぜSELECT RANKが重要なのか
数値計算の現場では、同じ物理計算アルゴリズムを「1次元（ベクトル）」「2次元（面）」「3次元（空間）」のそれぞれに対して実装することがよくあります。従来、これらは名前を変えた複数のルーチンを用意するか、あるいは単一の巨大な多次元配列に無理やり詰め込む必要がありました。しかし、Fortran 2018で導入されたSELECT RANK構文を使えば、配列の次元数を実行時に判定し、適切な処理へ自動的に分岐させることが可能です。これにより、コードの保守性が劇的に向上し、バグの混入リスクを最小限に抑えられます。
基礎知識：SELECT RANKの仕組み
SELECT RANKは、仮引数として渡された「想定ランク（Assumed Rank）配列」の次元数を調べるための制御構文です。想定ランク配列とは、引数の定義に `(..)` を用いることで、次元数を問わずあらゆる形状の配列を受け入れられる特別なインターフェースです。この構文を利用することで、プログラマは配列の次元を意識せず、計算対象の「中身」に集中した記述ができるようになります。
実装：想定ランク配列の活用手順
実装手順はシンプルです。まず、サブルーチンの引数定義で配列の次元を指定せず `(..)` と記述します。次に、サブルーチン内部で `SELECT RANK(a)` を呼び出し、各次元（rank）ごとの具体的な処理を記述します。引数の受け渡し時に次元情報が暗黙的に保持されるため、特別なメモリ管理を意識する必要はありません。
サンプルプログラム：次元を意識しない汎用加算ルーチン
以下は、1次元から3次元までの配列を受け取り、各要素に値を加算する汎用的なプログラム例です。
! 汎用加算サブルーチン
subroutine add_value(a, val)
! 想定ランク配列として宣言
real, intent(inout) :: a(..)
real, intent(in) :: val
! 次元数に応じて処理を分岐
select rank(a)
rank(1)
! 1次元配列の場合
a = a + val
rank(2)
! 2次元配列の場合
a = a + val
rank(3)
! 3次元配列の場合
a = a + val
rank()
! 想定外の次元数が渡された場合の例外処理
print , "エラー: 未対応の次元数です"
end select
end subroutine
! メイン処理（動作確認用）
program test_rank
real :: vec(3) = [1.0, 2.0, 3.0]
call add_value(vec, 10.0)
print , "結果:", vec
end program
応用・注意点：現場での活用と落とし穴
SELECT RANKを使用する際、最も注意すべきは「ランクごとの最適化」です。単にコードをまとめるだけでなく、各ランクのブロック内で、その次元に特化したループ最適化（SIMD化など）を記述できるのがこの構文の最大の利点です。
また、陥りやすいバグとして、`rank()`（未定義のランク）の書き忘れがあります。開発中の予期せぬ入力に対して安全に停止できるよう、必ずデフォルトケースを記述するようにしましょう。さらに、想定ランク配列は非常に強力ですが、配列の形状（Shape）そのものは `shape(a)` 関数で取得できることを覚えておくと、より柔軟な境界値処理が可能になります。
コメント