pbaylis

library(foreign)
library(data.table)
library(zoo)
library(Hmisc)
library(gstat)
library(sp)
rm(list=ls())

# Fill in missing data ----
sm.thresh <- 0.9 # Min proportion of obs that must be in a station-month

pbaylis

program define gen_relative_humidity
	* Calculate relative humidity given dry bulb temperature and dew point
	syntax, DRYbulbtemp(varname) DEWpoint(varname) [GENerate(string)]
	* Inputs in F
	local T `drybulbtemp'
	local TD `dewpoint'

	if "`generate'" == "" {
		local generate "rel_humidity"
	}

pbaylis

program define gen_heat_index
	* Calculat heat index
	* See https://en.wikipedia.org/wiki/Heat_index. Note that there is no exact heat index formula
	* since it is defined by the table. There are approximations, but this code uses the table itself.
	* The only modification is that I cap at 140.
	syntax, DRYbulbtemp(varname) RELativehumidity(varname) [GENerate(string)]
	
	if "`generate'" == "" {
		local generate "heat_index"
	}

pbaylis

program define makebins, rclass
    syntax varname, [PREfix(name) min(real 0) max(real 100) step(real 5) SEQuential DROPbin(real 70)]
    local var `varlist'

    if mod(`max' - `min',`step') != 0 {
        di as error "Warning: Step size `step' does not divide evenly into range. Top bin may extend beyond `max'."
    }

    * Declare prefix
    if "`prefix'" != "" {