program RelativeHumidity

! Written by Kate T-C
! February 26, 2009  For the Fortran Short Course
! This program estimates the relative humidity for different temperatures 
! at a common pressure and mixing ratio

implicit none

real :: Press = 1013., MR = 0.005
real, dimension(300) :: Temps = 0., RHs = 0.
real :: RelHum
integer :: i = 0

! Initialize the temperatures
do i = 1, 300
	Temps(i) = 273.15 + 100.0/300.0*(i-1)
enddo

do i = 1, 300
	RHs(i) = RelHum(Temps(i),Press,MR)
enddo

print *, RHs
print *, Temps

end program RelativeHumidity