femwind/femwind_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 addpath ../netcdf
  11
  12 % test same results for ndt_mult from matlab and fortran
  13 % if exist('fortran/f_assembly_test.exe')
  14     disp('testing if same result as fortran')
  15     p=params_defaults;
  16     p.run_fortran=0;
  17     p.run_matlab=1;
  18     p.femwind_fortran_test=0
  19     err=femwind_fortran(A,X,u0,p);
  20     if abs(err)<1e-6
  21     fprintf('error %g OK\n',err)
  22     else
  23     error(sprintf('error %g too large',err))
  24     end
  25 %else
  26 %    warning('fortran/f_assembly_test.exe not available')
  27 %end