Weather Software > Station Software Development

Wet Bulb Globe Temperature

(1/1)

**marshad**:

Wet Bulb Globe Temperature (WBGT) is an important indicator of heat stress and is used by the military, sports and occupational health authorities. More info here: https://en.wikipedia.org/wiki/Wet-bulb_globe_temperature.

It needs a black globe thermometer to be measured accurately as it takes into account solar radiation, sun position and wind speed as well as temperature and humidity. I have tried to create a program in python to simulate WBGT using standard weather station data. The main source for the maths is here: https://www.weather.gov/media/tsa/pdf/WBGTpaper2.pdf. (I don't pretend to understand all the maths!)

--- Code: ---import math

from math import log

from math import atan

from math import exp

from math import cos

from math import radians

# Temperature in °C

# Humidity in %

# Pressure Sea Level in mbars

# WindSpeed in Km/h

# SolarRad in w/m3

# SunAltitude in degrees above horizon

# Cloudcover in %

# Enter the data

Temperature = 29.7

Humidity = 27

Pressure = 1014

WindSpeed = 0

SolarRad = 780

SunAltitude = 45

Cloudcover = 0

T = Temperature

RH = Humidity

P = Pressure

W = WindSpeed

S = SolarRad

Alt = SunAltitude

f = Cloudcover

u = W * 1000

if u == 0:

u = 500

z = radians(90 - Alt)

a = 17.368

b = 238.88

g = log(RH/100) + ((a*T)/(b+T))

# Calculate Dew Point DP

DP = (b * g)/(a-g)

# Calculate Wet Bulb WB

WB = T * atan(0.151977 * (RH + 8.313659)**0.5) + atan(T + RH) - atan(RH - 1.676331) + 0.00391838 * RH**1.5 * atan(0.023101 * RH) - 4.686035

ea = exp(17.67 * (DP - T)/(DP + 243.5)) * (1.0007 + 0.00000346 * P) * 6.112 * exp(17.502 * T / (240.97 + T))

eb = 0.575 * ea**0.142857

fdif = (100 - f)/100

fdb = 1 - fdif

B = S * ((fdb / ((4*(5.67*10**-8)) * cos(z))) + (1.2/(5.67*10**-8)) * fdif) + eb * T**4

C = (0.315*u**0.58)/(5.3865*10**-8)

# Calculate Globe Temperature GT

GT = (B + C * T + 7680000)/(C+256000)

# Calculate Wet Bulb Globe Temperature WBGT

WBGT = 0.7*WB + 0.2*GT + 0.1*T

print(round(WBGT, 2))

--- End code ---

Has anybody else tried to do this?

Does anybody have access to real WBGT data to test it?

The program is running live on my site using my own data (in the Temperature box)

http://moonappreciationsociety.com/Meteo/VillarPerosa/index.php

**SteveFitz1**:

I created a WBGT page on my site a couple of years ago using PHP. You can see my page here:

https://www.tylertexasweather.com/wbgtforecast.htm

I use a page on the NWS site to compare my results to. They are very close.

https://www.weather.gov/tsa/wbgt

Steve

**marshad**:

Thanks Steve

What are you using to calculate WBGT?

Can you share the code?

**SteveFitz1**:

Check your messages. I sent you a couple of PMs.

Steve

**marshad**:

--- Quote from: SteveFitz1 on July 18, 2021, 11:39:51 AM ---Check your messages. I sent you a couple of PMs.

Steve

--- End quote ---

Thanks -have responded.

Navigation

[0] Message Index

Go to full version