使用MeteoInfoC.dll如何创建一个和已有结构一样的nc文件

看天吃饭

我有一个nc文件，目前需要做两个工作：
1）将原有的含60年资料的文件按月取平均值（此步骤已在楼主的指导下完成）。
2）要将计算好的平均值数据(12个月的数据最终得到了12个GridData数据)重新写一个nc文件，要求和原文件结构一样。

1. 
   
2. NetCDFDataInfo gridFile = new NetCDFDataInfo();
   
3.             gridFile.CreateNCFile(Application.StartupPath + "\\Sample\\result.nc");
   
4. 
   
5. //... 如何将GridData数据写入到相应的变量中？
   

<br /> NetCDFDataInfo gridFile = new NetCDFDataInfo();<br />             gridFile.CreateNCFile(Application.StartupPath + "<a href=\"http://bbs.06climate.com/file://\\Sample\\result.nc\" target=\"_blank\">\\Sample\\result.nc</a>");<br /> <br /> //... 如何将GridData数据写入到相应的变量中？<br />

mofangbao

貌似有python下脚本的例子，看看这个有帮助没：
http://bbs.06climate.com/forum.php?mod=viewthread&tid=5485

看天吃饭

十分感谢，这个贴子很有作用。做完了我把源码献上。
还有一个帖子也很有用：
http://hi.baidu.com/meteoinfo/item/369f9a2279ed661e2a0f1ce3

看天吃饭

请高人指点，目前我想要的效果已经基本可以实现了，即：
1）求历年12个月数据各个月的平均值，求完放在 mdatas数组中(mdatas 为 GridData[12])。
2）创建result.nc文件，希望将计算结果（即历年12个月数据各个月的平均值mdatas）保存到新的文件中。nc的格式需要同原来的数据格式相同。

注：目前我将每个月的数据都存为一个nc文件没有问题。可是不知道如何将12个月的数据都存在一个文件中。
先贴代码再解释我的做法。

1. 
   
2. NetCDFDataInfo aDataInfo = new NetCDFDataInfo();
   
3.             string aFile = Application.StartupPath + "\\Sample\\Hadisst_obs_1949-2010_T106_BCC_AGCM2.1.nc";
   
4.             //读取nc文件
   
5.             aDataInfo.ReadDataInfo(aFile);
   
6.             var t = aDataInfo.times;
   
7.             List<GridData>[] datas = new List<GridData>[12];
   
8.             int sYear = aDataInfo.times[0].Year;
   
9.             int eYear = aDataInfo.times[aDataInfo.times.Count - 1].Year;
   
10.             int yIndex = 0;
    
11.             GridData gData = new GridData();
    
12.             //循环年份
    
13.             for (int startYear = 1961; startYear < eYear; startYear++)
    
14.             {
    
15.                 yIndex = startYear - sYear;
    
16.                 //将每年12个月的数据都摆在相应月份的数组里
    
17.                 for (int mIndex = 0; mIndex < 12; mIndex++)
    
18.                 {
    
19.                     List<GridData> mList = datas[mIndex];
    
20.                     if (mList == null)
    
21.                     {
    
22.                         mList = new List<GridData>();
    
23.                     }
    
24.                     mList.Add(aDataInfo.GetGridData_LonLat(yIndex * 12 + mIndex, 4, 0));
    
25.                     datas[mIndex] = mList;
    
26.                 }
    
27.             }
    
28.             GridData[] mdatas = new GridData[12];
    
29.             //对每个月份的数据计算平均值
    
30.             for (int mIndex = 0; mIndex < 12; mIndex++)
    
31.             {
    
32.                 List<GridData> mList = datas[mIndex];
    
33.                 mdatas[mIndex] = DataMath.Average(mList, false);
    
34.             }
    
35.             NetCDFDataInfo gridFile = new NetCDFDataInfo();
    
36.             string outFilePath = Application.StartupPath + "\\Sample\\result.nc";
    
37.             gridFile.fileName = outFilePath;
    
38.             gridFile.UNDEF = aDataInfo.UNDEF;
    
39.             gridFile.unlimdimid = aDataInfo.unlimdimid;
    
40.             // Add dimensions:
    
41.             foreach (var dimension in aDataInfo.dimList)
    
42.             {
    
43.                 gridFile.AddDimension(dimension);
    
44.             }
    
45.             object[] lonDatas = new object[160];
    
46.             object[] latDatas = new object[320];
    
47.             // Add variables
    
48.             foreach (Variable variable in aDataInfo.varList)
    
49.             {
    
50.                 gridFile.AddVariable((Variable)variable.Clone());
    
51.                 //Add variable attributes
    
52.                 foreach (var attribute in variable.Attributes)
    
53.                 {
    
54.                     gridFile.AddVariableAttribute(variable.Name, attribute.attName, attribute.attValue.ToString());
    
55.                 }
    
56.                 if (variable.Name == "lon")
    
57.                 {
    
58.                     aDataInfo.GetVarData(variable, ref lonDatas);
    
59.                 }
    
60.                 if (variable.Name == "lat")
    
61.                 {
    
62.                     aDataInfo.GetVarData(variable, ref latDatas);
    
63.                 }
    
64.             }
    
65.             //Add global attributes
    
66.             foreach (var globalAttribute in aDataInfo.gAttList)
    
67.             {
    
68.                 gridFile.AddGlobalAttribute(globalAttribute.attName, globalAttribute.attValue.ToString());
    
69.             }
    
70.             //Create netCDF file
    
71.             gridFile.CreateNCFile(outFilePath);
    
72.             gridFile.WriteVar("lon", lonDatas);
    
73.             gridFile.WriteVar("lat", latDatas);
    
74.             //Add Time
    
75.             int[] startArray = new int[1];
    
76.             int[] countArray = new int[1];
    
77.             double[] dataArray = new double[12];
    
78.             for (int mIndex = 0; mIndex < 12; mIndex++)
    
79.             {
    
80.                 dataArray[mIndex] = mIndex;
    
81.             }
    
82.             startArray[0] = 0;
    
83.             countArray[0] = 1;
    
84.             gridFile.WriteVara("time", startArray, countArray, dataArray);
    
85. 
    
86.             //Add Data
    
87.             int[] startdataArray = new int[3];
    
88.             int[] countdataArray = new int[3];
    
89.             startdataArray[0] = 0;
    
90.             startdataArray[1] = 0;
    
91.             startdataArray[2] = 0;
    
92.             countdataArray[0] = 12;
    
93.             countdataArray[1] = aDataInfo.Y.Length;
    
94.             countdataArray[2] = aDataInfo.X.Length;
    
95.             
    
96.             //创建一个可以盛放12个月数据的数组（从这个地方开始就不知道对不对了）
    
97.             double[] dataResult = new double[aDataInfo.X.Length * aDataInfo.Y.Length * 12];
    
98.             int dataIndex = 0;
    
99.             //循环12个月的数据
    
100.             for (int mIndex = 0; mIndex < 12; mIndex++)
     
101.             {
     
102.                 //获取每个月的数据
     
103.                 var ddt = mdatas[mIndex].Data;
     
104.                 //将每个月的数据转换为数组赋值
     
105.                 for (int y = 0; y < aDataInfo.Y.Length; y++)
     
106.                 {
     
107.                     for (int x = 0; x < aDataInfo.X.Length; x++)
     
108.                     {
     
109.                         dataResult[dataIndex] = ddt[y, x];
     
110.                         dataIndex++;
     
111.                     }
     
112.                 }
     
113.             }
     
114.             // gridFile.w
     
115.             gridFile.WriteVara("SST_cpl", startdataArray, countdataArray, dataResult);
     
116.             gridFile.CloseNCFile();
     

<br /> NetCDFDataInfo aDataInfo = new NetCDFDataInfo();<br />             string aFile = Application.StartupPath + "<a href=\"http://bbs.06climate.com/file://\\Sample\\Hadisst_obs_1949-2010_T106_BCC_AGCM2.1.nc\" target=\"_blank\">\\Sample\\Hadisst_obs_1949-2010_T106_BCC_AGCM2.1.nc</a>";<br />             //读取nc文件<br />             aDataInfo.ReadDataInfo(aFile);<br />             var t = aDataInfo.times;<br />             List<GridData>[] datas = new List<GridData>[12];<br />             int sYear = aDataInfo.times[0].Year;<br />             int eYear = aDataInfo.times[aDataInfo.times.Count - 1].Year;<br />             int yIndex = 0;<br />             GridData gData = new GridData();<br />             //循环年份<br />             for (int startYear = 1961; startYear < eYear; startYear++)<br />             {<br />                 yIndex = startYear - sYear;<br />                 //将每年12个月的数据都摆在相应月份的数组里<br />                 for (int mIndex = 0; mIndex < 12; mIndex++)<br />                 {<br />                     List<GridData> mList = datas[mIndex];<br />                     if (mList == null)<br />                     {<br />                         mList = new List<GridData>();<br />                     }<br />                     mList.Add(aDataInfo.GetGridData_LonLat(yIndex * 12 + mIndex, 4, 0));<br />                     datas[mIndex] = mList;<br />                 }<br />             }<br />             GridData[] mdatas = new GridData[12];<br />             //对每个月份的数据计算平均值<br />             for (int mIndex = 0; mIndex < 12; mIndex++)<br />             {<br />                 List<GridData> mList = datas[mIndex];<br />                 mdatas[mIndex] = DataMath.Average(mList, false);<br />             }<br />             NetCDFDataInfo gridFile = new NetCDFDataInfo();<br />             string outFilePath = Application.StartupPath + "<a href=\"http://bbs.06climate.com/file://\\Sample\\result.nc\" target=\"_blank\">\\Sample\\result.nc</a>";<br />             gridFile.fileName = outFilePath;<br />             gridFile.UNDEF = aDataInfo.UNDEF;<br />             gridFile.unlimdimid = aDataInfo.unlimdimid;<br />             // Add dimensions:<br />             foreach (var dimension in aDataInfo.dimList)<br />             {<br />                 gridFile.AddDimension(dimension);<br />             }<br />             object[] lonDatas = new object[160];<br />             object[] latDatas = new object[320];<br />             // Add variables<br />             foreach (Variable variable in aDataInfo.varList)<br />             {<br />                 gridFile.AddVariable((Variable)variable.Clone());<br />                 //Add variable attributes<br />                 foreach (var attribute in variable.Attributes)<br />                 {<br />                     gridFile.AddVariableAttribute(variable.Name, attribute.attName, attribute.attValue.ToString());<br />                 }<br />                 if (variable.Name == "lon")<br />                 {<br />                     aDataInfo.GetVarData(variable, ref lonDatas);<br />                 }<br />                 if (variable.Name == "lat")<br />                 {<br />                     aDataInfo.GetVarData(variable, ref latDatas);<br />                 }<br />             }<br />             //Add global attributes<br />             foreach (var globalAttribute in aDataInfo.gAttList)<br />             {<br />                 gridFile.AddGlobalAttribute(globalAttribute.attName, globalAttribute.attValue.ToString());<br />             }<br />             //Create netCDF file<br />             gridFile.CreateNCFile(outFilePath);<br />             gridFile.WriteVar("lon", lonDatas);<br />             gridFile.WriteVar("lat", latDatas);<br />             //Add Time<br />             int[] startArray = new int[1];<br />             int[] countArray = new int[1];<br />             double[] dataArray = new double[12];<br />             for (int mIndex = 0; mIndex < 12; mIndex++)<br />             {<br />                 dataArray[mIndex] = mIndex;<br />             }<br />             startArray[0] = 0;<br />             countArray[0] = 1;<br />             gridFile.WriteVara("time", startArray, countArray, dataArray);<br /> <br />             //Add Data<br />             int[] startdataArray = new int[3];<br />             int[] countdataArray = new int[3];<br />             startdataArray[0] = 0;<br />             startdataArray[1] = 0;<br />             startdataArray[2] = 0;<br />             countdataArray[0] = 12;<br />             countdataArray[1] = aDataInfo.Y.Length;<br />             countdataArray[2] = aDataInfo.X.Length;<br />             <br />             //创建一个可以盛放12个月数据的数组（从这个地方开始就不知道对不对了）<br />             double[] dataResult = new double[aDataInfo.X.Length * aDataInfo.Y.Length * 12];<br />             int dataIndex = 0;<br />             //循环12个月的数据<br />             for (int mIndex = 0; mIndex < 12; mIndex++)<br />             {<br />                 //获取每个月的数据<br />                 var ddt = mdatas[mIndex].Data;<br />                 //将每个月的数据转换为数组赋值<br />                 for (int y = 0; y < aDataInfo.Y.Length; y++)<br />                 {<br />                     for (int x = 0; x < aDataInfo.X.Length; x++)<br />                     {<br />                         dataResult[dataIndex] = ddt[y, x];<br />                         dataIndex++;<br />                     }<br />                 }<br />             }<br />             // gridFile.w<br />             gridFile.WriteVara("SST_cpl", startdataArray, countdataArray, dataResult);<br />             gridFile.CloseNCFile();<br />

各位看官请看源码的最下面，我开始测试，如果我这么写：

1. 
   
2. 
   
3. //Add Data
   
4. int[] startdataArray = new int[3];
   
5. int[] countdataArray = new int[3];
   
6. startdataArray[0] = 0;
   
7. startdataArray[1] = 0;
   
8. startdataArray[2] = 0;
   
9. countdataArray[0] = 12;
   
10. countdataArray[1] = aDataInfo.Y.Length;
    
11. countdataArray[2] = aDataInfo.X.Length;
    
12. 
    
13.   //创建一个可以盛放数据的数组
    
14. double[] dataResult = new double[aDataInfo.X.Length * aDataInfo.Y.Length ];
    
15. int dataIndex = 0;
    
16. //获取1月份的数据
    
17. var ddt = mdatas[0].Data;
    
18. //将每个月的数据转换为数组赋值
    
19. for (int y = 0; y < aDataInfo.Y.Length; y++)
    
20. {
    
21. for (int x = 0; x < aDataInfo.X.Length; x++)
    
22. {
    
23. dataResult[dataIndex] = ddt[y, x];
    
24. dataIndex++;
    
25. }
    
26. }
    
27. // 写变量数据
    
28. gridFile.WriteVara("SST_cpl", startdataArray, countdataArray, dataResult);
    
29. 
    

<br /> <br /> //Add Data<br /> int[] startdataArray = new int[3];<br /> int[] countdataArray = new int[3];<br /> startdataArray[0] = 0;<br /> startdataArray[1] = 0;<br /> startdataArray[2] = 0;<br /> countdataArray[0] = 12;<br /> countdataArray[1] = aDataInfo.Y.Length;<br /> countdataArray[2] = aDataInfo.X.Length;<br /> <br />   //创建一个可以盛放数据的数组<br /> double[] dataResult = new double[aDataInfo.X.Length * aDataInfo.Y.Length ];<br /> int dataIndex = 0;<br /> //获取1月份的数据<br /> var ddt = mdatas[0].Data;<br /> //将每个月的数据转换为数组赋值<br /> for (int y = 0; y < aDataInfo.Y.Length; y++)<br /> {<br /> for (int x = 0; x < aDataInfo.X.Length; x++)<br /> {<br /> dataResult[dataIndex] = ddt[y, x];<br /> dataIndex++;<br /> }<br /> }<br /> // 写变量数据<br /> gridFile.WriteVara("SST_cpl", startdataArray, countdataArray, dataResult);<br /> <br />

目前最大的问题就是多个时间尺度的数据如何保存。
startdataArray, countdataArray, dataResult三个变量都是怎么设置的？
说实话我都是依葫芦画瓢，希望王老师或其他高手能再给点指导。

看天吃饭

很急，在线等。。。各位大侠出来帮忙看看

MeteoInfo

看天吃饭 发表于 2013-1-30 13:34

                               
登录/注册后可看大图

请高人指点，目前我想要的效果已经基本可以实现了，即：
1）求历年12个月数据各个月的平均值，求完放在 md ...

可以参考此帖：MeteoInfo脚本示例：MICAPS第4类数据转为NetCDF
http://bbs.06climate.com/forum.p ... 630&fromuid=106

看天吃饭

MeteoInfo 发表于 2013-1-30 14:50

                               
登录/注册后可看大图

可以参考此帖：MeteoInfo脚本示例：MICAPS第4类数据转为NetCDF
http://bbs.06climate.com/forum.php?mod ...

王老师，按照那个方式，我调整了一下关键代码，好像现在写的数据不对了。

1. 
   
2. //创建一个可以盛放12个月数据的数组（从这个地方开始就不知道对不对了）
   
3.             double[] dataResult = new double[aDataInfo.X.Length * aDataInfo.Y.Length];
   
4. 
   
5. 
   
6.             //Add Data
   
7.             int[] startdataArray = new int[4];
   
8.             int[] countdataArray = new int[4];
   
9. 
   
10.             //循环12个月的数据
    
11.             for (int mIndex = 0; mIndex < 12; mIndex++)
    
12.             {
    
13.                 dataResult = new double[aDataInfo.X.Length * aDataInfo.Y.Length];
    
14.                 int dataIndex = 0;
    
15.                 //获取每个月的数据
    
16.                 var ddt = mdatas[mIndex].Data;
    
17. 
    
18.                 //将每个月的数据转换为数组赋值
    
19.                 for (int y = 0; y < aDataInfo.Y.Length; y++)
    
20.                 {
    
21.                     for (int x = 0; x < aDataInfo.X.Length; x++)
    
22.                     {
    
23.                         dataResult[dataIndex] = ddt[y, x];
    
24. 
    
25.                         dataIndex++;
    
26.                     }
    
27.                 }
    
28. 
    
29.                 startdataArray[0] = mIndex;
    
30.                 startdataArray[1] = dataIndex;
    
31.                 startdataArray[2] = 0;
    
32.                 startdataArray[3] = 0;
    
33. 
    
34.                 countdataArray[0] = 1;
    
35.                 countdataArray[1] = 1;
    
36.                 countdataArray[2] = aDataInfo.Y.Length;
    
37.                 countdataArray[3] = aDataInfo.X.Length;
    
38.                 gridFile.WriteVara("SST_cpl", startdataArray, countdataArray, dataResult);
    
39.             }         

复制代码

看天吃饭

OK，搞定了。我这个数据不需要四维的数组。每个时间点的数据不一样就循环写12次。
每次写变量数据的时候把时间节点编码赋值给startArray[0].

以下附上本次内容的代码：

1. NetCDFDataInfo aDataInfo = new NetCDFDataInfo();
   
2.             string aFile = Application.StartupPath + "\\Sample\\Hadisst_obs_1949-2010_T106_BCC_AGCM2.1.nc";
   
3.             //读取nc文件
   
4.             aDataInfo.ReadDataInfo(aFile);
   
5. 
   
6.             var t = aDataInfo.times;
   
7. 
   
8.             List<GridData>[] datas = new List<GridData>[12];
   
9.             int sYear = aDataInfo.times[0].Year;
   
10.             int eYear = aDataInfo.times[aDataInfo.times.Count - 1].Year;
    
11.             int yIndex = 0;
    
12.             GridData gData = new GridData();
    
13.             //循环年份
    
14.             for (int startYear = 1961; startYear < eYear; startYear++)
    
15.             {
    
16.                 yIndex = startYear - sYear;
    
17. 
    
18.                 //将每年12个月的数据都摆在相应月份的数组里
    
19.                 for (int mIndex = 0; mIndex < 12; mIndex++)
    
20.                 {
    
21.                     List<GridData> mList = datas[mIndex];
    
22. 
    
23.                     if (mList == null)
    
24.                     {
    
25.                         mList = new List<GridData>();
    
26.                     }
    
27. 
    
28.                     mList.Add(aDataInfo.GetGridData_LonLat(yIndex * 12 + mIndex, 3, 0));
    
29. 
    
30.                     datas[mIndex] = mList;
    
31.                 }
    
32.             }
    
33. 
    
34.             GridData[] mdatas = new GridData[12];
    
35. 
    
36.             //对每个月份的数据计算平均值
    
37.             for (int mIndex = 0; mIndex < 12; mIndex++)
    
38.             {
    
39.                 List<GridData> mList = datas[mIndex];
    
40. 
    
41.                 mdatas[mIndex] = DataMath.Average(mList, false);
    
42.             }
    
43. 
    
44.             NetCDFDataInfo gridFile = new NetCDFDataInfo();
    
45.             string outFilePath = Application.StartupPath + "\\Sample\\result.nc";
    
46.             gridFile.fileName = outFilePath;
    
47.             gridFile.UNDEF = aDataInfo.UNDEF;
    
48.             gridFile.unlimdimid = aDataInfo.unlimdimid;
    
49. 
    
50.             // Add dimensions:
    
51.             foreach (var dimension in aDataInfo.dimList)
    
52.             {
    
53.                 gridFile.AddDimension(dimension);
    
54.             }
    
55. 
    
56.             object[] lonDatas = new object[160];
    
57.             object[] latDatas = new object[320];
    
58. 
    
59.             // Add variables
    
60.             foreach (Variable variable in aDataInfo.varList)
    
61.             {
    
62.                 gridFile.AddVariable((Variable)variable.Clone());
    
63. 
    
64.                 //Add variable attributes
    
65.                 foreach (var attribute in variable.Attributes)
    
66.                 {
    
67.                     gridFile.AddVariableAttribute(variable.Name, attribute.attName, attribute.attValue.ToString());
    
68.                 }
    
69. 
    
70.                 if (variable.Name == "lon")
    
71.                 {
    
72.                     aDataInfo.GetVarData(variable, ref lonDatas);
    
73.                 }
    
74. 
    
75.                 if (variable.Name == "lat")
    
76.                 {
    
77.                     aDataInfo.GetVarData(variable, ref latDatas);
    
78.                 }
    
79.             }
    
80. 
    
81.             //Add global attributes
    
82.             foreach (var globalAttribute in aDataInfo.gAttList)
    
83.             {
    
84.                 gridFile.AddGlobalAttribute(globalAttribute.attName, globalAttribute.attValue.ToString());
    
85.             }
    
86. 
    
87.             //Create netCDF file
    
88.             gridFile.CreateNCFile(outFilePath);
    
89. 
    
90.             gridFile.WriteVar("lon", lonDatas);
    
91.             gridFile.WriteVar("lat", latDatas);
    
92. 
    
93.             int[] startArray = new int[1];
    
94.             int[] countArray = new int[1];
    
95.             double[] dataArray = new double[12];
    
96. 
    
97.             for (int mIndex = 0; mIndex < 12; mIndex++)
    
98.             {
    
99.                 dataArray[mIndex] = mIndex;
    
100.             }
     
101. 
     
102.             startArray[0] = 0;
     
103.             countArray[0] = 12;
     
104.             gridFile.WriteVara("time", startArray, countArray, dataArray);
     
105. 
     
106. 
     
107.             //创建一个可以盛放12个月数据的数组（从这个地方开始就不知道对不对了）
     
108.             double[] dataResult = new double[aDataInfo.X.Length * aDataInfo.Y.Length];
     
109. 
     
110. 
     
111.             //Add Data
     
112.             int[] startdataArray = new int[3];
     
113.             int[] countdataArray = new int[3];
     
114. 
     
115.             //循环12个月的数据
     
116.             for (int mIndex = 0; mIndex < 12; mIndex++)
     
117.             {
     
118.                 dataResult = new double[aDataInfo.X.Length * aDataInfo.Y.Length];
     
119.                 int dataIndex = 0;
     
120.                 //获取每个月的数据
     
121.                 var ddt = mdatas[mIndex].Data;
     
122. 
     
123.                 //将每个月的数据转换为数组赋值
     
124.                 for (int y = 0; y < aDataInfo.Y.Length; y++)
     
125.                 {
     
126.                     for (int x = 0; x < aDataInfo.X.Length; x++)
     
127.                     {
     
128.                         dataResult[dataIndex] = ddt[y, x];
     
129. 
     
130.                         dataIndex++;
     
131.                     }
     
132.                 }
     
133. 
     
134.                 startdataArray[0] = mIndex;
     
135.                 startdataArray[1] = 0;
     
136.                 startdataArray[2] = 0;
     
137. 
     
138.                 countdataArray[0] = 1;
     
139.                 countdataArray[1] = aDataInfo.Y.Length;
     
140.                 countdataArray[2] = aDataInfo.X.Length;
     
141.                 gridFile.WriteVara("SST_cpl", startdataArray, countdataArray, dataResult);
     
142.             }         
     
143.             
     
144.             gridFile.CloseNCFile();

NetCDFDataInfo aDataInfo = new NetCDFDataInfo();<br />             string aFile = Application.StartupPath + "\\Sample\\Hadisst_obs_1949-2010_T106_BCC_AGCM2.1.nc";<br />             //读取nc文件<br />             aDataInfo.ReadDataInfo(aFile);<br /> <br />             var t = aDataInfo.times;<br /> <br />             List<GridData>[] datas = new List<GridData>[12];<br />             int sYear = aDataInfo.times[0].Year;<br />             int eYear = aDataInfo.times[aDataInfo.times.Count - 1].Year;<br />             int yIndex = 0;<br />             GridData gData = new GridData();<br />             //循环年份<br />             for (int startYear = 1961; startYear < eYear; startYear++)<br />             {<br />                 yIndex = startYear - sYear;<br /> <br />                 //将每年12个月的数据都摆在相应月份的数组里<br />                 for (int mIndex = 0; mIndex < 12; mIndex++)<br />                 {<br />                     List<GridData> mList = datas[mIndex];<br /> <br />                     if (mList == null)<br />                     {<br />                         mList = new List<GridData>();<br />                     }<br /> <br />                     mList.Add(aDataInfo.GetGridData_LonLat(yIndex * 12 + mIndex, 3, 0));<br /> <br />                     datas[mIndex] = mList;<br />                 }<br />             }<br /> <br />             GridData[] mdatas = new GridData[12];<br /> <br />             //对每个月份的数据计算平均值<br />             for (int mIndex = 0; mIndex < 12; mIndex++)<br />             {<br />                 List<GridData> mList = datas[mIndex];<br /> <br />                 mdatas[mIndex] = DataMath.Average(mList, false);<br />             }<br /> <br />             NetCDFDataInfo gridFile = new NetCDFDataInfo();<br />             string outFilePath = Application.StartupPath + "\\Sample\\result.nc";<br />             gridFile.fileName = outFilePath;<br />             gridFile.UNDEF = aDataInfo.UNDEF;<br />             gridFile.unlimdimid = aDataInfo.unlimdimid;<br /> <br />             // Add dimensions:<br />             foreach (var dimension in aDataInfo.dimList)<br />             {<br />                 gridFile.AddDimension(dimension);<br />             }<br /> <br />             object[] lonDatas = new object[160];<br />             object[] latDatas = new object[320];<br /> <br />             // Add variables<br />             foreach (Variable variable in aDataInfo.varList)<br />             {<br />                 gridFile.AddVariable((Variable)variable.Clone());<br /> <br />                 //Add variable attributes<br />                 foreach (var attribute in variable.Attributes)<br />                 {<br />                     gridFile.AddVariableAttribute(variable.Name, attribute.attName, attribute.attValue.ToString());<br />                 }<br /> <br />                 if (variable.Name == "lon")<br />                 {<br />                     aDataInfo.GetVarData(variable, ref lonDatas);<br />                 }<br /> <br />                 if (variable.Name == "lat")<br />                 {<br />                     aDataInfo.GetVarData(variable, ref latDatas);<br />                 }<br />             }<br /> <br />             //Add global attributes<br />             foreach (var globalAttribute in aDataInfo.gAttList)<br />             {<br />                 gridFile.AddGlobalAttribute(globalAttribute.attName, globalAttribute.attValue.ToString());<br />             }<br /> <br />             //Create netCDF file<br />             gridFile.CreateNCFile(outFilePath);<br /> <br />             gridFile.WriteVar("lon", lonDatas);<br />             gridFile.WriteVar("lat", latDatas);<br /> <br />             int[] startArray = new int[1];<br />             int[] countArray = new int[1];<br />             double[] dataArray = new double[12];<br /> <br />             for (int mIndex = 0; mIndex < 12; mIndex++)<br />             {<br />                 dataArray[mIndex] = mIndex;<br />             }<br /> <br />             startArray[0] = 0;<br />             countArray[0] = 12;<br />             gridFile.WriteVara("time", startArray, countArray, dataArray);<br /> <br /> <br />             //创建一个可以盛放12个月数据的数组（从这个地方开始就不知道对不对了）<br />             double[] dataResult = new double[aDataInfo.X.Length * aDataInfo.Y.Length];<br /> <br /> <br />             //Add Data<br />             int[] startdataArray = new int[3];<br />             int[] countdataArray = new int[3];<br /> <br />             //循环12个月的数据<br />             for (int mIndex = 0; mIndex < 12; mIndex++)<br />             {<br />                 dataResult = new double[aDataInfo.X.Length * aDataInfo.Y.Length];<br />                 int dataIndex = 0;<br />                 //获取每个月的数据<br />                 var ddt = mdatas[mIndex].Data;<br /> <br />                 //将每个月的数据转换为数组赋值<br />                 for (int y = 0; y < aDataInfo.Y.Length; y++)<br />                 {<br />                     for (int x = 0; x < aDataInfo.X.Length; x++)<br />                     {<br />                         dataResult[dataIndex] = ddt[y, x];<br /> <br />                         dataIndex++;<br />                     }<br />                 }<br /> <br />                 startdataArray[0] = mIndex;<br />                 startdataArray[1] = 0;<br />                 startdataArray[2] = 0;<br /> <br />                 countdataArray[0] = 1;<br />                 countdataArray[1] = aDataInfo.Y.Length;<br />                 countdataArray[2] = aDataInfo.X.Length;<br />                 gridFile.WriteVara("SST_cpl", startdataArray, countdataArray, dataResult);<br />             }          <br />             <br />             gridFile.CloseNCFile();

MeteoInfo

搂主还是下了功夫的，通过写这个小程序应该对NetCDF格式有了更深入的理解。

看天吃饭

呵呵，对于netCDF格式文件终于从听说转变成知道了。
以后还要继续努力，加紧学习啊。感谢王老师的指导。