def estratificacaoSolo(silencioso=0):
global sistemaResultados
if verificaVariaveisProfResi():
return -1
if not silencioso:
print u"Iniciando a estratificação do solo"
estratificacao.pho = resistividadeMedia
estratificacao.es = profundidade
if sistemaResultados["fEstratificacao"] == 0:
[p1, k, h] = estratificacao.estratifica2Camadas()
p2 = estratificacao.p2solo2Camadas(p1, k)
sistemaResultados["p1"] = p1
sistemaResultados["p2"] = p2
sistemaResultados["k"] = k
sistemaResultados["h"] = h
sistemaResultados["fEstratificacao"] = 1
if not silencioso:
print "Resistividade da primeira camada(ohm*m), ", sistemaResultados["p1"]
print "Resisitivdade da segunda camada(ohm*m), ", sistemaResultados["p2"]
print "Coeficiente de reflexao, ", sistemaResultados["k"]
print "Profundidade da primeira camada(m), ", sistemaResultados["h"]
return [sistemaResultados["p1"], sistemaResultados["k"], sistemaResultados["h"]]
def plotCurvaTeorica2Camadas():
global sistemaResultados
if verificaVariaveisProfResi(silencioso=1):
if lerTabelaExcel(silencioso=1):
return -1
if sistemaResultados["fEstratificacao"] == 0:
print u"Iniciar o processo de estratificação em 2 camadas[S/n]?",
if raw_input() != "n":
estratificacao.pho = resistividadeMedia
estratificacao.es = profundidade
[
sistemaResultados["p1"],
sistemaResultados["k"],
sistemaResultados["h"],
] = estratificacao.estratifica2Camadas()
sistemaResultados["fEstratificacao"] = 1
else:
print u"erro: é necessário iniciar a estratificação"
return -1
print sistemaResultados
a = arange(1e-3, profundidade[len(profundidade) - 1], 1e-2)
p = []
for i in a:
p.append(
estratificacao.funResistividade2Camadas(
i, sistemaResultados["p1"], sistemaResultados["k"], sistemaResultados["h"]
)
)
plt.xlabel("Profundidade [m]")
plt.ylabel("Resistividade Media [ohm*m]")
plt.title("Curva teoria e medida")
plt.grid(True)
plt.plot(a, p, profundidade, resistividadeMedia)
plt.show()
