7.4. Interannual pre-processingΒΆ

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


Scripts for the recovery of surface forcing data (based on CFSR reanalysis) for
interannual simulations


Scripts for the recovery of surface forcing data (based on ECMWF-ERAinterim simulations) for
interannual simulations


Scripts for the recovery of surface forcing data (based on ECMWF-ERA5 simulations) for
interannual simulations


Scripts for the recovery of surface forcing data (based on NCEP2 reanalysis) for
interannual simulations


Scripts for the recovery of wind stress from satellite scatterometer data (QuickSCAT)


Scripts for the generation of an operational oceanic forecast system


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
  1. 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


    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]
  2. 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



    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
  3. Check your generated files in CROCO_FILES

    You should have: