# 3 April 2002, SEY
# this file is /home/disk/gamera/kwajex/2001dualpol/INSTRUCT.cfad
# also on homesun at /home/disk/echo2/kwajex/dualpol

# instructions on making a time series of cfad data
# this example works uses a 1 day's worth of files,
# the final timeseries will be for all 6 months of files.
# for the surveillance scan a cfad is equivalent to making a time series
# of histogram data

## create the data file with a C program called cdf2d2cfad_list

on hydra

cd /home/disk/gamera/kwajex/dualpol

~yuter/src/cdf-src/cdf2d2cfad_list 1day.list maxdz 0 60 5 1day.output

# note for this example 1day.list was created with
ls -1 /zd/kwajex/kwajex_lowprf/*lowprf.20010601*cdf > 1day.list

# lists of any subset or all the 2001 data can be created with appropriate
# use of wild cards in the ls function

# 1day.output looks like:

maxdz: cfad min=0.00 max=60.00 bin=5.00 numint=12
inputfile x y z x_spacing y_spacing counts_in_each_bin 
20010601 0000 121 121 1 4 4     14 36 106 184 235 140 63 28 5 0 0 0 
20010601 0001 121 121 1 4 4     11 23 68 130 131 128 67 35 6 0 0 0 
20010601 0012 121 121 1 4 4     22 33 95 203 238 118 64 26 0 0 0 0 
20010601 0024 121 121 1 4 4     21 24 100 186 219 168 72 32 2 2 0 0 
20010601 0036 121 121 1 4 4     32 36 89 200 248 149 94 37 3 1 0 0 
20010601 0048 121 121 1 4 4     29 31 97 266 217 179 86 41 9 0 0 0 
...


## read in and plot the data with Splus

# note all the Splus functions are defined in ~yuter/.Data/
# to put them in your Splus folder just copy the filename corresponding
# to the function name to your .Data folder. There  will be functions
# called by a particular function that need to be copied as well.

# to read the data one can use the Splus function ~yuter/.Data/readcfadlist

test_readcfadlist("/home/disk/gamera/kwajex/dualpol/1day.output",12,122)

# the arguments are the number of intervals, in this case (60-0)/5 =12
# and the number of files, do 
 
hydra:yuter:907>wc 1day.output
     124    2331    7737 1day.output

# 124 is the number of lines in the file, there are two header lines so
# the number of files is 122

# now compute the area greater than a threshold of 20 dBZ, since
# the lower limits of the bins are 0, 5, 10, 15, 20, 25 etc. we want to
# start counting at 5

test.area_computecfadlistareas(test,5)

# now plot the timeseries

showmints(test.area$kmarea/(121*121*4*4)*100,test.area$julmin,atitle="KWAJ 2001 Precip Echo Area",minint=60*1,ylab="% Echo Area > 20 dBZ within 240 km range")

# note minint is the interval for the tics on the time axis, 60 is 60 minutes
# for daily use 60*24

# you can set a ylim as a parameter to showmints, for example
ylim=c(0,80), if ylim is not set it will scale the graph to fit the
maximum in the data.

# save the graph to a file
printgraph(file="/home/disk/gamera/kwajex/dualpol/1dayts.ps")