Script para cálculo do ponto de orvalho

 Para o primeiro post do blog vou compartilhar este script, que usará uma apróximação conhecida para calcular a temperatura do ponto de orvalho.

A dedução pode ser vista aqui

Para calcular a temperatura do ponto de orvalho, precisaremos apenas da temperatura instantânea (do bulbo seco, se preferir) e da umidade relativa (UR).

Este script é bem simples de escrever e pode ser usado pra ser executado no próprio prompt de comando inserindo os dados requisitados.

Atenção: Este script foi escrito em linguagem Python na versão 3+. Caso você utilize a 2 (versão descontinuada), favor fazer alterações necessárias.



import sys 

import numpy as np 

# Constantes utilizadas 

a = 17.271 

b = 237.7 

 

 #Aqui a temperatura será o argumento no 1  

T=float(sys.argv[1]) 

 

# Aqui a UR será o argumento no 2 

UR=float(sys.argv[2]) 

 

def ptorvalho_approx(T, UR):

     Td = (b * gamma(T,RH)) / (a - gamma(T,RH)) 

     return Td 

 

def gamma(T, UR): 

     g = (a * T / (b + T)) + np.log(UR/100.0) 

     return g 

 

Td = ptorvalho_approx(T, UR) 

print('T, UR', T, UR) 

print('Td = ', Td)


 

Como utilizar:

no terminal, você digita:

python3 script.py 29.2 98

Esse script calculará o ponto de orvalho com a temperatura de 29.2°C e Umidade Relativa de 98%

Qualquer dúvidas, só escrever nos comentários

Comentários

Postagens mais visitadas deste blog

Trabalhando com dados em NetCDF

Trabalhando com datas nos dados do INMET

Como ler dados em csv de estação do INMET