subroutine grid_get_gradient (record,igridgrad)
subroutine grid_set_gradient (record,igridgrad)
implicit double precision (a-h,o-z)
igridgrad=1 means grid derivatives will be taken into account
in computing gradients;
igridgrad=0 means they will not.