0
循环我有以下代码:写作值在csv文件
from math import *
from scipy.special import *
import matplotlib.pyplot as plt
import csv
## Definition der Parameter für Druckgleichung nach Rudnicki (1986) ##
q = 6.0/1000 # Fluidmenge pro Fläche und Zeit [l/s]
rhof = 1000 # Dichte Flüssigkeit [kg/m³]
lameu = 11.2*10**9 # Lamé-Parameter, undrained [Pa]
lame = 8.4*10**9 # Lamé-Parameter, drained [Pa]
pi # durch Pythonmodul "math" gegeben
alpha = 0.65 # Biot-Willis-Koeffizient
G = 8.4*10**9 # Schermodul [Pa]
k = 1.0e-15 # Permeabilität [m²] bzw. [Darcy]
eta = 0.001 # Viskosität des Fluids [Pa*s]
t = 1000*24*3600 # Zeit in [s]
## Beziehungen der Parameter untereinander ##
kappa = k/eta # Berechnung der Permeabilität nach Rudnicki (1986), [m³*s/kg]
print "kappa ist:",kappa # Ausgabe Permabilität
c = (kappa*(lameu-lame)*(lame+2*G))/((alpha**2)*(lameu+2*G)) # Berechnung der Diffusivität
print "c ist:",c # Ausgabe der Diffusivität
## Bereiche der Achsen in [m] ##
#def drange(start, stop, step):
# r = start
# while r <= stop:
# yield r
# r += step
xmin = 1
xmax = 20
ymin = 1
ymax = 20
## Druckberechnung um Bohrung ##
for x in range (xmin,xmax,1):
for y in range (ymin,ymax,1):
r = sqrt(x**2+y**2)
P = (q/(rhof*4*pi*kappa))*(expn(1,r**2/(4*c*t)))
z = P/1e6
#print x, " ", y, " ", z
对于我的输出中,我想写一个csv文件我的x,y和z值。它应该看起来像单行内的单元格中的每个x值和同一列内单元格中的每个y值,从而可以为每个x/y组合显示一个z值。它应该是这样的:
x x1 ...
YZ(X,Y)Z(X 1,Y)...
Y1 Z(X,Y)Z(X1,Y1)...
... ... ...
我是新来的csv模块,并没有真正习惯它,特别是如何设置值在正确的地方,以及何时使用它(例如里面或循环后)。这可能吗?或者可以将它只能是这样的(类似我outcommented打印):
XYZ(X,Y)
X Y1 Z(X,Y)
X Y2 Z(X,Y2) ......
对不起,这挺奇怪的和丑陋的格式,但我想你明白我的想法:)
任何帮助 - 一如既往 - 不胜感激!
什么是'scv'模块,什么是'pi'变量?它似乎没有初始化。 – Marcin 2014-12-19 01:27:15
抱歉Marcin,我忘了纠正这两个问题。感谢您的提示! – Alex 2014-12-19 12:00:43