%------------------------------------------------------------------
% STEP 2 - extraction
%------------------------------------------------------------------
% RUN EXTRACTION ROUTINES
% Note: these routine require a proper installation of
% RNT, RNC, RGRD, NETCDF and OPENDAP
% If you do not have this, just send the parameters of the selection
% to E. Di Lorenzo, globec_extract@o3d.org and we will run it for you.

% Load selection parameters
nepd_select_MAIN

% Load NEPD grid.
grd=rnt_gridload('nepd_conf.txt');  

% Grid Extraction Routine
[grdfile, grdn, gridconf] = nepd_ex_grid(sub_I, sub_J, grd, nameit, res)

% Forcing Extraction Routine
forcfiles=nepd_ex_forc(sub_I, sub_J, grd, grdn, nameit, res);

% Initial Condition Extraction Routine
[inifile, ctli]=nepd_ex_init(sub_I, sub_J, grd, grdn, nameit, res, inittime, npzd, np)

% Boundary Condition Extraction Routine (this function requires some time)
[bryfile, ctl]=nepd_ex_bry(sub_I, sub_J, grd, grdn, nameit, res, trange, npzd, np)

disp('done.');

return
% Additional stuff
% Determine DSTART
DSTART=num2str(ctli.time/24/60/60);
disp(['DSTART = ',DSTART]);

% Make figure of subdomain
clf; rnt_plcm(grd.h,grd); hold on; rnt_gridboxm(grdn);