macro wico
* 
* Plot histograms from the output of wico
* 12/17/07  ES
set *
*
* set options
*
hi/create/title_global 'Wico input and output'
option ndate
option nbox
*set stat 1111111
set stat 111
option nstat
set fit 111
option fit
option ngrid
*
* plotting options
*
*set xmgl 4.
*set ymgl 4.
*set asiz 0.4
*set xlab 2.
*set ylab 1.
*set xsiz 20.
*set xmgl 3.
*set ymgl 3.
*set ysiz 20.
*set gsiz 0.4
*
* set font definitions to bold roman 
*
set CFON -21
set GFON -21
set LFON -21
set TFON -21
set VFON -21
set txfp -21
set SMGU 0.02
set SMGR 0.02
set CSIZ 0.33
set VSIZ 0.3
set TSIZ 0.35
set YHTI 0.9
set HWID 3.0
set BWID 3.0
set asiz 0.38
*
* open metafile
*
for/file 66 wico.ps
meta 66 -111
*for/file 66 wico.eps
*meta 66 -113
*
npts1 = 101
*
* example vector creation
*
vector/create e([npts1]) R [npts1]*0.01
emin = 0.1
emax = 10
npts = [npts1] - 1
sigma j=array([npts1],0#[npts1]-1)
sigma e = [emin] + ([emax]-[emin])*j/[npts]
*
* input histograms
*
hi/file 1 wico.hbook
*
zone
hi/plot 1 'colz'
hi/plot 2 'colz'
zone 2 2
hi/plot 4
hi/plot 14
hi/plot 5
hi/plot 15
hi/plot 6
hi/plot 16
hi/plot 7
hi/plot 17
hi/plot 8
hi/plot 18
zone
hi/plot 9 'colz'
hi/plot 19 'colz'
*
close 66
*
exitm 
return