Making initial condition files for icosa_lmdz.exe from
an LMDZ start_archive.nc file.

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.

What you need to do/know before:
---------------------------
Run a dynamico simulation (with startphy_file=false for example) with the desired resolution to get
start.nc and startfi.nc. These files will be used as reference (rename as *_ref.nc) in the script.
You can use the files named *_no_start to have an example of files to run without startfiles

job_start_archive2icosa is the file you need to run start_archive2icosa on your platform (update 'arch' accordingly).

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.

This script is currently under development and still needs to be finalized.
Any suggestion for improvement is welcome.

In a nutshell:
--------------
1) input file start_archive.nc needs some minor adaptations:
- input file start_archive.nc need be converted to NetCDF4 format
nccopy -k 4 start_archive.nc start_archive_nc4.nc
2) script "compile" can be used to compile main tool start_archive2icosa.f90
   as well as rearrange_startphy.f90
3) to run, start_archive2icosa needs input files:
- "start_archive_nc4.nc" : from step 1
- "start_icosa_ref.nc" : a start.nc file from an icosa_lmdz.exe run (contains
  all the necessary information about the icosaedral grid at that resolution)
- "startphy_icosa_ref.nc" : a startfi.nc file from an icosa_lmdz.exe run (contains
  all the necessary information about colums from physics)
- "surface_nc4.nc" : topography file (also needs to be in NetCDF4 format)
NB: interpolation of HR topography requires lots of memory (see job file).
- "iodef.xml" : information about field names and grids for XIOS
4) start_archive2icosa outputs files "start_icosa_prefinalize.nc" and "startfi_prefinalize.nc"
  require a couple of extra modifications. Script "icosa_finalize.bash" does
  the required changes and generates files "start_icosa.nc" and "startfi.nc" that
  can be used as inputs by icosa_lmdz.exe
5) When using these start one must add in run_icosa.def the extra flags
etat0=start_file
etat0_start_file_colocated = true

Additional comments:
--------------------
- script "compile" assumes XIOS is to be found in "$LMDZ/XIOS",
  to be adapted to your own settings (set LMDZ variable or change script).
- "nccopy" is a standard netcdf tool available via e.g.:
module load intel/17.0 intelmpi/2017.0.098 hdf5/1.8.17 netcdf/4.4.0_fortran-4.4.2
- "nco" tools (ncrename, ncks, ...) are available via e.g.:
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
- in start_archive2icosa.f90, rearrange_startphy.f90, iodef.xml and
  icosa_finalize.bash input/output file names are hard-coded.
  Modify/adapt with caution!