質問編集履歴
1
コードの一部が見づらかったので一部書き直しました
title
CHANGED
|
File without changes
|
body
CHANGED
|
@@ -52,6 +52,7 @@
|
|
|
52
52
|
dz = sortz[line+5000,3]-sortz[line+25000,3]
|
|
53
53
|
dvorz = (surface2-surface1)/(dz)
|
|
54
54
|
print(dvorz)
|
|
55
|
+
|
|
55
56
|
'''以降の面を中心差分'''
|
|
56
57
|
for line in range(58):
|
|
57
58
|
surface1 = sortz[line:line+25000,3:7]
|
|
@@ -63,6 +64,7 @@
|
|
|
63
64
|
dvorz = np.r_(dvorz,dvorz1)
|
|
64
65
|
line += 25000
|
|
65
66
|
print(dvorz)
|
|
67
|
+
|
|
66
68
|
'''最後の面を後進差分'''
|
|
67
69
|
data1 = sortz[line:line+25000,3:7]
|
|
68
70
|
data2 = sortz[line+25000:line+50000,3:7]
|
|
@@ -70,6 +72,7 @@
|
|
|
70
72
|
dvorz1 = (data2-data1)/(dz)
|
|
71
73
|
dvorz = np.r_(dvorz,dvorz1)
|
|
72
74
|
print(dvorz)
|
|
75
|
+
|
|
73
76
|
np.savetxt('He_re1000_dvor-'+str(f),dvorz,delimiter=',',header="dvorx,dvory,dvorz")
|
|
74
77
|
f += 100
|
|
75
78
|
```
|