

gold = '/Users/thoar/git/DART_Manhattan/models/lorenz_96/work';
test = '../work';


cases = {'inf00','inf20','inf22','inf33','inf24','inf55'};

for icase = 1:length(cases)

   fprintf('Comparing inflation case %s \n',cases{icase})
   fnamegold = sprintf('%s/%s/filter_output.nc',gold,cases{icase});
   fnametest = sprintf('%s/%s/filter_output.nc',test,cases{icase});

   Compare_netCDF_files(fnamegold,fnametest)

   fprintf('\n')

end

