#$=head1 NAME
#$
#$causint - Causal integration
#$
#$=head1 SYNOPSIS
#$
#$C<ierr=causint_lop(adj,add,x,y)>
#$
#$=head1 PARAMETERS
#$
#$=over 4
#$
#$=item adj,add,x,y  - 
#$
#$      standard operator interface
#$
#$=back
#$
#$=head1 DESCRIPTION
#$
#$Causal integration
#$
#$=head1 SEE ALSO
#$
#$L<cdoubint>
#$
#$=head1 LIBRARY
#$
#$B<geef90>
#$
#$=cut
module causint {                                         # causal integration
#% _lop( xx,  yy) 
integer i, n       
real t
n = size (xx); t = 0
if( adj){ 
                 do i= n, 1, -1 {
                   t = t + yy(i)
                   xx(i) = xx(i) + t
		}
} else {
         	do i= 1, n {
                    t   = t + xx(i)
	            yy(i) = yy(i) + t
        	}
        }
}