
% run_get_model_and_obs_stats.m
%
% reichle, 17 Aug 2005
% reichle, 19 Feb 2008

clear

%for yr=2002:2011
  
  run_months = 1:12;
  
  exp_path = '/discover/nobackup/csdraper/LDASsa_m3_out/';
  
  exp_run    ='M09_ARS_fMERRA_openA-CPCU';
  
  domain   = 'SMAP_EASEv2_M09'; 
  
  start_year = 2002; % 2002;
  end_year   = 2011; % 2005;
%  start_year=yr; 
%  end_year=yr; 

  dt_assim = 10800;
 
  species = [ 1 ];  

  varname='sfmc'; 
  
  % radius of influence of each observation
  
  hscale = 0;          % degrees lat/lon
   
 % time_of_day_in_hours = 0:3:21;

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

   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(                             ...
	    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 ====================================================

