femwind/f_assembly_wrfout_test.m

   1 format compact
   2 % test for ndt_f_assembly from wrfout file
   3 path='./wrfout_d01_0001-01-01_00:00:00' % path to wrfout file
   4 [X,u0]=read_fmw_wrfout(path);
   5 A=diag([1 1 1]);
   6 lambda=[];
   7 params=[];
   8 iflags=[2 1 1];
   9 iflags=iflags(:);
  10
  11 % test same results for ndt_mult from matlab and fortran
  12 if exist('fortran/f_assembly_test.exe')
  13     disp('testing if same result as fortran')
  14     err=f_assembly_fortran(A,X,u0,[],[]);
  15     if abs(err)<1e-6
  16     fprintf('error %g OK\n',err)
  17     else
  18     error(sprintf('error %g too large',err))
  19     end
  20 else
  21     warning('fortran/f_assembly_test.exe not available')
  22 end