diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/001_pgd1/PRE_PGD1.nam b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/001_pgd1/PRE_PGD1.nam
new file mode 100644
index 0000000000000000000000000000000000000000..3ae48e8b0f44fc5ea746497eca9a0a0b2277c131
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/001_pgd1/PRE_PGD1.nam
@@ -0,0 +1,13 @@
+&NAM_CONFZ
+ NB_PROCIO_W=1 ,
+ !MPI_BUFFER_SIZE=40
+/
+&NAM_PGDFILE CPGDFILE="PGD00128" /
+&NAM_PGD_SCHEMES CNATURE='ISBA', CSEA='SEAFLX', CTOWN='NONE', CWATER='NONE' /
+&NAM_CONF_PROJ XLAT0=-11.5, XLON0=130.7, XRPK=0., XBETA=0. /
+&NAM_CONF_PROJ_GRID XLATCEN=-11.5, XLONCEN=130.7, NIMAX=128, NJMAX=128,
+ XDX=1600.00000, XDY=1600.00000 /
+&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' /
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/001_pgd1/get_pgd_files b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/001_pgd1/get_pgd_files
new file mode 100755
index 0000000000000000000000000000000000000000..56726933a0f7608e0c4918f42e80d5aea76b669d
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/001_pgd1/get_pgd_files
@@ -0,0 +1,30 @@
+#
+# Modif
+# J.Escobar 11/04/2014 get PGD files from 'dir_open' directory ( without psswd )
+# J.Escobar 25/04/2013 get LICENCE files
+#
+PREP_PGD_FILES=${PREP_PGD_FILES:-"$HOME/PREP_PGD_FILES_WWW"} ; export PREP_PGD_FILES
+echo "!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!"
+echo "!!!! WARNING !!!!"
+echo "!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!"
+echo
+echo you need 3GO to download this files in
+echo
+echo PREP_PGD_FILES=$PREP_PGD_FILES
+echo
+echo if OK press ENTER else CTRL-C
+read RIEN
+set -x
+mkdir -p $PREP_PGD_FILES
+cd $PREP_PGD_FILES
+PGD_URL="http://mesonh.aero.obs-mip.fr/mesonh/dir_open/dir_PGDFILES"
+WGET="wget"
+export PGD_URL
+for file in LICENSE_ECOCLIMAP.txt LICENSE_soil_data.txt \
+ gtopo30.dir gtopo30.hdr \
+ SAND_HWSD_MOY.hdr SAND_HWSD_MOY.dir CLAY_HWSD_MOY.hdr CLAY_HWSD_MOY.dir \
+ ECOCLIMAP_v2.0.hdr ECOCLIMAP_v2.0.dir
+do
+[ -f $file ] || ( ${WGET} -c -nd $PGD_URL/$file.gz ; gunzip $file.gz ; )
+done
+
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/001_pgd1/run_prep_pgd_xyz b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/001_pgd1/run_prep_pgd_xyz
new file mode 100755
index 0000000000000000000000000000000000000000..5c34e67e3f8ec036121d636d3d9c8cc5fadf197e
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/001_pgd1/run_prep_pgd_xyz
@@ -0,0 +1,35 @@
+#!/bin/bash
+
+export MPIRUN="Mpirun -np 4"
+
+
+PREP_PGD_FILES=${PREP_PGD_FILES:-"$HOME/PREP_PGD_FILES_WWW"} ; export PREP_PGD_FILES
+#
+if [ ! -d $PREP_PGD_FILES ]
+then
+cat << EOF
+
+Your directory PREP_PGD_FILES=$PREP_PGD_FILES
+
+containing the files gtopo30*, ECOCLIMAP_v2.0*
+doesn't exist ( or was not found !!! )
+use the script 'get_pgd_files' to download
+this files from the MESONH WEB PAGES !!!
+( or change the variable PREP_PGD_FILES ... )
+
+After fixing it , run this script again !!!
+
+EOF
+exit 1
+else
+set -x
+#
+rm -f OUTPUT_LISTING0 pipe_name
+rm -f gtopo30.??? sand_fao.??? clay_fao.???
+rm -f SAND_HWSD_MOY.??? ECOCLIMAP_v2.0.??? ecoclimats_v2.???
+rm -f PGD00128.*
+#
+ln -sf $PREP_PGD_FILES/*.dir $PREP_PGD_FILES/*.hdr .
+#
+time ${MPIRUN} PREP_PGD${XYZ}
+fi
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/002_prep_ideal_case/PRE_IDEA1.nam b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/002_prep_ideal_case/PRE_IDEA1.nam
new file mode 100644
index 0000000000000000000000000000000000000000..3baef58caeae016c444e5f7c447f821726afb861
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/002_prep_ideal_case/PRE_IDEA1.nam
@@ -0,0 +1,226 @@
+&NAM_CONFZ
+ !NB_PROCIO_R=1 ,
+ !NB_PROCIO_W=1 ,
+ !MPI_BUFFER_SIZE=40
+/
+
+&NAM_REAL_PGD
+ CPGD_FILE ="PGD00128" ,
+ LREAD_ZS =.TRUE., LREAD_GROUND_PARAM =.TRUE.
+/
+&NAM_DIMn_PRE NIMAX=40, NJMAX=32 /
+&NAM_CONF_PRE LCARTESIAN=.FALSE.
+ CIDEAL='RSOU' CZS='FLAT'
+ NVERB=0 /
+&NAM_CONFn LUSERV= T /
+&NAM_LUNITn CINIFILE = "DA0128" , CINIFILEPGD = "PGD00128" /
+&NAM_DYNn_PRE
+ CPRESOPT= 'ZRESI' ,
+ NITR=4 XRELAX=1. /
+&NAM_LBCn_PRE CLBCX= 2*'OPEN' CLBCY= 2*'OPEN' /
+&NAM_VPROF_PRE CTYPELOC='IJGRID' NILOC=1 NJLOC=1
+ CFUNU='ZZZ' CFUNV='ZZZ'
+ LGEOSBAL=.FALSE. /
+&NAM_VER_GRID NKMAX=126, YZGRID_TYPE='FUNCTN',
+ ZDZGRD=40., ZDZTOP=210., ZZMAX_STRGRD=2500.,
+ ZSTRGRD=7., ZSTRTOP=7. /
+&NAM_GRn_PRE
+ CSURF='EXTE'
+ ! CSURF='NONE'
+ /
+&NAM_PREP_ISBA XTG_SURF= 311., XTG_ROOT= 303., XTG_DEEP= 302.,
+ XHUG_SURF= 0.16, XHUG_ROOT= 0.16, XHUG_DEEP= 0.16 /
+&NAM_PREP_SEAFLUX XSST_UNIF= 304. /
+ RSOU
+ 2005 11 30 0
+ 'PUVTHDMR'
+ 0.0000000E+00
+ 100300.0
+ 303.3000
+ 1.9630000E-02
+ 92
+ 100000.0 3.637668 -3.637668
+ 99900.00 2.572220 -4.455216
+ 99100.00 -4.612377 0.4035313
+ 96600.00 5.124864 -0.4483674
+ 93000.00 5.836996 -2.009839
+ 92800.00 5.801031 -2.111402
+ 92500.00 5.466062 -1.464627
+ 85000.00 0.1345103 -1.537459
+ 81400.00 -0.9020693 -1.849517
+ 71800.00 -2.797467 1.304481
+ 70000.00 -4.472232 1.198331
+ 60500.00 -3.289776 -1.464702
+ 59900.00 -3.208611 -1.634869
+ 56900.00 -2.329535 -2.025034
+ 55600.00 -2.105100 -2.257443
+ 50100.00 -0.5325915 -1.987659
+ 50000.00 -0.2679967 -1.519885
+ 49500.00 0.2507798 -2.042438
+ 46700.00 5.466062 1.464627
+ 45300.00 5.594935 2.608961
+ 44800.00 5.450725 2.898202
+ 43700.00 4.994327 3.628591
+ 42200.00 4.561049 4.891127
+ 41900.00 4.387573 5.047325
+ 41800.00 4.298817 5.123130
+ 40000.00 6.305391 5.290851
+ 39000.00 5.970908 4.027428
+ 38600.00 5.671552 3.543979
+ 38300.00 5.346259 3.086664
+ 35500.00 2.976104 3.546782
+ 32700.00 4.455215 2.572220
+ 31400.00 3.712562 4.270814
+ 30700.00 3.086664 5.346259
+ 30000.00 3.540875 5.056894
+ 28800.00 3.786532 4.205370
+ 26700.00 3.329551 2.419061
+ 25000.00 2.797467 1.304481
+ 24400.00 2.819808 1.255459
+ 23800.00 2.841289 1.206056
+ 21400.00 3.424857 1.112804
+ 20000.00 3.478403 0.9320353
+ 19300.00 3.383935 1.231651
+ 17400.00 3.411947 2.301387
+ 15000.00 2.950729 4.214078
+ 12700.00 2.167201 2.875973
+ 10200.00 1.322713 1.576348
+ 10000.00 0.4466614 2.533142
+ 9120.000 -3.306782 3.940869
+ 8700.000 -7.950636 2.130367
+ 8200.000 -8.712269 0.7622259
+ 8000.000 -9.224754 -0.8070621
+ 7600.000 -4.009694 -2.314998
+ 7300.000 -5.066284 -0.8933228
+ 7000.000 -4.196201 -1.956721
+ 6700.000 -1.013257 0.1786647
+ 6400.000 -2.107039 1.475365
+ 6340.000 -2.558960 1.726041
+ 5480.000 -13.17234 2.322640
+ 5300.000 -15.37459 1.345105
+ 5000.000 -12.81216 -1.120920
+ 4900.000 -12.66571 -2.233307
+ 4860.000 -12.25463 -1.504677
+ 4700.000 -10.76221 0.9415731
+ 4500.000 -11.37801 7.966970
+ 3800.000 -25.11183 2.197004
+ 3600.000 -22.54940 -1.972818
+ 3200.000 -22.72465 -4.417223
+ 3000.000 -22.36116 -5.991656
+ 2490.000 -21.48829 -2.258509
+ 2300.000 -21.07936 -0.7361050
+ 2200.000 -20.57776 -1.7989648E-06
+ 2080.000 -19.00834 -0.9961855
+ 2000.000 -17.93702 -1.569287
+ 1600.000 -14.74927 -10.32755
+ 1500.000 -13.52109 -6.304989
+ 1430.000 -16.89510 -4.527030
+ 1400.000 -18.23862 -3.215962
+ 1300.000 -15.46941 -5.630407
+ 1200.000 -5.594935 2.608962
+ 1190.000 -5.671552 3.543979
+ 1100.000 1.786646 10.13257
+ 1050.000 -3.736842 7.333968
+ 1020.000 -5.435582 4.725079
+ 1000.000 -5.594935 2.608962
+ 950.0000 -4.347303 8.176087
+ 900.0000 2.143975 12.15908
+ 889.9999 1.031246 11.78719
+ 800.0000 -4.728969 4.728969
+ 760.0000 2.890712 7.154763
+ 710.0000 8.644947 -3.318484
+ 700.0000 7.487652 -6.282886
+ 600.0000 -4.214078 -2.950729
+ 91
+ 100000.0 303.1000 1.8050000E-02
+ 99900.00 303.1000 1.8040001E-02
+ 99100.00 303.1000 1.7969999E-02
+ 96600.00 303.1000 1.7729999E-02
+ 93000.00 303.0000 1.7370002E-02
+ 92800.00 303.1000 1.6980000E-02
+ 92500.00 303.2000 1.6400000E-02
+ 85000.00 307.1000 1.1190000E-02
+ 81400.00 309.6000 6.6700005E-03
+ 71800.00 312.3000 6.5600001E-03
+ 70000.00 312.9000 6.5400004E-03
+ 60500.00 316.2000 6.1800000E-03
+ 59900.00 318.1000 5.3400006E-03
+ 56900.00 318.7000 4.9899998E-03
+ 55600.00 322.2000 3.3100001E-03
+ 50100.00 325.5000 2.6000000E-03
+ 50000.00 325.5000 2.5800001E-03
+ 49500.00 325.7000 3.1700002E-03
+ 46700.00 328.6000 1.4500001E-03
+ 45300.00 330.1000 9.4000006E-04
+ 44800.00 330.1000 1.9600000E-03
+ 43700.00 331.0000 2.1500003E-03
+ 42200.00 332.2000 1.9400001E-03
+ 41900.00 332.9000 9.6000003E-04
+ 41800.00 333.0000 9.7000005E-04
+ 40000.00 334.2000 1.1600000E-03
+ 39000.00 334.6000 1.0300000E-03
+ 38600.00 335.6000 6.5000000E-04
+ 38300.00 335.8000 6.0000003E-04
+ 35500.00 338.2000 2.5000001E-04
+ 32700.00 340.8000 9.0000009E-05
+ 31400.00 342.0000 5.0000002E-05
+ 30700.00 342.1000 4.0000003E-05
+ 30000.00 342.3000 4.0000003E-05
+ 28800.00 342.3000 1.2000000E-04
+ 26700.00 344.8000 1.0000001E-05
+ 25000.00 345.4000 1.0000001E-05
+ 24400.00 346.3000 2.0000001E-05
+ 23800.00 346.4000 9.0000009E-05
+ 21400.00 348.1000 1.0000001E-05
+ 20000.00 348.8000 1.0000001E-05
+ 19300.00 348.9000 2.0000001E-05
+ 17400.00 351.5000 1.0000001E-05
+ 15000.00 354.0000 7.1596442E-06
+ 12700.00 355.7000 4.4376361E-06
+ 10200.00 364.6000 1.4789329E-06
+ 10000.00 365.4000 1.2422363E-06
+ 9120.000 371.6000 1.2919259E-06
+ 8700.000 379.9000 1.5496893E-06
+ 8200.000 390.6000 1.8819877E-06
+ 8000.000 395.1000 2.0217396E-06
+ 7600.000 404.7000 2.3198761E-06
+ 7300.000 412.3000 2.4844724E-06
+ 7000.000 420.4000 2.4844724E-06
+ 6700.000 431.1000 2.4844724E-06
+ 6400.000 442.4000 2.4844724E-06
+ 6340.000 444.8000 2.4844724E-06
+ 5480.000 467.8000 2.4844724E-06
+ 5300.000 471.0000 2.4844724E-06
+ 5000.000 476.5000 2.4844724E-06
+ 4900.000 477.6000 2.4844724E-06
+ 4860.000 478.0000 2.4844724E-06
+ 4700.000 485.1000 2.4844724E-06
+ 4500.000 494.6000 2.4844724E-06
+ 3800.000 532.9000 2.4844724E-06
+ 3600.000 545.6000 2.4844724E-06
+ 3200.000 574.4000 2.4844724E-06
+ 3000.000 585.1000 2.4844724E-06
+ 2490.000 624.0000 2.4844724E-06
+ 2300.000 648.9000 2.4844724E-06
+ 2200.000 655.4000 2.4844724E-06
+ 2080.000 663.6000 2.4844724E-06
+ 2000.000 673.5000 2.4844724E-06
+ 1600.000 726.4000 2.4844724E-06
+ 1500.000 742.5000 2.4844724E-06
+ 1430.000 750.7000 2.4844724E-06
+ 1400.000 759.0000 2.4844724E-06
+ 1300.000 788.6000 2.4844724E-06
+ 1200.000 821.7000 2.4844724E-06
+ 1190.000 825.2000 2.4844724E-06
+ 1100.000 833.9000 2.4844724E-06
+ 1050.000 839.0000 2.4844724E-06
+ 1020.000 845.3000 2.4844724E-06
+ 1000.000 851.6000 2.4844724E-06
+ 950.0000 877.0000 2.4844724E-06
+ 900.0000 885.6000 2.4844724E-06
+ 889.9999 887.3000 2.4844724E-06
+ 800.0000 926.6000 2.4844724E-06
+ 760.0000 946.0000 2.4844724E-06
+ 710.0000 958.0000 2.4844724E-06
+ 700.0000 963.1000 2.4844724E-06
+ 600.0000 20.30000 2.4844724E-06
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/002_prep_ideal_case/run_prep_ideal_xyz b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/002_prep_ideal_case/run_prep_ideal_xyz
new file mode 100755
index 0000000000000000000000000000000000000000..cc04c8357ad94625ece1e7492b5d68fdb8caad26
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/002_prep_ideal_case/run_prep_ideal_xyz
@@ -0,0 +1,20 @@
+#!/bin/bash
+
+NP=${NP:-16}
+NPS=${NPS:-4}
+
+export ACC_DEVICE_TYPE=HOST
+export MPIRUN="Mpirun -tag-output -report-bindings --mca mpi_cuda_support 0 -map-by ppr:${NPS}:socket -bind-to none -np ${NP} set_core_device_impair"
+
+#
+set -x
+set -e
+#
+rm -f OUTPUT_LISTING1 pipe_name file_for_xtransfer
+rm -f PGD00128.*
+rm -f DA0128.*
+#
+ln -sf ../001_pgd1/PGD00128.* .
+#
+time ${MPIRUN} PREP_IDEAL_CASE${XYZ}
+#
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/003_mesonh_step1/EXSEG1.nam b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/003_mesonh_step1/EXSEG1.nam
new file mode 100644
index 0000000000000000000000000000000000000000..cbfeca7d31195ad0c0dfbc8ba72acc279b1083f4
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/003_mesonh_step1/EXSEG1.nam
@@ -0,0 +1,53 @@
+&NAM_CONFZ
+ ! NB_PROCIO_R=1 ,
+ ! NB_PROCIO_W=1 ,
+ ! LMNH_MPI_BSEND = T ,
+ ! MPI_BUFFER_SIZE = 40
+/
+
+&NAM_LUNITn CINIFILE = "DA0128" , CINIFILEPGD = "PGD00128" /
+&NAM_CONFn LUSERV = T, LUSERC = T, LUSERR = T, LUSERI = T,
+ LUSERS = T, LUSERG = T, LUSERH = F, LUSECI = T /
+&NAM_DYNn XTSTEP = 15.00000 ,
+ CPRESOPT = "ZRESI", NITR = 12
+ !CPRESOPT = "ZSOLV", NITR = 12
+ XRELAX = 1., LHORELAX_UVWTH = T, LHORELAX_RV = T, LVE_RELAX = T,
+ NRIMX = 6, NRIMY = 6, XRIMKMAX = 0.0005, XT4DIFU = 4000 /
+&NAM_ADVn CUVW_ADV_SCHEME = "WENO_K", CMET_ADV_SCHEME = "PPM_01",
+ CSV_ADV_SCHEME = "PPM_01" /
+&NAM_PARAMn CTURB = "TKEL" ! "NONE" ! "TKEL",
+ CRAD = "ECMW",
+ CCLOUD = "ICE3" ! "NONE" ! "ICE3",
+ CDCONV = "NONE", CSCONV = "NONE" /
+&NAM_PARAM_RADn XDTRAD = 1800., XDTRAD_CLONLY = 1800., LCLEAR_SKY = F,
+ NRAD_COLNBR = 500 /
+&NAM_PARAM_MFSHALLn XIMPL_MF = 1, CMF_UPDRAFT = "EDKF", CMF_CLOUD = "NONE",
+ LMIXUV = T, LMF_FLX = F /
+&NAM_LBCn CLBCX = 2*"OPEN", CLBCY = 2*"OPEN", XCPHASE = 20. /
+&NAM_TURBn XIMPL = 1., CTURBLEN = "DEAR", CTURBDIM = "3DIM",
+ LTURB_FLX = F, LTURB_DIAG = F, CSUBG_AUCV = "NONE", LSIGMAS = F,
+ LSIG_CONV = F, LSUBG_COND = F /
+&NAM_CH_MNHCn /
+&NAM_CONF CCONF = "START", LFLAT = F, CEQNSYS = "DUR",
+ LLG=F, NMODEL = 1, NVERB = 0, CEXP = "D0128",
+ CSEG = "DAR01", CSPLIT = "BSPLITTING"
+ NHALO=1
+ /
+&NAM_DYN XSEGLEN = 14400.0000 ! 150. ! 14400.0000 , LCORIO = T,
+ LNUMDIFU = F, LNUMDIFTH = F,
+ XALKTOP = 0.001, XALZBOT = 22000. /
+&NAM_BLANK /
+&NAM_NESTING /
+&NAM_BACKUP
+ XBAK_TIME(1,1) = 3600.00000 ! 150. ! 3600.00000
+ XBAK_TIME(1,2) = 7200.00000
+ XBAK_TIME(1,3) = 10800.0000
+ XBAK_TIME(1,4) = 14400.0000
+ XBAK_TIME(1,5) = 18000.0000
+ XBAK_TIME(1,6) = 21600.0000
+ /
+&NAM_ISBAn /
+&NAM_SEAFLUXn /
+&NAM_DIAG_SURFn LSURF_BUDGET=T /
+&NAM_DIAG_ISBAn /
+&NAM_DIAG_SURF_ATMn /
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/003_mesonh_step1/parameters_mg.nam b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/003_mesonh_step1/parameters_mg.nam
new file mode 100644
index 0000000000000000000000000000000000000000..5d76f4dff351b9c8873c37072d3f1f0747877d70
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/003_mesonh_step1/parameters_mg.nam
@@ -0,0 +1,118 @@
+! *********************************************************************
+! *********************************************************************
+! ** **
+! ** Parameter file for geometric multigrid code **
+! ** **
+! *********************************************************************
+! *********************************************************************
+!
+! *********************************************************************
+! * General parameters
+! *********************************************************************
+¶meters_general
+ savefields = .F. ! Save fields to disk?
+/
+
+! *********************************************************************
+! * General solver parameters
+! *********************************************************************
+¶meters_solver
+ LUseO = F ,
+ LUseT = T ,
+ solvertype = 1, ! Solver type:
+ ! 1 : Richardson iteration
+ ! 2 : Conjugate gradient
+ resreduction = 1e-1 ! 1.0d-10 ! Required relative residual reduction
+ maxiter = 5 ! 50 ! Maximal number of iterations
+/
+
+! *********************************************************************
+! * Conjugate gradient parameters
+! *********************************************************************
+¶meters_conjugategradient
+ verbose = 10, ! Verbosity level
+ maxiter = 5, ! Maximal number of iterations
+ resreduction = 1.0e-1, ! Target residual reduction
+ n_prec = 1 ! Number of smoother applications in
+ ! preconditioner (N.B.: Using 0 is
+ ! inefficient, as the identity is used
+ ! for preconditioning, instead of using
+ ! unpreconditioned CG.)
+/
+
+! *********************************************************************
+! * Grid parameters
+! *********************************************************************
+¶meters_grid
+ n = 128 , ! 64, ! Number of horizontal grid cells
+ nz = 72, ! Number of vertical grid cells
+ L = 204800.0 ! 16000.0 metre * 256 pt grille ! 1.0, ! Size in horizontal direction
+ H = 19200.0 ! 150.0 metre * 128 niveaux ! 0.01, ! Size in vertical direction
+ vertbc = 2, ! Boundary conditions at top and bottom of the
+ ! atmosphere. 1 = DIRICHLET, 2 = NEUMANN
+ ! Note that Neumann boundary conditions only work
+ ! for coarsening in the horizontal only, as they are
+ ! not yet implemented in the prolongation operator.
+ graded = .F. ! Is the grid graded in the vertical direction?
+/
+
+! *********************************************************************
+! * Parallel communication parameters
+! *********************************************************************
+¶meters_communication
+ halo_size = 1 ! Size of halos (has to be 1 or 2)
+/
+
+! *********************************************************************
+! * Model parameters
+! *********************************************************************
+!
+! parameters of the Helmholtz operator
+!
+! -omega2*(d^2/dy^2 + d^2/dy^2 + lambda2*d^2/dz^2) u + delta u = RHS
+!
+¶meters_model
+ omega2 = 1.0,
+ lambda2 = 1.0 ! 100.0, ! Vertical coupling
+ delta = 0.0d0 ! Size of constant term
+/
+
+! *********************************************************************
+! * Smoother parameters
+! *********************************************************************
+!
+! parameters of the smoother
+!
+¶meters_smoother
+ smoother = 3, ! Smoother method
+ ! 3 = line SOR
+ ! 4 = line SSOR
+ ! 6 = line Jacobi
+ ordering = 2, ! Ordering of grid points (for smoother)
+
+ ! 1 = lexicographic
+ ! 2 = red-black ordering
+ !rho = 0.6666666666666666d0 ! Overrelaxation parameter
+ rho = 0.8d0 ! Overrelaxation parameter
+/
+
+! *********************************************************************
+! * Multigrid parameters
+! *********************************************************************
+¶meters_multigrid
+ verbose = 10, ! Verbosity level
+ n_lev = 8, ! Number of levels
+ lev_split = 3, ! First level where data is pulled together
+ n_presmooth = 1, ! Number of presmoothing steps
+ n_postsmooth = 1, ! Number of postsmoothing steps
+ n_coarsegridsmooth = 1, ! Number of smoothing steps on coarsest level
+ prolongation = 2 ! 2 best after modif !! ! Prologation method
+ ! 1 = constant interpolation
+ ! 2 = (tri-) linear interpolation
+ restriction = 1, ! Restriction method
+ ! 1 = cell average
+ ! 2 = Khalil
+ coarsegridsolver = 1 ! Solver on coarsest grid
+ ! 1 = use smoother
+ ! 2 = Conjugate gradient
+/
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/003_mesonh_step1/run_mesonh_step1_xyz b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/003_mesonh_step1/run_mesonh_step1_xyz
new file mode 100755
index 0000000000000000000000000000000000000000..d676734a931960542ae88ddc28a34cf7a11b509a
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/003_mesonh_step1/run_mesonh_step1_xyz
@@ -0,0 +1,23 @@
+
+
+#export PGI_ACC_POOL_ALLOC=0
+export PGI_ACC_SYNCHRONOUS=1
+
+NP=${NP:-16}
+NPS=${NPS:-4}
+
+export MPIRUN=${MPIRUN:-"Mpirun -tag-output -bind-to none -map-by ppr:${NPS}:socket -x PGI_ACC_POOL_ALLOC -x PGI_ACC_SYNCHRONOUS -np ${NP} set_core_device_impair "}
+
+
+#
+set -x
+set -e
+#
+rm -f DA0128.*
+rm -f PGD00128.*
+rm -f D0128.1.DAR01.* OUTPUT_LISTING* pipe_name
+#
+ln -sf ../001_pgd1/PGD00128.{des,nc} .
+ln -sf ../002_prep_ideal_case/DA0128.{des,nc} .
+#
+time ${MPIRUN} MESONH${XYZ}
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/004_mesonh_step2/EXSEG1.nam b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/004_mesonh_step2/EXSEG1.nam
new file mode 100644
index 0000000000000000000000000000000000000000..f2dcd82899ee955178d92a7bf90d484ac6c78920
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/004_mesonh_step2/EXSEG1.nam
@@ -0,0 +1,59 @@
+&NAM_CONFIO
+ LIO_ALLOW_NO_BACKUP = T , LIO_NO_WRITE = T
+/
+&NAM_CONFZ
+ ! NB_PROCIO_R=1 ,
+ ! NB_PROCIO_W=1 ,
+ LMNH_MPI_BSEND = F ,
+ ! MPI_BUFFER_SIZE = 40
+ !NZ_PROC=1
+/
+
+&NAM_LUNITn CINIFILE = "D0128.1.DAR01.004" , CINIFILEPGD = "PGD00128" /
+&NAM_CONFn LUSERV = T, LUSERC = T, LUSERR = T, LUSERI = T,
+ LUSERS = T, LUSERG = T, LUSERH = F, LUSECI = T /
+&NAM_DYNn XTSTEP = 15.00000 ,
+ CPRESOPT = "ZRESI",
+ !CPRESOPT = "ZSOLV",
+ NITR = 12,
+ XRELAX = 1., LHORELAX_UVWTH = T, LHORELAX_RV = T, LVE_RELAX = T,
+ NRIMX = 6, NRIMY = 6, XRIMKMAX = 0.0005, XT4DIFU = 4000 /
+&NAM_ADVn CUVW_ADV_SCHEME = "WENO_K", CMET_ADV_SCHEME = "PPM_01",
+ CSV_ADV_SCHEME = "PPM_01" /
+&NAM_PARAMn CTURB = "TKEL", ! "NONE" , ! "TKEL",
+ CRAD = "ECMWF",
+ CCLOUD = "ICE3", ! "NONE" , ! "ICE3",
+ CDCONV = "NONE", CSCONV = "NONE" /
+&NAM_PARAM_RADn XDTRAD = 1800., XDTRAD_CLONLY = 1800., LCLEAR_SKY = F,
+ NRAD_COLNBR = 500 /
+&NAM_PARAM_MFSHALLn XIMPL_MF = 1, CMF_UPDRAFT = "EDKF", CMF_CLOUD = "NONE",
+ LMIXUV = T, LMF_FLX = F /
+&NAM_LBCn CLBCX = 2*"OPEN", CLBCY = 2*"OPEN", XCPHASE = 20. /
+&NAM_TURBn XIMPL = 1., CTURBLEN = "DEAR", CTURBDIM = "3DIM",
+ LTURB_FLX = F, LTURB_DIAG = F, CSUBG_AUCV = "NONE", LSIGMAS = F,
+ LSIG_CONV = F, LSUBG_COND = F /
+&NAM_CH_MNHCn /
+&NAM_CONF CCONF = "RESTA", LFLAT = F, CEQNSYS = "DUR",
+ LLG=F, NMODEL = 1, NVERB = 0, CEXP = "D0128",
+ CSEG = "DAR02", CSPLIT = "BSPLITTING"
+ NHALO=1
+ LCHECK = F /
+&NAM_DYN XSEGLEN = 150.0 ! 150.0000 ,
+ LCORIO = T,
+ LNUMDIFU = F, LNUMDIFTH = F,
+ XALKTOP = 0.001, XALZBOT = 22000. /
+&NAM_BLANK /
+&NAM_NESTING /
+&NAM_BACKUP
+ XBAK_TIME(1,1) = 150.0 ! 3600.00000 ! 75. ! 3600.00000
+ XBAK_TIME(1,2) = 7200.00000
+ XBAK_TIME(1,3) = 10800.0000
+ XBAK_TIME(1,4) = 14400.0000
+ XBAK_TIME(1,5) = 18000.0000
+ XBAK_TIME(1,6) = 21600.0000
+ /
+&NAM_ISBAn /
+&NAM_SEAFLUXn /
+&NAM_DIAG_SURFn LSURF_BUDGET=T /
+&NAM_DIAG_ISBAn /
+&NAM_DIAG_SURF_ATMn /
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/004_mesonh_step2/MESONH_HOST2005 b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/004_mesonh_step2/MESONH_HOST2005
new file mode 100755
index 0000000000000000000000000000000000000000..07ad9fb226bfd7e7cc764749cb644e595ee25b13
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/004_mesonh_step2/MESONH_HOST2005
@@ -0,0 +1,5 @@
+#!/bin/bash
+export ACC_DEVICE=HOST
+ACC_DEVICE_TYPE=${ACC_DEVICE}
+exec /home/escj/DEV/MNH-55X-dev-OPENACC-juan-03/src/dir_obj-LXnvhpc2005-R8I4-MNH-V5-4-2-MPIAUTO-${OPTLEVEL}/ZSOLVER/MESONH
+
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/004_mesonh_step2/mppdb.nam b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/004_mesonh_step2/mppdb.nam
new file mode 100644
index 0000000000000000000000000000000000000000..39cb2adaff74600f63d7f252fc0a4429c9c2b67d
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/004_mesonh_step2/mppdb.nam
@@ -0,0 +1,7 @@
+&NAM_MPPDB
+MPPDB_DEBUG = .F. ,
+MPPDB_NBSON = 1 ,
+MPPDB_EXEC = "/home/escj/DEV/MNH-55X-dev-OPENACC-juan-03/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/004_mesonh_step2/MESONH_HOST2005" ,
+MPPDB_WDIR = "/home/escj/DEV/MNH-55X-dev-OPENACC-juan-03/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/004_mesonh_step2/dir_clone" ,
+MPPDB_CHECK_LB = .F.
+/
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/004_mesonh_step2/parameters_mg.nam b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/004_mesonh_step2/parameters_mg.nam
new file mode 100644
index 0000000000000000000000000000000000000000..5d76f4dff351b9c8873c37072d3f1f0747877d70
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/004_mesonh_step2/parameters_mg.nam
@@ -0,0 +1,118 @@
+! *********************************************************************
+! *********************************************************************
+! ** **
+! ** Parameter file for geometric multigrid code **
+! ** **
+! *********************************************************************
+! *********************************************************************
+!
+! *********************************************************************
+! * General parameters
+! *********************************************************************
+¶meters_general
+ savefields = .F. ! Save fields to disk?
+/
+
+! *********************************************************************
+! * General solver parameters
+! *********************************************************************
+¶meters_solver
+ LUseO = F ,
+ LUseT = T ,
+ solvertype = 1, ! Solver type:
+ ! 1 : Richardson iteration
+ ! 2 : Conjugate gradient
+ resreduction = 1e-1 ! 1.0d-10 ! Required relative residual reduction
+ maxiter = 5 ! 50 ! Maximal number of iterations
+/
+
+! *********************************************************************
+! * Conjugate gradient parameters
+! *********************************************************************
+¶meters_conjugategradient
+ verbose = 10, ! Verbosity level
+ maxiter = 5, ! Maximal number of iterations
+ resreduction = 1.0e-1, ! Target residual reduction
+ n_prec = 1 ! Number of smoother applications in
+ ! preconditioner (N.B.: Using 0 is
+ ! inefficient, as the identity is used
+ ! for preconditioning, instead of using
+ ! unpreconditioned CG.)
+/
+
+! *********************************************************************
+! * Grid parameters
+! *********************************************************************
+¶meters_grid
+ n = 128 , ! 64, ! Number of horizontal grid cells
+ nz = 72, ! Number of vertical grid cells
+ L = 204800.0 ! 16000.0 metre * 256 pt grille ! 1.0, ! Size in horizontal direction
+ H = 19200.0 ! 150.0 metre * 128 niveaux ! 0.01, ! Size in vertical direction
+ vertbc = 2, ! Boundary conditions at top and bottom of the
+ ! atmosphere. 1 = DIRICHLET, 2 = NEUMANN
+ ! Note that Neumann boundary conditions only work
+ ! for coarsening in the horizontal only, as they are
+ ! not yet implemented in the prolongation operator.
+ graded = .F. ! Is the grid graded in the vertical direction?
+/
+
+! *********************************************************************
+! * Parallel communication parameters
+! *********************************************************************
+¶meters_communication
+ halo_size = 1 ! Size of halos (has to be 1 or 2)
+/
+
+! *********************************************************************
+! * Model parameters
+! *********************************************************************
+!
+! parameters of the Helmholtz operator
+!
+! -omega2*(d^2/dy^2 + d^2/dy^2 + lambda2*d^2/dz^2) u + delta u = RHS
+!
+¶meters_model
+ omega2 = 1.0,
+ lambda2 = 1.0 ! 100.0, ! Vertical coupling
+ delta = 0.0d0 ! Size of constant term
+/
+
+! *********************************************************************
+! * Smoother parameters
+! *********************************************************************
+!
+! parameters of the smoother
+!
+¶meters_smoother
+ smoother = 3, ! Smoother method
+ ! 3 = line SOR
+ ! 4 = line SSOR
+ ! 6 = line Jacobi
+ ordering = 2, ! Ordering of grid points (for smoother)
+
+ ! 1 = lexicographic
+ ! 2 = red-black ordering
+ !rho = 0.6666666666666666d0 ! Overrelaxation parameter
+ rho = 0.8d0 ! Overrelaxation parameter
+/
+
+! *********************************************************************
+! * Multigrid parameters
+! *********************************************************************
+¶meters_multigrid
+ verbose = 10, ! Verbosity level
+ n_lev = 8, ! Number of levels
+ lev_split = 3, ! First level where data is pulled together
+ n_presmooth = 1, ! Number of presmoothing steps
+ n_postsmooth = 1, ! Number of postsmoothing steps
+ n_coarsegridsmooth = 1, ! Number of smoothing steps on coarsest level
+ prolongation = 2 ! 2 best after modif !! ! Prologation method
+ ! 1 = constant interpolation
+ ! 2 = (tri-) linear interpolation
+ restriction = 1, ! Restriction method
+ ! 1 = cell average
+ ! 2 = Khalil
+ coarsegridsolver = 1 ! Solver on coarsest grid
+ ! 1 = use smoother
+ ! 2 = Conjugate gradient
+/
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/004_mesonh_step2/run_mesonh_step2_xyz b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/004_mesonh_step2/run_mesonh_step2_xyz
new file mode 100755
index 0000000000000000000000000000000000000000..03470ed7a719790e446b249af6d459944b6f46f4
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/004_mesonh_step2/run_mesonh_step2_xyz
@@ -0,0 +1,26 @@
+#!/bin/bash
+#
+#
+set -x
+
+#export PGI_ACC_POOL_ALLOC=0
+export PGI_ACC_SYNCHRONOUS=1
+#export NVCOMPILER_ACC_POOL_ALLOC=0
+
+NP=${NP:-16}
+NPS=${NPS:-4}
+
+export MPIRUN=${MPIRUN:-"Mpirun -tag-output -report-bindings -bind-to none -map-by ppr:${NPS}:socket -x NVCOMPILER_ACC_POOL_ALLOC_MAXSIZE -x NVCOMPILER_ACC_POOL_ALLOC -x PGI_ACC_POOL_ALLOC -x PGI_ACC_SYNCHRONOUS -np ${NP} set_core_device_impair "}
+
+rm -f DA0128.*
+rm -f PGD00128.*
+rm -f D0128.1.DAR0?.* OUTPUT_LISTING1 pipe_name
+#
+ln -sf ../001_pgd1/PGD00128.* .
+ln -sf ../002_prep_ideal_case/DA0128.* .
+ln -sf ../003_mesonh_step1/D0128.1.DAR0?.* .
+#
+
+time ${MPIRUN} MESONH${XYZ}
+
+
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/004_mesonh_step2/run_mesonh_step2_xyz_mppdb b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/004_mesonh_step2/run_mesonh_step2_xyz_mppdb
new file mode 100755
index 0000000000000000000000000000000000000000..cd0f51b0bbd9e31776fa4e1cdb6abfbdb977a7f9
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/004_mesonh_step2/run_mesonh_step2_xyz_mppdb
@@ -0,0 +1,48 @@
+#!/bin/bash
+#
+#
+
+set -x
+DIRCLONE=dir_clone
+mkdir -p ${DIRCLONE}
+
+export PGI_ACC_POOL_ALLOC=0
+#export PGI_ACC_POOL_SIZE=100MB
+
+export NP=${NP:-1}
+
+export MPIRUN=${MPIRUN:-Mpirun -tag-output -np ${NP}}
+
+#export PGI_ACC_NOTIFY=$(( 1+2+4+8+16+32 )) PGI_ACC_DEBUG=$(( 1+2+4+8+16+32 ))
+#export PGI_ACC_TIME=1
+
+export PGI_ACC_SYNCHRONOUS=1
+
+(
+cd ${DIRCLONE}
+#
+rm -f DA0128.*
+rm -f PGD00128.*
+rm -f D0128.1.DAR0?.* OUTPUT_LISTING1 pipe_name
+#
+cp ../EXSEG1.nam .
+ln -sf ../../001_pgd1/PGD00128.* .
+ln -sf ../../002_prep_ideal_case/DA0128.* .
+ln -sf ../../003_mesonh_step1/D0128.1.DAR01.* .
+#
+)
+
+#
+rm -f DA0128.*
+rm -f PGD00128.*
+rm -f D0128.1.DAR0?.* OUTPUT_LISTING1 pipe_name
+#
+ln -sf ../001_pgd1/PGD00128.* .
+ln -sf ../002_prep_ideal_case/DA0128.* .
+ln -sf ../003_mesonh_step1/D0128.1.DAR0?.* .
+#
+
+killall -r 'MESONH.*'
+time ${MPIRUN} MESONH${XYZ}
+
+
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/005_mg_main_mnh_all/parameters_mg.nam b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/005_mg_main_mnh_all/parameters_mg.nam
new file mode 120000
index 0000000000000000000000000000000000000000..adeab0e1859deb44804ccc14f23c2ccb35775079
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/005_mg_main_mnh_all/parameters_mg.nam
@@ -0,0 +1 @@
+../004_mesonh_step2/parameters_mg.nam
\ No newline at end of file
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/005_mg_main_mnh_all/run_mg_main_mnh_all_xyz b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/005_mg_main_mnh_all/run_mg_main_mnh_all_xyz
new file mode 100755
index 0000000000000000000000000000000000000000..c3694f2de5a11d8600af4cf5874b3248c26057f4
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_128x128_MNH-55X/005_mg_main_mnh_all/run_mg_main_mnh_all_xyz
@@ -0,0 +1,18 @@
+#!/bin/bash
+#
+#
+set -x
+
+#export PGI_ACC_POOL_ALLOC=0
+export PGI_ACC_SYNCHRONOUS=1
+#export NVCOMPILER_ACC_POOL_ALLOC=0
+
+NP=${NP:-16}
+NPS=${NPS:-4}
+
+export MPIRUN=${MPIRUN:-"Mpirun -tag-output -report-bindings -bind-to none -map-by ppr:${NPS}:socket -x NVCOMPILER_ACC_POOL_ALLOC_MAXSIZE -x NVCOMPILER_ACC_POOL_ALLOC -x PGI_ACC_POOL_ALLOC -x PGI_ACC_SYNCHRONOUS -np ${NP} set_core_device_impair "}
+
+
+time ${MPIRUN} MG_MAIN_MNH_ALL${XYZ}
+
+
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_256x256_MNH-55X/001_pgd1/PRE_PGD1.nam b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_256x256_MNH-55X/001_pgd1/PRE_PGD1.nam
new file mode 100644
index 0000000000000000000000000000000000000000..a7edb249a9cc07b565c80de0c7b74d6a9528e220
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_256x256_MNH-55X/001_pgd1/PRE_PGD1.nam
@@ -0,0 +1,13 @@
+&NAM_CONFZ
+ NB_PROCIO_W=1 ,
+ !MPI_BUFFER_SIZE=40
+/
+&NAM_PGDFILE CPGDFILE="PGD00256" /
+&NAM_PGD_SCHEMES CNATURE='ISBA', CSEA='SEAFLX', CTOWN='NONE', CWATER='NONE' /
+&NAM_CONF_PROJ XLAT0=-11.5, XLON0=130.7, XRPK=0., XBETA=0. /
+&NAM_CONF_PROJ_GRID XLATCEN=-11.5, XLONCEN=130.7, NIMAX=256, NJMAX=256,
+ XDX=800.00000, XDY=800.00000 /
+&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' /
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_256x256_MNH-55X/001_pgd1/get_pgd_files b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_256x256_MNH-55X/001_pgd1/get_pgd_files
new file mode 100755
index 0000000000000000000000000000000000000000..56726933a0f7608e0c4918f42e80d5aea76b669d
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_256x256_MNH-55X/001_pgd1/get_pgd_files
@@ -0,0 +1,30 @@
+#
+# Modif
+# J.Escobar 11/04/2014 get PGD files from 'dir_open' directory ( without psswd )
+# J.Escobar 25/04/2013 get LICENCE files
+#
+PREP_PGD_FILES=${PREP_PGD_FILES:-"$HOME/PREP_PGD_FILES_WWW"} ; export PREP_PGD_FILES
+echo "!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!"
+echo "!!!! WARNING !!!!"
+echo "!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!"
+echo
+echo you need 3GO to download this files in
+echo
+echo PREP_PGD_FILES=$PREP_PGD_FILES
+echo
+echo if OK press ENTER else CTRL-C
+read RIEN
+set -x
+mkdir -p $PREP_PGD_FILES
+cd $PREP_PGD_FILES
+PGD_URL="http://mesonh.aero.obs-mip.fr/mesonh/dir_open/dir_PGDFILES"
+WGET="wget"
+export PGD_URL
+for file in LICENSE_ECOCLIMAP.txt LICENSE_soil_data.txt \
+ gtopo30.dir gtopo30.hdr \
+ SAND_HWSD_MOY.hdr SAND_HWSD_MOY.dir CLAY_HWSD_MOY.hdr CLAY_HWSD_MOY.dir \
+ ECOCLIMAP_v2.0.hdr ECOCLIMAP_v2.0.dir
+do
+[ -f $file ] || ( ${WGET} -c -nd $PGD_URL/$file.gz ; gunzip $file.gz ; )
+done
+
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_256x256_MNH-55X/001_pgd1/run_prep_pgd_xyz b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_256x256_MNH-55X/001_pgd1/run_prep_pgd_xyz
new file mode 100755
index 0000000000000000000000000000000000000000..c4b8da5bf37afca5b11f863c02bd21ae92cd9ea0
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_256x256_MNH-55X/001_pgd1/run_prep_pgd_xyz
@@ -0,0 +1,36 @@
+#!/bin/bash
+
+export ACC_DEVICE_TYPE=HOST
+export MPIRUN="Mpirun -tag-output -report-bindings -bind-to none -np 4 set_core_device_impair "
+
+
+PREP_PGD_FILES=${PREP_PGD_FILES:-"$HOME/PREP_PGD_FILES_WWW"} ; export PREP_PGD_FILES
+#
+if [ ! -d $PREP_PGD_FILES ]
+then
+cat << EOF
+
+Your directory PREP_PGD_FILES=$PREP_PGD_FILES
+
+containing the files gtopo30*, ECOCLIMAP_v2.0*
+doesn't exist ( or was not found !!! )
+use the script 'get_pgd_files' to download
+this files from the MESONH WEB PAGES !!!
+( or change the variable PREP_PGD_FILES ... )
+
+After fixing it , run this script again !!!
+
+EOF
+exit 1
+else
+set -x
+#
+rm -f OUTPUT_LISTING0 pipe_name
+rm -f gtopo30.??? sand_fao.??? clay_fao.???
+rm -f SAND_HWSD_MOY.??? ECOCLIMAP_v2.0.??? ecoclimats_v2.???
+rm -f PGD00256.*
+#
+ln -sf $PREP_PGD_FILES/*.dir $PREP_PGD_FILES/*.hdr .
+#
+time ${MPIRUN} PREP_PGD${XYZ}
+fi
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_256x256_MNH-55X/002_prep_ideal_case/PRE_IDEA1.nam b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_256x256_MNH-55X/002_prep_ideal_case/PRE_IDEA1.nam
new file mode 100644
index 0000000000000000000000000000000000000000..0f9bb438c273171b5f5d94e36e1c2a51e42f29a0
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_256x256_MNH-55X/002_prep_ideal_case/PRE_IDEA1.nam
@@ -0,0 +1,226 @@
+&NAM_CONFZ
+ !NB_PROCIO_R=1 ,
+ !NB_PROCIO_W=1 ,
+ !MPI_BUFFER_SIZE=40
+/
+
+&NAM_REAL_PGD
+ CPGD_FILE ="PGD00256" ,
+ LREAD_ZS =.TRUE., LREAD_GROUND_PARAM =.TRUE.
+/
+&NAM_DIMn_PRE NIMAX=40, NJMAX=32 /
+&NAM_CONF_PRE LCARTESIAN=.FALSE.
+ CIDEAL='RSOU' CZS='FLAT'
+ NVERB=0 /
+&NAM_CONFn LUSERV= T /
+&NAM_LUNITn CINIFILE = "DA0256" , CINIFILEPGD = "PGD00256" /
+&NAM_DYNn_PRE
+ CPRESOPT= 'ZRESI' ,
+ NITR=4 XRELAX=1. /
+&NAM_LBCn_PRE CLBCX= 2*'OPEN' CLBCY= 2*'OPEN' /
+&NAM_VPROF_PRE CTYPELOC='IJGRID' NILOC=1 NJLOC=1
+ CFUNU='ZZZ' CFUNV='ZZZ'
+ LGEOSBAL=.FALSE. /
+&NAM_VER_GRID NKMAX=126, YZGRID_TYPE='FUNCTN',
+ ZDZGRD=40., ZDZTOP=210., ZZMAX_STRGRD=2500.,
+ ZSTRGRD=7., ZSTRTOP=7. /
+&NAM_GRn_PRE
+ CSURF='EXTE'
+ ! CSURF='NONE'
+ /
+&NAM_PREP_ISBA XTG_SURF= 311., XTG_ROOT= 303., XTG_DEEP= 302.,
+ XHUG_SURF= 0.16, XHUG_ROOT= 0.16, XHUG_DEEP= 0.16 /
+&NAM_PREP_SEAFLUX XSST_UNIF= 304. /
+ RSOU
+ 2005 11 30 0
+ 'PUVTHDMR'
+ 0.0000000E+00
+ 100300.0
+ 303.3000
+ 1.9630000E-02
+ 92
+ 100000.0 3.637668 -3.637668
+ 99900.00 2.572220 -4.455216
+ 99100.00 -4.612377 0.4035313
+ 96600.00 5.124864 -0.4483674
+ 93000.00 5.836996 -2.009839
+ 92800.00 5.801031 -2.111402
+ 92500.00 5.466062 -1.464627
+ 85000.00 0.1345103 -1.537459
+ 81400.00 -0.9020693 -1.849517
+ 71800.00 -2.797467 1.304481
+ 70000.00 -4.472232 1.198331
+ 60500.00 -3.289776 -1.464702
+ 59900.00 -3.208611 -1.634869
+ 56900.00 -2.329535 -2.025034
+ 55600.00 -2.105100 -2.257443
+ 50100.00 -0.5325915 -1.987659
+ 50000.00 -0.2679967 -1.519885
+ 49500.00 0.2507798 -2.042438
+ 46700.00 5.466062 1.464627
+ 45300.00 5.594935 2.608961
+ 44800.00 5.450725 2.898202
+ 43700.00 4.994327 3.628591
+ 42200.00 4.561049 4.891127
+ 41900.00 4.387573 5.047325
+ 41800.00 4.298817 5.123130
+ 40000.00 6.305391 5.290851
+ 39000.00 5.970908 4.027428
+ 38600.00 5.671552 3.543979
+ 38300.00 5.346259 3.086664
+ 35500.00 2.976104 3.546782
+ 32700.00 4.455215 2.572220
+ 31400.00 3.712562 4.270814
+ 30700.00 3.086664 5.346259
+ 30000.00 3.540875 5.056894
+ 28800.00 3.786532 4.205370
+ 26700.00 3.329551 2.419061
+ 25000.00 2.797467 1.304481
+ 24400.00 2.819808 1.255459
+ 23800.00 2.841289 1.206056
+ 21400.00 3.424857 1.112804
+ 20000.00 3.478403 0.9320353
+ 19300.00 3.383935 1.231651
+ 17400.00 3.411947 2.301387
+ 15000.00 2.950729 4.214078
+ 12700.00 2.167201 2.875973
+ 10200.00 1.322713 1.576348
+ 10000.00 0.4466614 2.533142
+ 9120.000 -3.306782 3.940869
+ 8700.000 -7.950636 2.130367
+ 8200.000 -8.712269 0.7622259
+ 8000.000 -9.224754 -0.8070621
+ 7600.000 -4.009694 -2.314998
+ 7300.000 -5.066284 -0.8933228
+ 7000.000 -4.196201 -1.956721
+ 6700.000 -1.013257 0.1786647
+ 6400.000 -2.107039 1.475365
+ 6340.000 -2.558960 1.726041
+ 5480.000 -13.17234 2.322640
+ 5300.000 -15.37459 1.345105
+ 5000.000 -12.81216 -1.120920
+ 4900.000 -12.66571 -2.233307
+ 4860.000 -12.25463 -1.504677
+ 4700.000 -10.76221 0.9415731
+ 4500.000 -11.37801 7.966970
+ 3800.000 -25.11183 2.197004
+ 3600.000 -22.54940 -1.972818
+ 3200.000 -22.72465 -4.417223
+ 3000.000 -22.36116 -5.991656
+ 2490.000 -21.48829 -2.258509
+ 2300.000 -21.07936 -0.7361050
+ 2200.000 -20.57776 -1.7989648E-06
+ 2080.000 -19.00834 -0.9961855
+ 2000.000 -17.93702 -1.569287
+ 1600.000 -14.74927 -10.32755
+ 1500.000 -13.52109 -6.304989
+ 1430.000 -16.89510 -4.527030
+ 1400.000 -18.23862 -3.215962
+ 1300.000 -15.46941 -5.630407
+ 1200.000 -5.594935 2.608962
+ 1190.000 -5.671552 3.543979
+ 1100.000 1.786646 10.13257
+ 1050.000 -3.736842 7.333968
+ 1020.000 -5.435582 4.725079
+ 1000.000 -5.594935 2.608962
+ 950.0000 -4.347303 8.176087
+ 900.0000 2.143975 12.15908
+ 889.9999 1.031246 11.78719
+ 800.0000 -4.728969 4.728969
+ 760.0000 2.890712 7.154763
+ 710.0000 8.644947 -3.318484
+ 700.0000 7.487652 -6.282886
+ 600.0000 -4.214078 -2.950729
+ 91
+ 100000.0 303.1000 1.8050000E-02
+ 99900.00 303.1000 1.8040001E-02
+ 99100.00 303.1000 1.7969999E-02
+ 96600.00 303.1000 1.7729999E-02
+ 93000.00 303.0000 1.7370002E-02
+ 92800.00 303.1000 1.6980000E-02
+ 92500.00 303.2000 1.6400000E-02
+ 85000.00 307.1000 1.1190000E-02
+ 81400.00 309.6000 6.6700005E-03
+ 71800.00 312.3000 6.5600001E-03
+ 70000.00 312.9000 6.5400004E-03
+ 60500.00 316.2000 6.1800000E-03
+ 59900.00 318.1000 5.3400006E-03
+ 56900.00 318.7000 4.9899998E-03
+ 55600.00 322.2000 3.3100001E-03
+ 50100.00 325.5000 2.6000000E-03
+ 50000.00 325.5000 2.5800001E-03
+ 49500.00 325.7000 3.1700002E-03
+ 46700.00 328.6000 1.4500001E-03
+ 45300.00 330.1000 9.4000006E-04
+ 44800.00 330.1000 1.9600000E-03
+ 43700.00 331.0000 2.1500003E-03
+ 42200.00 332.2000 1.9400001E-03
+ 41900.00 332.9000 9.6000003E-04
+ 41800.00 333.0000 9.7000005E-04
+ 40000.00 334.2000 1.1600000E-03
+ 39000.00 334.6000 1.0300000E-03
+ 38600.00 335.6000 6.5000000E-04
+ 38300.00 335.8000 6.0000003E-04
+ 35500.00 338.2000 2.5000001E-04
+ 32700.00 340.8000 9.0000009E-05
+ 31400.00 342.0000 5.0000002E-05
+ 30700.00 342.1000 4.0000003E-05
+ 30000.00 342.3000 4.0000003E-05
+ 28800.00 342.3000 1.2000000E-04
+ 26700.00 344.8000 1.0000001E-05
+ 25000.00 345.4000 1.0000001E-05
+ 24400.00 346.3000 2.0000001E-05
+ 23800.00 346.4000 9.0000009E-05
+ 21400.00 348.1000 1.0000001E-05
+ 20000.00 348.8000 1.0000001E-05
+ 19300.00 348.9000 2.0000001E-05
+ 17400.00 351.5000 1.0000001E-05
+ 15000.00 354.0000 7.1596442E-06
+ 12700.00 355.7000 4.4376361E-06
+ 10200.00 364.6000 1.4789329E-06
+ 10000.00 365.4000 1.2422363E-06
+ 9120.000 371.6000 1.2919259E-06
+ 8700.000 379.9000 1.5496893E-06
+ 8200.000 390.6000 1.8819877E-06
+ 8000.000 395.1000 2.0217396E-06
+ 7600.000 404.7000 2.3198761E-06
+ 7300.000 412.3000 2.4844724E-06
+ 7000.000 420.4000 2.4844724E-06
+ 6700.000 431.1000 2.4844724E-06
+ 6400.000 442.4000 2.4844724E-06
+ 6340.000 444.8000 2.4844724E-06
+ 5480.000 467.8000 2.4844724E-06
+ 5300.000 471.0000 2.4844724E-06
+ 5000.000 476.5000 2.4844724E-06
+ 4900.000 477.6000 2.4844724E-06
+ 4860.000 478.0000 2.4844724E-06
+ 4700.000 485.1000 2.4844724E-06
+ 4500.000 494.6000 2.4844724E-06
+ 3800.000 532.9000 2.4844724E-06
+ 3600.000 545.6000 2.4844724E-06
+ 3200.000 574.4000 2.4844724E-06
+ 3000.000 585.1000 2.4844724E-06
+ 2490.000 624.0000 2.4844724E-06
+ 2300.000 648.9000 2.4844724E-06
+ 2200.000 655.4000 2.4844724E-06
+ 2080.000 663.6000 2.4844724E-06
+ 2000.000 673.5000 2.4844724E-06
+ 1600.000 726.4000 2.4844724E-06
+ 1500.000 742.5000 2.4844724E-06
+ 1430.000 750.7000 2.4844724E-06
+ 1400.000 759.0000 2.4844724E-06
+ 1300.000 788.6000 2.4844724E-06
+ 1200.000 821.7000 2.4844724E-06
+ 1190.000 825.2000 2.4844724E-06
+ 1100.000 833.9000 2.4844724E-06
+ 1050.000 839.0000 2.4844724E-06
+ 1020.000 845.3000 2.4844724E-06
+ 1000.000 851.6000 2.4844724E-06
+ 950.0000 877.0000 2.4844724E-06
+ 900.0000 885.6000 2.4844724E-06
+ 889.9999 887.3000 2.4844724E-06
+ 800.0000 926.6000 2.4844724E-06
+ 760.0000 946.0000 2.4844724E-06
+ 710.0000 958.0000 2.4844724E-06
+ 700.0000 963.1000 2.4844724E-06
+ 600.0000 20.30000 2.4844724E-06
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_256x256_MNH-55X/002_prep_ideal_case/run_prep_ideal_xyz b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_256x256_MNH-55X/002_prep_ideal_case/run_prep_ideal_xyz
new file mode 100755
index 0000000000000000000000000000000000000000..6493fd00e6f6d69c0ec9de8b6588c82c816bd577
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_256x256_MNH-55X/002_prep_ideal_case/run_prep_ideal_xyz
@@ -0,0 +1,20 @@
+#!/bin/bash
+
+NP=${NP:-16}
+NPS=${NPS:-4}
+
+export ACC_DEVICE_TYPE=HOST
+export MPIRUN="Mpirun -tag-output -report-bindings --mca mpi_cuda_support 0 -map-by ppr:${NPS}:socket -bind-to none -np ${NP} set_core_device_impair"
+
+#
+set -x
+set -e
+#
+rm -f OUTPUT_LISTING1 pipe_name file_for_xtransfer
+rm -f PGD00256.*
+rm -f DA0256.*
+#
+ln -sf ../001_pgd1/PGD00256.* .
+#
+time ${MPIRUN} PREP_IDEAL_CASE${XYZ}
+#
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_256x256_MNH-55X/003_mesonh_step1/EXSEG1.nam b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_256x256_MNH-55X/003_mesonh_step1/EXSEG1.nam
new file mode 100644
index 0000000000000000000000000000000000000000..1f733b80a6013aaa15ffc2721e7b9076b87bd5b8
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_256x256_MNH-55X/003_mesonh_step1/EXSEG1.nam
@@ -0,0 +1,53 @@
+&NAM_CONFZ
+ ! NB_PROCIO_R=1 ,
+ ! NB_PROCIO_W=1 ,
+ ! LMNH_MPI_BSEND = T ,
+ MPI_BUFFER_SIZE = 200
+/
+
+&NAM_LUNITn CINIFILE = "DA0256" , CINIFILEPGD = "PGD00256" /
+&NAM_CONFn LUSERV = T, LUSERC = T, LUSERR = T, LUSERI = T,
+ LUSERS = T, LUSERG = T, LUSERH = F, LUSECI = T /
+&NAM_DYNn XTSTEP = 15.00000 ,
+ CPRESOPT = "ZRESI", NITR = 12
+ !CPRESOPT = "ZSOLV", NITR = 12
+ XRELAX = 1., LHORELAX_UVWTH = T, LHORELAX_RV = T, LVE_RELAX = T,
+ NRIMX = 6, NRIMY = 6, XRIMKMAX = 0.0005, XT4DIFU = 4000 /
+&NAM_ADVn CUVW_ADV_SCHEME = "WENO_K", CMET_ADV_SCHEME = "PPM_01",
+ CSV_ADV_SCHEME = "PPM_01" /
+&NAM_PARAMn CTURB = "TKEL" ! "NONE" ! "TKEL",
+ CRAD = "ECMW",
+ CCLOUD = "ICE3" ! "NONE" ! "ICE3",
+ CDCONV = "NONE", CSCONV = "NONE" /
+&NAM_PARAM_RADn XDTRAD = 1800., XDTRAD_CLONLY = 1800., LCLEAR_SKY = F,
+ NRAD_COLNBR = 500 /
+&NAM_PARAM_MFSHALLn XIMPL_MF = 1, CMF_UPDRAFT = "EDKF", CMF_CLOUD = "NONE",
+ LMIXUV = T, LMF_FLX = F /
+&NAM_LBCn CLBCX = 2*"OPEN", CLBCY = 2*"OPEN", XCPHASE = 20. /
+&NAM_TURBn XIMPL = 1., CTURBLEN = "DEAR", CTURBDIM = "3DIM",
+ LTURB_FLX = F, LTURB_DIAG = F, CSUBG_AUCV = "NONE", LSIGMAS = F,
+ LSIG_CONV = F, LSUBG_COND = F /
+&NAM_CH_MNHCn /
+&NAM_CONF CCONF = "START", LFLAT = F, CEQNSYS = "DUR",
+ LLG=F, NMODEL = 1, NVERB = 0, CEXP = "D0256",
+ CSEG = "DAR01", CSPLIT = "BSPLITTING"
+ NHALO=1
+ /
+&NAM_DYN XSEGLEN = 14400.0000 ! 21600.0000 ! 150. ! 14400.0000 , LCORIO = T,
+ LNUMDIFU = F, LNUMDIFTH = F,
+ XALKTOP = 0.001, XALZBOT = 22000. /
+&NAM_BLANK /
+&NAM_NESTING /
+&NAM_BACKUP
+ XBAK_TIME(1,1) = 3600.00000 ! 150. ! 3600.00000
+ XBAK_TIME(1,2) = 7200.00000
+ XBAK_TIME(1,3) = 10800.0000
+ XBAK_TIME(1,4) = 14400.0000
+ XBAK_TIME(1,5) = 18000.0000
+ XBAK_TIME(1,6) = 21600.0000
+ /
+&NAM_ISBAn /
+&NAM_SEAFLUXn /
+&NAM_DIAG_SURFn LSURF_BUDGET=T /
+&NAM_DIAG_ISBAn /
+&NAM_DIAG_SURF_ATMn /
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_256x256_MNH-55X/003_mesonh_step1/parameters_mg.nam b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_256x256_MNH-55X/003_mesonh_step1/parameters_mg.nam
new file mode 100644
index 0000000000000000000000000000000000000000..209e1d572e4d230aa3a8ad8fd4fc9450a5652197
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_256x256_MNH-55X/003_mesonh_step1/parameters_mg.nam
@@ -0,0 +1,120 @@
+! *********************************************************************
+! *********************************************************************
+! ** **
+! ** Parameter file for geometric multigrid code **
+! ** **
+! *********************************************************************
+! *********************************************************************
+!
+! *********************************************************************
+! * General parameters
+! *********************************************************************
+¶meters_general
+ savefields = .F. ! Save fields to disk?
+/
+
+! *********************************************************************
+! * General solver parameters
+! *********************************************************************
+¶meters_solver
+ LUseO = F ,
+ LUseT = T ,
+ solvertype = 1, ! Solver type:
+ ! 1 : Richardson iteration
+ ! 2 : Conjugate gradient
+ resreduction = 1e-1 ! 1.0d-15 ! Required relative residual reduction
+ maxiter = 5 ! Maximal number of iterations
+/
+
+! *********************************************************************
+! * Conjugate gradient parameters
+! *********************************************************************
+¶meters_conjugategradient
+ verbose = 10, ! Verbosity level
+ maxiter = 5, ! Maximal number of iterations
+ resreduction = 1.0e-1, ! Target residual reduction
+ n_prec = 1 ! Number of smoother applications in
+ ! preconditioner (N.B.: Using 0 is
+ ! inefficient, as the identity is used
+ ! for preconditioning, instead of using
+ ! unpreconditioned CG.)
+/
+
+! *********************************************************************
+! * Grid parameters
+! *********************************************************************
+¶meters_grid
+ n = 256 , ! 64, ! Number of horizontal grid cells
+ nz = 128, ! Number of vertical grid cells
+ L = 204800.0 ! 8000.0 metre * 256 pt grille ! 1.0, ! Size in horizontal direction
+ H = 19200.0 ! 150.0 metre * 128 niveaux ! 0.01, ! Size in vertical direction
+ vertbc = 2, ! Boundary conditions at top and bottom of the
+ ! atmosphere. 1 = DIRICHLET, 2 = NEUMANN
+ ! Note that Neumann boundary conditions only work
+ ! for coarsening in the horizontal only, as they are
+ ! not yet implemented in the prolongation operator.
+ graded = .F. ! Is the grid graded in the vertical direction?
+/
+
+! *********************************************************************
+! * Parallel communication parameters
+! *********************************************************************
+¶meters_communication
+ halo_size = 1 ! Size of halos (has to be 1 or 2)
+/
+
+! *********************************************************************
+! * Model parameters
+! *********************************************************************
+!
+! parameters of the Helmholtz operator
+!
+! -omega2*(d^2/dy^2 + d^2/dy^2 + lambda2*d^2/dz^2) u + delta u = RHS
+!
+¶meters_model
+ omega2 = 1.0,
+ lambda2 = 1.0 ! 100.0, ! Vertical coupling
+ delta = 0.0d0 ! Size of constant term
+/
+
+! *********************************************************************
+! * Smoother parameters
+! *********************************************************************
+!
+! parameters of the smoother
+!
+¶meters_smoother
+ smoother = 3, ! Smoother method
+ ! 3 = line SOR
+ ! 4 = line SSOR
+ ! 6 = line Jacobi
+ ordering = 2, ! Ordering of grid points (for smoother)
+
+ ! 1 = lexicographic
+ ! 2 = red-black ordering
+! rho = 0.5d0 ! Overrelaxation parameter
+ !rho = 0.6666666666666666d0 ! Overrelaxation parameter
+ rho = 0.8d0 ! Overrelaxation parameter
+! rho = 1.0d0 ! Overrelaxation parameter
+/
+
+! *********************************************************************
+! * Multigrid parameters
+! *********************************************************************
+¶meters_multigrid
+ verbose = 10, ! Verbosity level
+ n_lev = 9, ! Number of levels
+ lev_split = 5, ! First level where data is pulled together
+ n_presmooth = 1, ! Number of presmoothing steps
+ n_postsmooth = 1, ! Number of postsmoothing steps
+ n_coarsegridsmooth = 1, ! Number of smoothing steps on coarsest level
+ prolongation = 2 ! 2 best after modif !! ! Prologation method
+ ! 1 = constant interpolation
+ ! 2 = (tri-) linear interpolation
+ restriction = 1, ! Restriction method
+ ! 1 = cell average
+ ! 2 = Khalil
+ coarsegridsolver = 1 ! Solver on coarsest grid
+ ! 1 = use smoother
+ ! 2 = Conjugate gradient
+/
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_256x256_MNH-55X/003_mesonh_step1/run_mesonh_step1_xyz b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_256x256_MNH-55X/003_mesonh_step1/run_mesonh_step1_xyz
new file mode 100755
index 0000000000000000000000000000000000000000..2526cfb89e68ce8b7e9cb46f87f68aab266c35d4
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_256x256_MNH-55X/003_mesonh_step1/run_mesonh_step1_xyz
@@ -0,0 +1,23 @@
+
+
+#export PGI_ACC_POOL_ALLOC=0
+export PGI_ACC_SYNCHRONOUS=1
+
+NP=${NP:-16}
+NPS=${NPS:-4}
+
+export MPIRUN=${MPIRUN:-"Mpirun -tag-output -map-by ppr:${NPS}:socket -bind-to none -x PGI_ACC_POOL_ALLOC -x PGI_ACC_SYNCHRONOUS -np ${NP} set_core_device_impair "}
+
+
+#
+set -x
+set -e
+#
+rm -f DA0256.*
+rm -f PGD00256.*
+rm -f D0256.1.DAR01.* OUTPUT_LISTING* pipe_name
+#
+ln -sf ../001_pgd1/PGD00256.{des,nc} .
+ln -sf ../002_prep_ideal_case/DA0256.{des,nc} .
+#
+time ${MPIRUN} MESONH${XYZ}
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_256x256_MNH-55X/004_mesonh_step2/EXSEG1.nam b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_256x256_MNH-55X/004_mesonh_step2/EXSEG1.nam
new file mode 100644
index 0000000000000000000000000000000000000000..3f1d451f31d3737a664765827fb8bf51c1e0252a
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_256x256_MNH-55X/004_mesonh_step2/EXSEG1.nam
@@ -0,0 +1,59 @@
+&NAM_CONFIO
+ LIO_ALLOW_NO_BACKUP = T , LIO_NO_WRITE = T
+/
+&NAM_CONFZ
+ ! NB_PROCIO_R=1 ,
+ ! NB_PROCIO_W=1 ,
+ LMNH_MPI_BSEND = F ,
+ ! MPI_BUFFER_SIZE = 40
+ !NZ_PROC=1
+/
+
+&NAM_LUNITn CINIFILE = "D0256.1.DAR01.004" , CINIFILEPGD = "PGD00256" /
+&NAM_CONFn LUSERV = T, LUSERC = T, LUSERR = T, LUSERI = T,
+ LUSERS = T, LUSERG = T, LUSERH = F, LUSECI = T /
+&NAM_DYNn XTSTEP = 15.00000 ,
+ !CPRESOPT = "ZRESI",
+ CPRESOPT = "ZSOLV",
+ NITR = 12,
+ XRELAX = 1., LHORELAX_UVWTH = T, LHORELAX_RV = T, LVE_RELAX = T,
+ NRIMX = 6, NRIMY = 6, XRIMKMAX = 0.0005, XT4DIFU = 4000 /
+&NAM_ADVn CUVW_ADV_SCHEME = "WENO_K", CMET_ADV_SCHEME = "PPM_01",
+ CSV_ADV_SCHEME = "PPM_01" /
+&NAM_PARAMn CTURB = "TKEL", ! "NONE" , ! "TKEL",
+ CRAD = "ECMWF",
+ CCLOUD = "ICE3", ! "NONE" , ! "ICE3",
+ CDCONV = "NONE", CSCONV = "NONE" /
+&NAM_PARAM_RADn XDTRAD = 1800., XDTRAD_CLONLY = 1800., LCLEAR_SKY = F,
+ NRAD_COLNBR = 500 /
+&NAM_PARAM_MFSHALLn XIMPL_MF = 1, CMF_UPDRAFT = "EDKF", CMF_CLOUD = "NONE",
+ LMIXUV = T, LMF_FLX = F /
+&NAM_LBCn CLBCX = 2*"OPEN", CLBCY = 2*"OPEN", XCPHASE = 20. /
+&NAM_TURBn XIMPL = 1., CTURBLEN = "DEAR", CTURBDIM = "3DIM",
+ LTURB_FLX = F, LTURB_DIAG = F, CSUBG_AUCV = "NONE", LSIGMAS = F,
+ LSIG_CONV = F, LSUBG_COND = F /
+&NAM_CH_MNHCn /
+&NAM_CONF CCONF = "RESTA", LFLAT = F, CEQNSYS = "DUR",
+ LLG=F, NMODEL = 1, NVERB = 0, CEXP = "D0256",
+ CSEG = "DAR02", CSPLIT = "BSPLITTING"
+ NHALO=1
+ LCHECK = F /
+&NAM_DYN XSEGLEN = 150.0 ! 150.0000 ,
+ LCORIO = T,
+ LNUMDIFU = F, LNUMDIFTH = F,
+ XALKTOP = 0.001, XALZBOT = 22000. /
+&NAM_BLANK /
+&NAM_NESTING /
+&NAM_BACKUP
+ XBAK_TIME(1,1) = 150.0 ! 3600.00000 ! 75. ! 3600.00000
+ XBAK_TIME(1,2) = 7200.00000
+ XBAK_TIME(1,3) = 10800.0000
+ XBAK_TIME(1,4) = 14400.0000
+ XBAK_TIME(1,5) = 18000.0000
+ XBAK_TIME(1,6) = 21600.0000
+ /
+&NAM_ISBAn /
+&NAM_SEAFLUXn /
+&NAM_DIAG_SURFn LSURF_BUDGET=T /
+&NAM_DIAG_ISBAn /
+&NAM_DIAG_SURF_ATMn /
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_256x256_MNH-55X/004_mesonh_step2/parameters_mg.nam b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_256x256_MNH-55X/004_mesonh_step2/parameters_mg.nam
new file mode 100644
index 0000000000000000000000000000000000000000..209e1d572e4d230aa3a8ad8fd4fc9450a5652197
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_256x256_MNH-55X/004_mesonh_step2/parameters_mg.nam
@@ -0,0 +1,120 @@
+! *********************************************************************
+! *********************************************************************
+! ** **
+! ** Parameter file for geometric multigrid code **
+! ** **
+! *********************************************************************
+! *********************************************************************
+!
+! *********************************************************************
+! * General parameters
+! *********************************************************************
+¶meters_general
+ savefields = .F. ! Save fields to disk?
+/
+
+! *********************************************************************
+! * General solver parameters
+! *********************************************************************
+¶meters_solver
+ LUseO = F ,
+ LUseT = T ,
+ solvertype = 1, ! Solver type:
+ ! 1 : Richardson iteration
+ ! 2 : Conjugate gradient
+ resreduction = 1e-1 ! 1.0d-15 ! Required relative residual reduction
+ maxiter = 5 ! Maximal number of iterations
+/
+
+! *********************************************************************
+! * Conjugate gradient parameters
+! *********************************************************************
+¶meters_conjugategradient
+ verbose = 10, ! Verbosity level
+ maxiter = 5, ! Maximal number of iterations
+ resreduction = 1.0e-1, ! Target residual reduction
+ n_prec = 1 ! Number of smoother applications in
+ ! preconditioner (N.B.: Using 0 is
+ ! inefficient, as the identity is used
+ ! for preconditioning, instead of using
+ ! unpreconditioned CG.)
+/
+
+! *********************************************************************
+! * Grid parameters
+! *********************************************************************
+¶meters_grid
+ n = 256 , ! 64, ! Number of horizontal grid cells
+ nz = 128, ! Number of vertical grid cells
+ L = 204800.0 ! 8000.0 metre * 256 pt grille ! 1.0, ! Size in horizontal direction
+ H = 19200.0 ! 150.0 metre * 128 niveaux ! 0.01, ! Size in vertical direction
+ vertbc = 2, ! Boundary conditions at top and bottom of the
+ ! atmosphere. 1 = DIRICHLET, 2 = NEUMANN
+ ! Note that Neumann boundary conditions only work
+ ! for coarsening in the horizontal only, as they are
+ ! not yet implemented in the prolongation operator.
+ graded = .F. ! Is the grid graded in the vertical direction?
+/
+
+! *********************************************************************
+! * Parallel communication parameters
+! *********************************************************************
+¶meters_communication
+ halo_size = 1 ! Size of halos (has to be 1 or 2)
+/
+
+! *********************************************************************
+! * Model parameters
+! *********************************************************************
+!
+! parameters of the Helmholtz operator
+!
+! -omega2*(d^2/dy^2 + d^2/dy^2 + lambda2*d^2/dz^2) u + delta u = RHS
+!
+¶meters_model
+ omega2 = 1.0,
+ lambda2 = 1.0 ! 100.0, ! Vertical coupling
+ delta = 0.0d0 ! Size of constant term
+/
+
+! *********************************************************************
+! * Smoother parameters
+! *********************************************************************
+!
+! parameters of the smoother
+!
+¶meters_smoother
+ smoother = 3, ! Smoother method
+ ! 3 = line SOR
+ ! 4 = line SSOR
+ ! 6 = line Jacobi
+ ordering = 2, ! Ordering of grid points (for smoother)
+
+ ! 1 = lexicographic
+ ! 2 = red-black ordering
+! rho = 0.5d0 ! Overrelaxation parameter
+ !rho = 0.6666666666666666d0 ! Overrelaxation parameter
+ rho = 0.8d0 ! Overrelaxation parameter
+! rho = 1.0d0 ! Overrelaxation parameter
+/
+
+! *********************************************************************
+! * Multigrid parameters
+! *********************************************************************
+¶meters_multigrid
+ verbose = 10, ! Verbosity level
+ n_lev = 9, ! Number of levels
+ lev_split = 5, ! First level where data is pulled together
+ n_presmooth = 1, ! Number of presmoothing steps
+ n_postsmooth = 1, ! Number of postsmoothing steps
+ n_coarsegridsmooth = 1, ! Number of smoothing steps on coarsest level
+ prolongation = 2 ! 2 best after modif !! ! Prologation method
+ ! 1 = constant interpolation
+ ! 2 = (tri-) linear interpolation
+ restriction = 1, ! Restriction method
+ ! 1 = cell average
+ ! 2 = Khalil
+ coarsegridsolver = 1 ! Solver on coarsest grid
+ ! 1 = use smoother
+ ! 2 = Conjugate gradient
+/
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_256x256_MNH-55X/004_mesonh_step2/run_mesonh_step2_xyz b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_256x256_MNH-55X/004_mesonh_step2/run_mesonh_step2_xyz
new file mode 100755
index 0000000000000000000000000000000000000000..4ed5fb5e7a97f97e04e8be9b56d1e0fbdacb2558
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_256x256_MNH-55X/004_mesonh_step2/run_mesonh_step2_xyz
@@ -0,0 +1,26 @@
+#!/bin/bash
+#
+#
+set -x
+
+#export PGI_ACC_POOL_ALLOC=0
+export PGI_ACC_SYNCHRONOUS=1
+#export NVCOMPILER_ACC_POOL_ALLOC=0
+
+NP=${NP:-16}
+NPS=${NPS:-4}
+
+export MPIRUN=${MPIRUN:-"Mpirun -tag-output -report-bindings -bind-to none -map-by ppr:${NPS}:socket -x NVCOMPILER_ACC_POOL_ALLOC_MAXSIZE -x NVCOMPILER_ACC_POOL_ALLOC -x PGI_ACC_POOL_ALLOC -x PGI_ACC_SYNCHRONOUS -np ${NP} set_core_device_impair "}
+
+rm -f DA0256.*
+rm -f PGD00256.*
+rm -f D0256.1.DAR0?.* OUTPUT_LISTING1 pipe_name
+#
+ln -sf ../001_pgd1/PGD00256.* .
+ln -sf ../002_prep_ideal_case/DA0256.* .
+ln -sf ../003_mesonh_step1/D0256.1.DAR0?.* .
+#
+
+time ${MPIRUN} MESONH${XYZ}
+
+
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_256x256_MNH-55X/005_mg_main_mnh_all/parameters_mg.nam b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_256x256_MNH-55X/005_mg_main_mnh_all/parameters_mg.nam
new file mode 120000
index 0000000000000000000000000000000000000000..adeab0e1859deb44804ccc14f23c2ccb35775079
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_256x256_MNH-55X/005_mg_main_mnh_all/parameters_mg.nam
@@ -0,0 +1 @@
+../004_mesonh_step2/parameters_mg.nam
\ No newline at end of file
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_256x256_MNH-55X/005_mg_main_mnh_all/run_mg_main_mnh_all_xyz b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_256x256_MNH-55X/005_mg_main_mnh_all/run_mg_main_mnh_all_xyz
new file mode 100755
index 0000000000000000000000000000000000000000..c3694f2de5a11d8600af4cf5874b3248c26057f4
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_256x256_MNH-55X/005_mg_main_mnh_all/run_mg_main_mnh_all_xyz
@@ -0,0 +1,18 @@
+#!/bin/bash
+#
+#
+set -x
+
+#export PGI_ACC_POOL_ALLOC=0
+export PGI_ACC_SYNCHRONOUS=1
+#export NVCOMPILER_ACC_POOL_ALLOC=0
+
+NP=${NP:-16}
+NPS=${NPS:-4}
+
+export MPIRUN=${MPIRUN:-"Mpirun -tag-output -report-bindings -bind-to none -map-by ppr:${NPS}:socket -x NVCOMPILER_ACC_POOL_ALLOC_MAXSIZE -x NVCOMPILER_ACC_POOL_ALLOC -x PGI_ACC_POOL_ALLOC -x PGI_ACC_SYNCHRONOUS -np ${NP} set_core_device_impair "}
+
+
+time ${MPIRUN} MG_MAIN_MNH_ALL${XYZ}
+
+
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_512x512_MNH-55X/001_pgd1/PRE_PGD1.nam b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_512x512_MNH-55X/001_pgd1/PRE_PGD1.nam
new file mode 100644
index 0000000000000000000000000000000000000000..07f27620734c560d2eeede6cb11f8ca3ccb8770b
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_512x512_MNH-55X/001_pgd1/PRE_PGD1.nam
@@ -0,0 +1,13 @@
+&NAM_CONFZ
+ NB_PROCIO_W=1 ,
+ !MPI_BUFFER_SIZE=40
+/
+&NAM_PGDFILE CPGDFILE="PGD00512" NHALO=0 /
+&NAM_PGD_SCHEMES CNATURE='ISBA', CSEA='SEAFLX', CTOWN='NONE', CWATER='NONE' /
+&NAM_CONF_PROJ XLAT0=-11.5, XLON0=130.7, XRPK=0., XBETA=0. /
+&NAM_CONF_PROJ_GRID XLATCEN=-11.5, XLONCEN=130.7, NIMAX=512, NJMAX=512,
+ XDX=400.00000, XDY=400.00000 /
+&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' /
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_512x512_MNH-55X/001_pgd1/get_pgd_files b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_512x512_MNH-55X/001_pgd1/get_pgd_files
new file mode 100755
index 0000000000000000000000000000000000000000..56726933a0f7608e0c4918f42e80d5aea76b669d
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_512x512_MNH-55X/001_pgd1/get_pgd_files
@@ -0,0 +1,30 @@
+#
+# Modif
+# J.Escobar 11/04/2014 get PGD files from 'dir_open' directory ( without psswd )
+# J.Escobar 25/04/2013 get LICENCE files
+#
+PREP_PGD_FILES=${PREP_PGD_FILES:-"$HOME/PREP_PGD_FILES_WWW"} ; export PREP_PGD_FILES
+echo "!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!"
+echo "!!!! WARNING !!!!"
+echo "!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!"
+echo
+echo you need 3GO to download this files in
+echo
+echo PREP_PGD_FILES=$PREP_PGD_FILES
+echo
+echo if OK press ENTER else CTRL-C
+read RIEN
+set -x
+mkdir -p $PREP_PGD_FILES
+cd $PREP_PGD_FILES
+PGD_URL="http://mesonh.aero.obs-mip.fr/mesonh/dir_open/dir_PGDFILES"
+WGET="wget"
+export PGD_URL
+for file in LICENSE_ECOCLIMAP.txt LICENSE_soil_data.txt \
+ gtopo30.dir gtopo30.hdr \
+ SAND_HWSD_MOY.hdr SAND_HWSD_MOY.dir CLAY_HWSD_MOY.hdr CLAY_HWSD_MOY.dir \
+ ECOCLIMAP_v2.0.hdr ECOCLIMAP_v2.0.dir
+do
+[ -f $file ] || ( ${WGET} -c -nd $PGD_URL/$file.gz ; gunzip $file.gz ; )
+done
+
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_512x512_MNH-55X/001_pgd1/run_prep_pgd_xyz b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_512x512_MNH-55X/001_pgd1/run_prep_pgd_xyz
new file mode 100755
index 0000000000000000000000000000000000000000..342c7e22cf511ecaa4bb018bf4e51d45b47757e3
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_512x512_MNH-55X/001_pgd1/run_prep_pgd_xyz
@@ -0,0 +1,40 @@
+#!/bin/bash
+
+export ACC_DEVICE_TYPE=HOST
+
+NP=${NP:-16}
+NPS=${NPS:-1}
+
+export MPIRUN=${MPIRUN:-"Mpirun -tag-output --mca mpi_cuda_support 0 -map-by ppr:${NPS}:socket -bind-to none -np ${NP} set_core_device_impair "}
+
+
+PREP_PGD_FILES=${PREP_PGD_FILES:-"$HOME/PREP_PGD_FILES_WWW"} ; export PREP_PGD_FILES
+#
+if [ ! -d $PREP_PGD_FILES ]
+then
+cat << EOF
+
+Your directory PREP_PGD_FILES=$PREP_PGD_FILES
+
+containing the files gtopo30*, ECOCLIMAP_v2.0*
+doesn't exist ( or was not found !!! )
+use the script 'get_pgd_files' to download
+this files from the MESONH WEB PAGES !!!
+( or change the variable PREP_PGD_FILES ... )
+
+After fixing it , run this script again !!!
+
+EOF
+exit 1
+else
+set -x
+#
+rm -f OUTPUT_LISTING0 pipe_name
+rm -f gtopo30.??? sand_fao.??? clay_fao.???
+rm -f SAND_HWSD_MOY.??? ECOCLIMAP_v2.0.??? ecoclimats_v2.???
+rm -f PGD00512.*
+#
+ln -sf $PREP_PGD_FILES/*.dir $PREP_PGD_FILES/*.hdr .
+#
+time ${MPIRUN} PREP_PGD${XYZ}
+fi
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_512x512_MNH-55X/002_prep_ideal_case/PRE_IDEA1.nam b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_512x512_MNH-55X/002_prep_ideal_case/PRE_IDEA1.nam
new file mode 100644
index 0000000000000000000000000000000000000000..5b3f7253b15aba7334130cfad5fab26c738e3793
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_512x512_MNH-55X/002_prep_ideal_case/PRE_IDEA1.nam
@@ -0,0 +1,226 @@
+&NAM_CONFZ
+ !NB_PROCIO_R=1 ,
+ !NB_PROCIO_W=1 ,
+ MPI_BUFFER_SIZE=200
+/
+
+&NAM_REAL_PGD
+ CPGD_FILE ="PGD00512" ,
+ LREAD_ZS =.TRUE., LREAD_GROUND_PARAM =.TRUE.
+/
+&NAM_DIMn_PRE NIMAX=40, NJMAX=32 /
+&NAM_CONF_PRE LCARTESIAN=.FALSE.
+ CIDEAL='RSOU' CZS='FLAT'
+ NVERB=0 /
+&NAM_CONFn LUSERV= T /
+&NAM_LUNITn CINIFILE = "DA0512" , CINIFILEPGD = "PGD00512" /
+&NAM_DYNn_PRE
+ CPRESOPT= 'ZRESI' ,
+ NITR=4 XRELAX=1. /
+&NAM_LBCn_PRE CLBCX= 2*'OPEN' CLBCY= 2*'OPEN' /
+&NAM_VPROF_PRE CTYPELOC='IJGRID' NILOC=1 NJLOC=1
+ CFUNU='ZZZ' CFUNV='ZZZ'
+ LGEOSBAL=.FALSE. /
+&NAM_VER_GRID NKMAX=126, YZGRID_TYPE='FUNCTN',
+ ZDZGRD=40., ZDZTOP=210., ZZMAX_STRGRD=2500.,
+ ZSTRGRD=7., ZSTRTOP=7. /
+&NAM_GRn_PRE
+ CSURF='EXTE'
+ ! CSURF='NONE'
+ /
+&NAM_PREP_ISBA XTG_SURF= 311., XTG_ROOT= 303., XTG_DEEP= 302.,
+ XHUG_SURF= 0.16, XHUG_ROOT= 0.16, XHUG_DEEP= 0.16 /
+&NAM_PREP_SEAFLUX XSST_UNIF= 304. /
+ RSOU
+ 2005 11 30 0
+ 'PUVTHDMR'
+ 0.0000000E+00
+ 100300.0
+ 303.3000
+ 1.9630000E-02
+ 92
+ 100000.0 3.637668 -3.637668
+ 99900.00 2.572220 -4.455216
+ 99100.00 -4.612377 0.4035313
+ 96600.00 5.124864 -0.4483674
+ 93000.00 5.836996 -2.009839
+ 92800.00 5.801031 -2.111402
+ 92500.00 5.466062 -1.464627
+ 85000.00 0.1345103 -1.537459
+ 81400.00 -0.9020693 -1.849517
+ 71800.00 -2.797467 1.304481
+ 70000.00 -4.472232 1.198331
+ 60500.00 -3.289776 -1.464702
+ 59900.00 -3.208611 -1.634869
+ 56900.00 -2.329535 -2.025034
+ 55600.00 -2.105100 -2.257443
+ 50100.00 -0.5325915 -1.987659
+ 50000.00 -0.2679967 -1.519885
+ 49500.00 0.2507798 -2.042438
+ 46700.00 5.466062 1.464627
+ 45300.00 5.594935 2.608961
+ 44800.00 5.450725 2.898202
+ 43700.00 4.994327 3.628591
+ 42200.00 4.561049 4.891127
+ 41900.00 4.387573 5.047325
+ 41800.00 4.298817 5.123130
+ 40000.00 6.305391 5.290851
+ 39000.00 5.970908 4.027428
+ 38600.00 5.671552 3.543979
+ 38300.00 5.346259 3.086664
+ 35500.00 2.976104 3.546782
+ 32700.00 4.455215 2.572220
+ 31400.00 3.712562 4.270814
+ 30700.00 3.086664 5.346259
+ 30000.00 3.540875 5.056894
+ 28800.00 3.786532 4.205370
+ 26700.00 3.329551 2.419061
+ 25000.00 2.797467 1.304481
+ 24400.00 2.819808 1.255459
+ 23800.00 2.841289 1.206056
+ 21400.00 3.424857 1.112804
+ 20000.00 3.478403 0.9320353
+ 19300.00 3.383935 1.231651
+ 17400.00 3.411947 2.301387
+ 15000.00 2.950729 4.214078
+ 12700.00 2.167201 2.875973
+ 10200.00 1.322713 1.576348
+ 10000.00 0.4466614 2.533142
+ 9120.000 -3.306782 3.940869
+ 8700.000 -7.950636 2.130367
+ 8200.000 -8.712269 0.7622259
+ 8000.000 -9.224754 -0.8070621
+ 7600.000 -4.009694 -2.314998
+ 7300.000 -5.066284 -0.8933228
+ 7000.000 -4.196201 -1.956721
+ 6700.000 -1.013257 0.1786647
+ 6400.000 -2.107039 1.475365
+ 6340.000 -2.558960 1.726041
+ 5480.000 -13.17234 2.322640
+ 5300.000 -15.37459 1.345105
+ 5000.000 -12.81216 -1.120920
+ 4900.000 -12.66571 -2.233307
+ 4860.000 -12.25463 -1.504677
+ 4700.000 -10.76221 0.9415731
+ 4500.000 -11.37801 7.966970
+ 3800.000 -25.11183 2.197004
+ 3600.000 -22.54940 -1.972818
+ 3200.000 -22.72465 -4.417223
+ 3000.000 -22.36116 -5.991656
+ 2490.000 -21.48829 -2.258509
+ 2300.000 -21.07936 -0.7361050
+ 2200.000 -20.57776 -1.7989648E-06
+ 2080.000 -19.00834 -0.9961855
+ 2000.000 -17.93702 -1.569287
+ 1600.000 -14.74927 -10.32755
+ 1500.000 -13.52109 -6.304989
+ 1430.000 -16.89510 -4.527030
+ 1400.000 -18.23862 -3.215962
+ 1300.000 -15.46941 -5.630407
+ 1200.000 -5.594935 2.608962
+ 1190.000 -5.671552 3.543979
+ 1100.000 1.786646 10.13257
+ 1050.000 -3.736842 7.333968
+ 1020.000 -5.435582 4.725079
+ 1000.000 -5.594935 2.608962
+ 950.0000 -4.347303 8.176087
+ 900.0000 2.143975 12.15908
+ 889.9999 1.031246 11.78719
+ 800.0000 -4.728969 4.728969
+ 760.0000 2.890712 7.154763
+ 710.0000 8.644947 -3.318484
+ 700.0000 7.487652 -6.282886
+ 600.0000 -4.214078 -2.950729
+ 91
+ 100000.0 303.1000 1.8050000E-02
+ 99900.00 303.1000 1.8040001E-02
+ 99100.00 303.1000 1.7969999E-02
+ 96600.00 303.1000 1.7729999E-02
+ 93000.00 303.0000 1.7370002E-02
+ 92800.00 303.1000 1.6980000E-02
+ 92500.00 303.2000 1.6400000E-02
+ 85000.00 307.1000 1.1190000E-02
+ 81400.00 309.6000 6.6700005E-03
+ 71800.00 312.3000 6.5600001E-03
+ 70000.00 312.9000 6.5400004E-03
+ 60500.00 316.2000 6.1800000E-03
+ 59900.00 318.1000 5.3400006E-03
+ 56900.00 318.7000 4.9899998E-03
+ 55600.00 322.2000 3.3100001E-03
+ 50100.00 325.5000 2.6000000E-03
+ 50000.00 325.5000 2.5800001E-03
+ 49500.00 325.7000 3.1700002E-03
+ 46700.00 328.6000 1.4500001E-03
+ 45300.00 330.1000 9.4000006E-04
+ 44800.00 330.1000 1.9600000E-03
+ 43700.00 331.0000 2.1500003E-03
+ 42200.00 332.2000 1.9400001E-03
+ 41900.00 332.9000 9.6000003E-04
+ 41800.00 333.0000 9.7000005E-04
+ 40000.00 334.2000 1.1600000E-03
+ 39000.00 334.6000 1.0300000E-03
+ 38600.00 335.6000 6.5000000E-04
+ 38300.00 335.8000 6.0000003E-04
+ 35500.00 338.2000 2.5000001E-04
+ 32700.00 340.8000 9.0000009E-05
+ 31400.00 342.0000 5.0000002E-05
+ 30700.00 342.1000 4.0000003E-05
+ 30000.00 342.3000 4.0000003E-05
+ 28800.00 342.3000 1.2000000E-04
+ 26700.00 344.8000 1.0000001E-05
+ 25000.00 345.4000 1.0000001E-05
+ 24400.00 346.3000 2.0000001E-05
+ 23800.00 346.4000 9.0000009E-05
+ 21400.00 348.1000 1.0000001E-05
+ 20000.00 348.8000 1.0000001E-05
+ 19300.00 348.9000 2.0000001E-05
+ 17400.00 351.5000 1.0000001E-05
+ 15000.00 354.0000 7.1596442E-06
+ 12700.00 355.7000 4.4376361E-06
+ 10200.00 364.6000 1.4789329E-06
+ 10000.00 365.4000 1.2422363E-06
+ 9120.000 371.6000 1.2919259E-06
+ 8700.000 379.9000 1.5496893E-06
+ 8200.000 390.6000 1.8819877E-06
+ 8000.000 395.1000 2.0217396E-06
+ 7600.000 404.7000 2.3198761E-06
+ 7300.000 412.3000 2.4844724E-06
+ 7000.000 420.4000 2.4844724E-06
+ 6700.000 431.1000 2.4844724E-06
+ 6400.000 442.4000 2.4844724E-06
+ 6340.000 444.8000 2.4844724E-06
+ 5480.000 467.8000 2.4844724E-06
+ 5300.000 471.0000 2.4844724E-06
+ 5000.000 476.5000 2.4844724E-06
+ 4900.000 477.6000 2.4844724E-06
+ 4860.000 478.0000 2.4844724E-06
+ 4700.000 485.1000 2.4844724E-06
+ 4500.000 494.6000 2.4844724E-06
+ 3800.000 532.9000 2.4844724E-06
+ 3600.000 545.6000 2.4844724E-06
+ 3200.000 574.4000 2.4844724E-06
+ 3000.000 585.1000 2.4844724E-06
+ 2490.000 624.0000 2.4844724E-06
+ 2300.000 648.9000 2.4844724E-06
+ 2200.000 655.4000 2.4844724E-06
+ 2080.000 663.6000 2.4844724E-06
+ 2000.000 673.5000 2.4844724E-06
+ 1600.000 726.4000 2.4844724E-06
+ 1500.000 742.5000 2.4844724E-06
+ 1430.000 750.7000 2.4844724E-06
+ 1400.000 759.0000 2.4844724E-06
+ 1300.000 788.6000 2.4844724E-06
+ 1200.000 821.7000 2.4844724E-06
+ 1190.000 825.2000 2.4844724E-06
+ 1100.000 833.9000 2.4844724E-06
+ 1050.000 839.0000 2.4844724E-06
+ 1020.000 845.3000 2.4844724E-06
+ 1000.000 851.6000 2.4844724E-06
+ 950.0000 877.0000 2.4844724E-06
+ 900.0000 885.6000 2.4844724E-06
+ 889.9999 887.3000 2.4844724E-06
+ 800.0000 926.6000 2.4844724E-06
+ 760.0000 946.0000 2.4844724E-06
+ 710.0000 958.0000 2.4844724E-06
+ 700.0000 963.1000 2.4844724E-06
+ 600.0000 20.30000 2.4844724E-06
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_512x512_MNH-55X/002_prep_ideal_case/run_prep_ideal_xyz b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_512x512_MNH-55X/002_prep_ideal_case/run_prep_ideal_xyz
new file mode 100755
index 0000000000000000000000000000000000000000..1a85e665eccd99a35eb1c4498766beba6593b344
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_512x512_MNH-55X/002_prep_ideal_case/run_prep_ideal_xyz
@@ -0,0 +1,20 @@
+#!/bin/bash
+
+NP=${NP:-16}
+NPS=${NPS:-4}
+
+export ACC_DEVICE_TYPE=HOST
+export MPIRUN="Mpirun -tag-output -report-bindings --mca mpi_cuda_support 0 -map-by ppr:${NPS}:socket -bind-to none -np ${NP} set_core_device_impair"
+
+#
+set -x
+set -e
+#
+rm -f OUTPUT_LISTING1 pipe_name file_for_xtransfer
+rm -f PGD00512.*
+rm -f DA0512.*
+#
+ln -sf ../001_pgd1/PGD00512.* .
+#
+time ${MPIRUN} PREP_IDEAL_CASE${XYZ}
+#
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_512x512_MNH-55X/003_mesonh_step1/EXSEG1.nam b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_512x512_MNH-55X/003_mesonh_step1/EXSEG1.nam
new file mode 100644
index 0000000000000000000000000000000000000000..878afa281253ff16017a9c0177e961b3675533fc
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_512x512_MNH-55X/003_mesonh_step1/EXSEG1.nam
@@ -0,0 +1,54 @@
+&NAM_CONFZ
+ ! NB_PROCIO_R=1 ,
+ ! NB_PROCIO_W=1 ,
+ ! LMNH_MPI_BSEND = T ,
+ MPI_BUFFER_SIZE = 1000
+/
+
+&NAM_LUNITn CINIFILE = "DA0512" , CINIFILEPGD = "PGD00512" /
+&NAM_CONFn LUSERV = T, LUSERC = T, LUSERR = T, LUSERI = T,
+ LUSERS = T, LUSERG = T, LUSERH = F, LUSECI = T /
+&NAM_DYNn XTSTEP = 10.0 ,
+ CPRESOPT = "ZRESI", NITR = 12
+ !CPRESOPT = "ZSOLV", NITR = 12
+ XRELAX = 1., LHORELAX_UVWTH = T, LHORELAX_RV = T, LVE_RELAX = T,
+ NRIMX = 6, NRIMY = 6, XRIMKMAX = 0.0005, XT4DIFU = 4000 /
+&NAM_ADVn CUVW_ADV_SCHEME = "WENO_K", CMET_ADV_SCHEME = "PPM_01",
+ CSV_ADV_SCHEME = "PPM_01" /
+&NAM_PARAMn CTURB = "TKEL" ! "NONE" ! "TKEL",
+ CRAD = "ECMW",
+ CCLOUD = "ICE3" ! "NONE" ! "ICE3",
+ CDCONV = "NONE", CSCONV = "NONE" /
+&NAM_PARAM_RADn XDTRAD = 1800., XDTRAD_CLONLY = 1800., LCLEAR_SKY = F,
+ NRAD_COLNBR = 500 /
+&NAM_PARAM_MFSHALLn XIMPL_MF = 1, CMF_UPDRAFT = "EDKF", CMF_CLOUD = "NONE",
+ LMIXUV = T, LMF_FLX = F /
+&NAM_LBCn CLBCX = 2*"OPEN", CLBCY = 2*"OPEN", XCPHASE = 20. /
+&NAM_TURBn XIMPL = 1., CTURBLEN = "DEAR", CTURBDIM = "3DIM",
+ LTURB_FLX = F, LTURB_DIAG = F, CSUBG_AUCV = "NONE", LSIGMAS = F,
+ LSIG_CONV = F, LSUBG_COND = F /
+&NAM_CH_MNHCn /
+&NAM_CONF CCONF = "START", LFLAT = F, CEQNSYS = "DUR",
+ LLG=F, NMODEL = 1, NVERB = 0, CEXP = "D0512",
+ CSEG = "DAR01", CSPLIT = "BSPLITTING"
+ NHALO=1
+ /
+&NAM_DYN XSEGLEN = 14400.0000 ! 100. ! 14400.0000 ,
+ LCORIO = T,
+ LNUMDIFU = F, LNUMDIFTH = F,
+ XALKTOP = 0.001, XALZBOT = 22000. /
+&NAM_BLANK /
+&NAM_NESTING /
+&NAM_BACKUP
+ XBAK_TIME(1,1) = 3600.00000 ! 100. ! 3600.00000
+ XBAK_TIME(1,2) = 7200.00000
+ XBAK_TIME(1,3) = 10800.0000
+ XBAK_TIME(1,4) = 14400.0000
+ XBAK_TIME(1,5) = 18000.0000
+ XBAK_TIME(1,6) = 21600.0000
+ /
+&NAM_ISBAn /
+&NAM_SEAFLUXn /
+&NAM_DIAG_SURFn LSURF_BUDGET=T /
+&NAM_DIAG_ISBAn /
+&NAM_DIAG_SURF_ATMn /
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_512x512_MNH-55X/003_mesonh_step1/parameters_mg.nam b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_512x512_MNH-55X/003_mesonh_step1/parameters_mg.nam
new file mode 100644
index 0000000000000000000000000000000000000000..1a5d92f6f3e36ecd068964fd6fdf580cb1940fe6
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_512x512_MNH-55X/003_mesonh_step1/parameters_mg.nam
@@ -0,0 +1,118 @@
+! *********************************************************************
+! *********************************************************************
+! ** **
+! ** Parameter file for geometric multigrid code **
+! ** **
+! *********************************************************************
+! *********************************************************************
+!
+! *********************************************************************
+! * General parameters
+! *********************************************************************
+¶meters_general
+ savefields = .F. ! Save fields to disk?
+/
+
+! *********************************************************************
+! * General solver parameters
+! *********************************************************************
+¶meters_solver
+ LUseO = F ,
+ LUseT = T ,
+ solvertype = 1, ! Solver type:
+ ! 1 : Richardson iteration
+ ! 2 : Conjugate gradient
+ resreduction = 1e-1 ! 1.0d-15 ! Required relative residual reduction
+ maxiter = 5 ! Maximal number of iterations
+/
+
+! *********************************************************************
+! * Conjugate gradient parameters
+! *********************************************************************
+¶meters_conjugategradient
+ verbose = 10, ! Verbosity level
+ maxiter = 5, ! Maximal number of iterations
+ resreduction = 1.0e-1, ! Target residual reduction
+ n_prec = 1 ! Number of smoother applications in
+ ! preconditioner (N.B.: Using 0 is
+ ! inefficient, as the identity is used
+ ! for preconditioning, instead of using
+ ! unpreconditioned CG.)
+/
+
+! *********************************************************************
+! * Grid parameters
+! *********************************************************************
+¶meters_grid
+ n = 512 , ! 64, ! Number of horizontal grid cells
+ nz = 128, ! Number of vertical grid cells
+ L = 204800.0 ! 400.0 metre * 512 pt grille ! 1.0, ! Size in horizontal direction
+ H = 25600.0 ! 200.0 metre * 128 niveaux ! 0.01, ! Size in vertical direction
+ vertbc = 2, ! Boundary conditions at top and bottom of the
+ ! atmosphere. 1 = DIRICHLET, 2 = NEUMANN
+ ! Note that Neumann boundary conditions only work
+ ! for coarsening in the horizontal only, as they are
+ ! not yet implemented in the prolongation operator.
+ graded = .F. ! Is the grid graded in the vertical direction?
+/
+
+! *********************************************************************
+! * Parallel communication parameters
+! *********************************************************************
+¶meters_communication
+ halo_size = 1 ! Size of halos (has to be 1 or 2)
+/
+
+! *********************************************************************
+! * Model parameters
+! *********************************************************************
+!
+! parameters of the Helmholtz operator
+!
+! -omega2*(d^2/dy^2 + d^2/dy^2 + lambda2*d^2/dz^2) u + delta u = RHS
+!
+¶meters_model
+ omega2 = 1.0,
+ lambda2 = 1.0 ! 100.0, ! Vertical coupling
+ delta = 0.0d0 ! Size of constant term
+/
+
+! *********************************************************************
+! * Smoother parameters
+! *********************************************************************
+!
+! parameters of the smoother
+!
+¶meters_smoother
+ smoother = 3, ! Smoother method
+ ! 3 = line SOR
+ ! 4 = line SSOR
+ ! 6 = line Jacobi
+ ordering = 2, ! Ordering of grid points (for smoother)
+
+ ! 1 = lexicographic
+ ! 2 = red-black ordering
+!rho = 0.6666666666666666d0 ! Overrelaxation parameter
+ rho = 0.8d0 ! Overrelaxation parameter
+/
+
+! *********************************************************************
+! * Multigrid parameters
+! *********************************************************************
+¶meters_multigrid
+ verbose = 10, ! Verbosity level
+ n_lev = 10, ! Number of levels
+ lev_split = 5, ! First level where data is pulled together
+ n_presmooth = 1, ! Number of presmoothing steps
+ n_postsmooth = 1, ! Number of postsmoothing steps
+ n_coarsegridsmooth = 1, ! Number of smoothing steps on coarsest level
+ prolongation = 2 ! 2 best after modif !! ! Prologation method
+ ! 1 = constant interpolation
+ ! 2 = (tri-) linear interpolation
+ restriction = 1, ! Restriction method
+ ! 1 = cell average
+ ! 2 = Khalil
+ coarsegridsolver = 1 ! Solver on coarsest grid
+ ! 1 = use smoother
+ ! 2 = Conjugate gradient
+/
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_512x512_MNH-55X/003_mesonh_step1/run_mesonh_step1_xyz b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_512x512_MNH-55X/003_mesonh_step1/run_mesonh_step1_xyz
new file mode 100755
index 0000000000000000000000000000000000000000..1a3cfc09ce327c8c20cc728df069ac08e720c746
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_512x512_MNH-55X/003_mesonh_step1/run_mesonh_step1_xyz
@@ -0,0 +1,23 @@
+
+
+#export PGI_ACC_POOL_ALLOC=0
+export PGI_ACC_SYNCHRONOUS=1
+
+NP=${NP:-16}
+NPS=${NPS:-4}
+
+export MPIRUN=${MPIRUN:-"Mpirun -tag-output -report-bindings -bind-to none -map-by ppr:${NPS}:socket -x PGI_ACC_POOL_ALLOC -x PGI_ACC_SYNCHRONOUS -np ${NP} set_core_device_impair "}
+
+
+#
+set -x
+set -e
+#
+rm -f DA0512.*
+rm -f PGD00512.*
+rm -f D0512.1.DAR01.* OUTPUT_LISTING* pipe_name
+#
+ln -sf ../001_pgd1/PGD00512.{des,nc} .
+ln -sf ../002_prep_ideal_case/DA0512.{des,nc} .
+#
+time ${MPIRUN} MESONH${XYZ}
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_512x512_MNH-55X/004_mesonh_step2/EXSEG1.nam b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_512x512_MNH-55X/004_mesonh_step2/EXSEG1.nam
new file mode 100644
index 0000000000000000000000000000000000000000..bf779d1d5b84596c64a86daeb5a56ac1c8915a0b
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_512x512_MNH-55X/004_mesonh_step2/EXSEG1.nam
@@ -0,0 +1,59 @@
+&NAM_CONFIO
+ LIO_ALLOW_NO_BACKUP = T , LIO_NO_WRITE = T
+/
+&NAM_CONFZ
+ ! NB_PROCIO_R=1 ,
+ ! NB_PROCIO_W=1 ,
+ ! LMNH_MPI_BSEND = T ,
+ MPI_BUFFER_SIZE = 1000
+/
+
+&NAM_LUNITn CINIFILE = "D0512.1.DAR01.004" ! "D0512.1.DAR01.003" ,
+ CINIFILEPGD = "PGD00512" /
+&NAM_CONFn LUSERV = T, LUSERC = T, LUSERR = T, LUSERI = T,
+ LUSERS = T, LUSERG = T, LUSERH = F, LUSECI = T /
+&NAM_DYNn XTSTEP = 10.00000 ,
+ !CPRESOPT = "ZRESI",
+ CPRESOPT = "ZSOLV",
+ NITR = 12,
+ XRELAX = 1., LHORELAX_UVWTH = T, LHORELAX_RV = T, LVE_RELAX = T,
+ NRIMX = 6, NRIMY = 6, XRIMKMAX = 0.0005, XT4DIFU = 4000 /
+&NAM_ADVn CUVW_ADV_SCHEME = "WENO_K", CMET_ADV_SCHEME = "PPM_01",
+ CSV_ADV_SCHEME = "PPM_01" /
+&NAM_PARAMn CTURB = "TKEL", ! "NONE" , ! "TKEL",
+ CRAD = "ECMWF",
+ CCLOUD = "ICE3", ! "NONE" , ! "ICE3",
+ CDCONV = "NONE", CSCONV = "NONE" /
+&NAM_PARAM_RADn XDTRAD = 1800., XDTRAD_CLONLY = 1800., LCLEAR_SKY = F,
+ NRAD_COLNBR = 500 /
+&NAM_PARAM_MFSHALLn XIMPL_MF = 1, CMF_UPDRAFT = "EDKF", CMF_CLOUD = "NONE",
+ LMIXUV = T, LMF_FLX = F /
+&NAM_LBCn CLBCX = 2*"OPEN", CLBCY = 2*"OPEN", XCPHASE = 20. /
+&NAM_TURBn XIMPL = 1., CTURBLEN = "DEAR", CTURBDIM = "3DIM",
+ LTURB_FLX = F, LTURB_DIAG = F, CSUBG_AUCV = "NONE", LSIGMAS = F,
+ LSIG_CONV = F, LSUBG_COND = F /
+&NAM_CH_MNHCn /
+&NAM_CONF CCONF = "RESTA", LFLAT = F, CEQNSYS = "DUR",
+ LLG=F, NMODEL = 1, NVERB = 0, CEXP = "D0512",
+ CSEG = "DAR02", CSPLIT = "BSPLITTING"
+ NHALO=1
+ LCHECK = F /
+&NAM_DYN XSEGLEN = 100.0 ! 150.0000 ,
+ LCORIO = T,
+ LNUMDIFU = F, LNUMDIFTH = F,
+ XALKTOP = 0.001, XALZBOT = 22000. /
+&NAM_BLANK /
+&NAM_NESTING /
+&NAM_BACKUP
+ XBAK_TIME(1,1) = 100.0 ! 3600.00000 ! 3600.00000
+ XBAK_TIME(1,2) = 7200.00000
+ XBAK_TIME(1,3) = 10800.0000
+ XBAK_TIME(1,4) = 14400.0000
+ XBAK_TIME(1,5) = 18000.0000
+ XBAK_TIME(1,6) = 21600.0000
+ /
+&NAM_ISBAn /
+&NAM_SEAFLUXn /
+&NAM_DIAG_SURFn LSURF_BUDGET=T /
+&NAM_DIAG_ISBAn /
+&NAM_DIAG_SURF_ATMn /
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_512x512_MNH-55X/004_mesonh_step2/parameters_mg.nam b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_512x512_MNH-55X/004_mesonh_step2/parameters_mg.nam
new file mode 100644
index 0000000000000000000000000000000000000000..1a5d92f6f3e36ecd068964fd6fdf580cb1940fe6
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_512x512_MNH-55X/004_mesonh_step2/parameters_mg.nam
@@ -0,0 +1,118 @@
+! *********************************************************************
+! *********************************************************************
+! ** **
+! ** Parameter file for geometric multigrid code **
+! ** **
+! *********************************************************************
+! *********************************************************************
+!
+! *********************************************************************
+! * General parameters
+! *********************************************************************
+¶meters_general
+ savefields = .F. ! Save fields to disk?
+/
+
+! *********************************************************************
+! * General solver parameters
+! *********************************************************************
+¶meters_solver
+ LUseO = F ,
+ LUseT = T ,
+ solvertype = 1, ! Solver type:
+ ! 1 : Richardson iteration
+ ! 2 : Conjugate gradient
+ resreduction = 1e-1 ! 1.0d-15 ! Required relative residual reduction
+ maxiter = 5 ! Maximal number of iterations
+/
+
+! *********************************************************************
+! * Conjugate gradient parameters
+! *********************************************************************
+¶meters_conjugategradient
+ verbose = 10, ! Verbosity level
+ maxiter = 5, ! Maximal number of iterations
+ resreduction = 1.0e-1, ! Target residual reduction
+ n_prec = 1 ! Number of smoother applications in
+ ! preconditioner (N.B.: Using 0 is
+ ! inefficient, as the identity is used
+ ! for preconditioning, instead of using
+ ! unpreconditioned CG.)
+/
+
+! *********************************************************************
+! * Grid parameters
+! *********************************************************************
+¶meters_grid
+ n = 512 , ! 64, ! Number of horizontal grid cells
+ nz = 128, ! Number of vertical grid cells
+ L = 204800.0 ! 400.0 metre * 512 pt grille ! 1.0, ! Size in horizontal direction
+ H = 25600.0 ! 200.0 metre * 128 niveaux ! 0.01, ! Size in vertical direction
+ vertbc = 2, ! Boundary conditions at top and bottom of the
+ ! atmosphere. 1 = DIRICHLET, 2 = NEUMANN
+ ! Note that Neumann boundary conditions only work
+ ! for coarsening in the horizontal only, as they are
+ ! not yet implemented in the prolongation operator.
+ graded = .F. ! Is the grid graded in the vertical direction?
+/
+
+! *********************************************************************
+! * Parallel communication parameters
+! *********************************************************************
+¶meters_communication
+ halo_size = 1 ! Size of halos (has to be 1 or 2)
+/
+
+! *********************************************************************
+! * Model parameters
+! *********************************************************************
+!
+! parameters of the Helmholtz operator
+!
+! -omega2*(d^2/dy^2 + d^2/dy^2 + lambda2*d^2/dz^2) u + delta u = RHS
+!
+¶meters_model
+ omega2 = 1.0,
+ lambda2 = 1.0 ! 100.0, ! Vertical coupling
+ delta = 0.0d0 ! Size of constant term
+/
+
+! *********************************************************************
+! * Smoother parameters
+! *********************************************************************
+!
+! parameters of the smoother
+!
+¶meters_smoother
+ smoother = 3, ! Smoother method
+ ! 3 = line SOR
+ ! 4 = line SSOR
+ ! 6 = line Jacobi
+ ordering = 2, ! Ordering of grid points (for smoother)
+
+ ! 1 = lexicographic
+ ! 2 = red-black ordering
+!rho = 0.6666666666666666d0 ! Overrelaxation parameter
+ rho = 0.8d0 ! Overrelaxation parameter
+/
+
+! *********************************************************************
+! * Multigrid parameters
+! *********************************************************************
+¶meters_multigrid
+ verbose = 10, ! Verbosity level
+ n_lev = 10, ! Number of levels
+ lev_split = 5, ! First level where data is pulled together
+ n_presmooth = 1, ! Number of presmoothing steps
+ n_postsmooth = 1, ! Number of postsmoothing steps
+ n_coarsegridsmooth = 1, ! Number of smoothing steps on coarsest level
+ prolongation = 2 ! 2 best after modif !! ! Prologation method
+ ! 1 = constant interpolation
+ ! 2 = (tri-) linear interpolation
+ restriction = 1, ! Restriction method
+ ! 1 = cell average
+ ! 2 = Khalil
+ coarsegridsolver = 1 ! Solver on coarsest grid
+ ! 1 = use smoother
+ ! 2 = Conjugate gradient
+/
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_512x512_MNH-55X/004_mesonh_step2/run_mesonh_step2_xyz b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_512x512_MNH-55X/004_mesonh_step2/run_mesonh_step2_xyz
new file mode 100755
index 0000000000000000000000000000000000000000..c10adb44e1c80d1049209cbdf0b33c94c6f17127
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_512x512_MNH-55X/004_mesonh_step2/run_mesonh_step2_xyz
@@ -0,0 +1,26 @@
+#!/bin/bash
+#
+#
+set -x
+
+#export PGI_ACC_POOL_ALLOC=0
+export PGI_ACC_SYNCHRONOUS=1
+#export NVCOMPILER_ACC_POOL_ALLOC=0
+
+NP=${NP:-16}
+NPS=${NPS:-4}
+
+export MPIRUN=${MPIRUN:-"Mpirun -tag-output -report-bindings -bind-to none -map-by ppr:${NPS}:socket -x NVCOMPILER_ACC_POOL_ALLOC_MAXSIZE -x NVCOMPILER_ACC_POOL_ALLOC -x PGI_ACC_POOL_ALLOC -x PGI_ACC_SYNCHRONOUS -np ${NP} set_core_device_impair "}
+
+rm -f DA0512.*
+rm -f PGD00512.*
+rm -f D0512.1.DAR0?.* OUTPUT_LISTING1 pipe_name
+#
+ln -sf ../001_pgd1/PGD00512.* .
+ln -sf ../002_prep_ideal_case/DA0512.* .
+ln -sf ../003_mesonh_step1/D0512.1.DAR0?.* .
+#
+
+time ${MPIRUN} MESONH${XYZ}
+
+
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_512x512_MNH-55X/005_mg_main_mnh_all/parameters_mg.nam b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_512x512_MNH-55X/005_mg_main_mnh_all/parameters_mg.nam
new file mode 120000
index 0000000000000000000000000000000000000000..adeab0e1859deb44804ccc14f23c2ccb35775079
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_512x512_MNH-55X/005_mg_main_mnh_all/parameters_mg.nam
@@ -0,0 +1 @@
+../004_mesonh_step2/parameters_mg.nam
\ No newline at end of file
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_512x512_MNH-55X/005_mg_main_mnh_all/run_mg_main_mnh_all_xyz b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_512x512_MNH-55X/005_mg_main_mnh_all/run_mg_main_mnh_all_xyz
new file mode 100755
index 0000000000000000000000000000000000000000..c3694f2de5a11d8600af4cf5874b3248c26057f4
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_512x512_MNH-55X/005_mg_main_mnh_all/run_mg_main_mnh_all_xyz
@@ -0,0 +1,18 @@
+#!/bin/bash
+#
+#
+set -x
+
+#export PGI_ACC_POOL_ALLOC=0
+export PGI_ACC_SYNCHRONOUS=1
+#export NVCOMPILER_ACC_POOL_ALLOC=0
+
+NP=${NP:-16}
+NPS=${NPS:-4}
+
+export MPIRUN=${MPIRUN:-"Mpirun -tag-output -report-bindings -bind-to none -map-by ppr:${NPS}:socket -x NVCOMPILER_ACC_POOL_ALLOC_MAXSIZE -x NVCOMPILER_ACC_POOL_ALLOC -x PGI_ACC_POOL_ALLOC -x PGI_ACC_SYNCHRONOUS -np ${NP} set_core_device_impair "}
+
+
+time ${MPIRUN} MG_MAIN_MNH_ALL${XYZ}
+
+
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/001_pgd1/PRE_PGD1.nam b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/001_pgd1/PRE_PGD1.nam
new file mode 100644
index 0000000000000000000000000000000000000000..98c1fb75ddf67ed2c004f85e38e92503dcaf6303
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/001_pgd1/PRE_PGD1.nam
@@ -0,0 +1,13 @@
+&NAM_CONFZ
+ NB_PROCIO_W=1 ,
+ !MPI_BUFFER_SIZE=40
+/
+&NAM_PGDFILE CPGDFILE="PGD00064" /
+&NAM_PGD_SCHEMES CNATURE='ISBA', CSEA='SEAFLX', CTOWN='NONE', CWATER='NONE' /
+&NAM_CONF_PROJ XLAT0=-11.5, XLON0=130.7, XRPK=0., XBETA=0. /
+&NAM_CONF_PROJ_GRID XLATCEN=-11.5, XLONCEN=130.7, NIMAX=64, NJMAX=64,
+ XDX=3200.00000, XDY=3200.00000 /
+&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' /
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/001_pgd1/get_pgd_files b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/001_pgd1/get_pgd_files
new file mode 100755
index 0000000000000000000000000000000000000000..56726933a0f7608e0c4918f42e80d5aea76b669d
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/001_pgd1/get_pgd_files
@@ -0,0 +1,30 @@
+#
+# Modif
+# J.Escobar 11/04/2014 get PGD files from 'dir_open' directory ( without psswd )
+# J.Escobar 25/04/2013 get LICENCE files
+#
+PREP_PGD_FILES=${PREP_PGD_FILES:-"$HOME/PREP_PGD_FILES_WWW"} ; export PREP_PGD_FILES
+echo "!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!"
+echo "!!!! WARNING !!!!"
+echo "!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!"
+echo
+echo you need 3GO to download this files in
+echo
+echo PREP_PGD_FILES=$PREP_PGD_FILES
+echo
+echo if OK press ENTER else CTRL-C
+read RIEN
+set -x
+mkdir -p $PREP_PGD_FILES
+cd $PREP_PGD_FILES
+PGD_URL="http://mesonh.aero.obs-mip.fr/mesonh/dir_open/dir_PGDFILES"
+WGET="wget"
+export PGD_URL
+for file in LICENSE_ECOCLIMAP.txt LICENSE_soil_data.txt \
+ gtopo30.dir gtopo30.hdr \
+ SAND_HWSD_MOY.hdr SAND_HWSD_MOY.dir CLAY_HWSD_MOY.hdr CLAY_HWSD_MOY.dir \
+ ECOCLIMAP_v2.0.hdr ECOCLIMAP_v2.0.dir
+do
+[ -f $file ] || ( ${WGET} -c -nd $PGD_URL/$file.gz ; gunzip $file.gz ; )
+done
+
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/001_pgd1/run_prep_pgd_xyz b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/001_pgd1/run_prep_pgd_xyz
new file mode 100755
index 0000000000000000000000000000000000000000..827253ef48a7fd218717d6e6305133e583f428a6
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/001_pgd1/run_prep_pgd_xyz
@@ -0,0 +1,46 @@
+#!/bin/bash
+#BSUB -a openmpi
+#BSUB -J Hector_80x64
+#BSUB -n 4
+#BSUB -q normal
+#BSUB -R "span[ptile=4]"
+#BSUB -o Sortie_Hector_80x64%J
+#BSUB -W 01:00
+#
+
+#
+#. /linkhome/rech/yiu/ryiu001/MNH-V5-1-4/conf/profile_mesonh-LXgfortranI4-MNH-V5-1-4-OMPI1-10-MPIAUTO-O2
+
+export MPIRUN="Mpirun -np 1"
+#export MPIRUN="Exec totalview"
+
+PREP_PGD_FILES=${PREP_PGD_FILES:-"$HOME/PREP_PGD_FILES_WWW"} ; export PREP_PGD_FILES
+#
+if [ ! -d $PREP_PGD_FILES ]
+then
+cat << EOF
+
+Your directory PREP_PGD_FILES=$PREP_PGD_FILES
+
+containing the files gtopo30*, ECOCLIMAP_v2.0*
+doesn't exist ( or was not found !!! )
+use the script 'get_pgd_files' to download
+this files from the MESONH WEB PAGES !!!
+( or change the variable PREP_PGD_FILES ... )
+
+After fixing it , run this script again !!!
+
+EOF
+exit 1
+else
+set -x
+#
+rm -f OUTPUT_LISTING0 pipe_name
+rm -f gtopo30.??? sand_fao.??? clay_fao.???
+rm -f SAND_HWSD_MOY.??? ECOCLIMAP_v2.0.??? ecoclimats_v2.???
+rm -f PGD00064.*
+#
+ln -sf $PREP_PGD_FILES/*.dir $PREP_PGD_FILES/*.hdr .
+#
+time ${MPIRUN} PREP_PGD${XYZ}
+fi
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/002_prep_ideal_case/PRE_IDEA1.nam b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/002_prep_ideal_case/PRE_IDEA1.nam
new file mode 100644
index 0000000000000000000000000000000000000000..d5542562c0c89612f6eda02f5a35001fb7924fe3
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/002_prep_ideal_case/PRE_IDEA1.nam
@@ -0,0 +1,226 @@
+&NAM_CONFZ
+ !NB_PROCIO_R=1 ,
+ !NB_PROCIO_W=1 ,
+ !MPI_BUFFER_SIZE=40
+/
+
+&NAM_REAL_PGD
+ CPGD_FILE ="PGD00064" ,
+ LREAD_ZS =.TRUE., LREAD_GROUND_PARAM =.TRUE.
+/
+&NAM_DIMn_PRE NIMAX=40, NJMAX=32 /
+&NAM_CONF_PRE LCARTESIAN=.FALSE.
+ CIDEAL='RSOU' CZS='FLAT'
+ NVERB=0 /
+&NAM_CONFn LUSERV= T /
+&NAM_LUNITn CINIFILE = "DA0064" , CINIFILEPGD = "PGD00064" /
+&NAM_DYNn_PRE
+ CPRESOPT= 'ZRESI' ,
+ NITR=4 XRELAX=1. /
+&NAM_LBCn_PRE CLBCX= 2*'OPEN' CLBCY= 2*'OPEN' /
+&NAM_VPROF_PRE CTYPELOC='IJGRID' NILOC=1 NJLOC=1
+ CFUNU='ZZZ' CFUNV='ZZZ'
+ LGEOSBAL=.FALSE. /
+&NAM_VER_GRID NKMAX=126, YZGRID_TYPE='FUNCTN',
+ ZDZGRD=40., ZDZTOP=210., ZZMAX_STRGRD=2500.,
+ ZSTRGRD=7., ZSTRTOP=7. /
+&NAM_GRn_PRE
+ CSURF='EXTE'
+ ! CSURF='NONE'
+ /
+&NAM_PREP_ISBA XTG_SURF= 311., XTG_ROOT= 303., XTG_DEEP= 302.,
+ XHUG_SURF= 0.16, XHUG_ROOT= 0.16, XHUG_DEEP= 0.16 /
+&NAM_PREP_SEAFLUX XSST_UNIF= 304. /
+ RSOU
+ 2005 11 30 0
+ 'PUVTHDMR'
+ 0.0000000E+00
+ 100300.0
+ 303.3000
+ 1.9630000E-02
+ 92
+ 100000.0 3.637668 -3.637668
+ 99900.00 2.572220 -4.455216
+ 99100.00 -4.612377 0.4035313
+ 96600.00 5.124864 -0.4483674
+ 93000.00 5.836996 -2.009839
+ 92800.00 5.801031 -2.111402
+ 92500.00 5.466062 -1.464627
+ 85000.00 0.1345103 -1.537459
+ 81400.00 -0.9020693 -1.849517
+ 71800.00 -2.797467 1.304481
+ 70000.00 -4.472232 1.198331
+ 60500.00 -3.289776 -1.464702
+ 59900.00 -3.208611 -1.634869
+ 56900.00 -2.329535 -2.025034
+ 55600.00 -2.105100 -2.257443
+ 50100.00 -0.5325915 -1.987659
+ 50000.00 -0.2679967 -1.519885
+ 49500.00 0.2507798 -2.042438
+ 46700.00 5.466062 1.464627
+ 45300.00 5.594935 2.608961
+ 44800.00 5.450725 2.898202
+ 43700.00 4.994327 3.628591
+ 42200.00 4.561049 4.891127
+ 41900.00 4.387573 5.047325
+ 41800.00 4.298817 5.123130
+ 40000.00 6.305391 5.290851
+ 39000.00 5.970908 4.027428
+ 38600.00 5.671552 3.543979
+ 38300.00 5.346259 3.086664
+ 35500.00 2.976104 3.546782
+ 32700.00 4.455215 2.572220
+ 31400.00 3.712562 4.270814
+ 30700.00 3.086664 5.346259
+ 30000.00 3.540875 5.056894
+ 28800.00 3.786532 4.205370
+ 26700.00 3.329551 2.419061
+ 25000.00 2.797467 1.304481
+ 24400.00 2.819808 1.255459
+ 23800.00 2.841289 1.206056
+ 21400.00 3.424857 1.112804
+ 20000.00 3.478403 0.9320353
+ 19300.00 3.383935 1.231651
+ 17400.00 3.411947 2.301387
+ 15000.00 2.950729 4.214078
+ 12700.00 2.167201 2.875973
+ 10200.00 1.322713 1.576348
+ 10000.00 0.4466614 2.533142
+ 9120.000 -3.306782 3.940869
+ 8700.000 -7.950636 2.130367
+ 8200.000 -8.712269 0.7622259
+ 8000.000 -9.224754 -0.8070621
+ 7600.000 -4.009694 -2.314998
+ 7300.000 -5.066284 -0.8933228
+ 7000.000 -4.196201 -1.956721
+ 6700.000 -1.013257 0.1786647
+ 6400.000 -2.107039 1.475365
+ 6340.000 -2.558960 1.726041
+ 5480.000 -13.17234 2.322640
+ 5300.000 -15.37459 1.345105
+ 5000.000 -12.81216 -1.120920
+ 4900.000 -12.66571 -2.233307
+ 4860.000 -12.25463 -1.504677
+ 4700.000 -10.76221 0.9415731
+ 4500.000 -11.37801 7.966970
+ 3800.000 -25.11183 2.197004
+ 3600.000 -22.54940 -1.972818
+ 3200.000 -22.72465 -4.417223
+ 3000.000 -22.36116 -5.991656
+ 2490.000 -21.48829 -2.258509
+ 2300.000 -21.07936 -0.7361050
+ 2200.000 -20.57776 -1.7989648E-06
+ 2080.000 -19.00834 -0.9961855
+ 2000.000 -17.93702 -1.569287
+ 1600.000 -14.74927 -10.32755
+ 1500.000 -13.52109 -6.304989
+ 1430.000 -16.89510 -4.527030
+ 1400.000 -18.23862 -3.215962
+ 1300.000 -15.46941 -5.630407
+ 1200.000 -5.594935 2.608962
+ 1190.000 -5.671552 3.543979
+ 1100.000 1.786646 10.13257
+ 1050.000 -3.736842 7.333968
+ 1020.000 -5.435582 4.725079
+ 1000.000 -5.594935 2.608962
+ 950.0000 -4.347303 8.176087
+ 900.0000 2.143975 12.15908
+ 889.9999 1.031246 11.78719
+ 800.0000 -4.728969 4.728969
+ 760.0000 2.890712 7.154763
+ 710.0000 8.644947 -3.318484
+ 700.0000 7.487652 -6.282886
+ 600.0000 -4.214078 -2.950729
+ 91
+ 100000.0 303.1000 1.8050000E-02
+ 99900.00 303.1000 1.8040001E-02
+ 99100.00 303.1000 1.7969999E-02
+ 96600.00 303.1000 1.7729999E-02
+ 93000.00 303.0000 1.7370002E-02
+ 92800.00 303.1000 1.6980000E-02
+ 92500.00 303.2000 1.6400000E-02
+ 85000.00 307.1000 1.1190000E-02
+ 81400.00 309.6000 6.6700005E-03
+ 71800.00 312.3000 6.5600001E-03
+ 70000.00 312.9000 6.5400004E-03
+ 60500.00 316.2000 6.1800000E-03
+ 59900.00 318.1000 5.3400006E-03
+ 56900.00 318.7000 4.9899998E-03
+ 55600.00 322.2000 3.3100001E-03
+ 50100.00 325.5000 2.6000000E-03
+ 50000.00 325.5000 2.5800001E-03
+ 49500.00 325.7000 3.1700002E-03
+ 46700.00 328.6000 1.4500001E-03
+ 45300.00 330.1000 9.4000006E-04
+ 44800.00 330.1000 1.9600000E-03
+ 43700.00 331.0000 2.1500003E-03
+ 42200.00 332.2000 1.9400001E-03
+ 41900.00 332.9000 9.6000003E-04
+ 41800.00 333.0000 9.7000005E-04
+ 40000.00 334.2000 1.1600000E-03
+ 39000.00 334.6000 1.0300000E-03
+ 38600.00 335.6000 6.5000000E-04
+ 38300.00 335.8000 6.0000003E-04
+ 35500.00 338.2000 2.5000001E-04
+ 32700.00 340.8000 9.0000009E-05
+ 31400.00 342.0000 5.0000002E-05
+ 30700.00 342.1000 4.0000003E-05
+ 30000.00 342.3000 4.0000003E-05
+ 28800.00 342.3000 1.2000000E-04
+ 26700.00 344.8000 1.0000001E-05
+ 25000.00 345.4000 1.0000001E-05
+ 24400.00 346.3000 2.0000001E-05
+ 23800.00 346.4000 9.0000009E-05
+ 21400.00 348.1000 1.0000001E-05
+ 20000.00 348.8000 1.0000001E-05
+ 19300.00 348.9000 2.0000001E-05
+ 17400.00 351.5000 1.0000001E-05
+ 15000.00 354.0000 7.1596442E-06
+ 12700.00 355.7000 4.4376361E-06
+ 10200.00 364.6000 1.4789329E-06
+ 10000.00 365.4000 1.2422363E-06
+ 9120.000 371.6000 1.2919259E-06
+ 8700.000 379.9000 1.5496893E-06
+ 8200.000 390.6000 1.8819877E-06
+ 8000.000 395.1000 2.0217396E-06
+ 7600.000 404.7000 2.3198761E-06
+ 7300.000 412.3000 2.4844724E-06
+ 7000.000 420.4000 2.4844724E-06
+ 6700.000 431.1000 2.4844724E-06
+ 6400.000 442.4000 2.4844724E-06
+ 6340.000 444.8000 2.4844724E-06
+ 5480.000 467.8000 2.4844724E-06
+ 5300.000 471.0000 2.4844724E-06
+ 5000.000 476.5000 2.4844724E-06
+ 4900.000 477.6000 2.4844724E-06
+ 4860.000 478.0000 2.4844724E-06
+ 4700.000 485.1000 2.4844724E-06
+ 4500.000 494.6000 2.4844724E-06
+ 3800.000 532.9000 2.4844724E-06
+ 3600.000 545.6000 2.4844724E-06
+ 3200.000 574.4000 2.4844724E-06
+ 3000.000 585.1000 2.4844724E-06
+ 2490.000 624.0000 2.4844724E-06
+ 2300.000 648.9000 2.4844724E-06
+ 2200.000 655.4000 2.4844724E-06
+ 2080.000 663.6000 2.4844724E-06
+ 2000.000 673.5000 2.4844724E-06
+ 1600.000 726.4000 2.4844724E-06
+ 1500.000 742.5000 2.4844724E-06
+ 1430.000 750.7000 2.4844724E-06
+ 1400.000 759.0000 2.4844724E-06
+ 1300.000 788.6000 2.4844724E-06
+ 1200.000 821.7000 2.4844724E-06
+ 1190.000 825.2000 2.4844724E-06
+ 1100.000 833.9000 2.4844724E-06
+ 1050.000 839.0000 2.4844724E-06
+ 1020.000 845.3000 2.4844724E-06
+ 1000.000 851.6000 2.4844724E-06
+ 950.0000 877.0000 2.4844724E-06
+ 900.0000 885.6000 2.4844724E-06
+ 889.9999 887.3000 2.4844724E-06
+ 800.0000 926.6000 2.4844724E-06
+ 760.0000 946.0000 2.4844724E-06
+ 710.0000 958.0000 2.4844724E-06
+ 700.0000 963.1000 2.4844724E-06
+ 600.0000 20.30000 2.4844724E-06
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/002_prep_ideal_case/run_prep_ideal_xyz b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/002_prep_ideal_case/run_prep_ideal_xyz
new file mode 100755
index 0000000000000000000000000000000000000000..2cdc37f2a30123c4eb10d6e7fe0c9991ad1dc6c4
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/002_prep_ideal_case/run_prep_ideal_xyz
@@ -0,0 +1,20 @@
+#!/bin/bash
+
+NP=${NP:-16}
+NPS=${NPS:-4}
+
+export ACC_DEVICE_TYPE=HOST
+export MPIRUN="Mpirun -tag-output -report-bindings --mca mpi_cuda_support 0 -map-by ppr:${NPS}:socket -bind-to none -np ${NP} set_core_device_impair"
+
+#
+set -x
+set -e
+#
+rm -f OUTPUT_LISTING1 pipe_name file_for_xtransfer
+rm -f PGD00064.*
+rm -f DA0064.*
+#
+ln -sf ../001_pgd1/PGD00064.* .
+#
+time ${MPIRUN} PREP_IDEAL_CASE${XYZ}
+#
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/003_mesonh_step1/EXSEG1.nam b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/003_mesonh_step1/EXSEG1.nam
new file mode 100644
index 0000000000000000000000000000000000000000..0898b916f1f8f098f34bdeffa2a342f7296dcd8e
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/003_mesonh_step1/EXSEG1.nam
@@ -0,0 +1,55 @@
+&NAM_CONFZ
+ ! NB_PROCIO_R=1 ,
+ ! NB_PROCIO_W=1 ,
+ ! LMNH_MPI_BSEND = T ,
+ ! MPI_BUFFER_SIZE = 40
+/
+
+&NAM_LUNITn CINIFILE = "DA0064" , CINIFILEPGD = "PGD00064" /
+&NAM_CONFn LUSERV = T, LUSERC = T, LUSERR = T, LUSERI = T,
+ LUSERS = T, LUSERG = T, LUSERH = F, LUSECI = T /
+&NAM_DYNn XTSTEP = 15.00000 ,
+ CPRESOPT = "ZRESI", NITR = 12
+ !CPRESOPT = "ZSOLV", NITR = 12
+ XRELAX = 1., LHORELAX_UVWTH = T, LHORELAX_RV = T, LVE_RELAX = T,
+ NRIMX = 6, NRIMY = 6, XRIMKMAX = 0.0005, XT4DIFU = 4000 /
+&NAM_ADVn CUVW_ADV_SCHEME = "WENO_K", CMET_ADV_SCHEME = "PPM_01",
+ CSV_ADV_SCHEME = "PPM_01" /
+&NAM_PARAMn CTURB = "TKEL" ! "NONE" ! "TKEL",
+ CRAD = "ECMW" ! "NONE " ! "FIXE" ! "ECMW",
+ CCLOUD = "ICE3" ! "NONE" ! "ICE3",
+ CDCONV = "NONE", CSCONV = "NONE" /
+&NAM_PARAM_RADn XDTRAD = 1800., XDTRAD_CLONLY = 1800., LCLEAR_SKY = F,
+ NRAD_COLNBR = 500 /
+&NAM_PARAM_MFSHALLn XIMPL_MF = 1, CMF_UPDRAFT = "EDKF", CMF_CLOUD = "NONE",
+ LMIXUV = T, LMF_FLX = F /
+&NAM_LBCn CLBCX = 2*"OPEN", CLBCY = 2*"OPEN", XCPHASE = 20. /
+&NAM_TURBn XIMPL = 1., CTURBLEN = "DEAR", CTURBDIM = "3DIM",
+ LTURB_FLX = F, LTURB_DIAG = F, CSUBG_AUCV = "NONE", LSIGMAS = F,
+ LSIG_CONV = F, LSUBG_COND = F /
+&NAM_CH_MNHCn /
+&NAM_CONF CCONF = "START", LFLAT = F, CEQNSYS = "DUR",
+ LLG=F, NMODEL = 1, NVERB = 0, CEXP = "D0064",
+ CSEG = "DAR01", CSPLIT = "BSPLITTING"
+ NHALO=1
+ LCHECK=F
+ /
+&NAM_DYN XSEGLEN = 14400.0000 ! 150. ! 21600.0000
+ LCORIO = T,
+ LNUMDIFU = F, LNUMDIFTH = F,
+ XALKTOP = 0.001, XALZBOT = 22000. /
+&NAM_BLANK /
+&NAM_NESTING /
+&NAM_BACKUP
+ XBAK_TIME(1,1) = 3600.00000 ! 150.
+ XBAK_TIME(1,2) = 7200.00000
+ XBAK_TIME(1,3) = 10800.0000
+ XBAK_TIME(1,4) = 14400.0000
+ XBAK_TIME(1,5) = 18000.0000
+ XBAK_TIME(1,6) = 21600.0000
+ /
+&NAM_ISBAn /
+&NAM_SEAFLUXn /
+&NAM_DIAG_SURFn LSURF_BUDGET=T /
+&NAM_DIAG_ISBAn /
+&NAM_DIAG_SURF_ATMn /
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/003_mesonh_step1/MESONH_HOST2005 b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/003_mesonh_step1/MESONH_HOST2005
new file mode 100755
index 0000000000000000000000000000000000000000..f81d930cef2e920eca3ab0fa0b8c3f4a18ce5cb3
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/003_mesonh_step1/MESONH_HOST2005
@@ -0,0 +1,5 @@
+#!/bin/bash
+export ACC_DEVICE=HOST
+ACC_DEVICE_TYPE=${ACC_DEVICE}
+exec /home/escj/DEV/MNH-55X-dev-OPENACC-juan-03/exe/MESONH${XYZ}
+
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/003_mesonh_step1/mppdb.nam b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/003_mesonh_step1/mppdb.nam
new file mode 100644
index 0000000000000000000000000000000000000000..3c7a9f6aa3a9a8d6f7cbcee9dfc444d31a62d63e
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/003_mesonh_step1/mppdb.nam
@@ -0,0 +1,8 @@
+&NAM_MPPDB
+MPPDB_DEBUG = .F. ,
+MPPDB_NBSON = 1 ,
+!MPPDB_HOST = "localhost" ,
+MPPDB_EXEC = "/home/escj/DEV/MNH-55X-dev-OPENACC-juan-03/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/003_mesonh_step1/MESONH_HOST2005" ,
+MPPDB_WDIR = "/home/escj/DEV/MNH-55X-dev-OPENACC-juan-03/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/003_mesonh_step1/dir_clone" ,
+MPPDB_CHECK_LB = .F.
+/
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/003_mesonh_step1/parameters_mg.nam b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/003_mesonh_step1/parameters_mg.nam
new file mode 100644
index 0000000000000000000000000000000000000000..da6e1ac67e7e9cd235de8d628379403ef298c3ad
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/003_mesonh_step1/parameters_mg.nam
@@ -0,0 +1,117 @@
+! *********************************************************************
+! *********************************************************************
+! ** **
+! ** Parameter file for geometric multigrid code **
+! ** **
+! *********************************************************************
+! *********************************************************************
+!
+! *********************************************************************
+! * General parameters
+! *********************************************************************
+¶meters_general
+ savefields = .F. ! Save fields to disk?
+/
+
+! *********************************************************************
+! * General solver parameters
+! *********************************************************************
+¶meters_solver
+ LUseO = T ,
+ LUseT = F ,
+ solvertype = 1, ! Solver type:
+ ! 1 : Richardson iteration
+ ! 2 : Conjugate gradient
+ resreduction = 1e-1 ! 1.0d-10 ! Required relative residual reduction
+ maxiter = 5 ! 50 ! Maximal number of iterations
+/
+
+! *********************************************************************
+! * Conjugate gradient parameters
+! *********************************************************************
+¶meters_conjugategradient
+ verbose = 10 ! Verbosity level
+ maxiter = 5, ! Maximal number of iterations
+ resreduction = 1.0e-1, ! Target residual reduction
+ n_prec = 1 ! Number of smoother applications in
+ ! preconditioner (N.B.: Using 0 is
+ ! inefficient, as the identity is used
+ ! for preconditioning, instead of using
+ ! unpreconditioned CG.)
+/
+
+! *********************************************************************
+! * Grid parameters
+! *********************************************************************
+¶meters_grid
+ n = 64, ! Number of horizontal grid cells
+ nz = 128, ! Number of vertical grid cells
+ L = 204800.0 ! 32000.0 metre * 64 pt grille ! 1.0, ! Size in horizontal direction
+ H = 19200.0 ! 150.0 metre * 128 niveaux ! 0.01, ! Size in vertical direction
+ vertbc = 2, ! Boundary conditions at top and bottom of the
+ ! atmosphere. 1 = DIRICHLET, 2 = NEUMANN
+ ! Note that Neumann boundary conditions only work
+ ! for coarsening in the horizontal only, as they are
+ ! not yet implemented in the prolongation operator.
+ graded = .F. ! Is the grid graded in the vertical direction?
+/
+
+! *********************************************************************
+! * Parallel communication parameters
+! *********************************************************************
+¶meters_communication
+ halo_size = 1 ! Size of halos (has to be 1 or 2)
+/
+
+! *********************************************************************
+! * Model parameters
+! *********************************************************************
+!
+! parameters of the Helmholtz operator
+!
+! -omega2*(d^2/dy^2 + d^2/dy^2 + lambda2*d^2/dz^2) u + delta u = RHS
+!
+¶meters_model
+ omega2 = 1.0,
+ lambda2 = 1.0 ! 100.0, ! Vertical coupling
+ delta = 0.0d0 ! Size of constant term
+/
+
+! *********************************************************************
+! * Smoother parameters
+! *********************************************************************
+!
+! parameters of the smoother
+!
+¶meters_smoother
+ smoother = 3, ! Smoother method
+ ! 3 = line SOR
+ ! 4 = line SSOR
+ ! 6 = line Jacobi
+ ordering = 2, ! Ordering of grid points (for smoother)
+ ! 1 = lexicographic
+ ! 2 = red-black ordering
+ !rho = 0.6666666666666666d0 ! Overrelaxation parameter
+ rho = 0.8d0 ! Overrelaxation parameter
+/
+
+! *********************************************************************
+! * Multigrid parameters
+! *********************************************************************
+¶meters_multigrid
+ verbose = 10, ! Verbosity level
+ n_lev = 7, ! Number of levels
+ lev_split = 3, ! First level where data is pulled together
+ n_presmooth = 1, ! Number of presmoothing steps
+ n_postsmooth = 1, ! Number of postsmoothing steps
+ n_coarsegridsmooth = 1, ! Number of smoothing steps on coarsest level
+ prolongation = 2 ! 2 best after modif !! ! Prologation method
+ ! 1 = constant interpolation
+ ! 2 = (tri-) linear interpolation
+ restriction = 1, ! Restriction method
+ ! 1 = cell average
+ ! 2 = Khalil
+ coarsegridsolver = 1 ! Solver on coarsest grid
+ ! 1 = use smoother
+ ! 2 = Conjugate gradient
+/
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/003_mesonh_step1/run_mesonh_step1_xyz b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/003_mesonh_step1/run_mesonh_step1_xyz
new file mode 100755
index 0000000000000000000000000000000000000000..b139bee388f05c7a19e396a1aa2d6d5d020fbc3a
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/003_mesonh_step1/run_mesonh_step1_xyz
@@ -0,0 +1,23 @@
+
+
+#export PGI_ACC_POOL_ALLOC=0
+export PGI_ACC_SYNCHRONOUS=1
+
+NP=${NP:-16}
+NPS=${NPS:-4}
+
+export MPIRUN=${MPIRUN:-"Mpirun -tag-output -bind-to none -map-by ppr:${NPS}:socket -x PGI_ACC_POOL_ALLOC -x PGI_ACC_SYNCHRONOUS -np ${NP} set_core_device_impair "}
+
+
+#
+set -x
+set -e
+#
+rm -f DA0064.*
+rm -f PGD00064.*
+rm -f D0064.1.DAR01.* OUTPUT_LISTING* pipe_name
+#
+ln -sf ../001_pgd1/PGD00064.{des,nc} .
+ln -sf ../002_prep_ideal_case/DA0064.{des,nc} .
+#
+time ${MPIRUN} MESONH${XYZ}
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/003_mesonh_step1/run_mesonh_step1_xyz_mppdb b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/003_mesonh_step1/run_mesonh_step1_xyz_mppdb
new file mode 100755
index 0000000000000000000000000000000000000000..dc87ac29b2da19c0c0ddca9d7902e0e5d16c9f20
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/003_mesonh_step1/run_mesonh_step1_xyz_mppdb
@@ -0,0 +1,33 @@
+#!/bin/bash
+#
+set -x
+#
+
+export MPIRUN=${MPIRUN:-"Mpirun --mca btl ^openib -tag-output -report-bindings -bind-to core -x ACC_DEVICE_TYPE=HOST -np 1 "}
+
+DIR_CLONE=dir_clone
+
+mkdir -p ${DIR_CLONE}
+
+(
+
+cd ${DIR_CLONE}
+
+rm -f DA0064.*
+rm -f PGD00064.*
+rm -f D0064.1.DAR01.* OUTPUT_LISTING* pipe_name
+#
+ln -sf ../../001_pgd1/PGD00064.{des,nc} .
+ln -sf ../../002_prep_ideal_case/DA0064.{des,nc} .
+cp ../EXSEG1.nam .
+
+)
+
+rm -f DA0064.*
+rm -f PGD00064.*
+rm -f D0064.1.DAR01.* OUTPUT_LISTING* pipe_name
+#
+ln -sf ../001_pgd1/PGD00064.{des,nc} .
+ln -sf ../002_prep_ideal_case/DA0064.{des,nc} .
+#
+time ${MPIRUN} MESONH${XYZ}
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/004_mesonh_step2/EXSEG1.nam b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/004_mesonh_step2/EXSEG1.nam
new file mode 100644
index 0000000000000000000000000000000000000000..a46063ee7c62303314fc9fc60781fdbb0b17ed81
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/004_mesonh_step2/EXSEG1.nam
@@ -0,0 +1,60 @@
+&NAM_CONFIO
+ LIO_ALLOW_NO_BACKUP = T , LIO_NO_WRITE = T
+/
+&NAM_CONFZ
+ ! NB_PROCIO_R=1 ,
+ ! NB_PROCIO_W=1 ,
+ LMNH_MPI_BSEND = F ,
+ ! MPI_BUFFER_SIZE = 40
+ !NZ_PROC=1
+/
+
+&NAM_LUNITn CINIFILE = "D0064.1.DAR01.004" , CINIFILEPGD = "PGD00064" /
+&NAM_CONFn LUSERV = T, LUSERC = T, LUSERR = T, LUSERI = T,
+ LUSERS = T, LUSERG = T, LUSERH = F, LUSECI = T /
+&NAM_DYNn XTSTEP = 15.00000 ,
+ CPRESOPT = "ZRESI",
+ !CPRESOPT = "ZSOLV",
+ NITR = 12,
+ XRELAX = 1., LHORELAX_UVWTH = T, LHORELAX_RV = T, LVE_RELAX = T,
+ NRIMX = 6, NRIMY = 6, XRIMKMAX = 0.0005, XT4DIFU = 4000 /
+&NAM_ADVn CUVW_ADV_SCHEME = "WENO_K", CMET_ADV_SCHEME = "PPM_01",
+ CSV_ADV_SCHEME = "PPM_01" /
+&NAM_PARAMn CTURB = "TKEL", ! "NONE" , ! "TKEL",
+ CRAD = "ECMWF",
+ CCLOUD = "ICE3", ! "NONE" , ! "ICE3",
+ CDCONV = "NONE", CSCONV = "NONE" /
+&NAM_PARAM_RADn XDTRAD = 1800., XDTRAD_CLONLY = 1800., LCLEAR_SKY = F,
+ NRAD_COLNBR = 500 /
+&NAM_PARAM_MFSHALLn XIMPL_MF = 1, CMF_UPDRAFT = "EDKF", CMF_CLOUD = "NONE",
+ LMIXUV = T, LMF_FLX = F /
+&NAM_LBCn CLBCX = 2*"OPEN", CLBCY = 2*"OPEN", XCPHASE = 20. /
+&NAM_TURBn XIMPL = 1., CTURBLEN = "DEAR", CTURBDIM = "3DIM",
+ LTURB_FLX = F, LTURB_DIAG = F, CSUBG_AUCV = "NONE", LSIGMAS = F,
+ LSIG_CONV = F, LSUBG_COND = F /
+&NAM_CH_MNHCn /
+&NAM_CONF CCONF = "RESTA" ! "START" ! "RESTA",
+ LFLAT = F, CEQNSYS = "DUR",
+ LLG=F, NMODEL = 1, NVERB = 0, CEXP = "D0064",
+ CSEG = "DAR02", CSPLIT = "BSPLITTING"
+ NHALO=1
+ LCHECK = F /
+&NAM_DYN XSEGLEN = 150.0 ! 150.0000 ,
+ LCORIO = T,
+ LNUMDIFU = F, LNUMDIFTH = F,
+ XALKTOP = 0.001, XALZBOT = 22000. /
+&NAM_BLANK /
+&NAM_NESTING /
+&NAM_BACKUP
+ XBAK_TIME(1,1) = 150.0 ! 3600.00000 ! 75. ! 3600.00000
+ XBAK_TIME(1,2) = 7200.00000
+ XBAK_TIME(1,3) = 10800.0000
+ XBAK_TIME(1,4) = 14400.0000
+ XBAK_TIME(1,5) = 18000.0000
+ XBAK_TIME(1,6) = 21600.0000
+ /
+&NAM_ISBAn /
+&NAM_SEAFLUXn /
+&NAM_DIAG_SURFn LSURF_BUDGET=T /
+&NAM_DIAG_ISBAn /
+&NAM_DIAG_SURF_ATMn /
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/004_mesonh_step2/MESONH_HOST2005 b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/004_mesonh_step2/MESONH_HOST2005
new file mode 100755
index 0000000000000000000000000000000000000000..07ad9fb226bfd7e7cc764749cb644e595ee25b13
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/004_mesonh_step2/MESONH_HOST2005
@@ -0,0 +1,5 @@
+#!/bin/bash
+export ACC_DEVICE=HOST
+ACC_DEVICE_TYPE=${ACC_DEVICE}
+exec /home/escj/DEV/MNH-55X-dev-OPENACC-juan-03/src/dir_obj-LXnvhpc2005-R8I4-MNH-V5-4-2-MPIAUTO-${OPTLEVEL}/ZSOLVER/MESONH
+
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/004_mesonh_step2/mppdb.nam b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/004_mesonh_step2/mppdb.nam
new file mode 100644
index 0000000000000000000000000000000000000000..19b2bcc6e70fa54b8b1d84b2b6ba7e8c677fdcb9
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/004_mesonh_step2/mppdb.nam
@@ -0,0 +1,7 @@
+&NAM_MPPDB
+MPPDB_DEBUG = .F. ,
+MPPDB_NBSON = 1 ,
+MPPDB_EXEC = "/home/escj/DEV/MNH-55X-dev-OPENACC-juan-03/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/004_mesonh_step2/MESONH_HOST2005" ,
+MPPDB_WDIR = "/home/escj/DEV/MNH-55X-dev-OPENACC-juan-03/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/004_mesonh_step2/dir_clone" ,
+MPPDB_CHECK_LB = .F.
+/
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/004_mesonh_step2/parameters_mg.nam b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/004_mesonh_step2/parameters_mg.nam
new file mode 100644
index 0000000000000000000000000000000000000000..da6e1ac67e7e9cd235de8d628379403ef298c3ad
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/004_mesonh_step2/parameters_mg.nam
@@ -0,0 +1,117 @@
+! *********************************************************************
+! *********************************************************************
+! ** **
+! ** Parameter file for geometric multigrid code **
+! ** **
+! *********************************************************************
+! *********************************************************************
+!
+! *********************************************************************
+! * General parameters
+! *********************************************************************
+¶meters_general
+ savefields = .F. ! Save fields to disk?
+/
+
+! *********************************************************************
+! * General solver parameters
+! *********************************************************************
+¶meters_solver
+ LUseO = T ,
+ LUseT = F ,
+ solvertype = 1, ! Solver type:
+ ! 1 : Richardson iteration
+ ! 2 : Conjugate gradient
+ resreduction = 1e-1 ! 1.0d-10 ! Required relative residual reduction
+ maxiter = 5 ! 50 ! Maximal number of iterations
+/
+
+! *********************************************************************
+! * Conjugate gradient parameters
+! *********************************************************************
+¶meters_conjugategradient
+ verbose = 10 ! Verbosity level
+ maxiter = 5, ! Maximal number of iterations
+ resreduction = 1.0e-1, ! Target residual reduction
+ n_prec = 1 ! Number of smoother applications in
+ ! preconditioner (N.B.: Using 0 is
+ ! inefficient, as the identity is used
+ ! for preconditioning, instead of using
+ ! unpreconditioned CG.)
+/
+
+! *********************************************************************
+! * Grid parameters
+! *********************************************************************
+¶meters_grid
+ n = 64, ! Number of horizontal grid cells
+ nz = 128, ! Number of vertical grid cells
+ L = 204800.0 ! 32000.0 metre * 64 pt grille ! 1.0, ! Size in horizontal direction
+ H = 19200.0 ! 150.0 metre * 128 niveaux ! 0.01, ! Size in vertical direction
+ vertbc = 2, ! Boundary conditions at top and bottom of the
+ ! atmosphere. 1 = DIRICHLET, 2 = NEUMANN
+ ! Note that Neumann boundary conditions only work
+ ! for coarsening in the horizontal only, as they are
+ ! not yet implemented in the prolongation operator.
+ graded = .F. ! Is the grid graded in the vertical direction?
+/
+
+! *********************************************************************
+! * Parallel communication parameters
+! *********************************************************************
+¶meters_communication
+ halo_size = 1 ! Size of halos (has to be 1 or 2)
+/
+
+! *********************************************************************
+! * Model parameters
+! *********************************************************************
+!
+! parameters of the Helmholtz operator
+!
+! -omega2*(d^2/dy^2 + d^2/dy^2 + lambda2*d^2/dz^2) u + delta u = RHS
+!
+¶meters_model
+ omega2 = 1.0,
+ lambda2 = 1.0 ! 100.0, ! Vertical coupling
+ delta = 0.0d0 ! Size of constant term
+/
+
+! *********************************************************************
+! * Smoother parameters
+! *********************************************************************
+!
+! parameters of the smoother
+!
+¶meters_smoother
+ smoother = 3, ! Smoother method
+ ! 3 = line SOR
+ ! 4 = line SSOR
+ ! 6 = line Jacobi
+ ordering = 2, ! Ordering of grid points (for smoother)
+ ! 1 = lexicographic
+ ! 2 = red-black ordering
+ !rho = 0.6666666666666666d0 ! Overrelaxation parameter
+ rho = 0.8d0 ! Overrelaxation parameter
+/
+
+! *********************************************************************
+! * Multigrid parameters
+! *********************************************************************
+¶meters_multigrid
+ verbose = 10, ! Verbosity level
+ n_lev = 7, ! Number of levels
+ lev_split = 3, ! First level where data is pulled together
+ n_presmooth = 1, ! Number of presmoothing steps
+ n_postsmooth = 1, ! Number of postsmoothing steps
+ n_coarsegridsmooth = 1, ! Number of smoothing steps on coarsest level
+ prolongation = 2 ! 2 best after modif !! ! Prologation method
+ ! 1 = constant interpolation
+ ! 2 = (tri-) linear interpolation
+ restriction = 1, ! Restriction method
+ ! 1 = cell average
+ ! 2 = Khalil
+ coarsegridsolver = 1 ! Solver on coarsest grid
+ ! 1 = use smoother
+ ! 2 = Conjugate gradient
+/
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/004_mesonh_step2/run_mesonh_step2_xyz b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/004_mesonh_step2/run_mesonh_step2_xyz
new file mode 100755
index 0000000000000000000000000000000000000000..a193cf5fec9815f718964f29cfab06aeb9aaf670
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/004_mesonh_step2/run_mesonh_step2_xyz
@@ -0,0 +1,26 @@
+#!/bin/bash
+#
+#
+set -x
+
+#export PGI_ACC_POOL_ALLOC=0
+export PGI_ACC_SYNCHRONOUS=1
+#export NVCOMPILER_ACC_POOL_ALLOC=0
+
+NP=${NP:-16}
+NPS=${NPS:-4}
+
+export MPIRUN=${MPIRUN:-"Mpirun -tag-output -report-bindings -bind-to none -map-by ppr:${NPS}:socket -x NVCOMPILER_ACC_POOL_ALLOC_MAXSIZE -x NVCOMPILER_ACC_POOL_ALLOC -x PGI_ACC_POOL_ALLOC -x PGI_ACC_SYNCHRONOUS -np ${NP} set_core_device_impair "}
+
+rm -f DA0064.*
+rm -f PGD00064.*
+rm -f D0064.1.DAR0?.* OUTPUT_LISTING1 pipe_name
+#
+ln -sf ../001_pgd1/PGD00064.* .
+ln -sf ../002_prep_ideal_case/DA0064.* .
+ln -sf ../003_mesonh_step1/D0064.1.DAR0?.* .
+#
+
+time ${MPIRUN} MESONH${XYZ}
+
+
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/004_mesonh_step2/run_mesonh_step2_xyz_mppdb b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/004_mesonh_step2/run_mesonh_step2_xyz_mppdb
new file mode 100755
index 0000000000000000000000000000000000000000..6418882ea5c933e6a77676ac1f8380356c3a19c4
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/004_mesonh_step2/run_mesonh_step2_xyz_mppdb
@@ -0,0 +1,48 @@
+#!/bin/bash
+#
+#
+
+set -x
+DIRCLONE=dir_clone
+mkdir -p ${DIRCLONE}
+
+export PGI_ACC_POOL_ALLOC=0
+#export PGI_ACC_POOL_SIZE=100MB
+
+export NP=${NP:-1}
+
+export MPIRUN=${MPIRUN:-Mpirun -tag-output -np ${NP}}
+
+#export PGI_ACC_NOTIFY=$(( 1+2+4+8+16+32 )) PGI_ACC_DEBUG=$(( 1+2+4+8+16+32 ))
+#export PGI_ACC_TIME=1
+
+export PGI_ACC_SYNCHRONOUS=1
+
+(
+cd ${DIRCLONE}
+#
+rm -f DA0064.*
+rm -f PGD00064.*
+rm -f D0064.1.DAR0?.* OUTPUT_LISTING1 pipe_name
+#
+cp ../EXSEG1.nam .
+ln -sf ../../001_pgd1/PGD00064.* .
+ln -sf ../../002_prep_ideal_case/DA0064.* .
+ln -sf ../../003_mesonh_step1/D0064.1.DAR01.* .
+#
+)
+
+#
+rm -f DA0064.*
+rm -f PGD00064.*
+rm -f D0064.1.DAR0?.* OUTPUT_LISTING1 pipe_name
+#
+ln -sf ../001_pgd1/PGD00064.* .
+ln -sf ../002_prep_ideal_case/DA0064.* .
+ln -sf ../003_mesonh_step1/D0064.1.DAR0?.* .
+#
+
+killall -r 'MESONH.*'
+time ${MPIRUN} MESONH${XYZ}
+
+
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/005_mg_main_mnh_all/parameters_mg.nam b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/005_mg_main_mnh_all/parameters_mg.nam
new file mode 120000
index 0000000000000000000000000000000000000000..adeab0e1859deb44804ccc14f23c2ccb35775079
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/005_mg_main_mnh_all/parameters_mg.nam
@@ -0,0 +1 @@
+../004_mesonh_step2/parameters_mg.nam
\ No newline at end of file
diff --git a/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/005_mg_main_mnh_all/run_mg_main_mnh_all_xyz b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/005_mg_main_mnh_all/run_mg_main_mnh_all_xyz
new file mode 100755
index 0000000000000000000000000000000000000000..c3694f2de5a11d8600af4cf5874b3248c26057f4
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/004_Hector_64x64_MNH-55X/005_mg_main_mnh_all/run_mg_main_mnh_all_xyz
@@ -0,0 +1,18 @@
+#!/bin/bash
+#
+#
+set -x
+
+#export PGI_ACC_POOL_ALLOC=0
+export PGI_ACC_SYNCHRONOUS=1
+#export NVCOMPILER_ACC_POOL_ALLOC=0
+
+NP=${NP:-16}
+NPS=${NPS:-4}
+
+export MPIRUN=${MPIRUN:-"Mpirun -tag-output -report-bindings -bind-to none -map-by ppr:${NPS}:socket -x NVCOMPILER_ACC_POOL_ALLOC_MAXSIZE -x NVCOMPILER_ACC_POOL_ALLOC -x PGI_ACC_POOL_ALLOC -x PGI_ACC_SYNCHRONOUS -np ${NP} set_core_device_impair "}
+
+
+time ${MPIRUN} MG_MAIN_MNH_ALL${XYZ}
+
+
diff --git a/MY_RUN/BENCH/CINES_MNH55X/Clean_MNH_files b/MY_RUN/BENCH/CINES_MNH55X/Clean_MNH_files
new file mode 100755
index 0000000000000000000000000000000000000000..e866afc41566ac11d859b95cff80967138c6b4fa
--- /dev/null
+++ b/MY_RUN/BENCH/CINES_MNH55X/Clean_MNH_files
@@ -0,0 +1,5 @@
+#!/bin/bash
+find . -name '*.des' -o -name '*.nc' -o -name '*.lfi' -o -name '*OUT*' -o -iname 'sorti*' \
+ -o -name '*.dir' -o -name '*.hdr' -o -name 'core' -o -name file_for_xtransfer \
+ -o -name PRESSURE -o -name 'REMAP*' -o -name 'fort.*' -o -name 'timing.txt*' | xargs -n1 rm -f
+find . -type d -name 'dir_*' | xargs -n1 rm -fr