.gitignore

 conf/profile_mesonh
 conf/profile_mesonh-*
 exe/*
+MY_RUN/INTEGRATION_CASES/**/*.des
+MY_RUN/INTEGRATION_CASES/**/*.dir
+MY_RUN/INTEGRATION_CASES/**/*.hdr
+MY_RUN/INTEGRATION_CASES/**/*.nc
+MY_RUN/INTEGRATION_CASES/**/ecmwf.OD.????????.??
+MY_RUN/INTEGRATION_CASES/**/file_for_xtransfer
+MY_RUN/INTEGRATION_CASES/**/OUTPUT_LISTING*
+MY_RUN/INTEGRATION_CASES/**/PRESSURE
 MY_RUN/KTEST/*/*/*.des
 MY_RUN/KTEST/*/*/*.lfi
 MY_RUN/KTEST/*/*/*.nc
@@ -36,10 +44,25 @@ MY_RUN/KTEST/009_ICARTT/001_pgd1/*.asc
 MY_RUN/KTEST/009_ICARTT/002_arp2lfi/ecmwf.OD.20040810.18-V2
 MY_RUN/KTEST/009_ICARTT/002_arp2lfi/mocage.GLOB22.20040810.18
 MY_RUN/KTEST/9??_*
-pub/ncl-?.?.?
-pub/ncl_ncarg*/
+pub/FILEPP/filepp*
+!pub/FILEPP/filepp*.tar.gz
+pub/FILEPP/FILEPP*
+pub/FILEPP/MNH_Expand_*
+!pub/FILEPP/MNH_EXPAND_*.tar.gz
 src/dir_obj-*
+src/LIB/eccodes*
+!src/LIB/eccodes*.tar.gz
 src/LIB/grib_api*
+!src/LIB/grib_api*.tar.gz
 src/LIB/hdf5*
+!src/LIB/hdf5*.tar.gz
 src/LIB/libaec*
+!src/LIB/libaec*.tar.gz
 src/LIB/netcdf*
+!src/LIB/netcdf*.tar.gz
+src/LIB/MEGAN*
+src/LIB/oasis3-mct*
+!src/LIB/oasis3-mct*.tar.gz
+src/LIB/toy*
+!src/LIB/toy*.tar.gz
+vscode_mod

A-INSTALL

 #
 # Version of PACKAGE MESONH "Open distribution"
-#              PACK-MNH-V5-4-2
-#              DATE : 03/04/2019
-#              VERSION : MESONH MASDEV5_4 + BUG-2
+#              PACK-MNH-V5-7-1
+#              DATE : 04/09/2024
+#              VERSION : MESONH MASDEV5_7 + BUG-1
 #
 #  MAP
 #
@@ -13,7 +13,7 @@
 #      III) COMPILING/INSTALLING THE MESONH PACKAGE ON YOUR LINUX COMPUTER 
 #       IV) RUN SOME "SMALL KTEST" EXAMPLES
 #        V) RECOMPILING YOUR 'OWN' SOURCES ONLY 
-#       VI) COMPILING/INSTALLING MESONH ON GENCI & ECMWF & METEO COMPUTERS
+#       VI) COMPILING/INSTALLING MESONH ON GENCI & ECMWF & METEO & CALMIP COMPUTERS
 #      VII) "SCANDOLLAR" = SCRIPTING YOUR OWN PROCEDURES 
 #
 #     VIII) TROUBLE-SHOOTING
@@ -28,6 +28,10 @@
 #           e) cleaning previous compiled version
 #
 #
+# NEW since MNH-56X : For conda package for graphic output in test-case examples , read
+# the 'README_MNH_CONDA' file in this (root) directory
+#
+#
 ^L
 #  0) TWO WAYS OF DOWNLOADING MESONH
 #  =================================
@@ -80,14 +84,14 @@
 #
 # or directly
 #
-# http://mesonh.aero.obs-mip.fr/mesonh/dir_open/dir_MESONH/MNH-V5-4-2.tar.gz
+# http://mesonh.aero.obs-mip.fr/mesonh/dir_open/dir_MESONH/MNH-V5-7-1.tar.gz
 #
-# Then untar the file "MNH-V5-4-2.tar.gz" where you want to.
+# Then untar the file "MNH-V5-7-1.tar.gz" where you want to.
 # For example, in your home directory:
 #
 
 cd ~
-tar xvfz MNH-V5-4-2.tar.gz
+tar xvfz MNH-V5-7-1.tar.gz
 
 #
 #  Process now to the chapter to configure the MesoNH package.
@@ -130,7 +134,7 @@ git lfs install
 #    only access) by following the next link:
 #
 
-http://mesonh.aero.obs-mip.fr/mesonh54/GitSources?action=AttachFile&do=get&target=anongitmesonh.key
+http://mesonh.aero.obs-mip.fr/mesonh57/GitSources?action=AttachFile&do=get&target=anongitmesonh.key
 
 #
 #    and save the file in your $HOME/.ssh/ directory.
@@ -161,68 +165,68 @@ git config --global http.sslverify false
 # This is necessary to disable the certificate checks because a self-signed
 # certificate was used for the LFS server.
 #
-#    c) Cloning the Meso-NH Source repository on the developpement branch MNH-54-branch
+#    c) Cloning the Meso-NH Source repository on the developpement branch MNH-57-branch
 #    ----------------------------------------------------------------------------------
 #
 # Finally you can clone the Meso-NH Git repository with the following command:
 #
 
-git lfs clone anongit@anongit_mesonh:/gitrepos/MNH-git_open_source-lfs.git -b MNH-54-branch MNH-V5-4-2
+git clone anongit@anongit_mesonh:/gitrepos/MNH-git_open_source-lfs.git -b MNH-57-branch MNH-V5-7-1
 
 #
-# that will create the MNH-V5-4-2 directory containing a clone (copy) of the
-# Meso-NH package on the remote developpement branch MNH-54-branch
+# that will create the MNH-V5-7-1 directory containing a clone (copy) of the
+# Meso-NH package on the remote developpement branch MNH-57-branch
 #
 #
 #   d) Checking out a given version of MESONH
 #   -----------------------------------------
 #
 # Once the repository is cloned, it's better for you to checkout your own branch
-# (by default, you are on HEAD of the MNH-54-branch development branch  ).
+# (by default, you are on HEAD of the MNH-57-branch development branch  ).
 #
-# To create your local branch corresponding to the V5-4-2 version, type:
+# To create your local branch corresponding to the V5-7-1 version, type:
 #
 
-cd MNH-V5-4-2
-git checkout -b MYB-MNH-V5-4-2 PACK-MNH-V5-4-2
+cd MNH-V5-7-1
+git checkout -b MYB-MNH-V5-7-1 PACK-MNH-V5-7-1
 
 #
-# MYB-MNH-V5-4-2 is the name of the local branch you created
+# MYB-MNH-V5-7-1 is the name of the local branch you created
 # and
-# PACK-MNH-V5-4-2 is the remote/origin tag on which it is based.
+# PACK-MNH-V5-7-1 is the remote/origin tag on which it is based.
 #
 # The advantage of this way of downloading the package is that in the future
 # you could check/update quickly differences with the new version of the
 # package without having to download entirely the full package.
 #
-#  Suppose that a new version, for example "PACK-MNH-V5-4-2", is announced.
+#  Suppose that a new version, for example "PACK-MNH-V5-7-1", is announced.
 #
 #  To see the differences with your working copy, do:
 #
 
 git fetch
-git diff HEAD PACK-MNH-V5-4-2
+git diff HEAD PACK-MNH-V5-7-1
 
 #
 #  To go to the new version, you can, for example, create a new local branch:
 #
 
-git checkout -b MYB-MNH-V5-4-2 PACK-MNH-V5-4-2
+git checkout -b MYB-MNH-V5-7-1 PACK-MNH-V5-7-1
 
 #
 # At any time, you can also check for "uptodate" changes in the Git branch
-# dedicated to the MNH54 version before the official release of the "bugN+1"
+# dedicated to the MNH57 version before the official release of the "bugN+1"
 # bugfix version.
 #
 
 git fetch
-git diff HEAD MNH-54-branch
+git diff HEAD MNH-57-branch
 
 #
 # And, test this development (not yet official) version by going to this branch:
 #
 
-git checkout --track origin/MNH-54-branch
+git checkout --track origin/MNH-57-branch
 #
 #    e) Cloning the Meso-NH Documentation repository
 #
@@ -253,7 +257,7 @@ git clone anongit@anongit_mesonh:/gitrepos/MNH-DOC.git
 #    use the "./configure" script like this 
 #
 
-cd ~/MNH-V5-4-2/src
+cd ~/MNH-V5-7-1/src
 ./configure
 . ../conf/profile_mesonh
 
@@ -304,7 +308,7 @@ export OPTLEVEL=O2         # Compile in O2, 4 times faster then DEBUG, but less
 
 #    and then source/load the new generate file
 
-. ../conf/profile_mesonh.LXifort.MNH-V5-4-2.MPIAUTO.O2
+. ../conf/profile_mesonh.LXifort.MNH-V5-7-1.MPIAUTO.O2
 
 #
 #  REM: 
@@ -329,7 +333,7 @@ export OPTLEVEL=O2         # Compile in O2, 4 times faster then DEBUG, but less
 # go to the directory  "src"
 #
 
-cd  ~/MNH-V5-4-2/src
+cd  ~/MNH-V5-7-1/src
 
 #
 # if you have not already configured your MESONH environment
@@ -535,8 +539,8 @@ make examples
 
 ^L
 #
-#   VI) COMPILING/INSTALLING ON GENCI & ECMWF & METEO COMPUTERS
-#   ===========================================================
+#   VI) COMPILING/INSTALLING ON GENCI & ECMWF & METEO & CALMIP COMPUTERS
+#   ====================================================================
 #
 #   After downloading "exactly" like on "any standalone PC" 
 #   run the "./configure" command :
@@ -547,112 +551,166 @@ make examples
 #
 #  - If you do not have sufficient space in your "$HOME" directory,
 #    install the whole package directly on the ${WORKDIR}
+#    /!\  the name of the WORKDIR differ in the differents computer center
+#    read the doc , most of them manage disk space throw 'multi-projet' with only one unique login . 
 #
 #   WARNING :
 #   ======== 
 #    - Think to do a backup of your installation 
-#    - This space is no more "purged" but a "crash disk" could/will
+#    - This space is no "purged" but a "crash disk" could/will
 #      probably occur !!!
 # 
 
 cd $WORKDIR
-cd MNH-V5-4-2/src
+cd MNH-V5-7-1/src
 
 ./configure
  
 #   
 #   Due to limitation in time & memory on interactive connection 
-#   then compile the MESONH PACKAGE in batch mode with the different "job_make_mesonh*" files
+#   in some computer you have to compile the MESONH PACKAGE in batch mode with the different "job_make_mesonh*" files
 #   
 #   at IDRIS :
 #   ---------
 #
 
-#   - On ADA ( IBM/x3750 ) with :
+#   - On JEAN-ZAY ( HPE ) the compilation is in interactive :
 
-llsubmit job_make_mesonh_IBM_ada
+cd MNH-V5-7-1/src
+. ../conf/profile_mesonh-LXifort-R8I4-MNH-V5-7-1-MPIINTEL-O2
+make -j16 |& tee error$XYZ
+make installmaster
 
-#
-#
-#   - On TURING (IBM BG/Q ) : 
+# You could also use the 'compil' partition
 
-llsubmit job_make_mesonh_BGQ
+sbatch job_make_mesonh_HPE_jeanzay
+
+#   - to run the test case examples  run
+
+sbatch -A {your_projet}@cpu job_make_examples_BullX_jeanzay
 
-#
-#
-#   BACKUP : /!\ backup your work on "ergon"
 #
 #   
-#   at CINES on JADE( SGI/ICE ) , OCCIGEN not tested yet  :
+#   at CINES on ADASTRA (BULLX)  :
 #   -------------------------------------------------------
 # 
-#   - install the PACKAGE in your "/store/${USER}" directory 
-#   - you could compile in interactive mode
+#   - install the PACKAGE in your $HOME ( default 50Go  of quota )
+#   - Compile in interactive mode ( see IDRIS )
 #
-#    BACKUP : /!\ backup your work on "/data/${USER}"
-# 
+#   - to run the test case examples  run
 
+sbatch job_make_examples_BullX_occigen
 
+#
+#   
+#   at TGCC on IRENE (BULLX)  :
+#   -------------------------------------------------------
+#
+#   At TGCC , you have two architectures accessible throw 2 differents frontals
+#   but with a commun disk space , connect to : 
+#
+#   - ssh irene-fr  : for Intel SkyLake processors
+#      On Intel processors the MPI use is OPENMPI/4.1.4
+#      the configure will generate a
+#      profile_mesonh-LXifort-R8I4-MNH-V5-7-1-MPIAUTO-O2
+#
+#   - ssh irene-amd : for AMD , processors
+#      On AMD processors the MPI use is OPENMPI/4.1.4
+#      the configure will generate a
+#      profile_mesonh-LXifort-R8I4-MNH-V5-7-1-AMD-MPIAUTO-O2
 # 
-#   at ECMWF on cca ( CRAY/XC30 ) :
-#   ----------------------------
 #
-#   - to install MESONH go to your $PERM directory
-#     if you need more disk space than allowed for 'standard' user 
-#     Ask to Dominique Lucas, look this email :
+#   At TGCC they use 'one login' for multi-project allocation .
+#   This induce 'strange' problem with the installation of eccodes
+#   resulting in file with the wrong default group and 'disk quota excedeed' error
 #
-#    http://mesonh.aero.obs-mip.fr/cgi-bin/mesonh_interne/mail2html.cgi?file=2009_09_24_15:31:54
+#   You could get your different project info by 'ccc_myproject'
 #
-#   - for the compilation use :
-
-qsub job_make_mesonh_CRAY_cca
+#   I recommand you first, supposing that the "Genci Allocation" you want to use
+#   as for projet group "genXXXX"
+#
+#   1) to add  in your "~/.bash_profile"
+#
+module switch dfldatadir dfldatadir/genXXXX
+newgrp genXXXX
 
+#   2) "logout & login" again
+#      and check that you have now the good default group
 #
-# REM : On CRAY computer cca , the default compiler is the CRAY one = crayftn .
-#      I you have trouble with this one , you could compile with the Intel/ifort one .
-#      Use the ifort one do :
+id -ng
+--> genXXXX
+echo $CCCHOME
+--> /ccc/.../home/genXXXX/{your_login}
 
-export ARCH=LXifort
+#  If all is OK, you can install the code interactively as usual
+#
+#   - install the PACKAGE in your ${CCCHOME} ( default 20Go  of quota )
+#  for example on Intel frontal
+cd ${CCCHOME}
+cd MNH-V5-7-1/src
 ./configure
+. ../conf/profile_mesonh-LXifort-R8I4-MNH-V5-7-1-MPIAUTO-O2
+make -j16 |& tee error$XYZ
+make installmaster
 
-...
+# REM: The eccode lib will be installed in your ${CCCWORKDIR}
+#      also for the quota problem ...
+#    
+
+#   - to run the test case examples  run
+#
+# On intel Skylake
+ccc_msub job_make_examples_BullX_irene
+
+# On intel AMD 
+ccc_msub job_make_examples_BullX_irene_AMD
+
+# 
+#   at ECMWF on hpc-login ( ATos/HPCF ) :
+#   ------------------------------------------
+#
+#   - to install MESONH go to your $HPCPERM directory
+#
+#   - for the compilation :
+#   after the ./configure in the login node
+#
+#   connecte to an "interactive compute node"  via ( 16 core & 16GO of memory 
 
-creation du fichier -->  ../conf/profile_mesonh-LXifortI4-MNH-V5-4-2-MPICRAY-O2
+ecinteractive -c16 -m 16G -t 12:00:00
+
+etc ...
 
-# And for the compilation & example job , switch the ARCH variable to LXifort :
 
-vi  job_make_mesonh_CRAY_cca(job_make_examples_CRAY_cca)
+#   - to run the test case examples  run
 
-ARCH=LXifort
-#ARCH=LXcray # this is the default one
-. ../conf/profile_mesonh-${ARCH}I4-MNH-V5-4-2-MPICRAY-O2
+sbatch  job_make_examples_Atos_HPCF
 
 #
-#   - At Meteo-France DSI on beaufix (or prolix)
+#   - At Meteo-France DSI on belenos
 #
 #  to install the whole package on your "$HOME" directory
-#  untar the file "PACK-MNH-V5-4-2.tar.gz" from its location :
+#  untar the file "MNH-V5-7-1.tar.gz" from its location :
 
 cd ~
-tar xvf $MESONH/pack_MNH-V5-4-2.tar.gz
+tar xvf $MESONH/MNH-V5-7-1.tar.gz
 
 #  run the "./configure" command :
 
-cd MNH-V5-4-2/src
+cd MNH-V5-7-1/src
 ./configure
 
 #
 #   Due to limitation in time & memory on interactive connection
-#   then compile the MESONH PACKAGE in batch mode with the job_make_mesonh_BullX file :
+#   then compile the MESONH PACKAGE in batch mode with the job_make_mesonh_BullX_belenos file :
 
-sbatch   job_make_mesonh_BullX
+sbatch   job_make_mesonh_BullX_belenos
 
 #   This job does : gmake -j 4
 #   then : make installmaster
 
 #   To run basic KTEST examples :
 
-sbatch   job_make_examples_BullX
+sbatch   job_make_examples_BullX_belenos
 
 #   Step-2 : configure/compiling with VER_USER=...
 # ----------------------------------------------
@@ -669,46 +727,19 @@ export VER_USER=MY_MODIF
 #   in job_make_mesonh_user_BullX insert " export VER_USER=MY_MODIF "
 #   then submit in batch mode
 
-sbatch  job_make_mesonh_user_BullX
-
-#
-#   EXAMPLES ON GENCI & ECMWF PLATFORMS & METEO COMPUTERS
-#   =====================================================
-#
-#  - At IDRIS
-#
-
-# on ADA
-
-llsubmit job_make_examples_IBM_ada
-
-#on TURING
-
-llsubmit job_make_examples_BGQ
+sbatch  job_make_mesonh_user_BullX_belenos
 
 #
-#   - At CINES on JADE( SGI/ICE ) :
-#
-# load ".profile_mesonh..." & set the variables :
- 
-export MONORUN="mpirun -np 1"
-export MPIRUN="mpirun -np 4 " 
-  
-# and run with 
-
-make examples
-
-#
-#   - At ECMWF on cca
-#
-
-llsubmit  job_make_examples_CRAY_cca
-
-#
-#   - At Meteo-France DSI on beaufix or prolix
+#   
+#   at CALMIP on OLYMPE (BULLX)  :
+#   -------------------------------------------------------
+# 
+#   - install the PACKAGE in your /tmpdir/$USER
+#   - Compile in interactive mode 
 #
+#   - to run the test case examples  run
 
-sbatch   job_make_examples_BullX
+sbatch job_make_examples_BullX_olympe
 
 #
 # That's all for the basic INSTALLATION of the "MESONH PACKAGE"
@@ -718,7 +749,7 @@ sbatch   job_make_examples_BullX
 #  VII) "SCANDOLLAR" = SCRIPTING YOUR OWN PROCEDURES 
 #  ====================================================
 
-REM : not tested on 5-4-X version ...
+REM : not tested on 5-5-X version ...
 
 # A small script 'scandollar' is in test in this version of MESONH.
 # It will help you to build a complete chaining of an experiment. 
@@ -764,7 +795,7 @@ scandollar
 
 ## OUTPUT ::
 
->#  read default config file :: --->  CONF_DOLLAR=/home/escj/DEV64/PACK-MNH-V5-4-2/conf/post/confdollar_aeropc_default
+>#  read default config file :: --->  CONF_DOLLAR=/home/escj/DEV64/PACK-MNH-V5-7-1/conf/post/confdollar_aeropc_default
 >#
 >#  read    user config file :: --->  CONFIG=confdollar
 >#
@@ -786,7 +817,7 @@ scandollar 0*
 
 ## OUTPUT ::
 >#
->#  read default config file :: --->  CONF_DOLLAR=/home/escj/DEV64/PACK-MNH-V5-4-2/conf/post/confdollar_aeropc_default
+>#  read default config file :: --->  CONF_DOLLAR=/home/escj/DEV64/PACK-MNH-V5-7-1/conf/post/confdollar_aeropc_default
 >#
 >#  read    user config file :: --->  CONFIG=confdollar
 >#
@@ -860,22 +891,22 @@ cp -R  007_16janvier_scandollar /.../your_directory
 #
 # use this "profile_mesonh" :
 
-.  /home/rech/mnh/rmnh007/DEV/MNH-V5-4-2/conf/profile_mesonh-SX8-MNH-V5-4-2-MPIAUTO-O4
+.  /home/rech/mnh/rmnh007/DEV/MNH-V5-7-1/conf/profile_mesonh-SX8-MNH-V5-7-1-MPIAUTO-O4
 
 # And the examples are here ( link to my $WORKDIR in actually )  
 
-/home/rech/mnh/rmnh007/DEV/MNH-V5-4-2/MY_RUN/KTEST/007_16janvier_scandollar
+/home/rech/mnh/rmnh007/DEV/MNH-V5-7-1/MY_RUN/KTEST/007_16janvier_scandollar
 
 #
 #  On vargas
 #  ---------
 #  use this "profile_mesonh" :
 
-. /workgpfs/rech/mnh/rmnh007/DEV/MNH-V5-4-2/conf/profile_mesonh-AIX64-MNH-V5-4-2-MPIAUTO-O2
+. /workgpfs/rech/mnh/rmnh007/DEV/MNH-V5-7-1/conf/profile_mesonh-AIX64-MNH-V5-7-1-MPIAUTO-O2
 
 # and examples here :
 
-/workgpfs/rech/mnh/rmnh007/DEV/MNH-V5-4-2/MY_RUN/KTEST/007_16janvier_scandollar
+/workgpfs/rech/mnh/rmnh007/DEV/MNH-V5-7-1/MY_RUN/KTEST/007_16janvier_scandollar
 
 # 
 #  - At CINES on JADE :
@@ -883,11 +914,11 @@ cp -R  007_16janvier_scandollar /.../your_directory
 #
 # use 
 
-. /work/escobar/DEV/MNH-V5-4-2/conf/profile_mesonh-LXifort-MNH-V5-4-2-MPIICE-O2
+. /work/escobar/DEV/MNH-V5-7-1/conf/profile_mesonh-LXifort-MNH-V5-7-1-MPIICE-O2
 
 # and the exemples
 
-/work/escobar/DEV/MNH-V5-4-2/MY_RUN/KTEST/007_16janvier_scandollar
+/work/escobar/DEV/MNH-V5-7-1/MY_RUN/KTEST/007_16janvier_scandollar
 
 #
 #  - At ECMWF on cxa :
@@ -895,11 +926,11 @@ cp -R  007_16janvier_scandollar /.../your_directory
 #
 # use
 
-. /c1a/ms_perm/au5/MNH-V5-4-2/conf/profile_mesonh-AIX64-MNH-V5-4-2-MPIAUTO-O2
+. /c1a/ms_perm/au5/MNH-V5-7-1/conf/profile_mesonh-AIX64-MNH-V5-7-1-MPIAUTO-O2
 
 # and the examples
 
-/c1a/ms_perm/au5/MNH-V5-4-2/MY_RUN/KTEST/007_16janvier_scandollar
+/c1a/ms_perm/au5/MNH-V5-7-1/MY_RUN/KTEST/007_16janvier_scandollar
 
 
 #
@@ -991,70 +1022,58 @@ git clone -b 2014.01 https://github.com/forefireAPI/firefront.git
 # b) MNH_RTTOV for optional radiative computation
 # --------------------------------------
 #
-# The RTTOV package was not included into the open source version of Meso-NH
+# The RTTOV 13.2 package was not included into the open source version of Meso-NH
 # because it needs a licence agrement.
 #
-# ----------------------------------
-# OPTION 1: Use version 11.3 of RTTOV 
-# -----------------------------------
-# Download the RTTOV package rttov113.tar.gz by following the instructions given on http://nwpsaf.eu/site/software/rttov/
-#
-# Install the RTTOV package rttov113.tar.gz
-cd MNH.../src/LIB
-mkdir RTTOV-11.3
-cd RTTOV-11.3
-tar xvfz rttov113.tar.gz
-cd src
-make ARCH=ifort        # Use Intel "ifort" compiler; other options: gfortran, NAG, pgf90, IBM
+# Run the 'configure' script preceded with the setting of the MNH_RTTOV variable:
 #
-# And then for the compilation, run the 'configure' script preceded with the setting of the MNH_RTTOV variable:
-#
-
-cd MNH.../src/
+cd $SRC_MESONH/src/
 export MNH_RTTOV=1
-export VER_RTTOV=11.3
-./configure
-
-etc ...
-
-# ----------------------------------
-# OPTION 2: Use version 8.7 of RTTOV 
-# ----------------------------------
-# For already(old) licencied MesoNH users (MNH-4-X version with research licence see here: http://mesonh.aero.obs-mip.fr/mesonh410/UserInformation)
-#
-# the package could be reloaded in this way  
-#
-# - With cvs access
-#
-
-cd MNH.../src/LIB
-cvs up -rPACK-MNH-V4-10-3 -d -P RTTOV
-
+export VER_RTTOV=13.2
+#
+# Compile the HDF5 library
+#
+make cdf
+#
+# Download the RTTOV package rttov132.tar.xz by following the instructions given on https://nwpsaf.eu/site/software/rttov/
+#
+# Install the RTTOV package rttov132.tar.xz
+cd $SRC_MESONH/src/LIB
+mkdir RTTOV-13.2
+cd RTTOV-13.2
+tar xJf rttov132.tar.xz
+cd build
+edit Makefile.local and set HDF5_PREFIX, FFLAGS_HDF5 and LDFLAGS_HDF5 as shown below
+"
+HDF5_PREFIX  = $(SRC_MESONH)/src/dir_obj${XYZ}/MASTER/NETCDF-${VERSION_CDFF}
+FFLAGS_HDF5  = -D_RTTOV_HDF $(FFLAG_MOD)$(HDF5_PREFIX)/include
+LDFLAGS_HDF5 = -L$(HDF5_PREFIX)/lib64 -lhdf5hl_fortran -lhdf5_hl -lhdf5_fortran -lhdf5 -lsz -laec -lz -ldl
+"
+cd ../src
+../build/Makefile.PL RTTOV_HDF=1
+make ARCH=ifort        # Use Intel "ifort" compiler; other options: gfortran, NAG, pgf90, IBM
 #
-# - With WEB access (with WEB login/pass as usually) the RTTOV package could also be retrieve in tarball with wget like this:
+# And then for the Meso-NH compilation, do
 #
+cd $SRC_MESONH/src/
+make
 
-cd MNH.../src/LIB
-wget --http-user=USER  --http-password=PASS 'http://mesonh.aero.obs-mip.fr/cgi-bin/mesonh_interne/viewcvs.cgi/MNH-VX-Y-Z/src/LIB/RTTOV/?view=tar&pathrev=PACK-MNH-V4-10-3' -O RTTOV.tar.gz
-tar xvfz RTTOV.tar.gz
-
-# /!\ don't forget the ''.  Otherwise, this will not work!
-
+# c) MNH_ECRAD for optional compilation of new ECRAD radiative library from ECMWF 
+# --------------------------------------
 #
-# And then for the compilation, run the 'configure' script preceded with the setting of the MNH_RTTOV variable:
+# The default version of ECRAD is 1.4.0 (open-source)
 #
-
-cd MNH.../src/
-export MNH_RTTOV=1
-export VER_RTTOV=8.7
+# Configure & Compilation
+export MNH_ECRAD=1
 ./configure
 
 etc ...
-
-# c) MNH_ECRAD for optional compilation of new ECRAD radiative library from ECMWF 
-# --------------------------------------
+# The version of ECRAD is set by (by default):
+# export VER_ECRAD=140
+#
+# To use the previous version 1.0.1:
 #
-# The full ECRAD package was not included into the open source version of Meso-NH
+# The full ECRAD package 1.0.1  was not included into the open source version of Meso-NH
 # because it needs a licence agrement.
 #
 # See here to get the licence & full sources : https://software.ecmwf.int/wiki/display/ECRAD/ECMWF+Radiation+Scheme+Home
@@ -1068,6 +1087,7 @@ tar xvfz ecrad-1.0.1.tar.gz
 
 # Configure & Compilation
 export MNH_ECRAD=1
+export VER_ECRAD=101
 ./configure
 
 etc ...
@@ -1076,9 +1096,13 @@ etc ...
 # REM : the 'profile_mesonh...' file & the 'dir_obj...' directory will be suffixed with an ECRAD extension 
 #
 #  Usage :
-#   1) In namelist replace RAD='ECMW' be RAD='ECRA' 
-#   2) Add link to 'ecrad-1.0.1/data' files 
-#      see 007_16janvier/008_run2 test case for example
+#   1) In namelist replace RAD='ECMW' by RAD='ECRA' 
+#   2) Add link to all 'ecrad-1.X.X/data' files in your mesonh run directory
+ln -sf ${SRC_MESONH}/src/LIB/RAD/ecrad-1.X.X/data/* .
+#
+# REM : you can replace CDATADIR = "." by CDATADIR = "data" of ini_radiations_ecrad.f90 to link only the data folder instead of all the files one by one
+#
+# See 007_16janvier/008_run2 test case for example
 #
 # d) MNH_MEGAN for optional compilation of MEGAN code 
 # --------------------------------------

LIBTOOLS/tools/lfi2cdf/scripts/lfi2cdfregex.sh

-#!/bin/sh
+#!/bin/bash
 #
 #
 usage(){

LIBTOOLS/tools/lfi2cdf/src/lfi2cdf.f90

-!MNH_LIC Copyright 1994-2019 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1994-2024 CNRS, Meteo-France and Universite Paul Sabatier
 !MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
 !MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
 !MNH_LIC for details. version 1.
 !-----------------------------------------------------------------
+! Modifications:
+!  P. Wautelet 19/09/2019: add possibility to provide a fallback file if some information are not found in the input file
+!-----------------------------------------------------------------
 program LFI2CDF
   USE MODD_CONF,          ONLY: CPROGRAM
   USE MODD_CONFZ,         ONLY: NB_PROCIO_R
   USE MODD_DIM_n,         ONLY: NIMAX_ll, NJMAX_ll, NKMAX
+  USE MODD_GRID_n,        ONLY: XXHAT, XXHATM, XYHAT, XYHATM
   USE MODD_IO,            ONLY: LVERB_OUTLST, LVERB_STDOUT, NIO_ABORT_LEVEL, NIO_VERB, NGEN_ABORT_LEVEL, NGEN_VERB
   USE MODD_PARAMETERS,    ONLY: JPHEXT, JPVEXT
   USE MODD_TIMEZ,         ONLY: TIMEZ
 
-  USE MODE_FIELD,         ONLY: INI_FIELD_LIST
+  use mode_field,         only: Ini_field_list
   USE MODE_IO,            ONLY: IO_Init, IO_Config_set
+  use mode_ll
+  use mode_modeln_handler, only: Goto_model
   USE mode_options
+  USE MODE_SET_GRID,      ONLY: INTERP_HORGRID_TO_MASSPOINTS
   USE MODE_SPLITTINGZ_ll, ONLY: INI_PARAZ_ll
   USE mode_util
 
@@ -32,7 +39,7 @@ program LFI2CDF
   INTEGER :: IINFO_ll         ! return code of // routines
   INTEGER :: nfiles_out   = 0 ! number of output files
   CHARACTER(LEN=:),allocatable :: hvarlist
-  TYPE(TFILE_ELT),DIMENSION(1)        :: infiles
+  TYPE(TFILE_ELT),DIMENSION(2)        :: infiles
   TYPE(TFILE_ELT),DIMENSION(MAXFILES) :: outfiles
 
   TYPE(workfield), DIMENSION(:), POINTER :: tzreclist
@@ -44,14 +51,16 @@ program LFI2CDF
 
   CPROGRAM = 'LFICDF'
 
+  CALL Goto_model(1)
+
   CALL IO_Init()
   CALL VERSION
   CALL INI_CST
 
   ALLOCATE(TIMEZ) !Used by IO_WRITE_FIELD
 
-  NIO_VERB  = NVERB_WARNING
-  NGEN_VERB = NVERB_WARNING
+  NIO_VERB  = NVERB_INFO
+  NGEN_VERB = NVERB_INFO
   NIO_ABORT_LEVEL = NVERB_FATAL
   NGEN_ABORT_LEVEL = NVERB_FATAL
   LVERB_OUTLST = .FALSE.
@@ -87,12 +96,12 @@ program LFI2CDF
      CALL IO_Config_set()
   END IF
 
-  CALL INI_FIELD_LIST(1)
+  CALL INI_FIELD_LIST()
 
   CALL OPEN_FILES(infiles, outfiles, nfiles_out, hinfile, houtfile, nbvar_infile, options, runmode)
   IF (options(OPTLIST)%set) STOP
 
-  !Set and initialize parallel variables (necessary to read splitted files)
+  !Set and initialize parallel variables (necessary to read split files)
   CALL SET_JP_ll(1,JPHEXT,JPVEXT,JPHEXT)
   CALL SET_DAD0_ll()
   CALL SET_DIM_ll(NIMAX_ll, NJMAX_ll, NKMAX)
@@ -104,6 +113,12 @@ program LFI2CDF
   CALL SET_YEND_ll(NJMAX_ll+2*JPHEXT, 1)
   CALL INI_PARAZ_ll(IINFO_ll)
 
+  ! This has to be done after INI_PARAZ_ll and after reading of XXHAT and XYHAT (in OPEN_FILES)
+  ALLOCATE(XXHATM(NIMAX_ll+2*JPHEXT))
+  ALLOCATE(XYHATM(NJMAX_ll+2*JPHEXT))
+  ! Interpolations of positions to mass points
+  CALL INTERP_HORGRID_TO_MASSPOINTS( XXHAT, XYHAT, XXHATM, XYHATM )
+
   IF (runmode == MODELFI2CDF .OR. runmode == MODECDF2CDF) THEN
      IF (options(OPTVAR)%set) THEN
         ! nbvar_tbr is computed from number of requested variables
@@ -151,7 +166,11 @@ program LFI2CDF
      CALL fill_files(infiles,outfiles,tzreclist,nbvar,options)
   END IF
 
-  CALL CLOSE_FILES(infiles, 1)
+  if ( options( OPTFALLBACK )%set ) then
+    CALL CLOSE_FILES(infiles, 2)
+  else
+    CALL CLOSE_FILES(infiles, 1)
+  end if
   CALL CLOSE_FILES(outfiles,nfiles_out)
-  
+
 end program LFI2CDF

LIBTOOLS/tools/lfi2cdf/src/mode_options.f90

-!MNH_LIC Copyright 2015-2018 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2015-2020 CNRS, Meteo-France and Universite Paul Sabatier
 !MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
-!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
+!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
 !MNH_LIC for details. version 1.
 !-----------------------------------------------------------------
+! Modifications:
+!  P. Wautelet 19/09/2019: add possibility to provide a fallback file if some information are not found in the input file
+!  P. Wautelet 21/10/2019: add OPTDIR option to set directory for writing outfiles
+!-----------------------------------------------------------------
 module mode_options
-  USE MODE_FIELD, ONLY: TYPEUNDEF, TYPEINT, TYPELOG, TYPEREAL, TYPECHAR, TYPEDATE
+  use modd_field, only: TYPEUNDEF, TYPEINT, TYPELOG, TYPEREAL, TYPECHAR, TYPEDATE
 
   implicit none
 
-  integer,parameter :: nbavailoptions = 10
+  integer,parameter :: NBAVAILOPTIONS = 12
   integer,parameter :: MODEUNDEF = -11, MODECDF2CDF = 11, MODELFI2CDF = 12, MODECDF2LFI = 13
 
-  integer,parameter :: OPTCOMPRESS = 1, OPTHELP   = 2, OPTLIST   = 3
-  integer,parameter :: OPTMERGE    = 4, OPTOUTPUT = 5, OPTREDUCE = 6
-  integer,parameter :: OPTMODE     = 7, OPTSPLIT  = 8, OPTVAR    = 9
-  integer,parameter :: OPTVERBOSE  = 10
+  integer,parameter :: OPTCOMPRESS = 1,  OPTHELP     = 2,  OPTLIST   = 3
+  integer,parameter :: OPTMERGE    = 4,  OPTOUTPUT   = 5,  OPTREDUCE = 6
+  integer,parameter :: OPTMODE     = 7,  OPTSPLIT    = 8,  OPTVAR    = 9
+  integer,parameter :: OPTVERBOSE  = 10, OPTFALLBACK = 11, OPTDIR    = 12
 
   type option
     logical :: set = .false.
@@ -153,6 +157,17 @@ subroutine init_options(options)
   options(OPTVERBOSE)%long_name    = "verbose"
   options(OPTVERBOSE)%short_name   = 'V'
   options(OPTVERBOSE)%has_argument = .false.
+
+  options(OPTFALLBACK)%long_name    = "fallback-file"
+  options(OPTFALLBACK)%short_name   = 'f'
+  options(OPTFALLBACK)%has_argument = .true.
+  options(OPTFALLBACK)%type         = TYPECHAR
+
+  options(OPTDIR)%long_name    = "outdir"
+  options(OPTDIR)%short_name   = 'd'
+  options(OPTDIR)%has_argument = .true.
+  options(OPTDIR)%type         = TYPECHAR
+
 end subroutine init_options
 
 subroutine get_option(options,finished)
@@ -278,12 +293,6 @@ subroutine check_options(options,infile,runmode)
     call help()
   end if
 
-  !Merge flag only supported if -v is set
-  if (options(OPTMERGE)%set .AND. .NOT.options(OPTVAR)%set) then
-    print *,'Error: merge option must be used with var option'
-    call help()
-  end if
-
   !Split flag only supported if -v is set
   if (options(OPTSPLIT)%set .AND. .NOT.options(OPTVAR)%set) then
       options(OPTSPLIT)%set = .false.
@@ -329,18 +338,21 @@ subroutine help()
 !TODO: -l option for cdf2cdf and cdf2lfi
   print *,"Usage : lfi2cdf [-h --help] [-l] [-v --var var1[,...]] [-r --reduce-precision]"
   print *,"                [-m --merge number_of_z_levels] [-s --split] [-o --output output-file.nc]"
-  print *,"                [-R --runmode mode] [-V --verbose]"
+  print *,"                [-R --runmode mode] [-V --verbose] [-f --fallback-file fallback-file]"
   print *,"                [-c --compress compression_level] input-file.lfi"
   print *,"        cdf2cdf [-h --help] [-v --var var1[,...]] [-r --reduce-precision]"
   print *,"                [-m --merge number_of_split_files] [-s --split] [-o --output output-file.nc]"
-  print *,"                [-R --runmode mode] [-V --verbose]"
+  print *,"                [-R --runmode mode] [-V --verbose] [-f --fallback-file fallback-file]"
   print *,"                [-c --compress compression_level] input-file.nc"
-  print *,"        cdf2lfi [-o --output output-file.lfi] [-R --runmode mode]  [-V --verbose] input-file.nc"
+  print *,"        cdf2lfi [-o --output output-file.lfi] [-R --runmode mode]  [-V --verbose]"
+  print *,"                [-f --fallback-file fallback-file] input-file.nc"
   print *,""
   print *,"Options:"
   print *,"  --compress, -c compression_level"
   print *,"     Compress data. The compression level should be in the 1 to 9 interval."
   print *,"     Only supported with the netCDF format (cdf2cdf and lfi2cdf only)"
+  print *,"  -f --fallback-file fallback-file"
+  print *,"     File to use to read some grid information if not found in input-file"
   print *,"  --help, -h"
   print *,"     Print this text"
   print *,"  --list, -l"

LIBTOOLS/tools/lfi2cdf/src/mode_util.f90

-!MNH_LIC Copyright 1994-2019 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1994-2024 CNRS, Meteo-France and Universite Paul Sabatier
 !MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
 !MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
 !MNH_LIC for details. version 1.
@@ -6,20 +6,29 @@
 ! Modifications:
 !  P. Wautelet 07/02/2019: force TYPE to a known value for IO_FILE_ADD2LIST
 !  P. Wautelet 10/04/2019: use IO_Err_handle_nc4 to handle netCDF errors
+!  P. Wautelet 25/06/2019: add support for 3D integer arrays
+!  P. Wautelet 01/08/2019: allow merge of entire Z-split files
+!  P. Wautelet 18/09/2019: correct support of 64bit integers (MNH_INT=8)
+!  P. Wautelet 19/09/2019: add possibility to provide a fallback file if some information are not found in the input file
+!  P. Wautelet 21/10/2019: add OPTDIR option to set directory for writing outfiles
+!  P. Wautelet 21/10/2019: if DTMOD and DTCUR not found, try to read the time coordinate
+!  P. Wautelet 10/11/2020: new data structures for netCDF dimensions
 !-----------------------------------------------------------------
 MODULE mode_util
-  USE MODD_IO,         ONLY: TFILE_ELT
-  USE MODD_NETCDF,     ONLY: DIMCDF, CDFINT
+  use modd_field,      only: tfieldmetadata, tfieldlist
+  USE MODD_IO,         ONLY: TFILEDATA, TFILE_ELT
+  USE MODD_NETCDF,     ONLY: CDFINT, tdimnc
   USE MODD_PARAMETERS, ONLY: NLFIMAXCOMMENTLENGTH, NMNHNAMELGTMAX
   use modd_precision,  only: LFIINT
 
-  USE MODE_FIELD
+  use mode_field,          only: Find_field_id_from_mnhname
   USE MODE_IO_FIELD_READ
   USE MODE_IO_FIELD_WRITE
   use mode_io_tools_nc4,   only: IO_Err_handle_nc4
-
+  use mode_msg
   USE mode_options
 
+
   USE NETCDF
 
   IMPLICIT NONE 
@@ -39,16 +48,16 @@ MODULE mode_util
      LOGICAL                               :: LSPLIT = .FALSE. ! TRUE if variable is split by vertical level
      INTEGER                               :: NSIZE = 0 ! Size of the variable (in number of elements)
      INTEGER                               :: NSRC = 0 ! Number of variables used to compute the variable (needed only if calc=.true.)
-     INTEGER                               :: NDIMS_FILE  ! Number of dims (as present in input file)
-     INTEGER,DIMENSION(:),ALLOCATABLE      :: NDIMSIZES_FILE  ! Dimensions sizes (as present in input file)
+     INTEGER(kind=CDFINT)                  :: NDIMS_FILE  ! Number of dims (as present in input file)
+     INTEGER(kind=CDFINT),    DIMENSION(:),ALLOCATABLE :: NDIMSIZES_FILE  ! Dimensions sizes (as present in input file)
      CHARACTER(LEN=NF90_MAX_NAME),DIMENSION(:),ALLOCATABLE :: CDIMNAMES_FILE  ! Dimensions names (as present in input file)
      CHARACTER(LEN=40)                     :: CUNITS_FILE = '' ! Units (as present in input file)
      INTEGER                               :: NGRID_FILE  ! Grid  number (as present in input file)
-     INTEGER                               :: NTYPE_FILE  ! netCDF datatype (NF90_CHAR, NF90_INT...) (as present in input file)
+     INTEGER(kind=CDFINT)                  :: NTYPE_FILE  ! netCDF datatype (NF90_CHAR, NF90_INT...) (as present in input file)
      INTEGER,DIMENSION(MAXRAW)             :: src    ! List of variables used to compute the variable (needed only if calc=.true.)
      INTEGER                               :: tgt    ! Target: id of the variable that use it (calc variable)
-     TYPE(TFIELDDATA)                      :: TFIELD ! Metadata about the field
-     TYPE(DIMCDF),DIMENSION(:),ALLOCATABLE :: TDIMS  ! Dimensions of the field
+     TYPE(TFIELDMETADATA)                  :: TFIELD ! Metadata about the field
+     TYPE(tdimnc),DIMENSION(:),ALLOCATABLE :: TDIMS  ! Dimensions of the field
   END TYPE workfield
 
   LOGICAL(KIND=LFIINT), PARAMETER :: ltrue  = .TRUE.
@@ -60,6 +69,7 @@ CONTAINS
   SUBROUTINE parse_infiles(infiles, outfiles, KNFILES_OUT, nbvar_infile, nbvar_tbr, nbvar_calc, nbvar_tbw, &
                            tpreclist, options, runmode)
     USE MODD_DIM_n,      ONLY: NIMAX_ll, NJMAX_ll, NKMAX
+    use modd_io,         only: nio_verb
     USE MODD_PARAMETERS, ONLY: JPHEXT, JPVEXT, NGRIDUNKNOWN
 
     use mode_io_tools_nc4, only: IO_Dimids_guess_nc4
@@ -94,6 +104,7 @@ CONTAINS
     LOGICAL                                  :: ladvan
     LOGICAL                                  :: GOK
     TYPE(TLFIDATE),DIMENSION(MAXDATES)       :: TLFIDATES
+    type(TFILEDATA)                          :: tzfile
 
     CALL PRINT_MSG(NVERB_DEBUG,'IO','parse_infiles','called')
 
@@ -116,10 +127,11 @@ CONTAINS
       kcdf_id = INFILES(1)%TFILE%NNCID
     END IF
 
-    PRINT *,'MESONH 3D, 2D articles DIMENSIONS used :'
-    PRINT *,'DIMX =',NIMAX_ll+2*JPHEXT
-    PRINT *,'DIMY =',NJMAX_ll+2*JPHEXT
-    PRINT *,'DIMZ =',NKMAX   +2*JPVEXT
+    WRITE( cmnhmsg(1), '( A )' ) 'MESONH 3D, 2D articles DIMENSIONS used :'
+    WRITE( cmnhmsg(2), '( "DIMX = ", I0 )' ) NIMAX_ll + 2 * JPHEXT
+    WRITE( cmnhmsg(3), '( "DIMY = ", I0 )' ) NJMAX_ll + 2 * JPHEXT
+    WRITE( cmnhmsg(4), '( "DIMZ = ", I0 )' ) NKMAX    + 2 * JPVEXT
+    call Print_msg( NVERB_INFO, 'IO', 'parse_infiles' )
 
     ! Phase 1 : build articles list to convert.
     !
@@ -228,6 +240,7 @@ CONTAINS
               !If we are merging, maybe it is one of the split variable
               !In that case, the 1st part of the variable is in the 1st split file with a 0001 suffix
               kcdf_id2 = INFILES(1)%TFILE%TFILES_IOZ(1)%TFILE%NNCID
+              tzfile = INFILES(1)%TFILE%TFILES_IOZ(1)%TFILE
               status = NF90_INQ_VARID(kcdf_id2,trim(yrecfm)//'0001',var_id)
               IF (status == NF90_NOERR) THEN
                 tpreclist(ji)%LSPLIT = .true.
@@ -242,11 +255,12 @@ CONTAINS
               call IO_Err_handle_nc4( status, 'parse_infiles', 'NF90_INQ_VARID', trim(yrecfm) )
             ELSE
               kcdf_id2 = kcdf_id
+              tzfile = INFILES(1)%TFILE
             ENDIF
             !
             IF (status == NF90_NOERR) THEN
               tpreclist(ji)%found = .true.
-              CALL IO_Metadata_get_nc4(kcdf_id2,var_id,tpreclist(ji))
+              CALL IO_Metadata_get_nc4(tzfile,var_id,tpreclist(ji))
             END IF
           END IF
 
@@ -258,10 +272,14 @@ CONTAINS
        END DO
 
        maxvar = nbvar_tbr+nbvar_calc
-DO ji=1,nbvar_tbr+nbvar_calc
-  print *,ji,'name=',trim(tpreclist(ji)%name),' calc=',tpreclist(ji)%calc,' tbw=',tpreclist(ji)%tbw,&
-          ' tbr=',tpreclist(ji)%tbr,' found=',tpreclist(ji)%found
-END DO
+
+       if ( nio_verb >= NVERB_INFO ) then
+        do ji = 1, nbvar_tbr + nbvar_calc
+          write ( cmnhmsg(1), '( I0, " name=", A, "calc=", L1, " tbw=", L1, " tbr=", L1, " found=", L1 )' )  &
+                ji, tpreclist(ji)%name, tpreclist(ji)%calc, tpreclist(ji)%tbw, tpreclist(ji)%tbr, tpreclist(ji)%found
+          call Print_msg( NVERB_INFO, 'IO', 'parse_infiles' )
+        end do
+      end if
 
     ELSE
        ! Entire file is converted
@@ -339,7 +357,7 @@ END DO
            if ( status /= NF90_NOERR ) &
              call IO_Err_handle_nc4( status, 'parse_infiles', 'NF90_INQUIRE_VARIABLE', tpreclist(ji)%name )
            tpreclist(ji)%found  = .TRUE.
-           CALL IO_Metadata_get_nc4(kcdf_id,var_id,tpreclist(ji))
+           CALL IO_Metadata_get_nc4(INFILES(1)%TFILE,var_id,tpreclist(ji))
          END DO
        END IF
 
@@ -376,7 +394,7 @@ END DO
       ELSE
         CALL FIND_FIELD_ID_FROM_MNHNAME(tpreclist(ji)%name,IID,IRESP,ONOWARNING=.TRUE.)
         IF (IRESP==0) THEN
-          tpreclist(ji)%TFIELD = TFIELDLIST(IID)
+          tpreclist(ji)%TFIELD = TFIELDMETADATA( TFIELDLIST(IID) )
           ! Determine TDIMS
           IF (runmode==MODELFI2CDF) THEN
             ALLOCATE(tpreclist(ji)%TDIMS(tpreclist(ji)%TFIELD%NDIMS))
@@ -399,11 +417,11 @@ END DO
         ELSE !Field not found in list, try to determine characteristics
           tpreclist(ji)%TFIELD%CMNHNAME = TRIM(tpreclist(ji)%name)
           tpreclist(ji)%TFIELD%CSTDNAME = ''
-          tpreclist(ji)%TFIELD%CLONGNAME = TRIM(tpreclist(ji)%name)
-          tpreclist(ji)%TFIELD%CUNITS = ''
+          !Set in IO_Metadata_get_nc4 (and not used for LFI) tpreclist(ji)%TFIELD%CLONGNAME = TRIM(tpreclist(ji)%name)
+          !Set in IO_Metadata_get_nc4 (and not used for LFI) tpreclist(ji)%TFIELD%CUNITS = ''
           tpreclist(ji)%TFIELD%CDIR = 'XY' !Assumption...
           tpreclist(ji)%TFIELD%CLBTYPE = 'NONE'
-          tpreclist(ji)%TFIELD%CCOMMENT = ''
+          !Set in IO_Metadata_get_nc4 (and not used for LFI) tpreclist(ji)%TFIELD%CCOMMENT = ''
           !
           IF (runmode==MODELFI2CDF) THEN
             tpreclist(ji)%TFIELD%NGRID = 1 !Assumption
@@ -417,7 +435,7 @@ END DO
               CALL IO_Dimids_guess_nc4(outfiles(idx_out)%TFILE,tpreclist(ji)%TFIELD,&
                                       tpreclist(ji)%NSIZE,tpreclist(ji)%TDIMS,IRESP)
               !
-              IF (IRESP/=0 .OR. tpreclist(ji)%TDIMS(3)%LEN==1) THEN
+              IF (IRESP/=0 .OR. tpreclist(ji)%TDIMS(3)%nlen==1) THEN
                 CALL PRINT_MSG(NVERB_DEBUG,'IO','parse_infiles',tpreclist(ji)%TFIELD%CMNHNAME//': try 2D')
                 !Try again with 2D
                 tpreclist(ji)%TFIELD%NDIMS = 2
@@ -425,7 +443,7 @@ END DO
                                         tpreclist(ji)%NSIZE,tpreclist(ji)%TDIMS,IRESP)
               END IF
               !
-              IF (IRESP/=0 .OR. tpreclist(ji)%TDIMS(2)%LEN==1) THEN
+              IF (IRESP/=0 .OR. tpreclist(ji)%TDIMS(2)%nlen==1) THEN
                 CALL PRINT_MSG(NVERB_DEBUG,'IO','parse_infiles',tpreclist(ji)%TFIELD%CMNHNAME//': try 1D')
                 !Try again with 1D
                 tpreclist(ji)%TFIELD%NDIMS = 1
@@ -558,7 +576,6 @@ END DO
         tpreclist(ji)%TFIELD%NTYPE    = tpreclist(idx_var)%TFIELD%NTYPE
         tpreclist(ji)%TFIELD%NDIMS    = tpreclist(idx_var)%TFIELD%NDIMS
 #if 0
-!PW: TODO?
         tpreclist(ji)%TFIELD%NFILLVALUE
         tpreclist(ji)%TFIELD%XFILLVALUE
         tpreclist(ji)%TFIELD%NVALIDMIN
@@ -629,7 +646,8 @@ END DO
 
     CHARACTER(LEN=16) :: YMNHVERSION
     CHARACTER(LEN=:),ALLOCATABLE :: YHISTORY
-    INTEGER :: ilen, ji
+    INTEGER              :: ji
+    INTEGER(KIND=CDFINT) :: ilen
     INTEGER(KIND=CDFINT) :: status
     INTEGER(KIND=CDFINT) :: kcdf_id
 
@@ -687,15 +705,17 @@ END DO
     INTEGER(KIND=CDFINT),DIMENSION(NF90_MAX_VAR_DIMS) :: IDIMLEN
     logical,dimension(knaf)                  :: gtimedep_in, gtimedep_out
 
-    CHARACTER(LEN=:),       ALLOCATABLE :: YTAB0D
-    INTEGER,DIMENSION(:),   ALLOCATABLE :: ITAB1D, ITAB1D2
-    INTEGER,DIMENSION(:,:), ALLOCATABLE :: ITAB2D, ITAB2D2
-    LOGICAL,DIMENSION(:),   ALLOCATABLE :: GTAB1D
-    REAL,DIMENSION(:),      ALLOCATABLE :: XTAB1D, XTAB1D2
-    REAL,DIMENSION(:,:),    ALLOCATABLE :: XTAB2D, XTAB2D2
-    REAL,DIMENSION(:,:,:),  ALLOCATABLE :: XTAB3D, XTAB3D2
-    REAL,DIMENSION(:,:,:,:),ALLOCATABLE :: XTAB4D, XTAB4D2
+    CHARACTER(LEN=:),            ALLOCATABLE :: YTAB0D
+    INTEGER, DIMENSION(:),       ALLOCATABLE :: ITAB1D, ITAB1D2
+    INTEGER, DIMENSION(:,:),     ALLOCATABLE :: ITAB2D, ITAB2D2
+    INTEGER, DIMENSION(:,:,:),   ALLOCATABLE :: ITAB3D, ITAB3D2
+    LOGICAL, DIMENSION(:),       ALLOCATABLE :: GTAB1D
+    REAL,    DIMENSION(:),       ALLOCATABLE :: XTAB1D, XTAB1D2
+    REAL,    DIMENSION(:,:),     ALLOCATABLE :: XTAB2D, XTAB2D2
+    REAL,    DIMENSION(:,:,:),   ALLOCATABLE :: XTAB3D, XTAB3D2
+    REAL,    DIMENSION(:,:,:,:), ALLOCATABLE :: XTAB4D, XTAB4D2
     TYPE(DATE_TIME)                     :: TZDATE
+    TYPE(TFILEDATA)                     :: TZFILE
 
 
     CALL PRINT_MSG(NVERB_DEBUG,'IO','fill_files','called')
@@ -718,7 +738,7 @@ END DO
 
       SELECT CASE(tpreclist(ji)%TFIELD%NTYPE)
       CASE (TYPEINT)
-        IDIMLEN(1:IDIMS) = tpreclist(ji)%TDIMS(1:IDIMS)%LEN
+        IDIMLEN(1:IDIMS) = tpreclist(ji)%TDIMS(1:IDIMS)%nlen
 
         IF (.NOT.tpreclist(ji)%calc) THEN
           INSRC = 1
@@ -742,6 +762,10 @@ END DO
           ALLOCATE(ITAB2D(IDIMLEN(1),IDIMLEN(2)))
           IF (tpreclist(ji)%calc) ALLOCATE(ITAB2D2(IDIMLEN(1),IDIMLEN(2)))
           CALL IO_Field_read(INFILES(1)%TFILE,tpreclist(ISRC)%TFIELD,ITAB2D)
+        CASE (3)
+          ALLOCATE(ITAB3D(IDIMLEN(1),IDIMLEN(2),IDIMLEN(3)))
+          IF (tpreclist(ji)%calc) ALLOCATE(ITAB3D2(IDIMLEN(1),IDIMLEN(2),IDIMLEN(3)))
+          CALL IO_Field_read(INFILES(1)%TFILE,tpreclist(ISRC)%TFIELD,ITAB3D)
         CASE DEFAULT
           CALL PRINT_MSG(NVERB_WARNING,'IO','fill_files','too many dimensions for ' &
                          //TRIM(tpreclist(ISRC)%name)//' => ignored')
@@ -762,6 +786,9 @@ END DO
           CASE (2)
             CALL IO_Field_read(INFILES(1)%TFILE,tpreclist(ISRC)%TFIELD,ITAB2D2)
             ITAB2D(:,:) = ITAB2D(:,:) + ITAB2D2(:,:)
+          CASE (3)
+            CALL IO_Field_read(INFILES(1)%TFILE,tpreclist(ISRC)%TFIELD,ITAB3D2)
+            ITAB3D(:,:,:) = ITAB3D(:,:,:) + ITAB3D2(:,:,:)
           END SELECT
         END DO
 
@@ -779,11 +806,15 @@ END DO
           CALL IO_Field_write(outfiles(idx)%TFILE,tpreclist(ji)%TFIELD,ITAB2D)
           DEALLOCATE(ITAB2D)
           IF (tpreclist(ji)%calc) DEALLOCATE(ITAB2D2)
+        CASE (3)
+          CALL IO_Field_write(outfiles(idx)%TFILE,tpreclist(ji)%TFIELD,ITAB3D)
+          DEALLOCATE(ITAB3D)
+          IF (tpreclist(ji)%calc) DEALLOCATE(ITAB3D2)
         END SELECT
 
 
       CASE (TYPELOG)
-        IDIMLEN(1:IDIMS) = tpreclist(ji)%TDIMS(1:IDIMS)%LEN
+        IDIMLEN(1:IDIMS) = tpreclist(ji)%TDIMS(1:IDIMS)%nlen
 
         tpreclist(ji)%TFIELD%LTIMEDEP = gtimedep_in(ji)
         SELECT CASE(IDIMS)
@@ -807,7 +838,7 @@ END DO
 
 
       CASE (TYPEREAL)
-        IDIMLEN(1:IDIMS) = tpreclist(ji)%TDIMS(1:IDIMS)%LEN
+        IDIMLEN(1:IDIMS) = tpreclist(ji)%TDIMS(1:IDIMS)%nlen
 
         IF (.NOT.tpreclist(ji)%calc) THEN
           INSRC = 1
@@ -834,7 +865,16 @@ END DO
         CASE (3)
           ALLOCATE(XTAB3D(IDIMLEN(1),IDIMLEN(2),IDIMLEN(3)))
           IF (tpreclist(ji)%calc) ALLOCATE(XTAB3D2(IDIMLEN(1),IDIMLEN(2),IDIMLEN(3)))
-          CALL IO_Field_read(INFILES(1)%TFILE,tpreclist(ISRC)%TFIELD,XTAB3D)
+          !Hack not very clean: 3D LB fields are not split
+          !If NSUBFILES_IOZ is set to 0, IO_Field_read will read it as a non-split field
+          !CAUTION: there are no guarantee the IO_Field_read will continue to use this information that way...
+          if ( tpreclist(ji)%tfield%clbtype /= 'NONE' .or. tpreclist(ji)%name(1:2) == 'LB' ) then
+            tzfile = infiles(1)%tfile
+            tzfile%nsubfiles_ioz=0
+            call IO_Field_read(tzfile,tpreclist(isrc)%tfield,xtab3d)
+          else
+            call IO_Field_read(infiles(1)%tfile,tpreclist(isrc)%tfield,xtab3d)
+          end if
         CASE (4)
           ALLOCATE(XTAB4D(IDIMLEN(1),IDIMLEN(2),IDIMLEN(3),IDIMLEN(4)))
           IF (tpreclist(ji)%calc) ALLOCATE(XTAB4D2(IDIMLEN(1),IDIMLEN(2),IDIMLEN(3),IDIMLEN(4)))
@@ -940,16 +980,19 @@ END DO
   SUBROUTINE OPEN_FILES(infiles,outfiles,KNFILES_OUT,hinfile,houtfile,nbvar_infile,options,runmode)
     USE MODD_CONF,          ONLY: LCARTESIAN
     USE MODD_CONF_n,        ONLY: CSTORAGE_TYPE
+    USE MODD_CONFZ,         ONLY: NB_PROCIO_R
     USE MODD_DIM_n,         ONLY: NIMAX_ll, NJMAX_ll, NKMAX
     USE MODD_GRID,          ONLY: XBETA, XRPK, XLAT0, XLON0, XLATORI, XLONORI
-    USE MODD_GRID_n,        ONLY: LSLEVE, XXHAT, XYHAT, XZHAT
+    USE MODD_GRID_n,        ONLY: LSLEVE, XXHAT, XXHATM, XYHAT, XYHATM, XZHAT, XZHATM, &
+                                  XXHAT_ll, XYHAT_ll, XXHATM_ll, XYHATM_ll
     USE MODD_IO,            ONLY: LIOCDF4
-    USE MODD_PARAMETERS,    ONLY: JPHEXT
+    USE MODD_PARAMETERS,    ONLY: JPHEXT, JPVEXT
     USE MODD_PARAMETERS_ll, ONLY: JPHEXT_ll=>JPHEXT, JPVEXT_ll=>JPVEXT
     USE MODD_TIME_n,        ONLY: TDTCUR, TDTMOD
 
     USE MODE_IO_FILE,          ONLY: IO_FILE_OPEN, IO_FILE_CLOSE
     USE MODE_IO_MANAGE_STRUCT, ONLY: IO_FILE_ADD2LIST
+    USE MODE_SET_GRID,         ONLY: INTERP_VERGRID_TO_MASSPOINTS
 
     TYPE(TFILE_ELT),DIMENSION(:),INTENT(OUT) :: infiles
     TYPE(TFILE_ELT),DIMENSION(:),INTENT(OUT) :: outfiles
@@ -960,10 +1003,16 @@ END DO
     TYPE(option),DIMENSION(:),   INTENT(IN)  :: options
     INTEGER,                     INTENT(IN)  :: runmode
 
-    INTEGER              :: idx, IRESP2
-    INTEGER(KIND=CDFINT) :: omode
-    INTEGER(KIND=CDFINT) :: status
-    INTEGER(KIND=LFIINT) :: ilu,iresp
+    character(len=:), allocatable :: yunits
+    INTEGER                       :: idx, IRESP2
+    integer                       :: iiu, iju, iku
+    integer                       :: inb_procio_r_save
+    INTEGER(KIND=CDFINT)          :: ioldmode
+    INTEGER(KIND=CDFINT)          :: istatus
+    INTEGER(KIND=CDFINT)          :: ivar_id
+    integer(kind=CDFINT)          :: ilen
+    INTEGER(KIND=LFIINT)          :: ilu,iresp
+    logical                       :: gok
 
 
     CALL PRINT_MSG(NVERB_DEBUG,'IO','OPEN_FILES','called')
@@ -979,7 +1028,20 @@ END DO
        CALL IO_FILE_ADD2LIST(INFILES(1)%TFILE,HINFILE,'MNH','READ',HFORMAT='NETCDF4')
        CALL IO_FILE_OPEN(INFILES(1)%TFILE)
 
-       nbvar_infile = INFILES(1)%TFILE%NNCNAR
+       istatus = NF90_INQUIRE( infiles(1)%tfile%nncid, nvariables=nbvar_infile )
+       if (istatus /= NF90_NOERR) then
+         call Print_msg( NVERB_FATAL, 'IO', 'OPEN_FILES', 'NF90_INQUIRE for ' &
+                         // trim(infiles(1)%tfile%cname) // '.nc: ' // NF90_STRERROR(istatus) )
+       end if
+
+       !Open fallback file if provided
+       if ( options( OPTFALLBACK )%set ) then
+         inb_procio_r_save = NB_PROCIO_R
+         NB_PROCIO_R = 1
+         CALL IO_FILE_ADD2LIST(INFILES(2)%TFILE,options( OPTFALLBACK )%cvalue,'MNH','READ',HFORMAT='NETCDF4')
+         CALL IO_FILE_OPEN(INFILES(2)%TFILE)
+         NB_PROCIO_R = inb_procio_r_save
+       end if
    ELSE
        !
        ! LFI
@@ -997,52 +1059,177 @@ END DO
           CALL IO_FILE_CLOSE(INFILES(1)%TFILE)
           return
        END IF
+
+       !Open fallback file if provided
+       if ( options( OPTFALLBACK )%set ) then
+         inb_procio_r_save = NB_PROCIO_R
+         NB_PROCIO_R = 1
+         CALL IO_FILE_ADD2LIST(INFILES(2)%TFILE,options( OPTFALLBACK )%cvalue,'UNKNOWN','READ', &
+                               HFORMAT='LFI',KLFIVERB=0)
+         CALL IO_FILE_OPEN(INFILES(2)%TFILE)
+         NB_PROCIO_R = inb_procio_r_save
+       end if
    END IF
    !
    !Read problem dimensions and some grid variables (needed to determine domain size and also by IO_FILE_OPEN to create netCDF files)
-   CALL IO_Field_read(INFILES(1)%TFILE,'JPHEXT',JPHEXT)
+   JPHEXT = 1
+   CALL IO_Field_read(INFILES(1)%TFILE,'JPHEXT',JPHEXT,IRESP2)
+   !If not found in main file, try the fallback one
+   if ( options( OPTFALLBACK )%set .and. iresp2 /= 0 ) CALL IO_Field_read(INFILES(2)%TFILE,'JPHEXT',JPHEXT,IRESP2)
+   if ( iresp2 /= 0 ) call Print_msg( NVERB_ERROR, 'IO', 'OPEN_FILES', 'JPHEXT not found')
+
    JPHEXT_ll = JPHEXT
    JPVEXT_ll = JPVEXT
    !
    ALLOCATE(NIMAX_ll,NJMAX_ll,NKMAX)
-   CALL IO_Field_read(INFILES(1)%TFILE,'IMAX',NIMAX_ll)
-   CALL IO_Field_read(INFILES(1)%TFILE,'JMAX',NJMAX_ll)
+   CALL IO_Field_read(INFILES(1)%TFILE,'IMAX',NIMAX_ll,IRESP2)
+   !If not found in main file, try the fallback one
+   if ( options( OPTFALLBACK )%set .and. iresp2 /= 0 ) CALL IO_Field_read(INFILES(2)%TFILE,'IMAX',NIMAX_ll,IRESP2)
+   if ( iresp2 /= 0 ) call Print_msg( NVERB_ERROR, 'IO', 'OPEN_FILES', 'IMAX not found')
+
+   CALL IO_Field_read(INFILES(1)%TFILE,'JMAX',NJMAX_ll,IRESP2)
+   !If not found in main file, try the fallback one
+   if ( options( OPTFALLBACK )%set .and. iresp2 /= 0 ) CALL IO_Field_read(INFILES(2)%TFILE,'JMAX',NJMAX_ll,IRESP2)
+   if ( iresp2 /= 0 ) call Print_msg( NVERB_ERROR, 'IO', 'OPEN_FILES', 'JMAX not found')
+
    CALL IO_Field_read(INFILES(1)%TFILE,'KMAX',NKMAX,IRESP2)
+   !If not found in main file, try the fallback one
+   if ( options( OPTFALLBACK )%set .and. iresp2 /= 0 ) CALL IO_Field_read(INFILES(2)%TFILE,'KMAX',NKMAX,IRESP2)
    IF (IRESP2/=0) NKMAX = 0
    !
-   CALL IO_Field_read(INFILES(1)%TFILE,'PROGRAM',CPROGRAM_ORIG)
+   CALL IO_Field_read(INFILES(1)%TFILE,'PROGRAM',CPROGRAM_ORIG,IRESP2)
+   !If not found in main file, try the fallback one
+   if ( options( OPTFALLBACK )%set .and. iresp2 /= 0 ) CALL IO_Field_read(INFILES(2)%TFILE,'PROGRAM',CPROGRAM_ORIG,IRESP2)
+   if ( iresp2 /= 0 ) call Print_msg( NVERB_ERROR, 'IO', 'OPEN_FILES', 'PROGRAM not found')
    !
    ALLOCATE(CSTORAGE_TYPE)
-   CALL IO_Field_read(INFILES(1)%TFILE,'STORAGE_TYPE',CSTORAGE_TYPE)
-   !
-   ALLOCATE(XXHAT(NIMAX_ll+2*JPHEXT))
-   CALL IO_Field_read(INFILES(1)%TFILE,'XHAT',XXHAT)
-   ALLOCATE(XYHAT(NJMAX_ll+2*JPHEXT))
-   CALL IO_Field_read(INFILES(1)%TFILE,'YHAT',XYHAT)
-   CALL IO_Field_read(INFILES(1)%TFILE,'CARTESIAN',LCARTESIAN)
-   !
-   CALL IO_Field_read(INFILES(1)%TFILE,'LAT0',XLAT0)
-   CALL IO_Field_read(INFILES(1)%TFILE,'LON0',XLON0)
-   CALL IO_Field_read(INFILES(1)%TFILE,'BETA',XBETA)
+   CALL IO_Field_read(INFILES(1)%TFILE,'STORAGE_TYPE',CSTORAGE_TYPE,IRESP2)
+   !If not found in main file, try the fallback one
+   if ( options( OPTFALLBACK )%set .and. iresp2 /= 0 ) CALL IO_Field_read(INFILES(2)%TFILE,'STORAGE_TYPE',CSTORAGE_TYPE,IRESP2)
+   if ( iresp2 /= 0 ) call Print_msg( NVERB_ERROR, 'IO', 'OPEN_FILES', 'STORAGE_TYPE not found')
    !
+   iiu = nimax_ll + 2 * JPHEXT
+   ALLOCATE( XXHAT(iiu), XXHATM(iiu) )
+   CALL IO_Field_read(INFILES(1)%TFILE,'XHAT',XXHAT,IRESP2)
+   !If not found in main file, try the fallback one
+   if ( options( OPTFALLBACK )%set .and. iresp2 /= 0 ) CALL IO_Field_read(INFILES(2)%TFILE,'XHAT',XXHAT,IRESP2)
+   if ( iresp2 /= 0 ) call Print_msg( NVERB_ERROR, 'IO', 'OPEN_FILES', 'XHAT not found')
+
+   iju = njmax_ll + 2 * JPHEXT
+   ALLOCATE( XYHAT(iju), XYHATM(iju) )
+   CALL IO_Field_read(INFILES(1)%TFILE,'YHAT',XYHAT,IRESP2)
+   !If not found in main file, try the fallback one
+   if ( options( OPTFALLBACK )%set .and. iresp2 /= 0 ) CALL IO_Field_read(INFILES(2)%TFILE,'YHAT',XYHAT,IRESP2)
+   if ( iresp2 /= 0 ) call Print_msg( NVERB_ERROR, 'IO', 'OPEN_FILES', 'YHAT not found')
+
+   ! Interpolations of positions to mass points
+   XXHATM(1:IIU-1) = 0.5 * XXHAT(1:IIU-1) + 0.5 * XXHAT(2:IIU)
+   XXHATM(  IIU)   = 1.5 * XXHAT(  IIU)   - 0.5 * XXHAT(IIU-1)
+
+   XYHATM(1:IJU-1) = 0.5 * XYHAT(1:IJU-1) + 0.5 * XYHAT(2:IJU)
+   XYHATM(  IJU)   = 1.5 * XYHAT(  IJU)   - 0.5 * XYHAT(IJU-1)
+
+   ! Set global domain boundaries (hypothesis: only 1 process)
+   XXHAT_ll => XXHAT
+   XYHAT_ll => XYHAT
+   XXHATM_ll => XXHATM
+   XYHATM_ll => XYHATM
+
+   CALL IO_Field_read(INFILES(1)%TFILE,'CARTESIAN',LCARTESIAN,IRESP2)
+   !If not found in main file, try the fallback one
+   if ( options( OPTFALLBACK )%set .and. iresp2 /= 0 ) CALL IO_Field_read(INFILES(2)%TFILE,'CARTESIAN',LCARTESIAN,IRESP2)
+   if ( iresp2 /= 0 ) call Print_msg( NVERB_ERROR, 'IO', 'OPEN_FILES', 'CARTESIAN not found')
+
+   CALL IO_Field_read(INFILES(1)%TFILE,'LAT0',XLAT0,IRESP2)
+   !If not found in main file, try the fallback one
+   if ( options( OPTFALLBACK )%set .and. iresp2 /= 0 ) CALL IO_Field_read(INFILES(2)%TFILE,'LAT0',XLAT0,IRESP2)
+   if ( iresp2 /= 0 ) call Print_msg( NVERB_ERROR, 'IO', 'OPEN_FILES', 'LAT0 not found')
+
+   CALL IO_Field_read(INFILES(1)%TFILE,'LON0',XLON0,IRESP2)
+   !If not found in main file, try the fallback one
+   if ( options( OPTFALLBACK )%set .and. iresp2 /= 0 ) CALL IO_Field_read(INFILES(2)%TFILE,'LON0',XLON0,IRESP2)
+   if ( iresp2 /= 0 ) call Print_msg( NVERB_ERROR, 'IO', 'OPEN_FILES', 'LON0 not found')
+
+   CALL IO_Field_read(INFILES(1)%TFILE,'BETA',XBETA,IRESP2)
+   !If not found in main file, try the fallback one
+   if ( options( OPTFALLBACK )%set .and. iresp2 /= 0 ) CALL IO_Field_read(INFILES(2)%TFILE,'BETA',XBETA,IRESP2)
+   if ( iresp2 /= 0 ) call Print_msg( NVERB_ERROR, 'IO', 'OPEN_FILES', 'BETA not found')
+
    IF (.NOT.LCARTESIAN) THEN
-     CALL IO_Field_read(INFILES(1)%TFILE,'RPK',   XRPK)
-     CALL IO_Field_read(INFILES(1)%TFILE,'LATORI',XLATORI)
-     CALL IO_Field_read(INFILES(1)%TFILE,'LONORI',XLONORI)
+     CALL IO_Field_read(INFILES(1)%TFILE,'RPK',   XRPK,   IRESP2)
+    !If not found in main file, try the fallback one
+    if ( options( OPTFALLBACK )%set .and. iresp2 /= 0 ) CALL IO_Field_read(INFILES(2)%TFILE,'RPK',   XRPK,IRESP2)
+    if ( iresp2 /= 0 ) call Print_msg( NVERB_ERROR, 'IO', 'OPEN_FILES', 'RPK not found')
+
+     CALL IO_Field_read(INFILES(1)%TFILE,'LATORI',XLATORI,IRESP2)
+    !If not found in main file, try the fallback one
+    if ( options( OPTFALLBACK )%set .and. iresp2 /= 0 ) CALL IO_Field_read(INFILES(2)%TFILE,'LATORI',XLATORI,IRESP2)
+    if ( iresp2 /= 0 ) call Print_msg( NVERB_ERROR, 'IO', 'OPEN_FILES', 'LATORI not found')
+
+     CALL IO_Field_read(INFILES(1)%TFILE,'LONORI',XLONORI,IRESP2)
+    !If not found in main file, try the fallback one
+    if ( options( OPTFALLBACK )%set .and. iresp2 /= 0 ) CALL IO_Field_read(INFILES(2)%TFILE,'LONORI',XLONORI,IRESP2)
+    if ( iresp2 /= 0 ) call Print_msg( NVERB_ERROR, 'IO', 'OPEN_FILES', 'LONORI not found')
    ENDIF
    !
    IF (TRIM(CPROGRAM_ORIG)/='PGD' .AND. TRIM(CPROGRAM_ORIG)/='NESPGD' .AND. TRIM(CPROGRAM_ORIG)/='ZOOMPG' &
        .AND. .NOT.(TRIM(CPROGRAM_ORIG)=='REAL' .AND. CSTORAGE_TYPE=='SU') ) THEN !condition to detect PREP_SURFEX
-     ALLOCATE(XZHAT(NKMAX+2*JPVEXT))
-     CALL IO_Field_read(INFILES(1)%TFILE,'ZHAT',XZHAT)
+     iku = nkmax + 2 * JPVEXT
+     ALLOCATE( XZHAT(iku), XZHATM(iku) )
+     CALL IO_Field_read(INFILES(1)%TFILE,'ZHAT',XZHAT,IRESP2)
+     !If not found in main file, try the fallback one
+     if ( options( OPTFALLBACK )%set .and. iresp2 /= 0 ) CALL IO_Field_read(INFILES(2)%TFILE,'ZHAT',XZHAT,IRESP2)
+     if ( iresp2 /= 0 ) call Print_msg( NVERB_ERROR, 'IO', 'OPEN_FILES', 'ZHAT not found')
+
+     ! Interpolations of vertical positions to mass points
+     CALL INTERP_VERGRID_TO_MASSPOINTS( XZHAT, XZHATM )
+
      ALLOCATE(LSLEVE)
-     CALL IO_Field_read(INFILES(1)%TFILE,'SLEVE',LSLEVE)
+     CALL IO_Field_read(INFILES(1)%TFILE,'SLEVE',LSLEVE,IRESP2)
+     !If not found in main file, try the fallback one
+     if ( options( OPTFALLBACK )%set .and. iresp2 /= 0 ) CALL IO_Field_read(INFILES(2)%TFILE,'SLEVE',LSLEVE,IRESP2)
+     if ( iresp2 /= 0 ) call Print_msg( NVERB_ERROR, 'IO', 'OPEN_FILES', 'SLEVE not found')
+
      ALLOCATE(TDTMOD)
      CALL IO_Field_read(INFILES(1)%TFILE,'DTMOD',TDTMOD,IRESP2)
      IF(IRESP2/=0) DEALLOCATE(TDTMOD)
+
      ALLOCATE(TDTCUR)
      CALL IO_Field_read(INFILES(1)%TFILE,'DTCUR',TDTCUR,IRESP2)
      IF(IRESP2/=0) DEALLOCATE(TDTCUR)
+
+     !If time values were not found, try to get it from the time coordinate
+     if ( .not. associated( tdtcur ) .and. infiles(1)%tfile%cformat == 'NETCDF4' ) then
+       gok = .false.
+
+       istatus = NF90_INQ_VARID( infiles(1)%tfile%nncid, 'time', ivar_id )
+       if ( istatus == NF90_NOERR ) then
+         allocate( tdtcur )
+         istatus = NF90_GET_VAR( infiles(1)%tfile%nncid, ivar_id, tdtcur%xtime )
+         if ( istatus == NF90_NOERR ) then
+           istatus = NF90_INQUIRE_ATTRIBUTE( infiles(1)%tfile%nncid, ivar_id, 'units', len = ilen )
+           if ( istatus == NF90_NOERR ) then
+             allocate( character(len = ilen ) :: yunits )
+             istatus = NF90_GET_ATT( infiles(1)%tfile%nncid, ivar_id, 'units', yunits )
+             ! Extract date from yunits
+             idx =  INDEX( yunits, 'since ' )
+             Read( yunits(idx+6 :idx+9 ) , '( I4.4 )' ) tdtcur%nyear
+             Read( yunits(idx+11:idx+12 ), '( I2.2 )' ) tdtcur%nmonth
+             Read( yunits(idx+14:idx+15 ), '( I2.2 )' ) tdtcur%nday
+
+             if ( .not. associated( tdtmod ) ) then
+               allocate( tdtmod )
+               tdtmod = tdtcur
+               tdtmod%xtime = 0.
+             end if
+
+             gok = .true.
+           end if
+         end if
+       end if
+
+       if ( .not. gok ) deallocate( tdtcur )
+     end if
    END IF
    !
    ! Outfiles
@@ -1055,8 +1242,13 @@ END DO
          KNFILES_OUT = KNFILES_OUT + 1
 
          idx = KNFILES_OUT
-         CALL IO_FILE_ADD2LIST(outfiles(idx)%TFILE,HOUTFILE,'MNH','WRITE', &
-                               HFORMAT='NETCDF4',OOLD=.TRUE.)
+         if ( options(OPTDIR)%set ) then
+           CALL IO_FILE_ADD2LIST(outfiles(idx)%TFILE,HOUTFILE,'MNH','WRITE', &
+                                 HFORMAT='NETCDF4',OOLD=.TRUE., hdirname = options(OPTDIR)%cvalue )
+         else
+           CALL IO_FILE_ADD2LIST(outfiles(idx)%TFILE,HOUTFILE,'MNH','WRITE', &
+                                 HFORMAT='NETCDF4',OOLD=.TRUE.)
+         end if
          CALL IO_FILE_OPEN(outfiles(idx)%TFILE,HPROGRAM_ORIG=CPROGRAM_ORIG)
 
          IF (options(OPTCOMPRESS)%set) THEN
@@ -1068,8 +1260,8 @@ END DO
            outfiles(idx)%tfile%LNCREDUCE_FLOAT_PRECISION = .TRUE.
          END IF
 
-         status = NF90_SET_FILL(outfiles(idx)%TFILE%NNCID,NF90_NOFILL,omode)
-         if ( status /= NF90_NOERR ) call IO_Err_handle_nc4( status, 'OPEN_FILES', 'NF90_SET_FILL', '' )
+         istatus = NF90_SET_FILL(outfiles(idx)%TFILE%NNCID,NF90_NOFILL,ioldmode)
+         if ( istatus /= NF90_NOERR ) call IO_Err_handle_nc4( istatus, 'OPEN_FILES', 'NF90_SET_FILL', '' )
        END IF ! .NOT.osplit
     ELSE
        !
@@ -1077,8 +1269,13 @@ END DO
        !
        KNFILES_OUT = KNFILES_OUT + 1
        idx = KNFILES_OUT
-       CALL IO_FILE_ADD2LIST(outfiles(idx)%TFILE,houtfile,'MNH','WRITE', &
-                             HFORMAT='LFI',KLFIVERB=0,OOLD=.TRUE.)
+       if ( options(OPTDIR)%set ) then
+         CALL IO_FILE_ADD2LIST(outfiles(idx)%TFILE,houtfile,'MNH','WRITE', &
+                               HFORMAT='LFI',KLFIVERB=0,OOLD=.TRUE., hdirname = options(OPTDIR)%cvalue )
+       else
+         CALL IO_FILE_ADD2LIST(outfiles(idx)%TFILE,houtfile,'MNH','WRITE', &
+                               HFORMAT='LFI',KLFIVERB=0,OOLD=.TRUE.)
+       end if
        LIOCDF4 = .FALSE. !Necessary to open correctly the LFI file
        CALL IO_FILE_OPEN(outfiles(idx)%TFILE,HPROGRAM_ORIG=CPROGRAM_ORIG)
        LIOCDF4 = .TRUE.
@@ -1089,12 +1286,17 @@ END DO
      KNFILES_OUT = KNFILES_OUT + 1
 
      idx = KNFILES_OUT
-     CALL IO_FILE_ADD2LIST(outfiles(idx)%TFILE,'dummy_file','MNH','WRITE', &
-                           HFORMAT='NETCDF4',OOLD=.TRUE.)
+     if ( options(OPTDIR)%set ) then
+       CALL IO_FILE_ADD2LIST(outfiles(idx)%TFILE,'dummy_file','MNH','WRITE', &
+                             HFORMAT='NETCDF4',OOLD=.TRUE., hdirname = options(OPTDIR)%cvalue )
+     else
+       CALL IO_FILE_ADD2LIST(outfiles(idx)%TFILE,'dummy_file','MNH','WRITE', &
+                             HFORMAT='NETCDF4',OOLD=.TRUE.)
+     end if
      CALL IO_FILE_OPEN(outfiles(idx)%TFILE,HPROGRAM_ORIG=CPROGRAM_ORIG)
    END IF
 
-   PRINT *,'--> Converted to file: ', TRIM(houtfile)
+   call Print_msg( NVERB_INFO, 'IO', 'parse_infiles', '--> Converted to file: ' // trim(houtfile) )
 
   END SUBROUTINE OPEN_FILES
 
@@ -1114,7 +1316,7 @@ END DO
     INTEGER                  :: ji
     INTEGER                  :: idx1, idx2
     INTEGER(KIND=CDFINT)     :: status
-    INTEGER(KIND=CDFINT)     :: omode
+    INTEGER(KIND=CDFINT)     :: ioldmode
 
     CALL PRINT_MSG(NVERB_DEBUG,'IO','OPEN_SPLIT_NCFILES_OUT','called')
 
@@ -1144,8 +1346,13 @@ END DO
 
     DO ji = 1,nbvar
       filename = trim(houtfile)//'.'//TRIM(YVARS(ji))
-      CALL IO_FILE_ADD2LIST(outfiles(ji)%TFILE,filename,'MNH','WRITE', &
-                            HFORMAT='NETCDF4')
+      if ( options(OPTDIR)%set ) then
+        CALL IO_FILE_ADD2LIST(outfiles(ji)%TFILE,filename,'MNH','WRITE', &
+                              HFORMAT='NETCDF4', hdirname = options(OPTDIR)%cvalue )
+      else
+        CALL IO_FILE_ADD2LIST(outfiles(ji)%TFILE,filename,'MNH','WRITE', &
+                              HFORMAT='NETCDF4')
+      end if
       CALL IO_FILE_OPEN(outfiles(ji)%TFILE,HPROGRAM_ORIG=CPROGRAM_ORIG)
 
       IF (options(OPTCOMPRESS)%set) THEN
@@ -1157,7 +1364,7 @@ END DO
         outfiles(ji)%tfile%LNCREDUCE_FLOAT_PRECISION = .TRUE.
       END IF
 
-      status = NF90_SET_FILL(outfiles(ji)%TFILE%NNCID,NF90_NOFILL,omode)
+      status = NF90_SET_FILL(outfiles(ji)%TFILE%NNCID,NF90_NOFILL,ioldmode)
       if ( status /= NF90_NOERR ) call IO_Err_handle_nc4( status, 'OPEN_SPLIT_NCFILES_OUT', 'NF90_SET_FILL', '' )
     END DO
 
@@ -1181,31 +1388,126 @@ END DO
   END SUBROUTINE CLOSE_FILES
 
 
-  SUBROUTINE IO_Metadata_get_nc4(KFILE_ID,KVAR_ID,TPREC)
+  SUBROUTINE IO_Metadata_get_nc4(TPFILE,KVAR_ID,TPREC)
     USE MODD_DIM_n,      ONLY: NKMAX
     USE MODD_PARAMETERS, ONLY: JPVEXT
 
-    INTEGER(KIND=CDFINT), INTENT(IN)    :: KFILE_ID
+    TYPE(TFILEDATA),      INTENT(IN)    :: TPFILE
     INTEGER(KIND=CDFINT), INTENT(IN)    :: KVAR_ID
     TYPE(workfield),      INTENT(INOUT) :: TPREC
 
+    character(len=:), allocatable            :: YSPLIT
+    character(len=:), allocatable            :: YTIMEDEP
+    integer                                  :: iblocks
     INTEGER                                  :: ILENG
     INTEGER                                  :: JDIM
     INTEGER(KIND=CDFINT)                     :: ISTATUS
+    INTEGER(KIND=CDFINT)                     :: IFILE_ID
+    INTEGER(KIND=CDFINT)                     :: IVAR_ID
     INTEGER(KIND=CDFINT),DIMENSION(NF90_MAX_VAR_DIMS) :: IDIMS_ID
+    LOGICAL                                  :: GSPLIT_AT_ENTRY
+    LOGICAL                                  :: GSPLIT_INFO_AVAILABLE
 
     CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Metadata_get_nc4','called')
 
-    ISTATUS = NF90_INQUIRE_VARIABLE(KFILE_ID,KVAR_ID,NDIMS = TPREC%NDIMS_FILE, &
+    !Necessary to know if we already are in a split file for determining correct number of dimensions
+    GSPLIT_AT_ENTRY = TPREC%LSPLIT
+
+    IFILE_ID = TPFILE%NNCID
+
+    iblocks = -1
+
+    ISTATUS = NF90_INQUIRE_VARIABLE(IFILE_ID, KVAR_ID, NDIMS = TPREC%NDIMS_FILE, &
                                     XTYPE = TPREC%NTYPE_FILE, DIMIDS = IDIMS_ID)
     if ( istatus /= NF90_NOERR ) call IO_Err_handle_nc4( istatus, 'IO_Metadata_get_nc4', 'NF90_INQUIRE_VARIABLE', '' )
 
+    !split_variable and other attributes were added in MesoNH > 5.4.2
+    GSPLIT_INFO_AVAILABLE = .FALSE.
+    ISTATUS = NF90_INQUIRE_ATTRIBUTE(IFILE_ID, KVAR_ID, 'split_variable', LEN=ILENG)
+    IF (ISTATUS == NF90_NOERR) THEN
+      GSPLIT_INFO_AVAILABLE = .TRUE.
+
+      IF (GSPLIT_AT_ENTRY) CALL PRINT_MSG(NVERB_ERROR,'IO','IO_Metadata_get_nc4','split variable declaration inside a split file')
+
+      ALLOCATE(CHARACTER(LEN=ILENG) :: YSPLIT)
+      ISTATUS = NF90_GET_ATT(IFILE_ID, KVAR_ID, 'split_variable', YSPLIT)
+      IF (istatus /= NF90_NOERR) CALL IO_Err_handle_nc4( istatus, 'IO_Metadata_get_nc4', 'NF90_GET_ATT', 'split_variable' )
+      IF ( YSPLIT == 'yes' ) then
+        TPREC%LSPLIT = .true.
+
+        ISTATUS = NF90_GET_ATT(IFILE_ID, KVAR_ID, 'ndims', TPREC%NDIMS_FILE)
+        IF (istatus /= NF90_NOERR) CALL IO_Err_handle_nc4( istatus, 'IO_Metadata_get_nc4', 'NF90_GET_ATT', 'ndims' )
+        IF ( TPREC%NDIMS_FILE/=3 ) CALL PRINT_MSG(NVERB_ERROR,'IO','IO_Metadata_get_nc4', &
+                                                  'split variable with ndims/=3 not supported')
+
+        ISTATUS = NF90_INQUIRE_ATTRIBUTE(IFILE_ID, KVAR_ID, 'time_dependent', LEN=ILENG)
+        IF (istatus /= NF90_NOERR) CALL IO_Err_handle_nc4( istatus, 'IO_Metadata_get_nc4', 'NF90_INQUIRE_ATTRIBUTE', &
+                                                           'time_dependent' )
+        ALLOCATE(CHARACTER(LEN=ILENG) :: YTIMEDEP)
+        ISTATUS = NF90_GET_ATT(IFILE_ID, KVAR_ID, 'time_dependent', YTIMEDEP)
+        IF (istatus /= NF90_NOERR) CALL IO_Err_handle_nc4( istatus, 'IO_Metadata_get_nc4', 'NF90_GET_ATT', 'time_dependent' )
+        IF ( YTIMEDEP == 'yes' ) then
+          TPREC%TFIELD%LTIMEDEP = .TRUE.
+        ELSE IF ( YTIMEDEP == 'no' ) THEN
+          TPREC%TFIELD%LTIMEDEP = .FALSE.
+        ELSE
+          CALL PRINT_MSG(NVERB_WARNING,'IO','IO_Metadata_get_nc4','unknown value '//trim(YTIMEDEP)// &
+                                                                  ' for time_dependent attribute' )
+        END IF
+
+        ISTATUS = NF90_GET_ATT(IFILE_ID, KVAR_ID, 'split_nblocks', iblocks)
+        IF (istatus /= NF90_NOERR) CALL IO_Err_handle_nc4( istatus, 'IO_Metadata_get_nc4', 'NF90_GET_ATT', 'split_nblocks' )
+
+        IFILE_ID = TPFILE%TFILES_IOZ(1)%TFILE%NNCID
+
+        istatus = NF90_INQ_VARID(IFILE_ID,trim(TPREC%NAME)//'0001',ivar_id)
+        IF (ISTATUS /= NF90_NOERR) CALL IO_Err_handle_nc4( istatus, 'IO_Metadata_get_nc4', 'NF90_INQ_VARID', &
+                                                           trim(TPREC%NAME)//'0001' )
+        ISTATUS = NF90_INQUIRE_VARIABLE(IFILE_ID, IVAR_ID, DIMIDS = IDIMS_ID)
+        IF (ISTATUS /= NF90_NOERR) CALL IO_Err_handle_nc4( istatus, 'IO_Metadata_get_nc4', 'NF90_INQUIRE_VARIABLE',&
+                                                           trim(TPREC%NAME)//'0001' )
+
+        DEALLOCATE(YTIMEDEP)
+      ELSE IF ( YSPLIT /= 'no' ) THEN
+        CALL PRINT_MSG(NVERB_WARNING,'IO','IO_Metadata_get_nc4','unknown value '//trim(YSPLIT)//' for split_variable attribute' )
+      END IF
+
+      DEALLOCATE(YSPLIT)
+    END IF
+
+    !Reset IFILE_ID to master file (if split files)
+    IFILE_ID = TPFILE%NNCID
+
+    ISTATUS = NF90_GET_ATT(IFILE_ID,KVAR_ID,'long_name',TPREC%TFIELD%CLONGNAME)
+    IF (ISTATUS /= NF90_NOERR) TPREC%TFIELD%CLONGNAME = TRIM( TPREC%TFIELD%CMNHNAME )
+
+    ISTATUS = NF90_GET_ATT(IFILE_ID,KVAR_ID,'comment',TPREC%TFIELD%CCOMMENT)
+    IF (ISTATUS /= NF90_NOERR) TPREC%TFIELD%CCOMMENT = ''
+
+    ISTATUS = NF90_GET_ATT(IFILE_ID,KVAR_ID,'grid',TPREC%NGRID_FILE)
+    !On MesoNH versions < 5.4.0, the grid number was stored in 'GRID' instead of 'grid'
+    IF (ISTATUS /= NF90_NOERR) ISTATUS = NF90_GET_ATT(IFILE_ID,KVAR_ID,'GRID',TPREC%NGRID_FILE)
+    IF (ISTATUS /= NF90_NOERR) TPREC%NGRID_FILE = 0
+
+    ISTATUS = NF90_GET_ATT(IFILE_ID,KVAR_ID,'units',TPREC%CUNITS_FILE)
+    IF (ISTATUS /= NF90_NOERR) TPREC%CUNITS_FILE = ''
+
     IF (.NOT.TPREC%LSPLIT) THEN
       ALLOCATE(TPREC%NDIMSIZES_FILE(TPREC%NDIMS_FILE))
       ALLOCATE(TPREC%CDIMNAMES_FILE(TPREC%NDIMS_FILE))
     ELSE
-      ALLOCATE(TPREC%NDIMSIZES_FILE(TPREC%NDIMS_FILE+1))
-      ALLOCATE(TPREC%CDIMNAMES_FILE(TPREC%NDIMS_FILE+1))
+      IF ( GSPLIT_AT_ENTRY ) THEN
+        ALLOCATE(TPREC%NDIMSIZES_FILE(TPREC%NDIMS_FILE+1))
+        ALLOCATE(TPREC%CDIMNAMES_FILE(TPREC%NDIMS_FILE+1))
+      ELSE
+        IF (TPREC%TFIELD%LTIMEDEP) THEN
+          ALLOCATE(TPREC%NDIMSIZES_FILE(TPREC%NDIMS_FILE+1))
+          ALLOCATE(TPREC%CDIMNAMES_FILE(TPREC%NDIMS_FILE+1))
+        ELSE
+          ALLOCATE(TPREC%NDIMSIZES_FILE(TPREC%NDIMS_FILE))
+          ALLOCATE(TPREC%CDIMNAMES_FILE(TPREC%NDIMS_FILE))
+        END IF
+      END IF
     END IF
 
     IF (TPREC%NDIMS_FILE == 0) THEN
@@ -1215,7 +1517,7 @@ END DO
       ! Fill dimensions info
       ILENG = 1
       DO JDIM=1,TPREC%NDIMS_FILE
-        ISTATUS = NF90_INQUIRE_DIMENSION(KFILE_ID,IDIMS_ID(JDIM),                    &
+        ISTATUS = NF90_INQUIRE_DIMENSION(IFILE_ID,IDIMS_ID(JDIM),                    &
                                                    len =  TPREC%NDIMSIZES_FILE(JDIM), &
                                                    name = TPREC%CDIMNAMES_FILE(JDIM)  )
         if ( istatus /= NF90_NOERR ) call IO_Err_handle_nc4( istatus, 'IO_Metadata_get_nc4', 'NF90_INQUIRE_DIMENSION', '' )
@@ -1233,44 +1535,63 @@ END DO
       END IF
 
       IF (TPREC%LSPLIT) THEN
+#if 0
         IF(     (.NOT.TPREC%TFIELD%LTIMEDEP .AND.  TPREC%NDIMS_FILE/=2)   &
             .OR. (     TPREC%TFIELD%LTIMEDEP .AND.  TPREC%NDIMS_FILE/=3) ) &
-          CALL PRINT_MSG(NVERB_FATAL,'IO','parse_infiles','split variables can only be 3D')
-          !Split variables are Z-split
-          ILENG = ILENG * (NKMAX+2*JPVEXT)
-          !Move time dimension to last (4th) position
-          IF (TPREC%TFIELD%LTIMEDEP) THEN
-            TPREC%NDIMSIZES_FILE(4) = TPREC%NDIMSIZES_FILE(3)
-            TPREC%CDIMNAMES_FILE(4) = TPREC%CDIMNAMES_FILE(3)
-          END IF
-          !Add vertical dimension
-          TPREC%NDIMSIZES_FILE(3) = NKMAX+2*JPVEXT
-          TPREC%CDIMNAMES_FILE(3) = 'level' !Could also be 'level_w'
+#else
+        IF(      (      GSPLIT_INFO_AVAILABLE .AND.  TPREC%NDIMS_FILE/=3 )               &
+            .OR. ( .NOT.GSPLIT_INFO_AVAILABLE .AND.                                      &
+                            ( (.NOT.TPREC%TFIELD%LTIMEDEP .AND.  TPREC%NDIMS_FILE/=2)    &
+                      .OR. (     TPREC%TFIELD%LTIMEDEP .AND.  TPREC%NDIMS_FILE/=3) ) ) ) &
+#endif
+          CALL PRINT_MSG(NVERB_FATAL,'IO','IO_Metadata_get_nc4',trim(TPREC%NAME)//': split variables can only be 3D')
+        !Split variables are Z-split
+        !Move time dimension to last (4th) position
+        IF (TPREC%TFIELD%LTIMEDEP) THEN
+          TPREC%NDIMSIZES_FILE(4) = TPREC%NDIMSIZES_FILE(3)
+          TPREC%CDIMNAMES_FILE(4) = TPREC%CDIMNAMES_FILE(3)
         END IF
-      END IF
 
-      TPREC%NSIZE = ILENG
+        !Add vertical/3rd dimension
+        SELECT CASE(TPREC%NGRID_FILE)
+          CASE (1, 2, 3, 5)
+            TPREC%CDIMNAMES_FILE(3) = 'level'
+          CASE (4, 6, 7, 8)
+            TPREC%CDIMNAMES_FILE(3) = 'level_w'
+          CASE DEFAULT
+            TPREC%CDIMNAMES_FILE(3) = 'unknown'
+        END SELECT
+
+        IF (iblocks == -1 ) then
+          TPREC%NDIMSIZES_FILE(3) = NKMAX+2*JPVEXT
+        else
 
-      ISTATUS = NF90_GET_ATT(KFILE_ID,KVAR_ID,'grid',TPREC%NGRID_FILE)
-      !On MesoNH versions < 5.4.0, the grid number was stored in 'GRID' instead of 'grid'
-      IF (ISTATUS /= NF90_NOERR) ISTATUS = NF90_GET_ATT(KFILE_ID,KVAR_ID,'GRID',TPREC%NGRID_FILE)
-      IF (ISTATUS /= NF90_NOERR) TPREC%NGRID_FILE = 0
+          if (TPREC%NGRID_FILE/=0 .and. iblocks/=NKMAX+2*JPVEXT) THEN
+            !If size is not as expected, reset its name
+            CALL PRINT_MSG(NVERB_WARNING,'IO','IO_Metadata_get_nc4',trim(TPREC%NAME)//': strange nblocks size')
+            TPREC%CDIMNAMES_FILE(3) = 'unknown'
+          end if
+          TPREC%NDIMSIZES_FILE(3) = iblocks
+        end if
+
+        ILENG = ILENG * TPREC%NDIMSIZES_FILE(3)
+      END IF
+    END IF
 
-      ISTATUS = NF90_GET_ATT(KFILE_ID,KVAR_ID,'units',TPREC%CUNITS_FILE)
-      IF (ISTATUS /= NF90_NOERR) TPREC%CUNITS_FILE = ''
+    TPREC%NSIZE = ILENG
   END SUBROUTINE IO_Metadata_get_nc4
 
 
   SUBROUTINE IO_Dims_fill_nc4(TPFILE,TPREC,KRESP)
     USE MODD_IO,           ONLY: TFILEDATA
-    use mode_io_tools_nc4, only: IO_Dimcdf_get_nc4, IO_Dim_find_byname_nc4
+    use mode_io_tools_nc4, only: IO_Dim_find_create_nc4, IO_Dim_find_byname_nc4
 
     TYPE(TFILEDATA),INTENT(IN)    :: TPFILE
     TYPE(workfield),INTENT(INOUT) :: TPREC
     INTEGER,        INTENT(OUT)   :: KRESP
 
+    integer              :: iidx
     INTEGER              :: JJ
-    TYPE(DIMCDF),POINTER :: TZDIMPTR
 
     CALL PRINT_MSG(NVERB_DEBUG,'IO','IO_Dims_fill_nc4','called')
 
@@ -1292,13 +1613,13 @@ END DO
       CALL IO_Dim_find_byname_nc4(TPFILE,TPREC%CDIMNAMES_FILE(JJ),TPREC%TDIMS(JJ),KRESP)
       !If dimension not found => create it
       IF (KRESP/=0)  THEN
-        TZDIMPTR => IO_Dimcdf_get_nc4(TPFILE,TPREC%NDIMSIZES_FILE(JJ))
-        TPREC%TDIMS(JJ) = TZDIMPTR
+        call IO_Dim_find_create_nc4( tpfile, tprec%ndimsizes_file(jj), iidx )
+        tprec%tdims(jj) = tpfile%tncdims%tdims(iidx)
         KRESP = 0
       END IF
-      IF (TRIM(TPREC%TDIMS(JJ)%name)/='time' .AND. &
-        TPREC%TDIMS(JJ)%len /= TPREC%NDIMSIZES_FILE(JJ)) THEN
-        CALL PRINT_MSG(NVERB_WARNING,'IO','parse_infiles','problem with dimensions for '//TPREC%TFIELD%CMNHNAME)
+      IF (TRIM(TPREC%TDIMS(JJ)%cname)/='time' .AND. &
+        TPREC%TDIMS(JJ)%nlen /= TPREC%NDIMSIZES_FILE(JJ)) THEN
+        CALL PRINT_MSG(NVERB_WARNING,'IO','IO_Dims_fill_nc4','problem with dimensions for '//TPREC%TFIELD%CMNHNAME)
         KRESP = -3
         EXIT
       END IF

MY_RUN/BENCH/003_MEDATL_FULLHD/003_mesonh/run_mesonh_xyz_SEGN.ihm

-#MNH_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
+#MNH_LIC Copyright 1994-2024 CNRS, Meteo-France and Universite Paul Sabatier
 #MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
 #MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
 #MNH_LIC for details. version 1.
@@ -54,7 +54,6 @@ cat > EXSEG1.nam  << EOF
 &NAM_DYNn 
  XTSTEP = ${XTSTEP} ,
  CPRESOPT = "CRESI",
- NITR = 12, 
  LITRADJ = T,
  XRELAX = 1.,
  LHORELAX_UVWTH = T,

MY_RUN/INTEGRATION_CASES/HPC/16JAN/PYTHON/plot_16JAN.py

+#!/usr/bin/env python3
+"""
+@author: Quentin Rodier
+Creation : 07/01/2021
+
+Last modifications
+"""
+import matplotlib as mpl
+mpl.use('Agg')
+import cartopy.crs as ccrs
+from read_MNHfile import read_netcdf
+from Panel_Plot import PanelPlot
+import os
+
+os.system('rm -f tempgraph*')
+#
+#  User's parameter / Namelist
+#
+#
+path="../RUN/"
+LnameFiles = ['16JA1.1.WENO5.004dg.nc', '16JA2.2.WENO5.002dg.nc']
+
+Dvar_input = {
+'f1':['MRV700HPA','THT850HPA','UT850HPA','VT850HPA','UT700HPA','VT700HPA', 'ALT_PRESSURE','ALT_U','ALT_V', 'ZS', 'latitude', 'longitude'],
+'f2':['MRV700HPA','THT850HPA','UT850HPA','VT850HPA','UT700HPA','VT700HPA', 'ALT_PRESSURE', 'ZS', 'ALT_U','ALT_V','latitude', 'longitude']
+}
+
+#  Read the variables in the files
+Dvar = {}
+Dvar = read_netcdf(LnameFiles, Dvar_input, path=path, removeHALO=True)
+
+################################################################
+#########          PANEL 1  
+###############################################################
+Panel1 = PanelPlot(2,2, [15,10],'16 JAN domaine 1 ', minmaxpad=1.05, minmaxTextSize=7.7, colorbaraspect=14)
+
+Lplot = [ Dvar['f1']['ZS'],Dvar['f1']['THT850HPA'], Dvar['f1']['MRV700HPA'],Dvar['f1']['ALT_PRESSURE']]
+lon = [Dvar['f1']['longitude']]*len(Lplot)
+lat = [Dvar['f1']['latitude']]*len(Lplot)
+Ltitle = ['Orography', 'Potential Temperature at 850hPa', 'Water vapor mixing at 700hPa','Pressure at z = 9000m']
+Lcbarlabel = ['m','K', 'g/kg', 'hPa']
+Lxlab = ['longitude']*len(Lplot)
+Lylab = ['latitude']*len(Lplot)
+Lminval = [0, 280, 0.5, 270]
+Lmaxval = [1500, 293, 4.2, 303]
+Lstep = [10, 0.25, 0.1, 1]
+Lstepticks = [500, 2, 0.5, 5]
+Lfacconv = [1.0, 1.0, 1.0, 1./100.0]
+Lcolormap = ['terrain', 'gist_rainbow_r', 'gist_rainbow_r', 'gist_rainbow_r']
+Lprojection = [ccrs.PlateCarree()]*len(Lplot)
+Llvl = [0, 0, 0, 0]
+fig1 = Panel1.psectionH(lon=lon, lat=lat, Lvar=Lplot, Lcarte=[], Llevel=Llvl, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, 
+                                Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lproj=Lprojection, Lfacconv=Lfacconv)
+
+Lplot1 = [ Dvar['f1']['UT850HPA'], Dvar['f1']['UT700HPA'], Dvar['f1']['ALT_U']]
+Lplot2 = [ Dvar['f1']['VT850HPA'], Dvar['f1']['VT700HPA'], Dvar['f1']['ALT_V']]
+Ltitle = ['Wind at 850hPa', 'Wind at 700hPa', 'Wind at 9000m']
+Lxlab = ['longitude']*len(Lplot1)
+Lylab = ['latitude']*len(Lplot1)
+Llegendval = [20,20,40]
+Llegendlabel = ['(m/s)']*len(Lplot1)
+Larrowstep = [2]*len(Lplot1)
+Lwidth = [0.002]*len(Lplot1)
+Lcolor = ['black']*len(Lplot1)
+Lprojection = [ccrs.PlateCarree()]*len(Lplot1)
+Llvl = [0]*len(Lplot1)
+fig2 = Panel1.pvector(Lxx=lon, Lyy=lat, Lvar1=Lplot1, Lvar2=Lplot2, Lcarte=[], Llevel=Llvl, Lxlab=Lxlab, Lylab=Lylab, 
+                      Ltitle=Ltitle, Lwidth=Lwidth, Larrowstep=Larrowstep, Lproj=Lprojection,
+                      Lcolor=Lcolor, Llegendval=Llegendval, Llegendlabel=Llegendlabel, Lid_overlap=[2,4,6], ax=fig1.axes)
+
+Panel1.save_graph(1,fig2)
+
+################################################################
+#########          PANEL 2 
+###############################################################
+Panel2 = PanelPlot(2,2, [15,10],'16 JAN domaine 2 ', minmaxpad=1.05, minmaxTextSize=8, colorbaraspect=14)
+
+Lplot = [ Dvar['f2']['ZS'],Dvar['f2']['THT850HPA'], Dvar['f2']['MRV700HPA'],Dvar['f2']['ALT_PRESSURE']]
+lon = [Dvar['f2']['longitude']]*len(Lplot)
+lat = [Dvar['f2']['latitude']]*len(Lplot)
+Ltitle = ['Orography', 'Potential Temperature at 850hPa', 'Water vapor mixing at 700hPa','Pressure at z = 9000m']
+Lcbarlabel = ['m','K', 'g/kg', 'hPa']
+Lxlab = ['longitude']*len(Lplot)
+Lylab = ['latitude']*len(Lplot)
+Lminval = [0, 285, 0.9, 280]
+Lmaxval = [2600, 293, 4, 294]
+Lstep = [10, 0.5, 0.1, 0.5]
+Lstepticks = [500, 1, 0.2, 5]
+Lfacconv = [1.0, 1.0, 1.0, 1./100.0]
+Lcolormap = ['terrain', 'gist_rainbow_r', 'gist_rainbow_r', 'gist_rainbow_r']
+Lprojection = [ccrs.PlateCarree()]*len(Lplot)
+Llvl = [0]*len(Lplot)
+fig1 = Panel2.psectionH(lon=lon, lat=lat, Lvar=Lplot, Lcarte=[], Llevel=Llvl, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, 
+                                Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lproj=Lprojection, Lfacconv=Lfacconv)
+
+Lplot1 = [ Dvar['f2']['UT850HPA'], Dvar['f2']['UT700HPA'], Dvar['f2']['ALT_U']]
+Lplot2 = [ Dvar['f2']['VT850HPA'], Dvar['f2']['VT700HPA'], Dvar['f2']['ALT_V']]
+Ltitle = ['Wind at 850hPa', 'Wind at 700hPa', 'Wind at 9000m']
+Llegendval = [20,20,40]
+Lxlab = ['longitude']*len(Lplot1)
+Lylab = ['latitude']*len(Lplot1)
+Llegendlabel = ['(m/s)']*len(Lplot1)
+Larrowstep = [2]*len(Lplot1)
+Lwidth = [0.002]*len(Lplot1)
+Lcolor = ['black']*len(Lplot1)
+Lprojection = [ccrs.PlateCarree()]*len(Lplot1)
+Llvl = [0]*len(Lplot1)
+fig2 = Panel2.pvector(Lxx=lon, Lyy=lat, Lvar1=Lplot1, Lvar2=Lplot2, Lcarte=[], Llevel=Llvl, Lxlab=Lxlab, Lylab=Lylab, 
+                      Ltitle=Ltitle, Lwidth=Lwidth, Larrowstep=Larrowstep, Lproj=Lprojection,
+                      Lcolor=Lcolor, Llegendval=Llegendval, Llegendlabel=Llegendlabel, Lid_overlap=[2,4,6], ax=fig1.axes)
+
+Panel2.save_graph(2,fig2)

MY_RUN/INTEGRATION_CASES/HPC/16JAN/PYTHON/run_python

+#!/bin/sh
+#SBATCH -J tools
+#SBATCH -N 1            # nodes number
+#SBATCH -n 1            # CPUs number (on all nodes) 
+#SBATCH -o tools.eo%j   #
+#SBATCH -e tools.eo%j   #
+#SBATCH -t 01:00:00    # time limit
+# Echo des commandes
+ulimit -c 0
+ulimit -s unlimited
+# Arrete du job des la premiere erreur
+set -e
+
+. ~rodierq/DEV_57/MNH-PHYEX070-b95d84d7/conf/profile_mesonh-LXifort-R8I4-MNH-V5-6-2-ECRAD140-MPIAUTO-O2
+
+ln -sf ${SRC_MESONH}/src/LIB/Python/* .
+
+module purge
+module load python/3.7.6
+
+python3 plot_16JAN.py
+convert *.png 16JAN.pdf

MY_RUN/INTEGRATION_CASES/HPC/16JAN/RUN/DIAG1.nam_1

+&NAM_CONFIO  LCDF4=T, LLFIOUT=F, LLFIREAD=F /
+&NAM_DIAG             
+  CISO='TKPREV',
+  LVAR_RS=T, 
+  LVAR_TURB=F,
+  NCONV_KF=1, 
+  LVAR_MRW=T,
+  LVAR_MRSV=F,
+  LTRAJ=F,
+  LTPZH=T,
+  CBLTOP="RICHA",
+  LMSLP=T,
+  LAGEO=T, 
+  LTHW=T, 
+  LCLD_COV=T, 
+  NCAPE=1,
+  LRADAR=T,
+  LDIAG(:)=.FALSE.,
+  LISOAL=T, XISOAL(1)=9000,
+  LISOPR=T, XISOPR(1)=850, XISOPR(2)=700 /
+
+&NAM_DIAG_FILE      YINIFILE(1) = "16JA1.1.WENO5.004" ,
+                    YINIFILEPGD(1) = "16JAN98_36km.neste1",
+                    YSUFFIX='dg' /
+&NAM_DIAG_SURFn  N2M=2 LSURF_BUDGET=T /
+&NAM_DIAG_ISBAn LPGD=F LSURF_EVAP_BUDGET=T /
+

MY_RUN/INTEGRATION_CASES/HPC/16JAN/RUN/DIAG1.nam_2

+&NAM_CONFIO  LCDF4=T, LLFIOUT=F, LLFIREAD=F /
+&NAM_DIAG CISO='', LVAR_RS=F, LVAR_LS=F, LVAR_PR=T, LTOTAL_PR=T,
+  LISOAL=T, XISOAL(1)=9000,
+  LISOPR=T, XISOPR(1)=850, XISOPR(2)=700 /
+&NAM_BLANK /
+&NAM_DIAG_FILE  YINIFILE(1) = "16JA2.2.WENO5.002",
+                YINIFILEPGD(1) = "16JAN98_9km.neste1",
+                YSUFFIX = "dg" /

MY_RUN/INTEGRATION_CASES/HPC/16JAN/RUN/EXSEG1.nam_run1

+&NAM_CONFIO  LCDF4=T, LLFIOUT=F, LLFIREAD=F /
+&NAM_LUNITn CINIFILE = "15JAN_12_MNH" ,
+            CINIFILEPGD = "16JAN98_36km.neste1",
+            CCPLFILE(1) = "15JAN_18_MNH",
+            CCPLFILE(2) = "16JAN_00_MNH" /
+&NAM_CONFZ
+   NZ_VERB=5 , NZ_PROC=0 , NB_PROCIO_R=1 , NB_PROCIO_W=1
+/
+&NAM_DYNn XTSTEP = 120.,  CPRESOPT = "ZRESI", NITR = 8,
+           LHORELAX_UVWTH = T, LHORELAX_RV = T, LVE_RELAX = T,
+           NRIMX = 5, NRIMY = 5, XRIMKMAX = 0.0083, XT4DIFU = 5000.,XT4DIFTH = 5000. /
+&NAM_ADVn  CUVW_ADV_SCHEME = "WENO_K" ,CMET_ADV_SCHEME = "PPM_01",NWENO_ORDER=5,CTEMP_SCHEME='RK53' /
+&NAM_PARAMn CTURB = "TKEL", CRAD = "ECMW",
+           CDCONV = "KAFR", CCLOUD= "KESS",CSCONV="KAFR" /
+&NAM_PARAM_RADn XDTRAD = 3600., XDTRAD_CLONLY = 3600., NRAD_COLNBR= 400 /
+&NAM_PARAM_KAFRn XDTCONV = 300., NICE = 1, LREFRESH_ALL = T, LDOWN = T /
+&NAM_LBCn CLBCX = 2*"OPEN", CLBCY = 2*"OPEN" /
+&NAM_TURBn CTURBLEN = "BL89", CTURBDIM = "1DIM",
+           LTURB_DIAG = F, LTURB_FLX = F /
+&NAM_CONF CCONF = "START", NMODEL = 1, NVERB = 5,
+           CEXP = "16JA1 ", CSEG = "WENO5",NHALO=3 /
+&NAM_DYN  XSEGLEN = 43200., LCORIO = T, LNUMDIFU = F,LNUMDIFTH = F,
+          XALKTOP = 0.001, XALZBOT = 15000. /
+&NAM_BACKUP XBAK_TIME(1,1)=10800., XBAK_TIME(1,2)=21600.,XBAK_TIME(1,3)=32400.,XBAK_TIME(1,4)=43200./
+&NAM_DIAG_SURFn  /
+&NAM_ISBA CALBEDO = "DRY" /
+&NAM_ISBAn CSCOND  = "NP89", 
+           CC1DRY='DEF', CSOILFRZ='DEF', CDIFSFCOND='DEF', CSNOWRES='DEF' /
+&NAM_SSOn CROUGH="Z01D"/
+&NAM_DIAG_ISBAn /
+&NAM_SEAFLUXn CSEA_ALB="UNIF" /
+&NAM_BLANK /
+&NAM_SGH_ISBAn CRUNOFF = "WSAT"/
+&NAM_NEBn LSUBG_COND= F   /

MY_RUN/INTEGRATION_CASES/HPC/16JAN/RUN/EXSEG1.nam_run2

+&NAM_CONFIO  LCDF4=T, LLFIOUT=F, LLFIREAD=F /
+&NAM_LUNITn  CINIFILE = "16JA1.1.WENO5.004",
+             CINIFILEPGD = "16JAN98_36km.neste1",
+             CCPLFILE(1) = "16JAN_06_MNH" /
+&NAM_DYNn  XTSTEP = 120.,  CPRESOPT = "ZRESI", NITR = 8,
+           LHORELAX_UVWTH = T, LHORELAX_RV = T, LVE_RELAX = T,
+           NRIMX = 5, NRIMY = 5, XRIMKMAX = 0.0083, XT4DIFU = 5000. /
+&NAM_ADVn  CUVW_ADV_SCHEME = "WENO_K", CMET_ADV_SCHEME = "PPM_01",NWENO_ORDER=5,CTEMP_SCHEME='RK53' /
+&NAM_PARAMn CTURB = "TKEL", CRAD = "ECMW",
+           CDCONV = "KAFR", CSCONV="KAFR",CCLOUD= "KESS" /
+&NAM_PARAM_RADn  XDTRAD = 3600., XDTRAD_CLONLY = 3600., NRAD_COLNBR = 400 /
+&NAM_LBCn  CLBCX = 2*"OPEN", CLBCY = 2*"OPEN" /
+&NAM_PARAM_KAFRn XDTCONV = 300., NICE = 1, LREFRESH_ALL = T,
+           LDOWN = T /
+&NAM_TURBn  CTURBLEN = "BL89", CTURBDIM = "1DIM" /
+&NAM_CH_MNHCn LUSECHEM = F /
+&NAM_CONF  CCONF = "START", NVERB=5, NMODEL = 2, 
+            CEXP = "16JA2", CSEG = "WENO5",NHALO=3 /
+&NAM_DYN  XSEGLEN = 21600., LCORIO = T, LNUMDIFU = F,
+           XALKTOP = 0.001, XALZBOT = 14500. /
+&NAM_NESTING NDAD(2) = 1, NDTRATIO(2) = 2, XWAY(2) = 2. /
+&NAM_BACKUP XBAK_TIME(1,1) = 10800., XBAK_TIME(1,2) = 21600. ,
+            XBAK_TIME(2,1) = 10800., XBAK_TIME(2,2) = 21600. /
+&NAM_DIAG_SURFn  /
+&NAM_ISBA CALBEDO = "DRY" /
+&NAM_ISBAn CSCOND  = "NP89", 
+           CC1DRY='DEF', CSOILFRZ='DEF', CDIFSFCOND='DEF', CSNOWRES='DEF' /
+&NAM_SSOn CROUGH="Z01D"/
+&NAM_DIAG_ISBAn /
+&NAM_SGH_ISBAn
+       CRUNOFF = "WSAT" /
+&NAM_SEAFLUXn CSEA_ALB="UNIF" /
+&NAM_BLANK /
+&NAM_NEBn LSUBG_COND = F   /

MY_RUN/INTEGRATION_CASES/HPC/16JAN/RUN/EXSEG2.nam_run2

+&NAM_LUNITn  CINIFILE = "16JAN_00_MNH2",
+             CINIFILEPGD = "16JAN98_9km.neste1" /
+&NAM_DYNn CPRESOPT = "ZRESI", NITR = 8, XRELAX = 1., 
+           LHORELAX_UVWTH = F, LHORELAX_RV = F, LHORELAX_RC= F, LHORELAX_RR= F, 
+           LHORELAX_RS= F, LHORELAX_RI= F, LHORELAX_RG= F, LHORELAX_TKE= F, 
+           NRIMX = 0, NRIMY = 0, LVE_RELAX = T, XT4DIFU = 1500. /
+&NAM_ADVn  CUVW_ADV_SCHEME = "WENO_K", CMET_ADV_SCHEME = "PPM_01",NWENO_ORDER=5,CTEMP_SCHEME='RK53' /
+&NAM_PARAMn  CCLOUD = "KESS", CTURB = "TKEL", CRAD = "ECMW", CDCONV = "KAFR" ,
+             CSCONV="KAFR" /
+&NAM_PARAM_RADn  XDTRAD = 1800., XDTRAD_CLONLY = 1800., LCLEAR_SKY = F,
+             NRAD_COLNBR = 400 /
+&NAM_PARAM_KAFRn XDTCONV = 300., NICE = 1, LREFRESH_ALL = T,
+              LDOWN = T /
+&NAM_LBCn  CLBCX = 2*"OPEN", CLBCY = 2*"OPEN", XCPHASE = 20. /
+&NAM_TURBn  XIMPL = 1., CTURBLEN = "BL89", CTURBDIM = "1DIM" ,
+            LTURB_FLX = F, LTURB_DIAG = F /
+&NAM_CH_MNHCn LUSECHEM = F /
+&NAM_DIAG_SURFn  /
+&NAM_ISBA CALBEDO = "DRY" /
+&NAM_ISBAn CSCOND  = "NP89", 
+           CC1DRY='DEF', CSOILFRZ='DEF', CDIFSFCOND='DEF', CSNOWRES='DEF' /
+&NAM_SSOn CROUGH="Z01D"/
+&NAM_DIAG_ISBAn /
+&NAM_SGH_ISBAn
+       CRUNOFF = "WSAT" /
+&NAM_SEAFLUXn CSEA_ALB="UNIF" /
+&NAM_BLANK /
+&NAM_NEBn LSUBG_COND = F   /

MY_RUN/INTEGRATION_CASES/HPC/16JAN/RUN/PRE_NEST_PGD1.nam

+&NAM_CONFIO  LCDF4=T, LLFIOUT=F, LLFIREAD=F /
+&NAM_PGD1 YPGD1 = '16JAN98_36km' /
+&NAM_PGD2 YPGD2 = '16JAN98_9km', IDAD = 1 /
+&NAM_NEST_PGD YNEST = 'e1' /
+" 

MY_RUN/INTEGRATION_CASES/HPC/16JAN/RUN/PRE_PGD1.nam_1

+&NAM_CONFIO  LCDF4=T, LLFIOUT=F, LLFIREAD=F /
+&NAM_PGDFILE CPGDFILE='16JAN98_36km' /
+&NAM_CONF_PROJ XLAT0=50., XLON0=-2.8, XRPK=0.58, XBETA=0. /
+&NAM_CONF_PROJ_GRID XLATCEN=48.,  XLONCEN=-1.,
+               NIMAX=80, NJMAX=60, XDX=36000., XDY=36000. /
+&NAM_COVER     YCOVER='ECOCLIMAP_v2.0', YCOVERFILETYPE='DIRECT' /
+&NAM_ZS        YZS='gtopo30',    YZSFILETYPE='DIRECT' /
+&NAM_ISBA      YCLAY='CLAY_HWSD_MOY', YCLAYFILETYPE='DIRECT' ,
+               YSAND='SAND_HWSD_MOY', YSANDFILETYPE='DIRECT' /

MY_RUN/INTEGRATION_CASES/HPC/16JAN/RUN/PRE_PGD1.nam_2

+&NAM_CONFIO  LCDF4=T, LLFIOUT=F, LLFIREAD=F /
+&NAM_PGDFILE CPGDFILE='16JAN98_9km' ,NHALO=20/
+&NAM_PGD_GRID YINIFILE='16JAN98_36km' YINIFILETYPE='MESONH' /
+&NAM_INIFILE_CONF_PROJ IXOR=40, IYOR=16, IXSIZE=18, IYSIZE=15, IDXRATIO=4, IDYRATIO=4 /
+&NAM_COVER     YCOVER='ECOCLIMAP_v2.0', YCOVERFILETYPE='DIRECT' /
+&NAM_ZS        YZS='gtopo30',    YZSFILETYPE='DIRECT' /
+&NAM_ISBA      YCLAY='CLAY_HWSD_MOY', YCLAYFILETYPE='DIRECT' ,
+               YSAND='SAND_HWSD_MOY', YSANDFILETYPE='DIRECT' /

MY_RUN/INTEGRATION_CASES/HPC/16JAN/RUN/PRE_REAL1.nam_1

+&NAM_CONFIO  LCDF4=T, LLFIOUT=F, LLFIREAD=F /
+&NAM_FILE_NAMES HATMFILE ='ecmwf.EI.19980115.12' , HATMFILETYPE='GRIBEX',
+                HPGDFILE ='16JAN98_36km.neste1' , CINIFILE='15JAN_12_MNH' /
+&NAM_VER_GRID NKMAX=50, YZGRID_TYPE='FUNCTN',
+       ZDZGRD=60., ZDZTOP=700., 
+       ZZMAX_STRGRD=2500., ZSTRGRD=9., ZSTRTOP=7. /

MY_RUN/INTEGRATION_CASES/HPC/16JAN/RUN/PRE_REAL1.nam_2

+&NAM_CONFIO  LCDF4=T, LLFIOUT=F, LLFIREAD=F /
+&NAM_FILE_NAMES HATMFILE ='ecmwf.EI.19980115.18' , HATMFILETYPE='GRIBEX',
+                HPGDFILE ='16JAN98_36km.neste1' , CINIFILE='15JAN_18_MNH' /
+&NAM_VER_GRID NKMAX=50, YZGRID_TYPE='FUNCTN',
+       ZDZGRD=60., ZDZTOP=700., 
+       ZZMAX_STRGRD=2500., ZSTRGRD=9., ZSTRTOP=7. /

MY_RUN/INTEGRATION_CASES/HPC/16JAN/RUN/PRE_REAL1.nam_3

+&NAM_CONFIO  LCDF4=T, LLFIOUT=F, LLFIREAD=F /
+&NAM_FILE_NAMES HATMFILE ='ecmwf.EI.19980116.00' , HATMFILETYPE='GRIBEX',
+                HPGDFILE ='16JAN98_36km.neste1' , CINIFILE='16JAN_00_MNH' /
+&NAM_VER_GRID NKMAX=50, YZGRID_TYPE='FUNCTN',
+       ZDZGRD=60., ZDZTOP=700., 
+       ZZMAX_STRGRD=2500., ZSTRGRD=9., ZSTRTOP=7. /