Skip to content

Commit

Permalink
Add Png for Readme
Browse files Browse the repository at this point in the history
  • Loading branch information
Yuhan-xue committed Jun 18, 2023
1 parent d77bdf1 commit ad7adb4
Show file tree
Hide file tree
Showing 11 changed files with 206 additions and 43 deletions.
154 changes: 149 additions & 5 deletions README.md
Original file line number Diff line number Diff line change
@@ -1,37 +1,181 @@
# HAN_Marine_element

- 2023年6月18日已全部完成

汉の2023春OUC海洋要素计算作业库

- 此仓库为中国海洋大学课程:海洋要素计算*(选课号:[070103101251])*的作业而建立
- 此仓库为中国海洋大学课程:海洋要素计算(选课号:[*070103101251*])的作业而建立
- 所有者:Hanxue Yu 余汉学
- 指导教师:王颖颖&杨俊超
- 参考:
- Github上可以找到的学长多年前写的[代码仓库][先輩のRepo]
- 海洋要素数据的质量控制中的格林布斯准则参考:[c-bata/outlier-utils][glbs]
- WORK2中的EOF计算来自于[Python EOF module][eofs]
- WORK2中的PDO index来自于[PDO index][pdo]
- WORK4中计算比容异常[GSW-Python][gsw]
- WORK5中cmap "Rocket"来自于[seaborn: statistical data visualization][seaborn]
- 刘永玲.杜凌.李静凯.翟方国 海洋要素计算上机实验指导书[M] 青岛:中国海洋大学出版社,2021
- GB/T 14914.6-2021, 海洋观测规范 第6部分:数据处理与质量控制[S].
- [GSW-Python][gsw]
- 字体[LXGW Bright][lxgw]

## 作业一 [海洋要素数据的质量控制][def1]

- 1.**奇异值判定与处理,生成新序列**
- 2.**分析两序列的平均值、标准差等统计特征的变化**
- 拓展作业见[拓展作业-空间二位插值][def1_d]

![WORK1代表图片](WORK1\results.png)

## 作业二 [利用EOF方法研究北太平洋海表面温度总体时空变化特征][def2]

- 选取1900-2020年为研究时间段,利用EOF方法对北太平洋海表面温度场数据进行分析。
- 讨论北太平洋SST在1900-2020年总体时空变化特征

![WORK2代表图片](WORK2\results.png)

## 作业三 [长期水位资料的调和分析][def3]

<table border=\"1\" class=\"dataframe\">
<thead>
<tr style=\"text-align: right;\">
<th></th>
<th>a</th>
<th>b</th>
<th>R</th>
<th>θ</th>
<th>H</th>
<th>g</th>
</tr>
<tr>
<th>分潮</th>
<th></th>
<th></th>
<th></th>
<th></th>
<th></th>
<th></th>
</tr>
</thead>
<tbody>
<tr>
<th>M2</th>
<td>-1460.102843</td>
<td>1239.338306</td>
<td>1915.165724</td>
<td>82.648415</td>
<td>1845.483605</td>
<td>8.659484</td>
</tr>
<tr>
<th>O1</th>
<td>-53.780013</td>
<td>-216.599202</td>
<td>223.175949</td>
<td>-59.769138</td>
<td>276.439846</td>
<td>64.557064</td>
</tr>
<tr>
<th>S2</th>
<td>368.429984</td>
<td>-377.902054</td>
<td>527.778946</td>
<td>-63.498182</td>
<td>527.778946</td>
<td>239.201911</td>
</tr>
<tr>
<th>K1</th>
<td>197.306115</td>
<td>228.400426</td>
<td>301.821898</td>
<td>49.177560</td>
<td>341.978206</td>
<td>307.634942</td>
</tr>
<tr>
<th>P1</th>
<td>51.051973</td>
<td>97.105078</td>
<td>109.707339</td>
<td>62.267323</td>
<td>109.707339</td>
<td>293.144287</td>
</tr>
<tr>
<th>N2</th>
<td>-349.550854</td>
<td>64.514394</td>
<td>355.454506</td>
<td>80.799918</td>
<td>342.521514</td>
<td>317.393238</td>
</tr>
<tr>
<th>Q1</th>
<td>13.719769</td>
<td>-43.172999</td>
<td>45.300551</td>
<td>-72.370479</td>
<td>56.112128</td>
<td>10.179002</td>
</tr>
<tr>
<th>K2</th>
<td>-51.515602</td>
<td>102.913995</td>
<td>115.087565</td>
<td>38.694182</td>
<td>153.633137</td>
<td>72.628583</td>
</tr>
<tr>
<th>M4</th>
<td>34.167625</td>
<td>-87.967445</td>
<td>94.370006</td>
<td>-68.773225</td>
<td>87.627747</td>
<td>17.526245</td>
</tr>
<tr>
<th>MS4</th>
<td>-21.067135</td>
<td>53.356617</td>
<td>57.365083</td>
<td>18.557827</td>
<td>53.266638</td>
<td>248.168533</td>
</tr>
</tbody>
</table>

![WORK3代表图片](WORK3\results.png)

## 作业四 [地转流计算][def4]

## 作业五 [海浪资料统计分析-施工ing 06\07\2023][def5]
![WORK4代表图片](WORK4\2018-01%2010%20dbar.png)

- 1.利用两月的月平均资料计算北太平洋*6°N-35°N*范围内的地转流,选取**1500db**作为参考零面画出*10db**100 db**250 db**500 db*等四个深度层上的流场和流速
- 2.利用上面计算结果,计算北赤道流水体输运(如130°E,8°N-18°N断面),比较讨论两月结果差异,也可进一步比较不同断面的差异

## 作业五 [海浪资料统计分析][def5]

- 1. 画出波面高度时间序列,利用上跨零点或者下跨零点的方法读取数据资料中的波高、周期等波浪要素,并画图展示。
- 2. 画出波高的概率密度分布,并计算给出有效波高和有效波周期

![WORK5代表图片](WORK5\result.png)

[先輩のRepo]: https://github.com/FinalTheory/oceanography-numerical-calculations
[glbs]: https://github.com/c-bata/outlier-utils
[gsw]: https://teos-10.github.io/GSW-Python/gsw_flat.html
[eofs]: https://ajdawson.github.io/eofs/latest/api/eofs.standard.html
[pdo]: https://www.ncei.noaa.gov/access/monitoring/pdo/
[seaborn]: http:https://seaborn.pydata.org/
[lxgw]: https://github.com/lxgw/LxgwBright
[def1]: https://nbviewer.org/github/Yuhan-xue/HAN_Marine_element/blob/main/WORK1/WORK1.ipynb
[def1_d]: https://nbviewer.org/github/Yuhan-xue/HAN_Marine_element/blob/main/WORK1/%E6%8B%93%E5%B1%95/%E6%8B%93%E5%B1%951.ipynb
[def2]: https://nbviewer.org/github/Yuhan-xue/HAN_Marine_element/blob/main/WORK2/work2.ipynb
[def2]: https://nbviewer.org/github/Yuhan-xue/HAN_Marine_element/blob/main/WORK2/WORK2.ipynb
[def3]: https://nbviewer.org/github/Yuhan-xue/HAN_Marine_element/blob/main/WORK3/WORK3.ipynb
[def4]: https://nbviewer.org/github/Yuhan-xue/HAN_Marine_element/blob/main/WORK4/WORK4.ipynb
[def5]: https://nbviewer.org/github/Yuhan-xue/HAN_Marine_element/blob/main/WORK5/WORK5.ipynb
[def5]: https://nbviewer.org/github/Yuhan-xue/HAN_Marine_element/blob/main/WORK5/WORK5.ipynb
30 changes: 15 additions & 15 deletions WORK1/WORK1.ipynb

Large diffs are not rendered by default.

Binary file added WORK1/__pycache__/grubbs.cpython-310.pyc
Binary file not shown.
Binary file added WORK1/results.png
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file added WORK2/results.png
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
23 changes: 12 additions & 11 deletions WORK2/work2.ipynb

Large diffs are not rendered by default.

3 changes: 2 additions & 1 deletion WORK3/WORK3.ipynb
Original file line number Diff line number Diff line change
Expand Up @@ -826,7 +826,8 @@
"plt.legend()\n",
"plt.title('Tide')\n",
"del i,n,n0,S0,σ\n",
"#plt.xlim(0,4000)"
"#plt.xlim(0,4000)\n",
"plt.savefig('results.png')"
]
},
{
Expand Down
Binary file added WORK3/results.png
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
14 changes: 7 additions & 7 deletions WORK3/潮差.csv
Original file line number Diff line number Diff line change
Expand Up @@ -5,7 +5,7 @@
3,1996-12-28 17:00:00,5268.4807233128195,1996-12-28 10:00:00,1355.662646022257,3912.8180772905625
4,1996-12-29 05:00:00,5954.533195204371,1996-12-28 22:00:00,2024.0669341634778,3930.466261040893
5,1996-12-29 17:00:00,5546.477971891545,1996-12-29 11:00:00,1233.0706549389301,4313.407316952615
6,1996-12-30 05:00:00,6043.107710558454,1996-12-29 23:00:00,1812.9462126929097,4230.161497865543
6,1996-12-30 05:00:00,6043.107710558454,1996-12-29 23:00:00,1812.9462126929095,4230.161497865544
7,1996-12-30 18:00:00,5780.578272203231,1996-12-30 11:00:00,1162.7332670091555,4617.845005194075
8,1996-12-31 06:00:00,6097.920039082281,1996-12-31 00:00:00,1658.354979321264,4439.5650597610165
9,1996-12-31 19:00:00,5836.999703996469,1996-12-31 12:00:00,1088.7238130131814,4748.275890983288
Expand All @@ -14,17 +14,17 @@
12,1997-01-02 07:00:00,5826.211136858037,1997-01-02 01:00:00,1435.5621069844046,4390.649029873633
13,1997-01-02 20:00:00,5882.690648827307,1997-01-02 13:00:00,1228.6520138690507,4654.038634958257
14,1997-01-03 08:00:00,5605.993594139584,1997-01-03 02:00:00,1508.5236057854236,4097.469988354161
15,1997-01-03 20:00:00,5854.236878431358,1997-01-03 14:00:00,1388.1639459917924,4466.072932439565
15,1997-01-03 20:00:00,5854.236878431358,1997-01-03 14:00:00,1388.163945991792,4466.072932439565
16,1997-01-04 08:00:00,5344.570359143128,1997-01-04 02:00:00,1533.778771997198,3810.7915871459304
17,1997-01-04 21:00:00,5689.629165182086,1997-01-04 14:00:00,1567.4151720654195,4122.213993116666
18,1997-01-05 09:00:00,5077.937043848184,1997-01-05 03:00:00,1655.451322626803,3422.485721221381
17,1997-01-04 21:00:00,5689.629165182087,1997-01-04 14:00:00,1567.4151720654195,4122.213993116668
18,1997-01-05 09:00:00,5077.937043848184,1997-01-05 03:00:00,1655.4513226268027,3422.485721221381
19,1997-01-05 21:00:00,5595.680500260027,1997-01-05 15:00:00,1775.222197360745,3820.458302899282
20,1997-01-06 09:00:00,4787.927325562969,1997-01-06 04:00:00,1850.9277118598382,2936.9996137031303
21,1997-01-06 22:00:00,5393.931646593164,1997-01-06 15:00:00,2068.2340661714065,3325.6975804217577
22,1997-01-07 10:00:00,4526.896368606392,1997-01-07 04:00:00,1962.707354754568,2564.189013851824
23,1997-01-07 22:00:00,5231.261035042546,1997-01-07 16:00:00,2276.076235550154,2955.184799492392
24,1997-01-08 11:00:00,4274.607670215859,1997-01-08 05:00:00,2065.6658535408515,2208.9418166750074
25,1997-01-08 23:00:00,5056.3139446246905,1997-01-08 17:00:00,2479.6882964258307,2576.62564819886
25,1997-01-08 23:00:00,5056.3139446246905,1997-01-08 17:00:00,2479.688296425831,2576.6256481988594
26,1997-01-09 13:00:00,4165.330681883333,1997-01-09 06:00:00,2099.509479130666,2065.8212027526674
27,1997-01-10 01:00:00,4980.376126190547,1997-01-09 18:00:00,2595.1546972543306,2385.221428936217
28,1997-01-10 14:00:00,4345.108163211693,1997-01-10 07:00:00,2020.405529386535,2324.702633825158
Expand All @@ -45,7 +45,7 @@
43,1997-01-18 08:00:00,6225.696810778043,1997-01-18 02:00:00,716.7410242705382,5508.955786507505
44,1997-01-18 20:00:00,6517.0234075302815,1997-01-18 14:00:00,905.5205377813445,5611.502869748937
45,1997-01-19 09:00:00,5867.744308320934,1997-01-19 02:00:00,738.1271578051369,5129.617150515796
46,1997-01-19 21:00:00,6350.4109515019145,1997-01-19 15:00:00,1215.8157550902683,5134.595196411647
46,1997-01-19 21:00:00,6350.4109515019145,1997-01-19 15:00:00,1215.8157550902688,5134.595196411646
47,1997-01-20 09:00:00,5530.124393897435,1997-01-20 03:00:00,902.8878324983402,4627.236561399095
48,1997-01-20 22:00:00,6056.240144068899,1997-01-20 15:00:00,1472.7437597222315,4583.496384346668
49,1997-01-21 10:00:00,5127.07223143284,1997-01-21 04:00:00,1163.6892333291307,3963.3829981037097
Expand All @@ -55,7 +55,7 @@
53,1997-01-23 13:00:00,4616.485231325356,1997-01-23 06:00:00,1666.8147468779837,2949.6704844473725
54,1997-01-24 01:00:00,5348.889981574496,1997-01-23 18:00:00,2290.2233605844776,3058.6666209900186
55,1997-01-24 14:00:00,4713.895915156884,1997-01-24 07:00:00,1784.4596777597494,2929.436237397135
56,1997-01-25 02:00:00,5373.983273279984,1997-01-24 19:00:00,2352.455919102961,3021.527354177023
56,1997-01-25 02:00:00,5373.983273279984,1997-01-24 19:00:00,2352.4559191029616,3021.5273541770225
57,1997-01-25 15:00:00,4964.152477410089,1997-01-25 08:00:00,1781.739304058788,3182.413173351301
58,1997-01-26 03:00:00,5516.068188334377,1997-01-25 21:00:00,2231.3877293258674,3284.6804590085094
59,1997-01-26 16:00:00,5276.983075250683,1997-01-26 09:00:00,1684.255137471869,3592.7279377788145
Expand Down
25 changes: 21 additions & 4 deletions WORK5/WORK5.ipynb
Original file line number Diff line number Diff line change
@@ -1,11 +1,10 @@
{
"cells": [
{
"attachments": {},
"cell_type": "markdown",
"metadata": {},
"source": [
"# 海洋要素计算 编程作业5\n",
"# <a id='toc1_'></a>[海洋要素计算 编程作业5](#toc0_)\n",
"\n",
"\n",
"- 1. 画出波面高度时间序列,利用上跨零点或者下跨零点的方法读取数据资料中的波高、周期等波浪要素,并画图展示。 \n",
Expand All @@ -19,6 +18,23 @@
"**Student ID: 20010006082**"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"**Table of contents**<a id='toc0_'></a> \n",
"- [海洋要素计算 编程作业5](#toc1_) \n",
"\n",
"<!-- vscode-jupyter-toc-config\n",
"\tnumbering=false\n",
"\tanchor=true\n",
"\tflat=false\n",
"\tminLevel=1\n",
"\tmaxLevel=6\n",
"\t/vscode-jupyter-toc-config -->\n",
"<!-- THIS CELL WILL BE REPLACED ON TOC UPDATE. DO NOT WRITE YOUR TEXT IN THIS CELL -->"
]
},
{
"cell_type": "code",
"execution_count": 24,
Expand Down Expand Up @@ -210,7 +226,7 @@
},
{
"cell_type": "code",
"execution_count": 34,
"execution_count": 36,
"metadata": {},
"outputs": [
{
Expand All @@ -233,6 +249,7 @@
}
],
"source": [
"\n",
"\n",
"def cal_H(arr, indices):\n",
" dists = [arr[indices[i]:indices[i+1]] for i in range(len(indices)-1)]\n",
Expand Down Expand Up @@ -295,7 +312,7 @@
"ax3.tick_params(right=True,left=True,labelright=True,labelleft=False,labelbottom=False)\n",
"ax.tick_params(right=True,left=True,top=True,bottom=True)\n",
"plt.colorbar(c,orientation='horizontal',label='Percent (%)')\n",
"plt.show()"
"plt.savefig('result.png')"
]
},
{
Expand Down
Binary file added WORK5/result.png
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.

0 comments on commit ad7adb4

Please sign in to comment.