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

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

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