def plotPhoH():
if verificaVariaveisProfResi():
return -1
# plt.subplot(211)
plt.plot(profundidade, resistividadeMedia)
plt.xlabel("Profundidade [m]")
plt.ylabel("Resistividade Media [ohm*m]")
plt.title("Curva de Resistividade")
plt.grid(True)
plt.savefig(dirCurvas + ajudante.separador() + "curvadeResistividade_" + idPlanilha + "_" + formataHora() + "_.png")
print "aviso: arquivo png da figura salvo na pasta <curvas>"
# plt.show()
# s = UnivariateSpline(profundidade, resistividadeMedia, s=1)
# xs = linspace(0, profundidade[len(profundidade)-1], 1000)
# ys = s(xs)
# plt.subplot(212)
# plt.plot(xs, ys)
plt.show()
# 7 - Valor máximo da resitência de terra de modo a ser compatível com a sensibilidade da
# proteção.
#-------------------------------------------------------------------------------
#
from __future__ import division
from math import sqrt, log, pi
import estratificacao
import potenciais
import ConfigParser
from os import getcwd, mkdir
from os.path import isdir, isfile
import ajudante
pastaTrabalho = getcwd()+ajudante.separador()+'tabelas'
nomeArquivoProjeto = 'projetoMalha.cfg'
nomeArquivoProjetoCompleto = pastaTrabalho+ajudante.separador()+nomeArquivoProjeto
projetoMalha = {
# configurações do solo
'peq' : 0,
'pn1' : 0,
'deq' : 0,
# configurações da brita
'psBrita' : 0,
'hsBrita' : 0,
# configurações da malha de terra
'mLargura' : 0,
# GLOBAIS
# Qual sistema eu estou usando?
# if os() == 'Linux':
# separador = '/'
# cmdLimpaTela = 'clear'
# else:
# separador = '\\' # para o Windows
# cmdLimpaTela = 'cls'
# identificação do arquivo do excel ou arquivo qualquer, usado na criação de
# arquivos com plot ou no armazenamento de variaveis
idPlanilha = ""
# diretório para armazenamento das curvas
ndirCurvas = "curvas"
dirCurvas = getcwd() + ajudante.separador() + ndirCurvas
# diretório para armazenamento dos resultados
ndirResultados = "resultados"
dirResultados = getcwd() + ajudante.separador() + ndirResultados
# Dicionário com todas as principais variáveis de controle do sistema
sistemaVar = {
# se 0 a entrada para o diâmetro da haste é dada em m
# se 1 a entrada para o diâmetro da haste é dada em mm
# se 2 a entrada para o diâmetro da haste é dada em polegas
"unidadeHaste": "pol",
"usaValoresArquivo": 1,
"debugAterramento": 0,
"idPlanilha": "",
"arqMalha": "",
"arqTabela": "",
def testeEstratificacaoArquivos():
dirAtual = getcwd()
pastaTabelas = "tabelas"
planilhas = [
dirAtual + ajudante.separador() + pastaTabelas + ajudante.separador() + "tabelaExemplo2_12GeraldoKindermann.xlsx",
dirAtual + ajudante.separador() + pastaTabelas + ajudante.separador() + "subestacaoMamede.xlsx"]
#--------------------------------------------------------------------------
# testa a estratificação em duas camadas
print 'Inciando teste de estratificacao em 2 camadas,'
for i in range(9):
print
print 'caso, ', i
iniciaConstantes(i, debug=True)
[p1, k, h] = estratifica2Camadas(debug=True)
print 'valores,'
print 'p1, ', p1
print 'k, ', k
print 'h, ', abs(h)
print 'p2, ', p2solo2Camadas(p1, k)
print '_' * 80
print 'testando individualmente a funcao de estratificacao'
iniciaConstantes(0, debug=True)
print 'saida, ', funcaoEstratificacao([364.67, -0.43, -2.827])
print '-' * 40
iniciaConstantes(3, debug=True)
print 'saida, ', funcaoEstratificacao([12, 0, 0])
print '-' * 40
iniciaConstantes(3, debug=True)
print 'saida, ', funcaoEstratificacao([472, 0.088811995, 7.8])
print
print '**fim'
print '_' * 80
#--------------------------------------------------------------------------
# testa a leitura da planilha com os dados
print 'Iniciando teste leitura de uma planilha'
for i in planilhas:
print 'lendo planilha, ', i
m = lerPlanilha(i, debug=True)
[phoTabela, esTabela] = resistividadeMediaPlanilha(m, debug=True)
print 'retorno, '
print phoTabela
print esTabela
pho = esTabela
es = phoTabela
print pho, es
if debugPlot:
plt.plot(phoTabela, esTabela)
plt.show()
[p1, k, h] = estratifica2Camadas(debug=True)
print 'valores, '
print 'p1, ', p1
print 'k, ', k
print 'h, ', h
print 'p2, ', p2solo2Camadas(p1, k)
print '-' * 80
print
