
% script for running "compare_FP_FPP_HISTORY.m"

% reichle, 23 Jan 2020

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

% plot surface type area fractions

startyear   = 2019;
endyear     = 2019;

startmonth  =  [ 10 ];
endmonth    =  [ 12 ];

plot_time   = {'all'};

plot_fields = {'FRLAND',       'const_2d_asm_Nx';    ...
	       'FRLAKE',       'const_2d_asm_Nx';    ...
	       'FRLANDICE',    'const_2d_asm_Nx';    ...
	       'FROCEAN',      'const_2d_asm_Nx';    ...
};

for ii=1:length(startmonth)

  for tt=1:length(plot_time)

   % compare_FP_FPP_HISTORY( startyear, startmonth(ii), endyear, endmonth(ii), plot_time{tt}, plot_fields )

  end
end

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

startyear   = 2019;
endyear     = 2019;

startmonth  =  [ 10; 11; 12; 10];
endmonth    =  [ 10; 11; 12; 12];

%startmonth  =  [ 10 ];
%endmonth    =  [ 10 ];

%plot_time   = {'all'}
plot_time =    {'430am','530am','630am','730am',...
                '430pm','530pm','630pm','730pm'};


plot_fields = {'TSURF',        'tavg1_2d_lnd_Nx';    ...
	       'TSOIL1',       'tavg1_2d_lnd_Nx';    ...
	       'TSOIL2',       'tavg1_2d_lnd_Nx';    ...
	       'TSOIL3',       'tavg1_2d_lnd_Nx';    ...
	       'TSOIL4',       'tavg1_2d_lnd_Nx';    ...
	       'GWETTOP',      'tavg1_2d_lnd_Nx';    ...
	       'SFMC',         'tavg1_2d_lnd_Nx';    ...
	       'PRECTOT',      'tavg1_2d_lnd_Nx';    ...
	       'TS',           'tavg1_2d_slv_Nx';    ...
	       'T2M',          'tavg1_2d_slv_Nx';    ...
	       'T10M'          'tavg1_2d_slv_Nx';    ...
	       'QV2M',         'tavg1_2d_slv_Nx';    ...
	       'U10M',         'tavg1_2d_slv_Nx';    ...
	       'V10M',         'tavg1_2d_slv_Nx';    ...
	       'PS',           'tavg1_2d_slv_Nx';    ...
};

%plot_fields = {'TSURF',        'tavg1_2d_lnd_Nx'}


for ii=1:length(startmonth)

  for tt=1:length(plot_time)

    compare_FP_FPP_HISTORY( startyear, startmonth(ii), endyear, endmonth(ii), plot_time{tt}, plot_fields )

  end
end




