| 1 | Making initial condition files for icosa_lmdz.exe from |
|---|
| 2 | an LMDZ start_archive.nc file. |
|---|
| 3 | |
|---|
| 4 | This tutorial should follow https://lmdz-forge.lmd.jussieu.fr/mediawiki/Planets/index.php/LMDZ_to_DYNAMICO_start_files. If it does not, please update it. |
|---|
| 5 | |
|---|
| 6 | What you need to do/know before: |
|---|
| 7 | --------------------------- |
|---|
| 8 | Run a dynamico simulation (with startphy_file=false for example) with the desired resolution to get |
|---|
| 9 | start.nc and startfi.nc. These files will be used as reference (rename as *_ref.nc) in the script. |
|---|
| 10 | You can use the files named *_no_start to have an example of files to run without startfiles |
|---|
| 11 | |
|---|
| 12 | job_start_archive2icosa is the file you need to run start_archive2icosa on your platform (update 'arch' accordingly). |
|---|
| 13 | |
|---|
| 14 | Warning: the input start_archive.nc file must include SSO parameters (albedodat, ZMEA, etc) in this version, which is only available with a modified start2archive. Be sure to check if that has been updated or not. |
|---|
| 15 | |
|---|
| 16 | This script is currently under development and still needs to be finalized. |
|---|
| 17 | Any suggestion for improvement is welcome. |
|---|
| 18 | |
|---|
| 19 | In a nutshell: |
|---|
| 20 | -------------- |
|---|
| 21 | 1) input file start_archive.nc needs some minor adaptations: |
|---|
| 22 | - input file start_archive.nc need be converted to NetCDF4 format |
|---|
| 23 | nccopy -k 4 start_archive.nc start_archive_nc4.nc |
|---|
| 24 | 2) script "compile" can be used to compile main tool start_archive2icosa.f90 |
|---|
| 25 | as well as rearrange_startphy.f90 |
|---|
| 26 | 3) to run, start_archive2icosa needs input files: |
|---|
| 27 | - "start_archive_nc4.nc" : from step 1 |
|---|
| 28 | - "start_icosa_ref.nc" : a start.nc file from an icosa_lmdz.exe run (contains |
|---|
| 29 | all the necessary information about the icosaedral grid at that resolution) |
|---|
| 30 | - "startphy_icosa_ref.nc" : a startfi.nc file from an icosa_lmdz.exe run (contains |
|---|
| 31 | all the necessary information about colums from physics) |
|---|
| 32 | - "surface_nc4.nc" : topography file (also needs to be in NetCDF4 format) |
|---|
| 33 | NB: interpolation of HR topography requires lots of memory (see job file). |
|---|
| 34 | - "iodef.xml" : information about field names and grids for XIOS |
|---|
| 35 | 4) start_archive2icosa outputs files "start_icosa_prefinalize.nc" and "startfi_prefinalize.nc" |
|---|
| 36 | require a couple of extra modifications. Script "icosa_finalize.bash" does |
|---|
| 37 | the required changes and generates files "start_icosa.nc" and "startfi.nc" that |
|---|
| 38 | can be used as inputs by icosa_lmdz.exe |
|---|
| 39 | 5) When using these start one must add in run_icosa.def the extra flags |
|---|
| 40 | etat0=start_file |
|---|
| 41 | etat0_start_file_colocated = true |
|---|
| 42 | |
|---|
| 43 | Additional comments: |
|---|
| 44 | -------------------- |
|---|
| 45 | - script "compile" assumes XIOS is to be found in "$LMDZ/XIOS", |
|---|
| 46 | to be adapted to your own settings (set LMDZ variable or change script). |
|---|
| 47 | - "nccopy" is a standard netcdf tool available via e.g.: |
|---|
| 48 | module load intel/17.0 intelmpi/2017.0.098 hdf5/1.8.17 netcdf/4.4.0_fortran-4.4.2 |
|---|
| 49 | - "nco" tools (ncrename, ncks, ...) are available via e.g.: |
|---|
| 50 | module load intel/17.0 openmpi/intel/2.0.1 hdf5/1.8.17 netcdf/4.4.0_fortran-4.4.2 nco/4.4.6 |
|---|
| 51 | - in start_archive2icosa.f90, rearrange_startphy.f90, iodef.xml and |
|---|
| 52 | icosa_finalize.bash input/output file names are hard-coded. |
|---|
| 53 | Modify/adapt with caution! |
|---|
