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]
- Edit
- 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
- Check your generated files in