clc 
clear all
close all

  
  run_months = 1:12;
  
  exp_path = '/discover/nobackup/jkolassa/AMSRE_ASCAT_DA_experiment';
  
  exp_run    ='AMSRE_ASCAT_DA_noDA_ASCAT';
  
  domain   = 'SMAP_EASEv2_M09_CONUS'; 
  
  start_year = 2007;
  end_year   = 2011; 
  
  dt_assim = 10800;
 
  species = [ 1 ];  

  varname='sfmc'; 
  
  % radius of influence of each observation
  
  hscale = 1.25;          % degrees lat/lon
   


% ------------------

   for s=1
  
% this does each month separately - seems to be the default in GDL and QL code

      %for k=1:length(run_months)
	
%	get_model_and_obs_stats(                             ...
%	    varname, run_months(k), exp_path, exp_run, domain,     ...
%	    start_year, end_year, ...
%	    dt_assim, species, hscale, time_of_day_in_hours(j)  );
	
      %end


% run for all months
	get_model_and_obs_stats_mod(                             ...
	    varname, run_months, exp_path, exp_run, domain,     ...
	    start_year, end_year, ...
	    dt_assim, species(s), hscale  );

         get_cdf_match( ...
                varname, run_months, exp_path, exp_run, domain, start_year, end_year, ...
                species(s), hscale,0);


 end

% end  % yr

% ============= EOF ====================================================
