//============================================================================
// 侵食方向Wrangle
// 勾配とユーザ定義の閾値に基づいて砂が滑る方向を決定します。
// これは、侵食/堆積サイクルにおける第1工程です。
//============================================================================

float grid_resolution = chf("grid_resolution"); // 各グリッドセルのサイズ
float slope_threshold = chf("slope_threshold"); // 砂が滑るために必要な最小勾配
float scale_factor = chf("scale_factor"); // 極端な結果を回避するためのスケール係数

// --- 滑りのための有効な勾配閾値を計算します ---
float sliding_threshold = grid_resolution * slope_threshold;

// --- 現在のセルの高さとグリッド位置をキャッシュ化します。 ---
int x = @ix;
int y = @iy;

/*
   4つの標準方位(左、右、上、下)を定義します。
   これらの方位は、勾配を評価し、あり得る滑べり方向を決定するために使用されます。
*/

vector dirs[] = array(
    set(-1, 0, 0),  // 左
    set( 1, 0, 0),  // 右
    set( 0, -1, 0), // 下
    set( 0,  1, 0)  // 上
);

vector slide_dirs[] = array();

/*
   すべての隣接セルをループして、勾配を評価します。
   隣接セルに向かった勾配が調整済みの閾値を超えると、
   その方向は砂が滑る有効な候補方向と見なされます。
*/

foreach (vector dir; dirs) {
    int nx = x + int(dir.x);
    int ny = y + int(dir.y);
    vector neighbor_pos = set(nx, ny, 0); // グリッド位置

    float h_neighbor = volumeindex(0, "height", neighbor_pos); // 隣接セルの高さ
    int unstable = (@height - h_neighbor) > sliding_threshold;

    if (unstable) {
        append(slide_dirs, dir); // 方向を保存
    }
}

// --- 有効な滑り方向が存在した場合、そのどれかを選択します。

int direction_code = -1; // 安定状態

if (len(slide_dirs) > 0) {
    float seed_x = @Time + v@P.x * @height;
    float seed_y = @Time + v@P.y * @height;
    float rand = noise(set(seed_x, seed_y, 0));

    int index = int(fit01(rand, 0, len(slide_dirs) - 0.001));
    vector chosen_dir = slide_dirs[index];

    if (chosen_dir == set(-1, 0, 0)) direction_code = 0;
    else if (chosen_dir == set(1, 0, 0)) direction_code = 1;
    else if (chosen_dir == set(0, -1, 0)) direction_code = 2;
    else if (chosen_dir == set(0, 1, 0)) direction_code = 3;
}

@delta = float(direction_code); // 整数レイヤとしての方向。
@mask = direction_code != -1 ? 1.0 : 0.0; // 不安定なセルをマークします。

//--- 後で利用できるように、重要な値をDetailアトリビュートとして保存します。

setdetailattrib(geoself(), "sliding_threshold", sliding_threshold, "set");
setdetailattrib(geoself(), "scale_factor", scale_factor, "set");

// ============================================================================
// 浸食の適用
// 事前に計算された方向に基づいて砂の移動を適用します。
// これは、侵食/堆積サイクルにおける第2工程です。
// ============================================================================

float sliding_threshold = detail(0, "sliding_threshold") * detail(0, "scale_factor");

/*
   4つの標準方位(左、右、上、下)を定義します。
   これらの方位は、勾配を評価し、あり得る滑べり方向を決定するために使用されます。
*/

vector dirs[] = array(
    set(-1, 0, 0), // 左
    set( 1, 0, 0), // 右
    set( 0, -1, 0), // 下
    set( 0, 1, 0)  // 上
);

float height_increment = 0;
int x = @ix;
int y = @iy;

int opposite[] = array(1, 0, 3, 2); // 正しい結果を得るために方向を反転します: 左 -> 右, 下 -> 上。

// --- deltaレイヤの値を読み込み、それらの値を現在処理中のセルに適用します ---

for (int i = 0; i < len(dirs); i++) {
    int nx = x + int(dirs[i].x);
    int ny = y + int(dirs[i].y);
    float neighbor_delta = volumeindex(0, "delta", set(nx, ny, 0));

    if (int(neighbor_delta) == opposite[i]) {
        height_increment += sliding_threshold;
    }
}

if (@delta != -1) {
    height_increment -= sliding_threshold;
}

// --- 累積した高さの変化を適用します ---

@height += height_increment;