util1_jan/grad.m

   1 function [dfdx,dfdy]=grad(x,y,f)
   2 % compute directional central derivatives of 2D field
   3 % in:
   4 %    x, y   curvilinear rectangular grid of x,y nodal coordinates
   5 %    f      values of function on the nodes
   6 % out:
   7 %    dfdx,dfdy   derivatives of f along the x,y lines by central differences
   8
   9 dxx=  x(3:end, 2:end-1)-x(1:end-2, 2:end-1);
  10 dyx=  y(3:end, 2:end-1)-y(1:end-2, 2:end-1);
  11 dxn=  sqrt(dxx.^2+dyx.^2);
  12 dfdx=(f(3:end, 2:end-1)-f(1:end-2, 2:end-1))./dxn;
  13
  14 dxy=  x(2:end-1, 3:end)-x(2:end-1, 1:end-2);
  15 dyy=  y(2:end-1, 3:end)-y(2:end-1, 1:end-2);
  16 dyn=  sqrt(dxy.^2+dyy.^2);
  17 dfdy=(f(2:end-1, 3:end)-f(2:end-1, 1:end-2))./dyn;
  18
  19 end