Tools from NERSC to run the Ensemble Kalman Filter, adapted to NEMO.
|
5 months ago | |
---|---|---|
EnKF-MPI-TOPAZ | 5 months ago | |
conversion_uf | 5 months ago | |
interpolate_observations | 5 years ago | |
sanity_check | 3 years ago | |
README.md | 5 months ago | |
enkf_modules_to_load.txt | 5 months ago |
November 2016 Updated November 2017 Updated March 2019 at UCL Updated March 2024 A. Delhasse (alison.delhasse@uclouvain.be)
François Massonnet francois.massonnet@uclouvain.be/bsc.es
Everything available to set up the Ensemble Kalman filter (EnKF) for NEMO and EC-Earth
Go to the HPC
Clone the git repository (some HPCs don't allow communication with the web, in this case clone locally and rsync the git project to HPC)
git clone https://gogs.elic.ucl.ac.be/fmasson/EnKF.git
Compile the EnKF code
The folder EnKF-MPI-TOPAZ is a modified copy of the source code found at NERSC (http://enkf.nersc.no/Code/NERSC_implementation/V2/)
It has been modified as to work with NEMO and LIM, by F. Massonnet and C. König Beatty. The code itself should not be modified, unless there are good reasons to do so.
cd EnKF-MPI-TOPAZ
make clean
configure_ecearth
make
An executable EnKF should be produced. Try to run it to see if you get the following information message:
bsc32526@login2:/gpfs/projects/bsc32/bsc32526/enkf/EnKF-assimilation/EnKF-MPI-TOPAZ> ./EnKF Inne i start_mpi: qmpi_proc_num = 0 master = T MPI started with 1 processors Usage: EnKF
EnKF -h
Options:
-h -- describe parameter fie format
cd ../conversion_uf/ make clean source ../enkf_modules_to_load.txt make
Two executables named prep_obs_ORCA1 and prep_obs_ORCA25 should be created. Try them:
./prep_obs_ORCA1
0
(prep_obs) takes a real obs, extracts the desired variable and outputs it in a format that the EnKF can read & treat ('observations.uf').
A file named mask.nc containing the variables tmaskutil, nav_lon and nav_lat is expected to be in the current directory (ORCA-file)
Three command line arguments are expected:
Hope to see you again soon.
(prep_obs): Stopped.
The folder sanity_check is a set of routines that make sure that the restarts provided by the EnKF can be read by NEMO. These are home-made (François Massonnet, see also description in http://www.sciencedirect.com/science/article/pii/S1463500315000050
cd ../sanity_check/ make clean source ../enkf_modules_to_load.txt make
An executable named sanity_checked is created, try it:
./sanity_check
sanity_check_LIM3 needs arguments: -analysis_file_ice -forecast_file_ice -analysis_file_oce -forecast_file_oce Checks NEMO-LIM3 ice and ocean analyses restarst (netcdf) file for sanity and fixesthem if necessary.
Sanity means for now: Strongly follow limupdate.F90 Files mask.nc and mesh_hgr.nc need to be in the current directory
Hope to see you again soon.
Chris König Beatty Francois Massonnet -- francois.massonnet@uclouvain.be Last update: 2013 Last update: 2016 (to work with NEMO3.6) (sanity_check): Stopped.
cd ../conversion_uf/
open convert.bash, check it, and run it (always on HPC)