quicwind/mat_gen_wind_flux_div.m

   1 function DM=mat_gen_wind_flux_div(X,type)
   2 % get matrix of flux divergence
   3 % input:
   4 %     X   mesh
   5 %     type 'D' create divergence matrix D
   6 %          'M' create wind to flux matrix M
   7 %          'DM' create D*M (default)
   8 % output:
   9 %     DM  flux divergence matrix
  10 %
  11 % Error checking
  12 check_mesh(X);
  13
  14 if ~exist('type','var')
  15     type = 'DM'
  16 end
  17
  18 % Generate wind flux matrix
  19 Mmat = mat_flux(X);
  20
  21 % Generate divergence matrix
  22 Dmat = mat_div3(X);
  23
  24 % Generate flux divergence matrix
  25 if (type == 'M')
  26     DM = Mmat;
  27 elseif (type == 'D')
  28     DM = Dmat;
  29 else
  30     DM = Dmat * Mmat;
  31 end
  32
  33 end