The following programs
were used to edit the aes data for Prince George.
They can be run or copied from stratus: /home/lindsayv/aes_pg
Unless otherwise specified the output files are named: program.f- the fortran program file composed
in emacs. program.exe
- executable version of the program program.out
- output file from program
Before use check what
the actual input files are. Some programs were used to extract information
from several files by changing the input file names. Therefore input/output
may not match program names. Also, programs may have been modefied
to use a program.in
and program.out
file.
diff.f
This program compares two files and outputs the data which is found
in only one of them.
same.f
This program compares two files based on absolute dates and outputs
the data for all the dates which are the same (in both files).
the input variables are fix.in and fix2.in
joinday.f
This program merges the daily aes data for PG according to: yr, mo,
day, tmp_max, tmp_mean, ppt,gust_spd, gust_dir
joinhly.f
merges the hourly data for Prince George according to: yr, mo, day,
hr, tmp_dry, tmp_wet, prs, gust_spd, gust_dir, sol
The pressure data goes to 2004-02-03 but all the rest (temp. being
the longest) only goes until 2003-10-31 at the latest so I ended the file
at 2003-10-31.
jointwo.f
this program is the base for the two above program (joinhly.f and joinday.f)
and may be useful for re-writing of a similar program. It joins two files
into a third based on date comparison.
date.f
adds an absolute and Julian date to the file. The input file is currently
joinday.out but the program should work fine for the other join programs
if the variables are specified to meet the input variables.
dateJA.f
Isolates all the July and August data from date.out.
allhotdays.f
outputs the same data as fair_weather.f (including N values) for all
days with temperatures greater or equal to 20ºC. (Not limited to day
3 values)
fair_weather.f
output is hotday3.out.
extracts the day three information from a date.out file.
A variable N is also added which is the total number of days in the
heating cycle.
Each output string contains: dat, jd, yr, mo, dy, N, tmp, tmp3, ppt,
spd, dir
dat=absolute date
jd=julien day
N=number of days in cycle
tmp=maximum tmp.
tmp3=mean tmp.
The criteria for a fair weather heating cycle are daily maximum temperature
of more then 20ºC for a minimum of 4 consecutive days. The output file
is hotday.out
fair_weather_JA.f
same as fair_weather.f except further limited to July and August dates
only.
peakall.f
This program outputs all the days and cycle numbers for peak cycles.
peakall_fix.f was created to fix an error in this file but is no longer
necessary since the error was also fixed in the original program.
ben.f :-)
This program extracts specific dates from the input file. At the moment
the input file is an original (1096450.156.out) file and the joinhly.out
data file.
This may be a good base for extracting other dates of interest.peakall.f
Filed in hourly directory:
hlyhotday3.f hlyJA.f hlypeak.f
These three programs extract the corresponding data sets from the joinhly.out
file. Therefore the output file contains the matching hourly data for the
hotday3.out, all July and August dates, and peak.out files.
hly15.f
Limit hourly data input to July15 to August 15 days.
wind.f
This program was desined to create input files for the rose diagram
macro in excel. It converts the wind speed and direction from km/hr
to m/s and direction from 10's of degrees to degrees. It also removes -999
values (since these confuse the results and would have to be removed before
running macro or R. They are placed speraticaly somwhere in the 75888 original
values - not at begining or end.)
The input file is wind.in and
the output files are :
wind.out which outputs the
converted wind speed (as a real of format 00.0000) and direction
day_wind.out which limits
the above output to hours between 10 and 16 (10:00am to 4:00pm).
the desired data set copied into the wind.in file should be in the
date.out format (yr, mo, dy, hr, tmp, tmp2, prs, spd, dir, sol).
Filed in hothr directory:
hothr.f -data for fair weather day three (hotday3.out) hothrJA.f hothrPeak.f
These programs extract the data for the first hottest hour in each day of the
input file - corresponding to the above three programs.
Due to the differences in starting dates of these three data sets changing
the input file alone will not work. Therefore these three versions
of the program were created.
hothrall.f hothrJAall.f hothrallPeak.f
These programs extract the data for all the hottest hours for each day of the
input file . For example, the output for hothrall.f contains the data
for all the hottest hours of each day three of a fair weather
cycle.
These files use hothr.out, hothrJA.out, and hothrPeak.out as wall as
joinhly.out as input.