igrad2 - 2D gradient operator
ierr=igrad2_init(n1,n2)
ierr=igrad2_lop(adj,add,xx,yy)
dimensions of dataset
standard operator interface
Calculate 2-D gradient
the igrad1 manpage
geef90