// カスタムパラメータを定義します。

float base_angle_deg = chf("wind_angle");
float variation_pos_deg = chf("angle_variation_pos");
float variation_neg_deg = chf("angle_variation_neg");
float turbulence_strength = chf("turbulence_strength");
float wind_strength = chf("wind_strength");
float wind_shear_factor = chf("wind_shear_factor");

// 風向を差異ありで計算します。
// この部分では、角度の度単位をラジアンに変換します。
// オフセットは、時間、位置、任意のベクトルからランダム値を算出します。
// そして、その結果をvariation_neg_degからvariation_pos_degまでの範囲に再マッピングします。
// 最終的な風向を得るために、それらの値を加算します。

float base_angle_rad = radians(base_angle_deg + 180.0);
float angle_offset = fit01(rand(@Time + dot(@P, set(12.345, 67.89, 0))), -variation_neg_deg, variation_pos_deg);
float final_angle = base_angle_rad + radians(angle_offset);

// 乱流を含むメインの風ベクトル。
// final_angleのサイン値とコサイン値を組み合わせてベクトルを作成します。
// curlnoise関数は、その風に微小な乱流を加えます。
// 3つのベクトルのwind_dir_base、wind_turbulence、turbulence_strengthは、風速ではなく風向を取得するために正規化します。

vector wind_dir_base = set(cos(final_angle), 0.0, sin(final_angle));
vector wind_turbulence = curlnoise(@P * 0.1 + @Time);
vector wind_dir = normalize(wind_dir_base + wind_turbulence * turbulence_strength);

// Wind Shear(ウィンドシア:風向や風速が急激に変化する現象)。
// ウィンドシアは、XZ平面において90度の時計回りの回転をします。
// ウィンドシアの強度は、wind_strengthとwind_shear_factorの積です。

vector cross_wind_dir = normalize(set(-wind_dir.z, 0.0, wind_dir.x));
float cross_strength  = wind_strength * wind_shear_factor;

// target_posは、現在のボクセルの砂の目標位置を決めます。

vector target_pos = @P + wind_dir * wind_strength + cross_wind_dir * cross_strength;

// 確率係数。
// よりランダムな外観を得るために、このスクリプトは、ボクセル毎に輸送確率を計算します。

float prob_factor = fit01(rand(@Time + dot(@P, wind_dir)), 0.5, 1.0);
float transport_prob = wind_strength * prob_factor;

// この部分では、目標位置と現在のボクセル位置での地形の高さ値をサンプリングします。

float target_height = volumesample(0, "height", target_pos);
float current_height = volumesample(0, "height", @P);

// 現在のボクセルにおける地形の高さを設定します。
// このスクリプトは、ボクセルの砂が前のサンプルと計算に基づいて移動するかどうかを確認します。
// 砂が移動すれば、このスクリプトは、その砂を除去するか累積するか決めます。
// 砂が除去される場合、その砂の量が直接、地形の実際の高さに変換されます。

if (current_height > target_height) {
    if (rand(@P + @Time) < transport_prob) {
        float sand_amount = current_height - target_height;

        if (sand_amount > 0.0) {
            current_height -= sand_amount;
            target_height  += sand_amount;

            @height = current_height;
        }
    }
}