.gitignore

-bin_tools/cdf2cdf
-bin_tools/cdf2lfi
-bin_tools/lfi2cdf
 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
@@ -39,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

LIBTOOLS/lib/COMPRESS/src/compress.f90

-!-----------------------------------------------------------------
-!--------------- special set of characters for RCS information
-!-----------------------------------------------------------------
-! $Source$ $Revision$ $Date$
+!MNH_LIC Copyright 1994-2019 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.
 !-----------------------------------------------------------------
 SUBROUTINE COMPRESS_FIELD(XTAB,KX,KY,KNBTOT,KNBUSE)
 USE MODD_COMPPAR
 USE MODE_SEARCHGRP
 
-#ifdef NAGf95
 USE,INTRINSIC :: IEEE_ARITHMETIC
-#endif
 
 IMPLICIT NONE 
 
@@ -38,17 +36,21 @@ INTEGER :: IEXTCOD
 CHARACTER(LEN=8),PARAMETER :: KEYWORD='COMPRESS'
 REAL,DIMENSION(KNBTOT) :: XWORKTAB
 LOGICAL :: LUPREAL,LNAN
-#ifndef NAGf95
-LOGICAL, EXTERNAL :: IEEE_IS_NAN
-#endif
+logical :: gnansupport
 
 ILEVNBELT = KX*KY
 LUPREAL = .FALSE.
 LNAN    = .FALSE.
 
+if ( IEEE_SUPPORT_NAN( xtab(1)) ) then
+  gnansupport=.true.
+else
+  gnansupport=.false.
+end if
+
 ! Check for NAN and change Upper and Lower bound according to 32bits real limits.
 DO JI=1,KNBTOT
-  IF (IEEE_IS_NAN(XTAB(JI))) THEN 
+  IF ( gnansupport .and. IEEE_IS_NAN(XTAB(JI)) ) THEN
     XTAB(JI)=0.
     LNAN = .TRUE.
   ELSE IF (ABS(XTAB(JI)) > HUGE(1.0_4)) THEN

LIBTOOLS/lib/COMPRESS/src/ieee_is_nan.c

-#include <math.h>
-
-#ifdef NO_UNDERSCORE
-# define IEEE_IS_NAN ieee_is_nan
-#else
-# define IEEE_IS_NAN ieee_is_nan_
-#endif
-
-int IEEE_IS_NAN(double *x){
-    return isnan(*x);
-}

LIBTOOLS/tools/lfi2cdf/scripts/lfi2cdfregex.sh

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

LIBTOOLS/tools/lfi2cdf/src/lfi2cdf.f90

-!MNH_LIC Copyright 1994-2018 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 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_IO_ll,         ONLY: NIO_VERB, NGEN_VERB, LVERB_OUTLST, LVERB_STDOUT, NVERB_DEBUG
+  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_IO_ll,         ONLY: INITIO_ll, SET_CONFIO_ll
-  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
+
   USE MODI_VERSION
 
   USE MODN_CONFIO, ONLY: LCDF4, LLFIOUT, LLFIREAD
@@ -31,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
@@ -43,19 +51,28 @@ program LFI2CDF
 
   CPROGRAM = 'LFICDF'
 
-  CALL INITIO_ll()
+  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.
   LVERB_STDOUT = .TRUE.
 
   call read_commandline(options,hinfile,houtfile,runmode)
 
+  if (options(OPTVERBOSE)%set) then
+    NIO_VERB  = NVERB_DEBUG
+    NGEN_VERB = NVERB_DEBUG
+  end if
+
   IF (options(OPTMERGE)%set) THEN
     NB_PROCIO_R = options(OPTMERGE)%ivalue
   ELSE
@@ -66,25 +83,25 @@ program LFI2CDF
      LCDF4    = .TRUE.
      LLFIOUT  = .FALSE.
      LLFIREAD = .TRUE.
-     CALL SET_CONFIO_ll()
+     CALL IO_Config_set()
   ELSE IF (runmode == MODECDF2CDF) THEN
      LCDF4    = .TRUE.
      LLFIOUT  = .FALSE.
      LLFIREAD = .FALSE.
-     CALL SET_CONFIO_ll()
+     CALL IO_Config_set()
   ELSE
      LCDF4    = .TRUE.
      LLFIOUT  = .TRUE.
      LLFIREAD = .FALSE.
-     CALL SET_CONFIO_ll()
+     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)
@@ -96,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
@@ -127,14 +150,14 @@ program LFI2CDF
      ! Conversion LFI -> NetCDF
      IF (options(OPTSPLIT)%set) call open_split_ncfiles_out(outfiles,nfiles_out,houtfile,nbvar_tbw,options)
      CALL parse_infiles(infiles,outfiles,nfiles_out,nbvar_infile,nbvar_tbr,nbvar_calc,nbvar_tbw,tzreclist,options,runmode)
-     CALL def_ncdf(outfiles,nfiles_out)
+     CALL def_ncdf(infiles,outfiles,nfiles_out)
      CALL fill_files(infiles,outfiles,tzreclist,nbvar,options)
 
   ELSE IF (runmode == MODECDF2CDF) THEN
      ! Conversion netCDF -> netCDF
      IF (options(OPTSPLIT)%set) call open_split_ncfiles_out(outfiles,nfiles_out,houtfile,nbvar_tbw,options)
      CALL parse_infiles(infiles,outfiles,nfiles_out,nbvar_infile,nbvar_tbr,nbvar_calc,nbvar_tbw,tzreclist,options,runmode)
-     CALL def_ncdf(outfiles,nfiles_out)
+     CALL def_ncdf(infiles,outfiles,nfiles_out)
      CALL fill_files(infiles,outfiles,tzreclist,nbvar,options)
 
   ELSE
@@ -143,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 = 8
+  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 :: OPTSPLIT    = 7, OPTVAR    = 8
+  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.
@@ -58,8 +63,6 @@ subroutine read_commandline(options,hinfile,houtfile,runmode)
       runmode = MODELFI2CDF
     case default
       runmode = MODEUNDEF
-      print *,'Error: program started with unknown command: ',command
-      call help()
   end select
   deallocate(command,fullcommand)
 
@@ -137,6 +140,11 @@ subroutine init_options(options)
   options(OPTREDUCE)%short_name   = 'r'
   options(OPTREDUCE)%has_argument = .false.
 
+  options(OPTMODE)%long_name    = "runmode"
+  options(OPTMODE)%short_name   = 'R'
+  options(OPTMODE)%has_argument = .true.
+  options(OPTMODE)%type         = TYPECHAR
+
   options(OPTSPLIT)%long_name    = "split"
   options(OPTSPLIT)%short_name   = 's'
   options(OPTSPLIT)%has_argument = .false.
@@ -146,6 +154,20 @@ subroutine init_options(options)
   options(OPTVAR)%has_argument = .true.
   options(OPTVAR)%type         = TYPECHAR
 
+  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)
@@ -228,7 +250,7 @@ subroutine check_options(options,infile,runmode)
 
   type(option),dimension(:),intent(inout) :: options
   character(len=:),allocatable,intent(in) :: infile
-  integer,intent(in)                      :: runmode
+  integer,intent(inout)                   :: runmode
 
   integer :: idx1, idx2
 
@@ -237,6 +259,26 @@ subroutine check_options(options,infile,runmode)
     call help()
   end if
 
+  !Check runmode
+  if (options(OPTMODE)%set) then
+    select case (options(OPTMODE)%cvalue)
+      case ('cdf2cdf')
+        runmode = MODECDF2CDF
+      case ('lfi2cdf')
+        runmode = MODELFI2CDF
+      case ('cdf2lfi')
+        runmode = MODECDF2LFI
+      case default
+        print *,'Error: invalid runmode option'
+        call help()
+     end select
+  else
+    if(runmode==MODEUNDEF) then
+      print *,'Error: program started with unknown command'
+      call help()
+    end if
+  end if
+
   !Check compression level
   if (options(OPTCOMPRESS)%set) then
     if (options(OPTCOMPRESS)%ivalue < 1 .OR. options(OPTCOMPRESS)%ivalue > 9 ) then
@@ -251,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.
@@ -302,16 +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] [-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] [-f --fallback-file fallback-file]"
   print *,"                [-c --compress compression_level] input-file.nc"
-  print *,"        cdf2lfi [-o --output output-file.lfi] 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"
@@ -322,12 +363,16 @@ subroutine help()
   print *,"     Name of file for the output"
   print *,"  --reduce-precision, -r"
   print *,"     Reduce the precision of the floating point variables to single precision (cdf2cdf and lfi2cdf only)"
+  print *,"  --runmode, -R"
+  print *,"     Force runmode (lfi2cdf, cdf2cdf or cdf2lfi)"
   print *,"  --split, -s"
   print *,"     Split variables specified with the -v option (one per file) (cdf2cdf and lfi2cdf only)"
   print *,"  --var, -v var1[,...]"
   print *,"     List of the variable to write in the output file. Variables names have to be separated by commas (,)."
   print *,"     A variable can be computed from the sum of existing variables (format: new_var=var1+var2[+...])"
   print *,"     (cdf2cdf and lfi2cdf only)"
+  print *,"  --verbose, -V"
+  print *,"     Be verbose (for debugging purpose)"
   print *,""
   stop
 

LIBTOOLS/tools/lfiz/src/lfiz.f90

@@ -14,25 +14,27 @@ INTEGER :: arglen
 INTEGER :: inarg
 CHARACTER(LEN=50) :: yexe
 
+LOGICAL(KIND=LFI_INT),PARAMETER :: GTRUE  = .TRUE.
+LOGICAL(KIND=LFI_INT),PARAMETER :: GFALSE = .FALSE.
 
 INTEGER, PARAMETER :: FM_FIELD_SIZE = 16
-INTEGER, PARAMETER :: ISRCLU  = 11
-INTEGER, PARAMETER :: IDESTLU = 12
+INTEGER(KIND=LFI_INT), PARAMETER :: ISRCLU  = 11
+INTEGER(KIND=LFI_INT), PARAMETER :: IDESTLU = 12
 INTEGER :: JPHEXT
-INTEGER :: iverb
-INTEGER :: inap ! nb d'articles prevus (utile a la creation)
-INTEGER :: inaf ! nb d'articles presents dans un fichier existant
-INTEGER :: inafdest
+INTEGER(KIND=LFI_INT) :: iverb
+INTEGER(KIND=LFI_INT) :: inap ! nb d'articles prevus (utile a la creation)
+INTEGER(KIND=LFI_INT) :: inaf ! nb d'articles presents dans un fichier existant
+INTEGER(KIND=LFI_INT) :: inafdest
 
 CHARACTER(LEN=128) :: filename,DESTFNAME
 INTEGER :: JI,JJ
-INTEGER :: IRESP
+INTEGER(KIND=LFI_INT) :: IRESP
 CHARACTER(LEN=FM_FIELD_SIZE),DIMENSION(:),ALLOCATABLE :: yrecfm
-INTEGER,                     DIMENSION(:),ALLOCATABLE :: ileng
+INTEGER(KIND=LFI_INT),       DIMENSION(:),ALLOCATABLE :: ileng
 INTEGER(KIND=8),             DIMENSION(:),ALLOCATABLE :: iwork
 
-INTEGER :: ilengs
-INTEGER :: ipos
+INTEGER(KIND=LFI_INT) :: ilengs
+INTEGER(KIND=LFI_INT) :: ipos
 INTEGER :: sizemax
 
 INTEGER            :: IGRID
@@ -46,7 +48,9 @@ INTEGER :: LFICOMP
 INTEGER :: NEWSIZE
 INTEGER :: searchndx
 INTEGER :: INDDATIM
-INARG = IARGC()
+
+!OLD: INARG = IARGC()
+INARG = COMMAND_ARGUMENT_COUNT()
 
 #if defined(F90HP)
 #define HPINCR 1
@@ -54,6 +58,9 @@ INARG = IARGC()
 #define HPINCR 0
 #endif
 
+  CALL GET_COMMAND_ARGUMENT(0,yexe)
+#if 0
+!OLD:
 #if defined(FUJI) || defined(NAGf95) || defined(NEC) || defined(HP) || defined(pgf) || defined(G95) || defined(GFORTRAN)
   CALL GETARG(0+HPINCR,yexe)
   IF (LEN_TRIM(yexe) == 0) THEN
@@ -63,12 +70,17 @@ INARG = IARGC()
 #else
   CALL PXFGETARG(0,yexe,arglen,iresp)
 #endif
+#endif
 !  PRINT *,yexe, ' avec ',INARG,' arguments.'
   IF (INARG == 1) THEN 
+     CALL GET_COMMAND_ARGUMENT(1,filename)
+#if 0
+!OLD:
 #if defined(FUJI) || defined(NAGf95) || defined(NEC) || defined(HP) || defined(pgf) || defined(G95)|| defined(GFORTRAN)
      CALL GETARG(1+HPINCR,filename)
 #else
      CALL PXFGETARG(1,filename,arglen,iresp)
+#endif
 #endif
   ELSE 
      PRINT *,'Usage : ', TRIM(yexe), ' [fichier lfi]'
@@ -91,8 +103,8 @@ IDIMY = 0
 IDIMZ = 0
 GUSEDIM = .FALSE.
 
-CALL LFIOUV(IRESP,ISRCLU,.TRUE.,filename,'OLD',.FALSE.&
-            & ,.FALSE.,iverb,inap,inaf)
+CALL LFIOUV(IRESP,ISRCLU,GTRUE,filename,'OLD',GFALSE&
+            & ,GFALSE,iverb,inap,inaf)
 
 CALL FMREADLFIN1(ISRCLU,'LFI_COMPRESSED',LFICOMP,iresp)
 IF (iresp == 0) THEN
@@ -139,8 +151,8 @@ END IF
 
 
 PRINT *,'compressed file : ',DESTFNAME
-CALL LFIOUV(IRESP,IDESTLU,.TRUE.,DESTFNAME,'NEW'&
-     & ,.FALSE.,.FALSE.,iverb,inaf+1,inafdest)
+CALL LFIOUV(IRESP,IDESTLU,GTRUE,DESTFNAME,'NEW'&
+     & ,GFALSE,GFALSE,iverb,inaf+1,inafdest)
 
 CALL LFIPOS(IRESP,ISRCLU)
 ALLOCATE(yrecfm(inaf))
@@ -148,7 +160,7 @@ ALLOCATE(ileng(inaf))
 yrecfm(:) = ''
 sizemax=0
 DO ji=1,inaf
-  CALL LFICAS(IRESP,ISRCLU,yrecfm(ji),ileng(ji),ipos,.TRUE.)
+  CALL LFICAS(IRESP,ISRCLU,yrecfm(ji),ileng(ji),ipos,GTRUE)
   IF (ileng(ji) > sizemax) sizemax=ileng(ji)
 END DO
 PRINT *,' Nombre total d''articles dans fichier source :', inaf
@@ -218,13 +230,13 @@ CALL LFIFER(IRESP,IDESTLU,'KEEP')
 CONTAINS 
 
 SUBROUTINE FMREADLFIN1(klu,hrecfm,kval,kresp)
-INTEGER, INTENT(IN)         :: klu ! logical fortran unit au lfi file
-CHARACTER(LEN=*),INTENT(IN) :: hrecfm ! article name to be read
-INTEGER, INTENT(OUT)        :: kval ! integer value for hrecfm article
-INTEGER, INTENT(OUT)        :: kresp! return code null if OK
+INTEGER(KIND=LFI_INT), INTENT(IN)  :: klu ! logical fortran unit au lfi file
+CHARACTER(LEN=*),      INTENT(IN)  :: hrecfm ! article name to be read
+INTEGER,               INTENT(OUT) :: kval ! integer value for hrecfm article
+INTEGER(KIND=LFI_INT), INTENT(OUT) :: kresp! return code null if OK
 !
-INTEGER(KIND=8),DIMENSION(:),ALLOCATABLE::iwork
-INTEGER :: iresp,ilenga,iposex,icomlen
+INTEGER(KIND=8),DIMENSION(:),ALLOCATABLE :: iwork
+INTEGER(KIND=LFI_INT)                    :: iresp,ilenga,iposex,icomlen
 !
 CALL LFINFO(iresp,klu,hrecfm,ilenga,iposex)
 IF (iresp /=0 .OR. ilenga == 0) THEN

LIBTOOLS/tools/lfiz/src/unlfiz.f90

@@ -14,24 +14,26 @@ INTEGER :: arglen
 INTEGER :: inarg
 CHARACTER(LEN=50) :: yexe
 
+LOGICAL(KIND=LFI_INT),PARAMETER :: GTRUE  = .TRUE.
+LOGICAL(KIND=LFI_INT),PARAMETER :: GFALSE = .FALSE.
 
 INTEGER, PARAMETER :: FM_FIELD_SIZE = 16
-INTEGER, PARAMETER :: ISRCLU  = 11
-INTEGER, PARAMETER :: IDESTLU = 12
-INTEGER :: iverb
-INTEGER :: inap ! nb d'articles prevus (utile a la creation)
-INTEGER :: inaf ! nb d'articles presents dans un fichier existant
-INTEGER :: inafdest
+INTEGER(KIND=LFI_INT), PARAMETER :: ISRCLU  = 11
+INTEGER(KIND=LFI_INT), PARAMETER :: IDESTLU = 12
+INTEGER(KIND=LFI_INT) :: iverb
+INTEGER(KIND=LFI_INT) :: inap ! nb d'articles prevus (utile a la creation)
+INTEGER(KIND=LFI_INT) :: inaf ! nb d'articles presents dans un fichier existant
+INTEGER(KIND=LFI_INT) :: inafdest
 
 CHARACTER(LEN=128) :: filename,DESTFNAME
 INTEGER :: JI,JJ
-INTEGER :: IRESP
+INTEGER(KIND=LFI_INT) :: IRESP
 CHARACTER(LEN=FM_FIELD_SIZE),DIMENSION(:),ALLOCATABLE :: yrecfm
-INTEGER,                     DIMENSION(:),ALLOCATABLE :: ileng
+INTEGER(KIND=LFI_INT),       DIMENSION(:),ALLOCATABLE :: ileng
 INTEGER(KIND=8),             DIMENSION(:),ALLOCATABLE :: iwork,iworknew
 
-INTEGER :: ilengs
-INTEGER :: ipos
+INTEGER(KIND=LFI_INT) :: ilengs
+INTEGER(KIND=LFI_INT) :: ipos
 INTEGER :: sizemax
 
 INTEGER            :: ICOMLEN
@@ -43,9 +45,10 @@ INTEGER :: CPT
 INTEGER :: LFICOMP
 INTEGER :: searchndx
 INTEGER :: ITYPCOD
-INTEGER :: ITOTAL,ITOTALMAX
+INTEGER(KIND=LFI_INT) :: ITOTAL,ITOTALMAX
 
-INARG = IARGC()
+!OLD: INARG = IARGC()
+INARG = COMMAND_ARGUMENT_COUNT()
 
 #if defined(F90HP)
 #define HPINCR 1
@@ -53,6 +56,9 @@ INARG = IARGC()
 #define HPINCR 0
 #endif
 
+  CALL GET_COMMAND_ARGUMENT(0,yexe)
+#if 0
+!OLD:
 #if defined(FUJI) || defined(NAGf95) || defined(NEC) || defined(HP) || defined(pgf) || defined(G95) || defined(GFORTRAN)
   CALL GETARG(0+HPINCR,yexe)
   IF (LEN_TRIM(yexe) == 0) THEN
@@ -62,12 +68,17 @@ INARG = IARGC()
 #else
   CALL PXFGETARG(0,yexe,arglen,iresp)
 #endif
+#endif
 !  PRINT *,yexe, ' avec ',INARG,' arguments.'
   IF (INARG == 1) THEN 
+     CALL GET_COMMAND_ARGUMENT(1,filename)
+#if 0
+!OLD:
 #if defined(FUJI) || defined(NAGf95) || defined(NEC) || defined(HP) || defined(pgf) || defined(G95) || defined(GFORTRAN)
      CALL GETARG(1+HPINCR,filename)
 #else
      CALL PXFGETARG(1,filename,arglen,iresp)
+#endif
 #endif
   ELSE 
      PRINT *,'Usage : ', TRIM(yexe), ' [fichier lfi]'
@@ -93,8 +104,8 @@ IDIMY = 0
 IDIMZ = 0
 GUSEDIM = .FALSE.
 
-CALL LFIOUV(IRESP,ISRCLU,.TRUE.,filename,'OLD',.FALSE.&
-            & ,.FALSE.,iverb,inap,inaf)
+CALL LFIOUV(IRESP,ISRCLU,GTRUE,filename,'OLD',GFALSE&
+            & ,GFALSE,iverb,inap,inaf)
 
 CALL FMREADLFIN1(ISRCLU,'LFI_COMPRESSED',LFICOMP,iresp)
 IF (iresp /= 0 .OR. LFICOMP /= 1) THEN
@@ -104,8 +115,8 @@ IF (iresp /= 0 .OR. LFICOMP /= 1) THEN
 END IF  
 
 PRINT *,'Uncompressed (but 32 bits REAL precision) file : ',DESTFNAME
-CALL LFIOUV(IRESP,IDESTLU,.TRUE.,DESTFNAME,'NEW'&
-     & ,.FALSE.,.FALSE.,iverb,inaf,inafdest)
+CALL LFIOUV(IRESP,IDESTLU,GTRUE,DESTFNAME,'NEW'&
+     & ,GFALSE,GFALSE,iverb,inaf,inafdest)
 
 CALL LFIPOS(IRESP,ISRCLU)
 ALLOCATE(yrecfm(inaf))
@@ -113,7 +124,7 @@ ALLOCATE(ileng(inaf))
 yrecfm(:) = ''
 sizemax=0
 DO ji=1,inaf
-  CALL LFICAS(IRESP,ISRCLU,yrecfm(ji),ileng(ji),ipos,.TRUE.)
+  CALL LFICAS(IRESP,ISRCLU,yrecfm(ji),ileng(ji),ipos,GTRUE)
   IF (ileng(ji) > sizemax) sizemax=ileng(ji)
 END DO
 PRINT *,' Nombre total d''articles dans fichier source :', inaf
@@ -173,13 +184,13 @@ CALL LFIFER(IRESP,IDESTLU,'KEEP')
 CONTAINS 
 
 SUBROUTINE FMREADLFIN1(klu,hrecfm,kval,kresp)
-INTEGER, INTENT(IN)         :: klu ! logical fortran unit au lfi file
-CHARACTER(LEN=*),INTENT(IN) :: hrecfm ! article name to be read
-INTEGER, INTENT(OUT)        :: kval ! integer value for hrecfm article
-INTEGER, INTENT(OUT)        :: kresp! return code null if OK
+INTEGER(KIND=LFI_INT), INTENT(IN)  :: klu ! logical fortran unit au lfi file
+CHARACTER(LEN=*),      INTENT(IN)  :: hrecfm ! article name to be read
+INTEGER,               INTENT(OUT) :: kval ! integer value for hrecfm article
+INTEGER(KIND=LFI_INT), INTENT(OUT) :: kresp! return code null if OK
 !
-INTEGER(KIND=8),DIMENSION(:),ALLOCATABLE::iwork
-INTEGER :: iresp,ilenga,iposex,icomlen
+INTEGER(KIND=8),DIMENSION(:),ALLOCATABLE :: iwork
+INTEGER(KIND=LFI_INT)                    :: iresp,ilenga,iposex,icomlen
 !
 CALL LFINFO(iresp,klu,hrecfm,ilenga,iposex)
 IF (iresp /=0 .OR. ilenga == 0) THEN

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' /