igrad1 - 1D gradient operator
ierr=igrad1_lop(adj,add,xx,yy)
standard operator interface
Calculate 1-D gradient
the igrad2 manpage
geef90