cdoubint - Double causal integration
call cdoubint_init(nx)
ierr=cdoubint_lop(adj,add,x,y)
number of data points
standard operator interface
Double causal integration
the causint manpage
geef90