! compute the wind dir freq distribution for wind_gen station subroutine get_dirf(dirf) include 'wind_db.fi' ! summation of monthly frequencies, partial summation of monthly ! frequencies, and directional frequency for each month real sum(0:NMONTHS-1) real sumpar(0:NMONTHS-1) real dirf(0:NDIR,0:NMONTHS-1) integer i, m ! do the following computations to populate the dirf() array do 60 m = 0, NMONTHS-1 sum(m) = 0.0 do 60 i = 0, NDIR sum(m) = sum(m) + dir(i,m) 60 continue ! compute dirf(i,m) for each month do 70 m = 0, NMONTHS-1 sumpar(m) = 0.0 do 70 i = 0, NDIR sumpar(m) = sumpar(m) + dir(i,m) dirf(i,m) = sumpar(m) / sum(m) 70 continue return end