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()