#$=head1 NAME
#$
#$igrad1 - 1D gradient operator
#$
#$=head1 SYNOPSIS
#$
#$C<ierr=igrad1_lop(adj,add,xx,yy)>
#$
#$=head1 PARAMETERS
#$
#$=over 4
#$
#$=item  adj,add,xx,yy - 
#$
#$       standard operator interface
#$
#$=back
#$
#$=head1 DESCRIPTION
#$
#$ Calculate 1-D gradient
#$
#$=head1 SEE ALSO
#$
#$L<igrad2>
#$
#$=head1 LIBRARY
#$
#$B<geef90>
#$
#$=cut
#$
module igrad1 {	                                   # gradient in one dimension
#%  _lop( xx,  yy)
integer i
do i= 1, size(xx)-1 {
        if( adj) {
                xx(i+1) = xx(i+1) + yy(i)
                xx(i  ) = xx(i  ) - yy(i)
		}
	else
                yy(i) = yy(i) + xx(i+1) - xx(i)
	}
}