From 98b47016473640a20abbc21177e70d83378b36b1 Mon Sep 17 00:00:00 2001
From: Adam Frisby
Date: Sun, 22 Jul 2007 02:48:53 +0000
Subject: * Aerobic erosion now uses Navier Stokes algorithms for wind
 calculations.

---
 .../Region/Terrain.BasicTerrain/TerrainEngine.cs   |  2 +-
 .../Channel/Manipulators/AerobicErosion.cs         | 29 +++++++++++++++++-----
 .../Channel/Manipulators/NavierStokes.cs           |  6 ++---
 3 files changed, 27 insertions(+), 10 deletions(-)

(limited to 'OpenSim/Region')

diff --git a/OpenSim/Region/Terrain.BasicTerrain/TerrainEngine.cs b/OpenSim/Region/Terrain.BasicTerrain/TerrainEngine.cs
index 5ef2963..2e40531 100644
--- a/OpenSim/Region/Terrain.BasicTerrain/TerrainEngine.cs
+++ b/OpenSim/Region/Terrain.BasicTerrain/TerrainEngine.cs
@@ -409,7 +409,7 @@ namespace OpenSim.Region.Terrain
             {
                 case "aerobic":
                     // WindSpeed, PickupMinimum,DropMinimum,Carry,Rounds,Lowest
-                    heightmap.AerobicErosion(Convert.ToDouble(args[2]), Convert.ToDouble(args[3]), Convert.ToDouble(args[4]), Convert.ToDouble(args[5]), Convert.ToInt32(args[6]), Convert.ToBoolean(args[7]));
+                    heightmap.AerobicErosion(Convert.ToDouble(args[2]), Convert.ToDouble(args[3]), Convert.ToDouble(args[4]), Convert.ToDouble(args[5]), Convert.ToInt32(args[6]), Convert.ToBoolean(args[7]), true);
                     break;
                 case "thermal":
                     heightmap.thermalWeathering(Convert.ToDouble(args[2]), Convert.ToInt32(args[3]), Convert.ToDouble(args[4]));
diff --git a/OpenSim/Region/Terrain.BasicTerrain/libTerrainBSD/Channel/Manipulators/AerobicErosion.cs b/OpenSim/Region/Terrain.BasicTerrain/libTerrainBSD/Channel/Manipulators/AerobicErosion.cs
index 870b9b2..97c2762 100644
--- a/OpenSim/Region/Terrain.BasicTerrain/libTerrainBSD/Channel/Manipulators/AerobicErosion.cs
+++ b/OpenSim/Region/Terrain.BasicTerrain/libTerrainBSD/Channel/Manipulators/AerobicErosion.cs
@@ -74,7 +74,7 @@ namespace libTerrain
         /// <param name="carry">The percentage of rock which can be picked up to pickup 0..1</param>
         /// <param name="rounds">The number of erosion rounds (recommended: 25+)</param>
         /// <param name="lowest">Drop sediment at the lowest point?</param>
-        public void AerobicErosion(double windspeed, double pickup_talus_minimum, double drop_talus_minimum, double carry, int rounds, bool lowest)
+        public void AerobicErosion(double windspeed, double pickup_talus_minimum, double drop_talus_minimum, double carry, int rounds, bool lowest, bool usingFluidDynamics)
         {
             Channel wind = new Channel(w, h) ;
             Channel sediment = new Channel(w, h);
@@ -83,7 +83,15 @@ namespace libTerrain
             wind = this.copy();
             wind.normalise();   // Cheap wind calculations
             wind *= windspeed;
-            wind.pertubation(30);   // Can do better later
+
+            if (usingFluidDynamics)
+            {
+                wind.navierStokes(20, 0.1, 0.0, 0.0);
+            }
+            else
+            {
+                wind.pertubation(30);   // Can do better later
+            }
 
             for (i = 0; i < rounds; i++)
             {
@@ -115,10 +123,19 @@ namespace libTerrain
                         }
                     }
                 }
-                sediment.pertubation(10); // Sediment is blown around a bit
-                sediment.seed++;
-                wind.pertubation(15);   // So is the wind
-                wind.seed++;
+
+                if (usingFluidDynamics)
+                {
+                    sediment.navierStokes(7, 0.1, 0.0, 0.1);
+                    wind.navierStokes(10, 0.1, 0.0, 0.0);
+                }
+                else
+                {
+                    wind.pertubation(15);   // Can do better later
+                    wind.seed++;
+                    sediment.pertubation(10); // Sediment is blown around a bit
+                    sediment.seed++;
+                }
 
                 // Convert some sand to rock
                 for (x = 1; x < w - 1; x++)
diff --git a/OpenSim/Region/Terrain.BasicTerrain/libTerrainBSD/Channel/Manipulators/NavierStokes.cs b/OpenSim/Region/Terrain.BasicTerrain/libTerrainBSD/Channel/Manipulators/NavierStokes.cs
index f01b2860..43987d4 100644
--- a/OpenSim/Region/Terrain.BasicTerrain/libTerrainBSD/Channel/Manipulators/NavierStokes.cs
+++ b/OpenSim/Region/Terrain.BasicTerrain/libTerrainBSD/Channel/Manipulators/NavierStokes.cs
@@ -248,9 +248,9 @@ namespace libTerrain
         /// Performs computational fluid dynamics on a channel
         /// </summary>
         /// <param name="rounds">The number of steps to perform (Recommended: 20)</param>
-        /// <param name="dt">Delta Time - The time between steps</param>
-        /// <param name="diff">Fluid diffusion rate</param>
-        /// <param name="visc">Fluid viscosity</param>
+        /// <param name="dt">Delta Time - The time between steps (Recommended: 0.1)</param>
+        /// <param name="diff">Fluid diffusion rate (Recommended: 0.0)</param>
+        /// <param name="visc">Fluid viscosity (Recommended: 0.0)</param>
         public void navierStokes(int rounds, double dt, double diff, double visc)
         {
             nsSimulate(this.h, rounds, dt, diff, visc);
-- 
cgit v1.1