调试读取NC文件

Models/SimulationCommon/SeaSJ.cs

@@ -52,15 +52,15 @@ public class SeaSJ
     }
 
     // 交付类方法
-    public static List<double[]> GetDrift(NCread nCread, double[] initialPosition, double dt, double totalTime, int times, int latitudes, int longitudes, int days, int hour, double initlatitudes, double initlongitudes)
+    public static List<double[]> GetDrift(NCread nCread, double[] initialPosition, double dt, double totalTime, int times, int latitudes, int longitudes, int times1, int latitudes1, int longitudes1, int times2, int latitudes2, int longitudes2, int times3, int latitudes3, int longitudes3, int days, int hour, double initlatitudes, double initlongitudes)
     {
         // run 获取轨迹的函数
         // initialPosition --> 初始位置; dt --> 时间步长; totalTime --> 总时长
-        List<double[]> trajectory = CalculateDriftTrajectory(nCread, initialPosition, dt, totalTime, times, latitudes, longitudes, days, hour, initlatitudes, initlongitudes);
+        List<double[]> trajectory = CalculateDriftTrajectory(nCread, initialPosition, dt, totalTime, times, latitudes, longitudes, times1, latitudes1, longitudes1, times2, latitudes2, longitudes2, times3, latitudes3, longitudes3, days, hour, initlatitudes, initlongitudes);
         return trajectory;
     }
 
-    public static List<double[]> CalculateDriftTrajectory(NCread nCread, double[] initialPosition, double dt, double totalTime, int times, int latitudes, int longitudes, int days, int hour, double initlatitudes, double initlongitudes)
+    public static List<double[]> CalculateDriftTrajectory(NCread nCread, double[] initialPosition, double dt, double totalTime, int times, int latitudes, int longitudes, int times1, int latitudes1, int longitudes1, int times2, int latitudes2, int longitudes2, int times3, int latitudes3, int longitudes3, int days, int hour, double initlatitudes, double initlongitudes)
     {
         int timeSteps = (int)(totalTime / dt);
         List<double[]> trajectory = new List<double[]>();
@@ -72,20 +72,36 @@ public class SeaSJ
 
             int outputTime = 3600 * t;
             // 动态获取当前位置的风力和洋流数据
-            double[] windVelocity = GetWindVelocityFromAPI(nCread, currentPos[0], currentPos[1], outputTime, times, latitudes, longitudes, days, hour);
-            double[] currentVelocity = GetCurrentVelocityFromAPI(nCread, currentPos[0], currentPos[1], outputTime, times, latitudes, longitudes, days, hour);
-
+            double[] windVelocity = GetWindVelocityFromAPI(nCread, currentPos[0], currentPos[1], outputTime, times, latitudes, longitudes, times1, latitudes1, longitudes1, days, hour);
+            double[] currentVelocity = GetCurrentVelocityFromAPI(nCread, currentPos[0], currentPos[1], outputTime, times2, latitudes2, longitudes2, times3, latitudes3, longitudes3, days, hour);
+            //foreach (var item in windVelocity)
+            //{
+            //    Console.WriteLine("windVelocity:" + item);
+            //}
+            //foreach (var item in currentVelocity)
+            //{
+            //    Console.WriteLine("currentVelocity:" + item);
+            //}
             // 计算漂移速度(m/s)
             double[] driftVelocity = {
                 currentVelocity[0] + windVelocity[0],
                 currentVelocity[1] + windVelocity[1]
             };
-
+            //foreach (var item in driftVelocity)
+            //{
+            //    Console.WriteLine("driftVelocity:" + item);
+            //}
             // 更新位置（漂移速度单位是m/s，需要转换成经纬度的变化量）
             double[] newPosition = {
                 currentPos[0] + (driftVelocity[0] * dt * 3600) / ((6371000 * Math.Cos(currentPos[0] * Math.PI / 180)) * 180 / Math.PI), // 纬度变化量
                 currentPos[1] + (driftVelocity[1] * dt * 3600) / ((6371000 * Math.Cos(currentPos[1] * Math.PI / 180)) * 180 / Math.PI)  // 经度变化量
             };
+            //foreach (var item in newPosition)
+            //{
+            //    Console.WriteLine("newPosition:" + item);
+            //}
+            //Console.WriteLine("currentPos[0]:" + (driftVelocity[0] * dt * 3600) / ((6371000 * Math.Cos(currentPos[0] * Math.PI / 180)) * 180 / Math.PI));
+            //Console.WriteLine("currentPos[1]:" + (driftVelocity[1] * dt * 3600) / ((6371000 * Math.Cos(currentPos[1] * Math.PI / 180)) * 180 / Math.PI));
 
             trajectory.Add(newPosition);
         }
@@ -93,7 +109,7 @@ public class SeaSJ
         return trajectory;//（纬度，经度）
     }
 
-    public static double[] GetWindVelocityFromAPI(NCread windNCread, double latitude, double longitude, double temptime, int times, int latitudes, int longitudes, int days, int hour)
+    public static double[] GetWindVelocityFromAPI(NCread windNCread, double latitude, double longitude, double temptime, int times, int latitudes, int longitudes, int times1, int latitudes1, int longitudes1, int days, int hour)
     {
         float[] longitudeArray = windNCread.longitudeArray;//经度一维数组来源自文件'Text_readNC.cs'
         float[] latitudeArray = windNCread.latitudeArray;//纬度一维数组来源自文件'Text_readNC.cs'
@@ -141,13 +157,23 @@ public class SeaSJ
         return new double[] { windX, windY };
     }
 
-    public static double[] GetCurrentVelocityFromAPI(NCread CurrentNCread, double latitude, double longitude, double temptime, int times, int latitudes, int longitudes, int days, int hour)
+    public static double[] GetCurrentVelocityFromAPI(NCread CurrentNCread, double latitude, double longitude, double temptime, int times2, int latitudes2, int longitudes2, int times3, int latitudes3, int longitudes3, int days, int hour)
     {
         float[] longitudeArray = CurrentNCread.longitudeArray;//经度一维数组来源自文件'Text_readNC.cs'
         float[] latitudeArray = CurrentNCread.latitudeArray;//纬度一维数组来源自文件'Text_readNC.cs'
         float[][][] u10Array = CurrentNCread.u10Array;//风的10米U（向东）分量三维数组来源自文件'Text_readNC.cs'
         float[][][] v10Array = CurrentNCread.v10Array;//风的10米V（向北）分量三维数组来源自文件'Text_readNC.cs'
         float[][][] p140208Array = CurrentNCread.p140208Array;//海洋上空的自由对流速度三维数组来源自文件'Text_readNC.cs'
+        //foreach (var item in p140208Array)
+        //{
+        //    foreach (var item1 in item)
+        //    {
+        //        foreach (var item2 in item1)
+        //        {
+        //            Console.WriteLine("p140208Array:" + item2);
+        //        }
+        //    }
+        //}
         float[][][] mwdArray = CurrentNCread.mwdArray;//平均波向（单位：度；0度表示北方，90度表示东方）三维数组来源自文件'Text_readNC.cs'
 
         int longitudeNum = 0;
@@ -155,7 +181,7 @@ public class SeaSJ
         int temptimeNum = 0;
         //定义NC文件中读取不到的坐标点海洋数据信息（后续可采用插值法等）
         //经度连续化
-        for (int i = 0; i < longitudes - 1; i++)
+        for (int i = 0; i < longitudes2 - 1; i++)
         {
             if (longitude >= longitudeArray[i] && longitude < longitudeArray[i + 1])
             {
@@ -165,7 +191,7 @@ public class SeaSJ
         }
 
         //纬度连续化
-        for (int i = 0; i < latitudes - 1; i++)
+        for (int i = 0; i < latitudes2 - 1; i++)
         {
             if (latitude >= latitudeArray[i] && latitude < latitudeArray[i + 1])
             {
@@ -175,7 +201,7 @@ public class SeaSJ
         }
 
         //时间连续化
-        for (int i = 0; i < times - 1; i++)
+        for (int i = 0; i < times2 - 1; i++)
         {
             if (temptime >= 3600 * i && temptime < 3600 * (i + 1))
             {
@@ -183,9 +209,20 @@ public class SeaSJ
             }
         }
         temptimeNum += (days - 1) * 24 + hour;
+        //foreach (var item in p140208Array)
+        //{
+        //    Console.WriteLine("p140208Array:" + (double)p140208Array[0][1][2]);
+        //}
+        //Console.WriteLine("p140208Array:" + (double)p140208Array[2640][1][2]);
         double currentSpeed = (double)p140208Array[temptimeNum][latitudeNum][longitudeNum];
+        Random random = new Random();
+        if (currentSpeed > 10)
+            currentSpeed = random.NextDouble();
+        //Console.WriteLine("temptimeNum:" + temptimeNum + "_" + "latitudeNum:" + latitudeNum + "_" + "longitudeNum:" + longitudeNum);
         double currentDirection = (double)mwdArray[temptimeNum][latitudeNum][longitudeNum];
-
+        //double currentSpeed = (double)p140208Array[2401][0][0];
+        //Console.WriteLine("currentSpeed:" + currentSpeed);
+        //Console.WriteLine("currentDirection:" + currentDirection);
         double currentDirectionInRadians = currentDirection * (Math.PI / 180);
         double currentX = currentSpeed * Math.Cos(currentDirectionInRadians);
         double currentY = currentSpeed * Math.Sin(currentDirectionInRadians);
@@ -194,7 +231,7 @@ public class SeaSJ
     }
 
     //海浪高度获取
-    public static double GetWaveHeightFromAPI(NCread WaveNCread, double latitude, double longitude, double temptime, int times, int latitudes, int longitudes, int days, int hour)//temptime:仿真时间（秒）
+    public static double GetWaveHeightFromAPI(NCread WaveNCread, double latitude, double longitude, double temptime, int times4, int latitudes4, int longitudes4, int days, int hour)//temptime:仿真时间（秒）
     {
         float[] longitudeArray = WaveNCread.longitudeArray;//经度一维数组来源自文件'Text_readNC.cs'
         float[] latitudeArray = WaveNCread.latitudeArray;//纬度一维数组来源自文件'Text_readNC.cs'
@@ -206,7 +243,7 @@ public class SeaSJ
 
         //定义NC文件中读取不到的坐标点海洋数据信息（后续可采用插值法等）
         //经度连续化
-        for (int i = 0; i < longitudes - 1; i++)
+        for (int i = 0; i < longitudes4 - 1; i++)
         {
             if (longitude >= longitudeArray[i] && longitude < longitudeArray[i + 1])
             {
@@ -216,7 +253,7 @@ public class SeaSJ
         }
 
         //纬度连续化
-        for (int i = 0; i < latitudes - 1; i++)
+        for (int i = 0; i < latitudes4 - 1; i++)
         {
             if (latitude >= latitudeArray[i] && latitude < latitudeArray[i + 1])
             {
@@ -226,7 +263,7 @@ public class SeaSJ
         }
 
         //时间连续化
-        for (int i = 0; i < times - 1; i++)
+        for (int i = 0; i < times4 - 1; i++)
         {
             if (temptime >= 3600 * i && temptime < (3600 * (i + 1)))
             {

Models/SimulationCommon/Text_readNC.cs

@@ -9,18 +9,19 @@ public class Text_readNC
     public int times;
     public int latitudes;
     public int longitudes;
-    //public int timeSteps2;
-    //public int latitudes2;
-    //public int longitudes2;
-    //public int timeSteps3;
-    //public int latitudes3;
-    //public int longitudes3;
-    //public int timeSteps4;
-    //public int latitudes4;
-    //public int longitudes4;
-    //public int timeSteps1;
-    //public int latitudes1;
-    //public int longitudes1;
+    public int times1;
+    public int latitudes1;
+    public int longitudes1;
+    public int times2;
+    public int latitudes2;
+    public int longitudes2;
+    public int times3;
+    public int latitudes3;
+    public int longitudes3;
+    public int times4;
+    public int latitudes4;
+    public int longitudes4;
+
     public double initlatitudes;
     public double initlongitudes;
 
@@ -39,7 +40,7 @@ public class Text_readNC
             using (StreamReader reader1 = new StreamReader(response1.GetResponseStream()))
             {
                 string result1 = reader1.ReadToEnd();
-
+                //Console.WriteLine("result1:" + result1);
                 // Parse the JSON string into a JObject
                 JObject resultObject1 = JObject.Parse(result1);
 
@@ -125,23 +126,23 @@ public class Text_readNC
                 JArray dataArray2 = (JArray)resultObject2["data"];
 
                 // Determine the dimensions of the data array
-                times = dataArray2.Count;
-                latitudes = ((JArray)dataArray2[0]).Count;
-                longitudes = ((JArray)dataArray2[0][0]).Count;
+                times1 = dataArray2.Count;
+                latitudes1 = ((JArray)dataArray2[0]).Count;
+                longitudes1 = ((JArray)dataArray2[0][0]).Count;
 
                 // Initialize float[][][] array
-                float[][][] u10Array = new float[times][][];
+                float[][][] u10Array = new float[times1][][];
 
                 // Convert each element to float and populate v10Array
-                for (int t = 0; t < times; t++)
+                for (int t = 0; t < times1; t++)
                 {
-                    u10Array[t] = new float[latitudes][];
+                    u10Array[t] = new float[latitudes1][];
 
-                    for (int lat = 0; lat < latitudes; lat++)
+                    for (int lat = 0; lat < latitudes1; lat++)
                     {
-                        u10Array[t][lat] = new float[longitudes];
+                        u10Array[t][lat] = new float[longitudes1];
 
-                        for (int lon = 0; lon < longitudes; lon++)
+                        for (int lon = 0; lon < longitudes1; lon++)
                         {
                             u10Array[t][lat][lon] = (float)dataArray2[t][lat][lon];
                         }
@@ -174,23 +175,23 @@ public class Text_readNC
                 JArray dataArray3 = (JArray)resultObject3["data"];
 
                 // Determine the dimensions of the data array
-                times = dataArray3.Count;
-                latitudes = ((JArray)dataArray3[0]).Count;
-                longitudes = ((JArray)dataArray3[0][0]).Count;
+                times2 = dataArray3.Count;
+                latitudes2 = ((JArray)dataArray3[0]).Count;
+                longitudes2 = ((JArray)dataArray3[0][0]).Count;
 
                 // Initialize float[][][] array
-                float[][][] p140208Array = new float[times][][];
+                float[][][] p140208Array = new float[times2][][];
 
                 // Convert each element to float and populate v10Array
-                for (int t = 0; t < times; t++)
+                for (int t = 0; t < times2; t++)
                 {
-                    p140208Array[t] = new float[latitudes][];
+                    p140208Array[t] = new float[latitudes2][];
 
-                    for (int lat = 0; lat < latitudes; lat++)
+                    for (int lat = 0; lat < latitudes2; lat++)
                     {
-                        p140208Array[t][lat] = new float[longitudes];
+                        p140208Array[t][lat] = new float[longitudes2];
 
-                        for (int lon = 0; lon < longitudes; lon++)
+                        for (int lon = 0; lon < longitudes2; lon++)
                         {
                             p140208Array[t][lat][lon] = (float)dataArray3[t][lat][lon];
                         }
@@ -221,23 +222,23 @@ public class Text_readNC
                 JArray dataArray4 = (JArray)resultObject4["data"];
 
                 // Determine the dimensions of the data array
-                times = dataArray4.Count;
-                latitudes = ((JArray)dataArray4[0]).Count;
-                longitudes = ((JArray)dataArray4[0][0]).Count;
+                times3 = dataArray4.Count;
+                latitudes3 = ((JArray)dataArray4[0]).Count;
+                longitudes3 = ((JArray)dataArray4[0][0]).Count;
 
                 // Initialize float[][][] array
-                float[][][] mwdArray = new float[times][][];
+                float[][][] mwdArray = new float[times3][][];
 
                 // Convert each element to float and populate v10Array
-                for (int t = 0; t < times; t++)
+                for (int t = 0; t < times3; t++)
                 {
-                    mwdArray[t] = new float[latitudes][];
+                    mwdArray[t] = new float[latitudes3][];
 
-                    for (int lat = 0; lat < latitudes; lat++)
+                    for (int lat = 0; lat < latitudes3; lat++)
                     {
-                        mwdArray[t][lat] = new float[longitudes];
+                        mwdArray[t][lat] = new float[longitudes3];
 
-                        for (int lon = 0; lon < longitudes; lon++)
+                        for (int lon = 0; lon < longitudes3; lon++)
                         {
                             mwdArray[t][lat][lon] = (float)dataArray4[t][lat][lon];
                         }
@@ -272,22 +273,23 @@ public class Text_readNC
                 JArray dataArray1 = (JArray)resultObject1["data"];
 
                 // Determine the dimensions of the data array
-                times = dataArray1.Count;
-                latitudes = ((JArray)dataArray1[0]).Count;
-                longitudes = ((JArray)dataArray1[0][0]).Count;
+                times4 = dataArray1.Count;
+                latitudes4 = ((JArray)dataArray1[0]).Count;
+                Console.WriteLine("latitudes2:" + latitudes);
+                longitudes4 = ((JArray)dataArray1[0][0]).Count;
 
                 // Initialize float[][][] array
-                float[][][] hmaxArray = new float[times][][];
+                float[][][] hmaxArray = new float[times4][][];
 
-                for (int t = 0; t < times; t++)
+                for (int t = 0; t < times4; t++)
                 {
-                    hmaxArray[t] = new float[latitudes][];
+                    hmaxArray[t] = new float[latitudes4][];
 
-                    for (int lat = 0; lat < latitudes; lat++)
+                    for (int lat = 0; lat < latitudes4; lat++)
                     {
-                        hmaxArray[t][lat] = new float[longitudes];
+                        hmaxArray[t][lat] = new float[longitudes4];
 
-                        for (int lon = 0; lon < longitudes; lon++)
+                        for (int lon = 0; lon < longitudes4; lon++)
                         {
                             hmaxArray[t][lat][lon] = (float)dataArray1[t][lat][lon];
                         }

Models/SimulationCommon/temp_readNC.cs

@@ -37,22 +37,22 @@ public class GetNCData
                 JArray dataArray1 = (JArray)resultObject1["data"];
 
                 // Determine the dimensions of the data array
-                int timeSteps = dataArray1.Count;
-                int latitudes = ((JArray)dataArray1[0]).Count;
-                int longitudes = ((JArray)dataArray1[0][0]).Count;
+                int times5 = dataArray1.Count;
+                int latitudes5 = ((JArray)dataArray1[0]).Count;
+                int longitudes5 = ((JArray)dataArray1[0][0]).Count;
 
                 // Initialize float[][][] array
-                float[][][] tempArray = new float[timeSteps][][];
+                float[][][] tempArray = new float[times5][][];
 
-                for (int t = 0; t < timeSteps; t++)
+                for (int t = 0; t < times5; t++)
                 {
-                    tempArray[t] = new float[latitudes][];
+                    tempArray[t] = new float[latitudes5][];
 
-                    for (int lat = 0; lat < latitudes; lat++)
+                    for (int lat = 0; lat < latitudes5; lat++)
                     {
-                        tempArray[t][lat] = new float[longitudes];
+                        tempArray[t][lat] = new float[longitudes5];
 
-                        for (int lon = 0; lon < longitudes; lon++)
+                        for (int lon = 0; lon < longitudes5; lon++)
                         {
                             tempArray[t][lat][lon] = (float)dataArray1[t][lat][lon] - 273.15f;
                         }

SimulationServer/Entity/AircraftSJ.cs

@@ -76,7 +76,7 @@ public class AircraftSJ : AircraftEntity
         var nCread = text_ReadNC.windNCread;
 
         //漂移轨迹
-        List<double[]> trajectory = SeaSJ.CalculateDriftTrajectory(nCread, initialPosition, dt, totalTime, text_ReadNC.times, text_ReadNC.latitudes, text_ReadNC.longitudes, Days, Hour, text_ReadNC.initlatitudes, text_ReadNC.initlongitudes);
+        List<double[]> trajectory = SeaSJ.CalculateDriftTrajectory(nCread, initialPosition, dt, totalTime, text_ReadNC.times, text_ReadNC.latitudes, text_ReadNC.longitudes, text_ReadNC.times1, text_ReadNC.latitudes1, text_ReadNC.longitudes1, text_ReadNC.times2, text_ReadNC.latitudes2, text_ReadNC.longitudes2, text_ReadNC.times3, text_ReadNC.latitudes3, text_ReadNC.longitudes3, Days, Hour, text_ReadNC.initlatitudes, text_ReadNC.initlongitudes);
 
         // 生成任务终点
         MissionEndPoint = new MissionEndPoint
@@ -253,13 +253,13 @@ public class AircraftSJ : AircraftEntity
 
 
                 var wind = SeaSJ.GetWindVelocityFromAPI(nCread, currentLocation.CurrentLat, currentLocation.CurrentLon,
-                    temptime, text_ReadNC.times, text_ReadNC.latitudes, text_ReadNC.longitudes, Days, Hour);
+                    temptime, text_ReadNC.times, text_ReadNC.latitudes, text_ReadNC.longitudes, text_ReadNC.times1, text_ReadNC.latitudes1, text_ReadNC.longitudes1, Days, Hour);
 
 
                 var windSpeed = Math.Sqrt(wind[0] * wind[0] + wind[1] * wind[1]);
 
                 var waveHigh = SeaSJ.GetWaveHeightFromAPI(nCread, currentLocation.CurrentLon,
-                    currentLocation.CurrentLat, temptime, text_ReadNC.times, text_ReadNC.latitudes, text_ReadNC.longitudes, Days, Hour);
+                    currentLocation.CurrentLat, temptime, text_ReadNC.times4, text_ReadNC.latitudes4, text_ReadNC.longitudes4, Days, Hour);
 
                 var distance = Utils.Util.GetDistance(currentLocation.CurrentLon, targetPoint.x,
                     currentLocation.CurrentLat,