7.4. Interannual pre-processing¶

Dedicated scripts for interannual pre-processing can be found for the different forcing datasets in:

Aforc_CFSR	Scripts for the recovery of surface forcing data (based on CFSR reanalysis) for interannual simulations
Aforc_ECMWF	Scripts for the recovery of surface forcing data (based on ECMWF-ERAinterim simulations) for interannual simulations
Aforc_ERA5	Scripts for the recovery of surface forcing data (based on ECMWF-ERA5 simulations) for interannual simulations
Aforc_NCEP	Scripts for the recovery of surface forcing data (based on NCEP2 reanalysis) for interannual simulations
Aforc_QuikSCAT	Scripts for the recovery of wind stress from satellite scatterometer data (QuickSCAT)
Forecast_tools	Scripts for the generation of an operational oceanic forecast system
Oforc_OGCM	Scripts for the recovery of initial and lateral boundary conditions from global OGCMs (SODA (Carton et al., 2005), ECCO (Stammer et al., 1999) or CMEMS-GLORYS12) for inter-annual simulations

Edit crocotools_param.m

First section should already be set if you have completed the previous tutorial.

In the second section, check the path to forcing data directory:

% 2 - Generic file and directory names

%  Forcing data directory (ncep, quikscat, datasets download with opendap, etc..)
%
FORC_DATA_DIR = ['~/DATA/'];

In section 4, select only ini and bry (but no clim files, set: makeclim = 0;) to avoid too long pre-processing, and as it is the most usual set up:

%  initial/boundary data options (1 = process)
%  (used in make_clim, make_biol, make_bry,
%   make_OGCM.m and make_OGCM_frcst.m)
%
makeini    = 1;   % initial data
makeclim   = 0;   % climatological data (for boundaries and nudging layers)
makebry    = 1;  @ % lateral boundary data

Edit section 6 for running January to March 2005:

% 6 - Reference date and simulation times

Ymin          = 2005;          % first forcing year
Ymax          = 2005;          % last  forcing year
Mmin          = 1;             % first forcing month
Mmax          = 3;             % last  forcing month

Note

An important aspect is the definition of time and especially the choice of a time origin. The origin of time Yorig should be kept the same for all the preprocessing and postprocessing steps.

Edit section 7 for using CFSR and SODA forcing sets:

% 7 - Parameters for Interannual forcing (SODA, ECCO, CFSR, NCEP, ...)

Download_data = 0;   % Get data from OPENDAP sites
level         = 0;   % AGRIF level; 0 = parent grid
%
NCEP_version  = 3;   % NCEP version:
                     % [ CFSR up-to-date product are recommandated ]
                     %  1: NCEP/NCAR Reanalysis, 1/1/1948 - present
                     %  2: NCEP-DOE Reanalysis, 1/1/1979 - present
                     %  3: CFSR (Climate Forecast System Reanalysis),
                     %           1/1/1979 - 31/3/2011

NCEP_dir = [FORC_DATA_DIR,'CFSR_',CROCO_config,'/']; % CFSR data dir. [croco format]
makefrc      = 0;       % 1: create forcing files
makeblk      = 1;       % 1: create bulk files
QSCAT_blk    = 0;       % 1: a) correct NCEP frc/bulk files with
                        %        u,v,wspd fields from daily QSCAT data
                        %    b) download u,v,wspd in QSCAT frc file
add_tides    = 0;       % 1: add tides

% ...

OGCM        = 'SODA';        % Select the OGCM: SODA, ECCO
%
OGCM_dir    = [FORC_DATA_DIR,OGCM,'_',CROCO_config,'/'];  % OGCM data dir. [croco format]

Then you can run the Matlab pre-processing for these interannual forcing:
You should already have you grid set up. Otherwise, run make_grid

To build your interannual atmospheric forcing, the useful script is make_CFSR

To build your interannual ocean forcing, he useful script is make_OGCM
start make_CFSR make_OGCM
Warning
As this pluri-month preprocessing can be longer and uses more CPU ressources, you may need to submit it as a job. A few example scripts are provided:
cp ~/croco/croco_tools/example_job_prepro_matlab.pbs .
Launch your pre-processing job:
qsub example_job_prepro_matlab.pbs

Check your generated files in CROCO_FILES

You should have:

croco_blk_CFSR_Y????M?.nc
croco_bry_SODA_Y????M?.nc
croco_ini_SODA_Y????M?.nc