diff --git a/MY_RUN/INTEGRATION_CASES/HPC/PANAME/000_get_external_files/script_get_transfert b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/000_get_external_files/script_get_transfert
new file mode 100755
index 0000000000000000000000000000000000000000..90197900b63944a5084caf5383b102027651de0a
--- /dev/null
+++ b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/000_get_external_files/script_get_transfert
@@ -0,0 +1,24 @@
+#!/bin/bash
+#SBATCH -J transfert
+#SBATCH -N 1 # nodes number
+#SBATCH -n 1 # CPUs number (on all nodes)
+#SBATCH -o transf.eo%j #
+#SBATCH -e transf.eo%j #
+#SBATCH -t 01:00:00 # time limit
+#SBATCH -p transfert
+
+HERE=$(pwd)
+ftget -h hendrix -u rodierq /home/rodierq/PHYSIO_AND_AROME_PANAME.tar.gz $HERE/PHYSIO_AND_AROME_PANAME.tar.gz
+
+# For users without access to MF/hendrix, data are also available at
+# http://mesonh.aero.obs-mip.fr/mesonh/dir_doc/dir_namelist_examples/GRIB_KTEST/PHYSIO_AND_AROME_PANAME.tar.gz
+# wget -c -nd http://mesonh.aero.obs-mip.fr/mesonh/dir_doc/dir_namelist_examples/GRIB_KTEST/PHYSIO_AND_AROME_PANAME.tar.gz
+
+tar xvfz PHYSIO_AND_AROME_PANAME.tar.gz
+mv PHYSIO_AND_AROME_PANAME ../.
+rm -f PHYSIO_AND_AROME_PANAME.tar.gz
+
+cd ../001_pgd1
+sbatch run_pgd1
+
+ja
diff --git a/MY_RUN/INTEGRATION_CASES/HPC/PANAME/001_pgd1/PRE_PGD1.nam b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/001_pgd1/PRE_PGD1.nam
new file mode 100644
index 0000000000000000000000000000000000000000..d110188d4e4246fae48fb6cfdb8daa05f6917517
--- /dev/null
+++ b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/001_pgd1/PRE_PGD1.nam
@@ -0,0 +1,62 @@
+&NAM_PGDFILE CPGDFILE='PGD_1200m', NHALO=0/
+
+&NAM_CONFIO LCDF4 = .TRUE.,
+ LLFIOUT = .FALSE.,
+ LLFIREAD = .FALSE./
+
+&NAM_PGD_GRID CGRID='CONF PROJ'/
+
+&NAM_CONF_PROJ XLAT0=48.8562481409,
+ XLON0=2.34581195868,
+ XRPK=0.0,
+ XBETA=0.0 /
+
+&NAM_CONF_PROJ_GRID XLATCEN=48.8562481409,
+ XLONCEN=2.34581195868,
+ NIMAX=400,
+ NJMAX=400,
+ XDX=1200.,
+ XDY=1200. /
+
+&NAM_PGD_SCHEMES CNATURE='ISBA',
+ CSEA='SEAFLX',
+ CTOWN='TEB',
+ CWATER='WATFLX',
+ LGARDEN= .T./
+
+&NAM_FRAC LECOSG=.F.
+/
+
+&NAM_ZS YZS='srtm_france_90m',
+ YZSFILETYPE='DIRECT' /
+
+&NAM_ZSFILTER NZSFILTER=1 /
+
+&NAM_ISBA YSAND = 'sand_0-200cm_mean_int',
+ YSANDFILETYPE = 'DIRECT',
+ YCLAY = 'clay_0-200cm_mean_int',
+ YCLAYFILETYPE = 'DIRECT',
+ CISBA = 'DIF',
+ CPHOTO = 'AST',
+ NGROUND_LAYER = 12 ,
+ CALBEDO = 'CM13',
+ NPATCH = 12 /
+
+&NAM_TREEDRAG LTREEDRAG=.FALSE. /
+
+&NAM_ISBAn LCANOPY_DRAG = .FALSE.,
+ CDIFSFCOND='MLCH'/
+
+&NAM_TEB
+
+ LCHECK_TEB = .FALSE. /
+
+
+&NAM_TEBn CZ0H = "KAND07",
+ CCH_BEM = 'ROW30'/
+
+&NAM_SEAFLUXn CSEA_FLUX="ECUME6" /
+
+&NAM_COVER YCOVER='ECOCLIMAP_I_GLOBAL',
+ YCOVERFILETYPE='DIRECT',
+ LRM_RIVER = F/
diff --git a/MY_RUN/INTEGRATION_CASES/HPC/PANAME/001_pgd1/run_pgd1 b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/001_pgd1/run_pgd1
new file mode 100755
index 0000000000000000000000000000000000000000..b00b932a98458c7807a9ed26afb243e4a2c4eb3f
--- /dev/null
+++ b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/001_pgd1/run_pgd1
@@ -0,0 +1,29 @@
+#!/bin/sh
+#SBATCH -J PGD1
+#SBATCH -N 1 # nodes number
+#SBATCH -n 16 # CPUs number (on all nodes)
+#SBATCH -o PGD_PANAME.eo%j
+#SBATCH -e PGD_PANAME.eo%j
+#SBATCH -t 00:05:00
+#SBATCH -p normal256
+
+module purge
+
+# Arret du job des la premiere erreur (optionnel)
+set -e
+
+# Profil utilisateur MesoNH (obligatoire)
+. ~/MNH-V5-7-0/conf/profile_mesonh-LXifort-R8I4-MNH-V5-7-0-ECRAD140-MPIAUTO-O2-PANAME_FOR571
+
+ln -sf ~rodierq/SAVE/mesonh/PGD/sand_0-200cm* .
+ln -sf ~rodierq/SAVE/mesonh/PGD/clay_0-200cm* .
+ln -sf ~rodierq/SAVE/mesonh/PGD/srtm_france_90m* .
+ln -sf ~rodierq/SAVE/mesonh/PGD/ECOCLIMAP_I_GLOBAL* .
+
+# Exécution (obligatoire)
+export MPIRUN="Mpirun -np 16 "
+time ${MPIRUN} PREP_PGD${XYZ}
+
+### Job suivant
+cd ../002_pgd2
+sbatch run_pgd2
diff --git a/MY_RUN/INTEGRATION_CASES/HPC/PANAME/002_pgd2/PRE_PGD1.nam b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/002_pgd2/PRE_PGD1.nam
new file mode 100644
index 0000000000000000000000000000000000000000..b20cba420064a9cf4627e6a53e6fd530f6a8ef0f
--- /dev/null
+++ b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/002_pgd2/PRE_PGD1.nam
@@ -0,0 +1,150 @@
+&NAM_PGDFILE CPGDFILE='PGD_300m', NHALO=0/
+
+&NAM_CONFIO LCDF4 = .TRUE.,
+ LLFIOUT = .FALSE.,
+ LLFIREAD = .FALSE./
+
+&NAM_PGD_GRID YINIFILE='PGD_1200m',
+ YINIFILETYPE='MESONH' /
+
+&NAM_INIFILE_CONF_PROJ IXOR=150,
+ IYOR=150,
+ IXSIZE=100,
+ IYSIZE=100,
+ IDXRATIO=4,
+ IDYRATIO=4 /
+
+&NAM_PGD_SCHEMES CNATURE='ISBA',
+ CSEA='SEAFLX',
+ CTOWN='TEB',
+ CWATER='WATFLX',
+ LGARDEN= .T./
+
+&NAM_ISBA YSAND = 'sand_0-200cm_mean_int',
+ YSANDFILETYPE = 'DIRECT',
+ YCLAY = 'clay_0-200cm_mean_int',
+ YCLAYFILETYPE = 'DIRECT',
+ CISBA = 'DIF',
+ CPHOTO = 'AST',
+ NGROUND_LAYER = 12,
+ CALBEDO = 'CM13',
+ NPATCH = 12
+/
+&NAM_COVER YCOVER='ECOCLIMAP_I_GLOBAL',
+ YCOVERFILETYPE='DIRECT' ,
+ LRM_RIVER = F/
+
+&NAM_TREEDRAG LTREEDRAG=.FALSE. /
+
+&NAM_ISBAn LCANOPY_DRAG = .FALSE. ,
+ CDIFSFCOND='MLCH'/
+
+&NAM_TEBn CZ0H = "KAND07",
+ CCH_BEM = 'ROW30'/
+
+&NAM_SEAFLUXn CSEA_FLUX="DIRECT" /
+
+&NAM_ZS YZS='srtm_france_30m',
+ YZSFILETYPE='DIRECT' /
+
+&NAM_ZSFILTER NZSFILTER=1 /
+
+&NAM_FRAC LECOCLIMAP = TRUE,
+ LECOSG=.F.,
+ CFNAM_SEA = "./SEA.txt",
+ CFTYP_SEA = "ASCLLV",
+ CFNAM_WATER = "./WATER.txt",
+ CFTYP_WATER = "ASCLLV",
+ CFNAM_NATURE = "./NATURE.txt",
+ CFTYP_NATURE = "ASCLLV",
+ CFNAM_TOWN = "./TOWN.txt",
+ CFTYP_TOWN = "ASCLLV" /
+
+&NAM_TEB NROOF_LAYER = 6,
+ NWALL_LAYER = 6,
+ NFLOOR_LAYER = 6,
+ NMASS_LAYER = 6,
+ NBEMCOMP = 1, ! Pas d'usage de batiment donc 1 seul compartiment
+ NTIME_CHANGE = 0,
+ CBEM = "BEM",
+ LURBHYDRO = .F.,
+ CURBTREE = "TREE",
+ LSPARTACUS = .F.,
+ LEXPLW = .F.,
+ CROAD_GRID = 'MEDIUM',
+ CCOOL_COIL = "DXCOIL",
+ CHEAT_COIL = "FINCAP",
+ LAUTOSIZE = .TRUE.,
+ NTEB_PATCH=1,
+ LGREENROOF=.FALSE.,
+ LSOLAR_PANEL=.FALSE.,
+ LCHECK_TEB = .FALSE. /
+
+&NAM_DATA_TEB CFNAM_BLD = "./BLD.txt",
+ CFTYP_BLD = "ASCLLV",
+ CFNAM_BLD_HEIGHT = "./BLD_HEIGHT.txt",
+ CFTYP_BLD_HEIGHT = "ASCLLV",
+ CFNAM_ROAD = './ROAD.txt' ,
+ CFTYP_ROAD = 'ASCLLV' ,
+ CFNAM_WALL_O_HOR = "./WALL_O_HOR.txt",
+ CFTYP_WALL_O_HOR = "ASCLLV",
+ CFNAM_FRAC_HVEG = "./FRAC_VEG_HTE.txt",
+ CFTYP_FRAC_HVEG = "ASCLLV",
+ CFNAM_FRAC_LVEG = "./FRAC_VEG_BAS.txt",
+ CFTYP_FRAC_LVEG = "ASCLLV",
+ CFNAM_FRAC_NVEG = "./FRAC_NVEG.txt",
+ CFTYP_FRAC_NVEG = "ASCLLV"/
+
+
+&NAM_DATA_BEM
+ XUNIF_T_SIZE_MAX = 313. ,
+ XUNIF_T_SIZE_MIN = 258. ,
+ XUNIF_TCOOL_OCCD = 313. ,
+ XUNIF_TCOOL_OCCN = 313. ,
+ XUNIF_TCOOL_VCDD = 313. ,
+ XUNIF_TCOOL_VCDN = 313. ,
+ XUNIF_TCOOL_VCLD = 313. ,
+ XUNIF_THEAT_OCCD = 293.16 ,
+ XUNIF_THEAT_OCCN = 293.16 ,
+ XUNIF_THEAT_VCDD = 293.16 ,
+ XUNIF_THEAT_VCDN = 293.16 ,
+ XUNIF_THEAT_VCLD = 293.16 ,
+
+/
+
+&NAM_DATA_TEB_GARDEN CTYP_GARDEN_HVEG = 'TEBD',
+ CTYP_GARDEN_LVEG = 'GRAS',
+ CTYP_GARDEN_NVEG = 'NO ',
+ NTIME_GD = 12,
+ XUNIF_LAI_HVEG(01) = 1.0, ! LAI Diane ??
+ XUNIF_LAI_HVEG(02) = 1.0,
+ XUNIF_LAI_HVEG(03) = 1.0,
+ XUNIF_LAI_HVEG(04) = 1.7,
+ XUNIF_LAI_HVEG(05) = 2.3,
+ XUNIF_LAI_HVEG(06) = 3.0,
+ XUNIF_LAI_HVEG(07) = 3.0,
+ XUNIF_LAI_HVEG(08) = 2.7,
+ XUNIF_LAI_HVEG(09) = 2.4,
+ XUNIF_LAI_HVEG(10) = 2.1,
+ XUNIF_LAI_HVEG(11) = 1.8,
+ XUNIF_LAI_HVEG(12) = 1.5,
+ XUNIF_LAI_LVEG(01) = 1.5,
+ XUNIF_LAI_LVEG(02) = 1.5,
+ XUNIF_LAI_LVEG(03) = 1.5,
+ XUNIF_LAI_LVEG(04) = 1.5,
+ XUNIF_LAI_LVEG(05) = 1.5,
+ XUNIF_LAI_LVEG(06) = 1.5,
+ XUNIF_LAI_LVEG(07) = 1.5,
+ XUNIF_LAI_LVEG(08) = 1.5,
+ XUNIF_LAI_LVEG(09) = 1.5,
+ XUNIF_LAI_LVEG(10) = 1.5,
+ XUNIF_LAI_LVEG(11) = 1.5,
+ XUNIF_LAI_LVEG(12) = 1.5,
+ CSHAPE_GARDEN_HVEG = 'CYL' ,
+ CFNAM_H_HVEG = 'TREE_HEIGHT.txt' ,
+ CFTYP_H_HVEG = 'ASCLLV' ,
+ CFNAM_HTRUNK_HVEG = 'TREE_TRUNK.txt' ,
+ CFTYP_HTRUNK_HVEG = 'ASCLLV' ,
+ CFNAM_WCROWN_HVEG = 'TREE_CROWN.txt' ,
+ CFTYP_WCROWN_HVEG = 'ASCLLV' /
+
diff --git a/MY_RUN/INTEGRATION_CASES/HPC/PANAME/002_pgd2/run_pgd2 b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/002_pgd2/run_pgd2
new file mode 100755
index 0000000000000000000000000000000000000000..e6b81fc3ee1a547295f4b32951c20ef15b7282b9
--- /dev/null
+++ b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/002_pgd2/run_pgd2
@@ -0,0 +1,33 @@
+#!/bin/sh
+#SBATCH -J PGD2
+#SBATCH -N 1 # nodes number
+#SBATCH -n 16 # CPUs number (on all nodes)
+#SBATCH -o PGD2_PANAME.eo%j
+#SBATCH -e PGD2_PANAME.eo%j
+#SBATCH -t 00:15:00
+#SBATCH -p normal256
+module purge
+
+# Arret du job des la premiere erreur (optionnel)
+set -e
+
+# Profil utilisateur MesoNH (obligatoire)
+. ~/MNH-V5-7-0/conf/profile_mesonh-LXifort-R8I4-MNH-V5-7-0-ECRAD140-MPIAUTO-O2-PANAME_FOR571
+
+ln -sf ../001_pgd1/PGD_1200m.des .
+ln -sf ../001_pgd1/PGD_1200m.nc .
+
+ln -sf ~rodierq/SAVE/mesonh/PGD/sand_0-200cm* .
+ln -sf ~rodierq/SAVE/mesonh/PGD/clay_0-200cm* .
+ln -sf ~rodierq/SAVE/mesonh/PGD/srtm_france_30m* .
+ln -sf ~rodierq/SAVE/mesonh/PGD/ECOCLIMAP_I_GLOBAL* .
+ln -sf ../PHYSIO_AND_AROME_PANAME/OSM/* .
+
+
+# Exécution (obligatoire)
+export MPIRUN="Mpirun -np 16 "
+time ${MPIRUN} PREP_PGD${XYZ}
+
+### Job suivant
+cd ../003_pgd3
+sbatch run_pgd3
diff --git a/MY_RUN/INTEGRATION_CASES/HPC/PANAME/003_pgd3/PRE_PGD1.nam b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/003_pgd3/PRE_PGD1.nam
new file mode 100644
index 0000000000000000000000000000000000000000..23207e4f855235d0682eb7f5b95bdff626ea7a03
--- /dev/null
+++ b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/003_pgd3/PRE_PGD1.nam
@@ -0,0 +1,193 @@
+&NAM_PGDFILE CPGDFILE='PGD_100m', NHALO=0/
+
+&NAM_CONFIO LCDF4 = .TRUE.,
+ LLFIOUT = .FALSE.,
+ LLFIREAD = .FALSE./
+
+&NAM_PGD_GRID YINIFILE='PGD_300m',
+ YINIFILETYPE='MESONH' /
+
+&NAM_INIFILE_CONF_PROJ IXOR=150,
+ IYOR=150,
+ IXSIZE=100,
+ IYSIZE=100,
+ IDXRATIO=3,
+ IDYRATIO=3 /
+
+&NAM_PGD_SCHEMES CNATURE='ISBA',
+ CSEA ='SEAFLX',
+ CWATER ='WATFLX',
+ CTOWN ='TEB',
+ LGARDEN= .T./
+
+&NAM_ISBA YSAND = 'sand_0-200cm_mean_int',
+ YSANDFILETYPE = 'DIRECT',
+ YCLAY = 'clay_0-200cm_mean_int',
+ YCLAYFILETYPE = 'DIRECT',
+ CISBA = 'DIF',
+ CPHOTO = 'AST',
+ NGROUND_LAYER = 12 ,
+ CALBEDO = 'CM13',
+ NPATCH = 12
+/
+
+&NAM_COVER YCOVER='ECOCLIMAP_I_GLOBAL',
+ YCOVERFILETYPE='DIRECT',
+ LRM_RIVER = F/
+
+&NAM_TREEDRAG LTREEDRAG=.FALSE. /
+
+&NAM_ISBAn LCANOPY_DRAG = .FALSE. ,
+ CDIFSFCOND='MLCH'/
+
+&NAM_TEBn CZ0H = "KAND07",
+ CCH_BEM = 'ROW30'/
+
+&NAM_SEAFLUXn CSEA_FLUX="DIRECT" /
+
+&NAM_ZS YZS='srtm_france_30m',
+ YZSFILETYPE='DIRECT' /
+
+&NAM_ZSFILTER NZSFILTER=1 /
+
+&NAM_FRAC LECOCLIMAP = TRUE,
+ LECOSG=.F.,
+ CFNAM_SEA = "./SEA.txt",
+ CFTYP_SEA = "ASCLLV",
+ CFNAM_WATER = "./WATER.txt",
+ CFTYP_WATER = "ASCLLV",
+ CFNAM_NATURE = "./NATURE.txt",
+ CFTYP_NATURE = "ASCLLV",
+ CFNAM_TOWN = "./TOWN.txt",
+ CFTYP_TOWN = "ASCLLV" /
+
+&NAM_TEB NROOF_LAYER = 6,
+ NWALL_LAYER = 6,
+ NFLOOR_LAYER = 6,
+ NMASS_LAYER = 6,
+ NBEMCOMP = 6,
+ NTIME_CHANGE = 0,
+ CBEM = "BEM",
+ LURBHYDRO = .F.,
+ CURBTREE = "TREE",
+ LSPARTACUS = .F.,
+ LEXPLW = .F.,
+ CROAD_GRID = 'MEDIUM',
+ CCOOL_COIL = "DXCOIL",
+ CHEAT_COIL = "FINCAP",
+ LAUTOSIZE = .TRUE.,
+ NTEB_PATCH=1,
+ LGREENROOF=.FALSE.,
+ LSOLAR_PANEL=.FALSE.,
+ LCHECK_TEB = .FALSE. /
+
+&NAM_DATA_TEB CFNAM_BLD = "./BLD.txt",
+ CFTYP_BLD = "ASCLLV",
+ CFNAM_BLD_HEIGHT = "./BLD_HEIGHT.txt",
+ CFTYP_BLD_HEIGHT = "ASCLLV",
+ CFNAM_ROAD = './ROAD.txt' ,
+ CFTYP_ROAD = 'ASCLLV' ,
+ CFNAM_WALL_O_HOR = "./WALL_O_HOR.txt",
+ CFTYP_WALL_O_HOR = "ASCLLV",
+ CFNAM_FRAC_HVEG = "./FRAC_VEG_HTE.txt",
+ CFTYP_FRAC_HVEG = "ASCLLV",
+ CFNAM_FRAC_LVEG = "./FRAC_VEG_BAS.txt",
+ CFTYP_FRAC_LVEG = "ASCLLV",
+ CFNAM_FRAC_NVEG = "./FRAC_NVEG.txt",
+ CFTYP_FRAC_NVEG = "ASCLLV",
+!DATA MAPUCE
+ CFNAM_BLDTYPE = "./ISLETTYPE.txt",
+ CFTYP_BLDTYPE = "ASCLLV",
+ CFNAM_IND_BLD_AGE = "./PERMAJ_IND.txt",
+ CFTYP_IND_BLD_AGE = "ASCLLV",
+ CFNAM_COL_BLD_AGE = "./PERMAJ_COL.txt",
+ CFTYP_COL_BLD_AGE = "ASCLLV",
+ CFNAM_P1TERRITORY = "./P1_TERR.txt",
+ CFTYP_P1TERRITORY = "ASCLLV",
+ CFNAM_PXTERRITORY = "./PX_TERR.txt",
+ CFTYP_PXTERRITORY = "ASCLLV",
+ CFNAM_FRACIHS = "./FRAIHS.txt",
+ CFTYP_FRACIHS = "ASCLLV",
+ CFNAM_FRACCHS = "./FRACHS.txt",
+ CFTYP_FRACCHS = "ASCLLV",
+ CFNAM_FRACCOM = "./FRACOM.txt",
+ CFTYP_FRACCOM = "ASCLLV",
+ CFNAM_FRACTER = "./FRATER.txt",
+ CFTYP_FRACTER = "ASCLLV",
+ CFNAM_FRACIND = "./FRAIND.txt",
+ CFTYP_FRACIND = "ASCLLV",
+ CFNAM_FRACNHE = "./FRANHE.txt",
+ CFTYP_FRACNHE = "ASCLLV",
+ CFNAM_FRACPAV = "./FRAPAV.txt",
+ CFTYP_FRACPAV = "ASCLLV",
+ CFNAM_FRACMRI = "./FRAMRI.txt",
+ CFTYP_FRACMRI = "ASCLLV",
+ CFNAM_FRACHRI = "./FRAHRI.txt",
+ CFTYP_FRACHRI = "ASCLLV",
+ CFNAM_FRACATB = "./FRAATB.txt",
+ CFTYP_FRACATB = "ASCLLV",
+ CFNAM_FOEQI_MAIS = "./FORT_EQI_MAIS.txt",
+ CFTYP_FOEQI_MAIS = "ASCLLV",
+ CFNAM_FOEQI_APPT = "./FORT_EQI_APPT.txt",
+ CFTYP_FOEQI_APPT = "ASCLLV",
+ CFNAM_FAEQI_MAIS = "./FAIB_EQI_MAIS.txt",
+ CFTYP_FAEQI_MAIS = "ASCLLV",
+ CFNAM_FAEQI_APPT = "./FAIB_EQI_APPT.txt",
+ CFTYP_FAEQI_APPT = "ASCLLV",
+ CFNAM_CRE_MAIS = "./FORT_CRE_MAIS.txt",
+ CFTYP_CRE_MAIS = "ASCLLV",
+ CFNAM_CRE_APPT = "./FORT_CRE_APPT.txt",
+ CFTYP_CRE_APPT = "ASCLLV",
+ CFNAM_USETYPE = "./USAGE.txt",
+ CFTYP_USETYPE = "ASCLLV",
+ CCSVFILEARCHI = "TABLEAU_MAPUCE_ARCHI.csv",
+ CCSVFILECOMPO = "TABLEAU_MAPUCE_COMPO.csv"
+/
+
+&NAM_DATA_BEM
+ XUNIF_T_SIZE_MAX = 313. ,
+ XUNIF_T_SIZE_MIN = 258. ,
+ XUNIF_TCOOL_OCCD = 313. ,
+ XUNIF_TCOOL_OCCN = 313. ,
+ XUNIF_TCOOL_VCDD = 313. ,
+ XUNIF_TCOOL_VCDN = 313. ,
+ XUNIF_TCOOL_VCLD = 313. ,
+/
+
+&NAM_DATA_TEB_GARDEN CTYP_GARDEN_HVEG = 'TEBD',
+ CTYP_GARDEN_LVEG = 'GRAS',
+ CTYP_GARDEN_NVEG = 'NO ',
+ NTIME_GD = 12,
+ XUNIF_LAI_HVEG(01) = 1.0,
+ XUNIF_LAI_HVEG(02) = 1.0,
+ XUNIF_LAI_HVEG(03) = 1.0,
+ XUNIF_LAI_HVEG(04) = 1.7,
+ XUNIF_LAI_HVEG(05) = 2.3,
+ XUNIF_LAI_HVEG(06) = 3.0,
+ XUNIF_LAI_HVEG(07) = 3.0,
+ XUNIF_LAI_HVEG(08) = 2.7,
+ XUNIF_LAI_HVEG(09) = 2.4,
+ XUNIF_LAI_HVEG(10) = 2.1,
+ XUNIF_LAI_HVEG(11) = 1.8,
+ XUNIF_LAI_HVEG(12) = 1.5,
+ XUNIF_LAI_LVEG(01) = 1.5,
+ XUNIF_LAI_LVEG(02) = 1.5,
+ XUNIF_LAI_LVEG(03) = 1.5,
+ XUNIF_LAI_LVEG(04) = 1.5,
+ XUNIF_LAI_LVEG(05) = 1.5,
+ XUNIF_LAI_LVEG(06) = 1.5,
+ XUNIF_LAI_LVEG(07) = 1.5,
+ XUNIF_LAI_LVEG(08) = 1.5,
+ XUNIF_LAI_LVEG(09) = 1.5,
+ XUNIF_LAI_LVEG(10) = 1.5,
+ XUNIF_LAI_LVEG(11) = 1.5,
+ XUNIF_LAI_LVEG(12) = 1.5,
+ CSHAPE_GARDEN_HVEG = 'CYL' ,
+ CFNAM_H_HVEG = 'TREE_HEIGHT.txt' ,
+ CFTYP_H_HVEG = 'ASCLLV' ,
+ CFNAM_HTRUNK_HVEG = 'TREE_TRUNK.txt' ,
+ CFTYP_HTRUNK_HVEG = 'ASCLLV' ,
+ CFNAM_WCROWN_HVEG = 'TREE_CROWN.txt' ,
+ CFTYP_WCROWN_HVEG = 'ASCLLV' /
+
+
diff --git a/MY_RUN/INTEGRATION_CASES/HPC/PANAME/003_pgd3/run_pgd3 b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/003_pgd3/run_pgd3
new file mode 100755
index 0000000000000000000000000000000000000000..db24bd448a94def5ce053040141844c55419e384
--- /dev/null
+++ b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/003_pgd3/run_pgd3
@@ -0,0 +1,35 @@
+#!/bin/sh
+#SBATCH -J PGD3
+#SBATCH -N 1 # nodes number
+#SBATCH -n 16 # CPUs number (on all nodes)
+#SBATCH -o PGD3_PANAME.eo%j
+#SBATCH -e PGD3_PANAME.eo%j
+#SBATCH -t 00:15:00
+#SBATCH -p normal256
+module purge
+
+# Arret du job des la premiere erreur (optionnel)
+set -e
+
+# Profil utilisateur MesoNH (obligatoire)
+. ~/MNH-V5-7-0/conf/profile_mesonh-LXifort-R8I4-MNH-V5-7-0-ECRAD140-MPIAUTO-O2-PANAME_FOR571
+
+ln -sf ../001_pgd1/PGD_1200m.des .
+ln -sf ../001_pgd1/PGD_1200m.nc .
+ln -sf ../002_pgd2/PGD_300m.des .
+ln -sf ../002_pgd2/PGD_300m.nc .
+
+ln -sf ~rodierq/SAVE/mesonh/PGD/sand_0-200cm* .
+ln -sf ~rodierq/SAVE/mesonh/PGD/clay_0-200cm* .
+ln -sf ~rodierq/SAVE/mesonh/PGD/srtm_france_30m* .
+ln -sf ~rodierq/SAVE/mesonh/PGD/ECOCLIMAP_I_GLOBAL* .
+ln -sf ../PHYSIO_AND_AROME_PANAME/MAPUCE/* .
+ln -sf ../PHYSIO_AND_AROME_PANAME/OSM/* .
+
+# Exécution (obligatoire)
+export MPIRUN="Mpirun -np 16 "
+time ${MPIRUN} PREP_PGD${XYZ}
+
+### Job suivant
+cd ../004_nest_pgd
+sbatch run_nest_pgd
diff --git a/MY_RUN/INTEGRATION_CASES/HPC/PANAME/004_nest_pgd/PRE_NEST_PGD1.nam b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/004_nest_pgd/PRE_NEST_PGD1.nam
new file mode 100644
index 0000000000000000000000000000000000000000..d6e77dd0d4ec79a455221c55f875b002a6b42cbf
--- /dev/null
+++ b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/004_nest_pgd/PRE_NEST_PGD1.nam
@@ -0,0 +1,10 @@
+&NAM_CONFIO LCDF4 = .TRUE.,
+ LLFIOUT = .FALSE.,
+ LLFIREAD = .FALSE./
+
+&NAM_PGD1 YPGD1 = 'PGD_1200m' /
+&NAM_PGD2 YPGD2 = 'PGD_300m', IDAD = 1 /
+&NAM_PGD3 YPGD3 = 'PGD_100m', IDAD = 2 /
+
+&NAM_NEST_PGD YNEST='e1' /
+
diff --git a/MY_RUN/INTEGRATION_CASES/HPC/PANAME/004_nest_pgd/run_nest_pgd b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/004_nest_pgd/run_nest_pgd
new file mode 100755
index 0000000000000000000000000000000000000000..d00f769ed6ea8bd74620ae14243d72d9c8d24a30
--- /dev/null
+++ b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/004_nest_pgd/run_nest_pgd
@@ -0,0 +1,36 @@
+#!/bin/sh
+#SBATCH -J NEST
+#SBATCH -N 1 # nodes number
+#SBATCH -n 16 # CPUs number (on all nodes)
+#SBATCH -o NEST.eo%j
+#SBATCH -e NEST.eo%j
+#SBATCH -t 0:05:00
+#SBATCH -p normal256
+#SBATCH --no-requeue
+
+module purge
+
+
+# Arret du job des la premiere erreur (optionnel)
+set -e
+
+
+# Profil utilisateur MesoNH (obligatoire)
+#. ~/MNH-V5-6-2/conf/profile_mesonh-LXifort-R8I4-MNH-V5-6-2-MPIAUTO-O2-PANAME_FOR571
+. ~/MNH-V5-6-2/conf/profile_mesonh-LXifort-R8I4-MNH-V5-6-2-ECRAD140-MPIAUTO-O2-PANAME_FOR571
+
+
+
+ln -sf ../001_pgd1/PGD_1200m.des .
+ln -sf ../001_pgd1/PGD_1200m.nc .
+ln -sf ../002_pgd2/PGD_300m.des .
+ln -sf ../002_pgd2/PGD_300m.nc .
+ln -sf ../003_pgd3/PGD_100m.des .
+ln -sf ../003_pgd3/PGD_100m.nc .
+
+# Exécution (obligatoire)
+export MPIRUN="Mpirun -np 16 "
+time ${MPIRUN} PREP_NEST_PGD${XYZ}
+
+cd ../005_coupling1
+sbatch run_preal_xyz_1
diff --git a/MY_RUN/INTEGRATION_CASES/HPC/PANAME/005_coupling1/PRE_REAL1.nam_20230712_03H b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/005_coupling1/PRE_REAL1.nam_20230712_03H
new file mode 100644
index 0000000000000000000000000000000000000000..9b4fee46c6c41fd17daf825c2c6d3764be7288af
--- /dev/null
+++ b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/005_coupling1/PRE_REAL1.nam_20230712_03H
@@ -0,0 +1,23 @@
+&NAM_CONFIO LCDF4=T,
+ LLFIREAD=F,
+ NIO_ABORT_LEVEL=1,
+ NGEN_ABORT_LEVEL=1/
+
+&NAM_FILE_NAMES HATMFILE='arome.AN.20230712.03',
+ HATMFILETYPE='GRIBEX',
+ HPGDFILE='PGD_1200m.neste1',
+ CINIFILE='20230712_03H' /
+
+&NAM_REAL_CONF CEQNSYS='DUR',
+ CPRESOPT='ZRESI',
+ LCOUPLING=T /
+
+&NAM_CONFZ MPI_BUFFER_SIZE=500/
+
+&NAM_VER_GRID NKMAX=68,
+ YZGRID_TYPE='FUNCTN',
+ ZDZGRD=4,
+ ZDZTOP=500,
+ ZZMAX_STRGRD=1000,
+ ZSTRGRD=17,
+ ZSTRTOP=4 /
diff --git a/MY_RUN/INTEGRATION_CASES/HPC/PANAME/005_coupling1/PRE_REAL1.nam_20230712_06H b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/005_coupling1/PRE_REAL1.nam_20230712_06H
new file mode 100644
index 0000000000000000000000000000000000000000..c2d8a4c92393f9a19489aba6762318d93e8dc3a2
--- /dev/null
+++ b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/005_coupling1/PRE_REAL1.nam_20230712_06H
@@ -0,0 +1,23 @@
+&NAM_CONFIO LCDF4=T,
+ LLFIREAD=F,
+ NIO_ABORT_LEVEL=1,
+ NGEN_ABORT_LEVEL=1/
+
+&NAM_FILE_NAMES HATMFILE='arome.AN.20230712.06',
+ HATMFILETYPE='GRIBEX',
+ HPGDFILE='PGD_1200m.neste1',
+ CINIFILE='20230712_06H' /
+
+&NAM_REAL_CONF CEQNSYS='DUR',
+ CPRESOPT='ZRESI',
+ LCOUPLING=T /
+
+&NAM_CONFZ MPI_BUFFER_SIZE=500/
+
+&NAM_VER_GRID NKMAX=68,
+ YZGRID_TYPE='FUNCTN',
+ ZDZGRD=4,
+ ZDZTOP=500,
+ ZZMAX_STRGRD=1000,
+ ZSTRGRD=17,
+ ZSTRTOP=4 /
diff --git a/MY_RUN/INTEGRATION_CASES/HPC/PANAME/005_coupling1/PRE_REAL1.nam_20230712_09H b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/005_coupling1/PRE_REAL1.nam_20230712_09H
new file mode 100644
index 0000000000000000000000000000000000000000..8d2a92fbedfd5ba2b417f02263e672f988ed735f
--- /dev/null
+++ b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/005_coupling1/PRE_REAL1.nam_20230712_09H
@@ -0,0 +1,23 @@
+&NAM_CONFIO LCDF4=T,
+ LLFIREAD=F,
+ NIO_ABORT_LEVEL=1,
+ NGEN_ABORT_LEVEL=1/
+
+&NAM_FILE_NAMES HATMFILE='arome.AN.20230712.09',
+ HATMFILETYPE='GRIBEX',
+ HPGDFILE='PGD_1200m.neste1',
+ CINIFILE='20230712_09H' /
+
+&NAM_REAL_CONF CEQNSYS='DUR',
+ CPRESOPT='ZRESI',
+ LCOUPLING=T /
+
+&NAM_CONFZ MPI_BUFFER_SIZE=500/
+
+&NAM_VER_GRID NKMAX=68,
+ YZGRID_TYPE='FUNCTN',
+ ZDZGRD=4,
+ ZDZTOP=500,
+ ZZMAX_STRGRD=1000,
+ ZSTRGRD=17,
+ ZSTRTOP=4 /
diff --git a/MY_RUN/INTEGRATION_CASES/HPC/PANAME/005_coupling1/PRE_REAL1.nam_20230712_12H b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/005_coupling1/PRE_REAL1.nam_20230712_12H
new file mode 100644
index 0000000000000000000000000000000000000000..d4cd64d263d70a80335dfc1415f60cfbee696fc8
--- /dev/null
+++ b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/005_coupling1/PRE_REAL1.nam_20230712_12H
@@ -0,0 +1,23 @@
+&NAM_CONFIO LCDF4=T,
+ LLFIREAD=F,
+ NIO_ABORT_LEVEL=1,
+ NGEN_ABORT_LEVEL=1/
+
+&NAM_FILE_NAMES HATMFILE='arome.AN.20230712.12',
+ HATMFILETYPE='GRIBEX',
+ HPGDFILE='PGD_1200m.neste1',
+ CINIFILE='20230712_12H' /
+
+&NAM_REAL_CONF CEQNSYS='DUR',
+ CPRESOPT='ZRESI',
+ LCOUPLING=T /
+
+&NAM_CONFZ MPI_BUFFER_SIZE=500/
+
+&NAM_VER_GRID NKMAX=68,
+ YZGRID_TYPE='FUNCTN',
+ ZDZGRD=4,
+ ZDZTOP=500,
+ ZZMAX_STRGRD=1000,
+ ZSTRGRD=17,
+ ZSTRTOP=4 /
diff --git a/MY_RUN/INTEGRATION_CASES/HPC/PANAME/005_coupling1/PRE_REAL1.nam_20230712_15H b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/005_coupling1/PRE_REAL1.nam_20230712_15H
new file mode 100644
index 0000000000000000000000000000000000000000..954eecea8aabf7e517dd7cc33cc284655f57f307
--- /dev/null
+++ b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/005_coupling1/PRE_REAL1.nam_20230712_15H
@@ -0,0 +1,23 @@
+&NAM_CONFIO LCDF4=T,
+ LLFIREAD=F,
+ NIO_ABORT_LEVEL=1,
+ NGEN_ABORT_LEVEL=1/
+
+&NAM_FILE_NAMES HATMFILE='arome.AN.20230712.15',
+ HATMFILETYPE='GRIBEX',
+ HPGDFILE='PGD_1200m.neste1',
+ CINIFILE='20230712_15H' /
+
+&NAM_REAL_CONF CEQNSYS='DUR',
+ CPRESOPT='ZRESI',
+ LCOUPLING=T /
+
+&NAM_CONFZ MPI_BUFFER_SIZE=500/
+
+&NAM_VER_GRID NKMAX=68,
+ YZGRID_TYPE='FUNCTN',
+ ZDZGRD=4,
+ ZDZTOP=500,
+ ZZMAX_STRGRD=1000,
+ ZSTRGRD=17,
+ ZSTRTOP=4 /
diff --git a/MY_RUN/INTEGRATION_CASES/HPC/PANAME/005_coupling1/PRE_REAL1.nam_20230712_18H b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/005_coupling1/PRE_REAL1.nam_20230712_18H
new file mode 100644
index 0000000000000000000000000000000000000000..bfccc41f17f808a3701611fbfe9c1b0b07acd0ef
--- /dev/null
+++ b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/005_coupling1/PRE_REAL1.nam_20230712_18H
@@ -0,0 +1,23 @@
+&NAM_CONFIO LCDF4=T,
+ LLFIREAD=F,
+ NIO_ABORT_LEVEL=1,
+ NGEN_ABORT_LEVEL=1/
+
+&NAM_FILE_NAMES HATMFILE='arome.AN.20230712.18',
+ HATMFILETYPE='GRIBEX',
+ HPGDFILE='PGD_1200m.neste1',
+ CINIFILE='20230712_18H' /
+
+&NAM_REAL_CONF CEQNSYS='DUR',
+ CPRESOPT='ZRESI',
+ LCOUPLING=T /
+
+&NAM_CONFZ MPI_BUFFER_SIZE=500/
+
+&NAM_VER_GRID NKMAX=68,
+ YZGRID_TYPE='FUNCTN',
+ ZDZGRD=4,
+ ZDZTOP=500,
+ ZZMAX_STRGRD=1000,
+ ZSTRGRD=17,
+ ZSTRTOP=4 /
diff --git a/MY_RUN/INTEGRATION_CASES/HPC/PANAME/005_coupling1/PRE_REAL1.nam_20230712_21H b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/005_coupling1/PRE_REAL1.nam_20230712_21H
new file mode 100644
index 0000000000000000000000000000000000000000..489bc0666abb8da75c96357723400e9e13cd2d33
--- /dev/null
+++ b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/005_coupling1/PRE_REAL1.nam_20230712_21H
@@ -0,0 +1,23 @@
+&NAM_CONFIO LCDF4=T,
+ LLFIREAD=F,
+ NIO_ABORT_LEVEL=1,
+ NGEN_ABORT_LEVEL=1/
+
+&NAM_FILE_NAMES HATMFILE='arome.AN.20230712.21',
+ HATMFILETYPE='GRIBEX',
+ HPGDFILE='PGD_1200m.neste1',
+ CINIFILE='20230712_21H' /
+
+&NAM_REAL_CONF CEQNSYS='DUR',
+ CPRESOPT='ZRESI',
+ LCOUPLING=T /
+
+&NAM_CONFZ MPI_BUFFER_SIZE=500/
+
+&NAM_VER_GRID NKMAX=68,
+ YZGRID_TYPE='FUNCTN',
+ ZDZGRD=4,
+ ZDZTOP=500,
+ ZZMAX_STRGRD=1000,
+ ZSTRGRD=17,
+ ZSTRTOP=4 /
diff --git a/MY_RUN/INTEGRATION_CASES/HPC/PANAME/005_coupling1/PRE_REAL1.nam_20230713_00H b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/005_coupling1/PRE_REAL1.nam_20230713_00H
new file mode 100644
index 0000000000000000000000000000000000000000..341fb86695fc387553347c759415a94298523139
--- /dev/null
+++ b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/005_coupling1/PRE_REAL1.nam_20230713_00H
@@ -0,0 +1,23 @@
+&NAM_CONFIO LCDF4=T,
+ LLFIREAD=F,
+ NIO_ABORT_LEVEL=1,
+ NGEN_ABORT_LEVEL=1/
+
+&NAM_FILE_NAMES HATMFILE='arome.AN.20230713.00',
+ HATMFILETYPE='GRIBEX',
+ HPGDFILE='PGD_1200m.neste1',
+ CINIFILE='20230713_00H' /
+
+&NAM_REAL_CONF CEQNSYS='DUR',
+ CPRESOPT='ZRESI',
+ LCOUPLING=T /
+
+&NAM_CONFZ MPI_BUFFER_SIZE=500/
+
+&NAM_VER_GRID NKMAX=68,
+ YZGRID_TYPE='FUNCTN',
+ ZDZGRD=4,
+ ZDZTOP=500,
+ ZZMAX_STRGRD=1000,
+ ZSTRGRD=17,
+ ZSTRTOP=4 /
diff --git a/MY_RUN/INTEGRATION_CASES/HPC/PANAME/005_coupling1/PRE_REAL1_INIT.nam_20230712_00H b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/005_coupling1/PRE_REAL1_INIT.nam_20230712_00H
new file mode 100644
index 0000000000000000000000000000000000000000..d5250a36b7f24aa5a7cf0aaafdf936dc83dd88d7
--- /dev/null
+++ b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/005_coupling1/PRE_REAL1_INIT.nam_20230712_00H
@@ -0,0 +1,36 @@
+&NAM_CONFIO LCDF4=T,
+ LLFIREAD=F,
+ NIO_ABORT_LEVEL=1,
+ NGEN_ABORT_LEVEL=1/
+
+&NAM_FILE_NAMES HATMFILE='arome.AN.20230712.00',
+ HATMFILETYPE='GRIBEX',
+ HPGDFILE='PGD_1200m.neste1',
+ CINIFILE='20230712_00H' /
+
+&NAM_REAL_CONF CEQNSYS='DUR',
+ CPRESOPT='ZRESI',
+ LCOUPLING=F /
+
+&NAM_PREP_SURF_ATM CFILE='INIT_SURF.20230712.00',
+ CFILETYPE='MESONH',
+ CFILEPGD='PGD_oper_41t1.01km30',
+ CFILEPGDTYPE='MESONH',
+ NHALO_PREP=100/
+
+&NAM_CONFZ MPI_BUFFER_SIZE=500/
+
+&NAM_PREP_ISBA LISBA_CANOPY=F /
+
+&NAM_PREP_TEB LTEB_CANOPY = T / !Booleen qui ne sert plus par la suite, mais par soucis de clarte je mets a .T.
+
+&NAM_PREP_ISBA_SNOW CSNOW = '3-L' / !pas le choix avec DIF
+
+&NAM_VER_GRID NKMAX=68,
+ YZGRID_TYPE='FUNCTN',
+ ZDZGRD=4,
+ ZDZTOP=500,
+ ZZMAX_STRGRD=1000,
+ ZSTRGRD=17,
+ ZSTRTOP=4 /
+
diff --git a/MY_RUN/INTEGRATION_CASES/HPC/PANAME/005_coupling1/run_preal_xyz_1 b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/005_coupling1/run_preal_xyz_1
new file mode 100755
index 0000000000000000000000000000000000000000..f4c8e17889471d3bd487bfd7c88ed31f67248eed
--- /dev/null
+++ b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/005_coupling1/run_preal_xyz_1
@@ -0,0 +1,85 @@
+#!/bin/sh
+#SBATCH -J PRE_REAL1
+#SBATCH -N 2 # nodes number
+#SBATCH -n 16 # CPUs number (on all nodes)
+#SBATCH -o output_PRE_REAL.eo%j
+#SBATCH -e output_PRE_REAL.eo%j
+#SBATCH -t 0:45:00 # time limit
+#SBATCH --no-requeue # afin d'éviter qu'un job coupable d'un dépassement mémoire soit remis en queue infiniment
+#SBATCH -p normal256
+
+#N=7 n=7 00:45:00
+
+module purge
+
+ulimit -s unlimited
+ulimit -c 0
+
+# Arret du job des la premiere erreur
+set -e
+#Nom de la machine
+hostname
+
+
+#Date du jour
+AN=$(date +%Y)
+MM=$(date +%m)
+DD=$(date +%d)
+
+
+DATE=20230712
+DATEP1=$((DATE+1))
+# Profil utilisateur MesoNH
+. ~/MNH-V5-6-2/conf/profile_mesonh-LXifort-R8I4-MNH-V5-6-2-ECRAD140-MPIAUTO-O2-PANAME_FOR571
+
+
+# Import fichiers PGD nestes
+ln -sf ../004_nest_pgd/PGD*.neste1.* .
+
+DIR_FILES_ARO=../PHYSIO_AND_AROME_PANAME/AROME
+# Import fichiers AROME
+ln -sf $DIR_FILES_ARO/arome.AN.${DATE}* .
+ln -sf $DIR_FILES_ARO/INIT_SURF.${DATE}*.lfi .
+
+ln -sf $DIR_FILES_ARO/arome.AN.${DATEP1}* .
+ln -sf $DIR_FILES_ARO/INIT_SURF.${DATEP1}*.lfi .
+
+PREP_PGD_FILES=/home/cnrm_other/ge/mrmh/rodierq/SAVE/mesonh/PGD/
+ln -sf $PREP_PGD_FILES/PGD_oper_41t1.01km30.lfi
+
+#On cree des fichiers vides .des
+>PGD_oper_41t1.01km30.des
+>INIT_SURF.${DATE}.00.des
+
+#Execution
+export MPIRUN="Mpirun -np 16 "
+
+HEURES="03 06 09 12 15 18 21 24" #couplage 3h
+
+
+cp PRE_REAL1_INIT.nam_${DATE}_00H PRE_REAL1.nam
+time ${MPIRUN} PREP_REAL_CASE${XYZ}
+mv OUTPUT_LISTING0 OUTPUT_LISTING0_INIT.${DATE}
+#Nettoyage
+rm -f PRE_REAL1.nam file_for_xtransfer
+
+ln -sf $DIR_FILES_ARO/arome.AN.${DATEP1}* .
+
+for HH in ${HEURES}
+do
+DATE_HH=${DATE}_${HH}H
+if (($HH > 21 )); then
+HHP1=$(($HH-24))
+DATE_HH=${DATEP1}_0${HHP1}H
+fi
+cp PRE_REAL1.nam_${DATE_HH} PRE_REAL1.nam
+time ${MPIRUN} PREP_REAL_CASE${XYZ}
+mv OUTPUT_LISTING0 OUTPUT_LISTING0_COUPLING.${DATE}
+#Nettoyage
+rm -f PRE_REAL1.nam file_for_xtransfer
+done
+
+
+#job suivant :
+cd ../006_spawning12/
+sbatch run_spawning_xyz_1_2 ${DATE}
diff --git a/MY_RUN/INTEGRATION_CASES/HPC/PANAME/006_spawning12/SPAWN1.nam_1_2 b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/006_spawning12/SPAWN1.nam_1_2
new file mode 100644
index 0000000000000000000000000000000000000000..fcd27db728dcb131bfd0772d2cb0a82e6c56e6fe
--- /dev/null
+++ b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/006_spawning12/SPAWN1.nam_1_2
@@ -0,0 +1,8 @@
+
+&NAM_LUNIT2_SPA CINIFILE = '20230712_00H',
+ CINIFILEPGD='PGD_1200m.neste1',
+ YDOMAIN = 'PGD_300m.neste1',
+ YSPANBR = '01' /
+
+&NAM_SPAWN_SURF LSPAWN_SURF=F /
+
diff --git a/MY_RUN/INTEGRATION_CASES/HPC/PANAME/006_spawning12/creation_namelist_SPAWN.sh b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/006_spawning12/creation_namelist_SPAWN.sh
new file mode 100755
index 0000000000000000000000000000000000000000..6abf67fa7e20d89cc83936ff63cd631540ad4d33
--- /dev/null
+++ b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/006_spawning12/creation_namelist_SPAWN.sh
@@ -0,0 +1,26 @@
+#!/bin/ksh
+#Programme qui copie la namelist PRE REAL MesoNH de façon automatique pour toutes les échéances voulues
+#adapte du programme de Marc Mandement Février 2020
+#dans /scratch/work/mandementm/SIMU_CAS_AUDE_2MODELS/004_arome2lfi_preal1
+# Clotilde Augros septembre 2020
+# adapte pour simu CEP
+
+AN=$(date +%Y)
+MM=$(date +%m)
+DD=$(date +%d)
+
+#DD=30
+
+DATE=$1
+
+cat <<EOT | tr '"' "'" > SPAWN1.nam_1_2
+
+&NAM_LUNIT2_SPA CINIFILE = '${DATE}_00H',
+ CINIFILEPGD="PGD_1200m.neste1",
+ YDOMAIN = 'PGD_300m.neste1',
+ YSPANBR = '01' /
+
+&NAM_SPAWN_SURF LSPAWN_SURF=F /
+
+EOT
+
diff --git a/MY_RUN/INTEGRATION_CASES/HPC/PANAME/006_spawning12/run_spawning_xyz_1_2 b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/006_spawning12/run_spawning_xyz_1_2
new file mode 100755
index 0000000000000000000000000000000000000000..d4e8027be35747b20bd4f091be4d7d46ba61573a
--- /dev/null
+++ b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/006_spawning12/run_spawning_xyz_1_2
@@ -0,0 +1,38 @@
+#!/bin/sh
+#SBATCH -J SPAWN1
+#SBATCH -N 1 # nodes number
+#SBATCH -n 1 # CPUs number (on all nodes)
+#SBATCH -o spawn_1.eo%j
+#SBATCH -e spawn_1.eo%j
+#SBATCH -t 0:05:00 # time limit
+#SBATCH -p normal256
+module purge
+
+ulimit -s unlimited
+ulimit -c 0
+
+# Arret du job des la premiere erreur
+set -e
+# Echo des commandes
+set -x
+
+DATE=$1
+
+. creation_namelist_SPAWN.sh $DATE
+
+. ~/MNH-V5-6-2/conf/profile_mesonh-LXifort-R8I4-MNH-V5-6-2-ECRAD140-MPIAUTO-O2-PANAME_FOR571
+
+# Import fichiers PGD nestes
+ln -sf ../004_nest_pgd/PGD*.neste1.* .
+ln -sf ../005_coupling1/${DATE}* .
+
+export MPIRUN="Mpirun -np 1"
+
+cp SPAWN1.nam_1_2 SPAWN1.nam
+time ${MPIRUN} SPAWNING${XYZ}
+mv OUTPUT_LISTING0 OUTPUT_LISTING0_SPAWN_1_2
+rm SPAWN1.nam
+
+#job suivant
+cd ../007_coupling2/
+sbatch run_preal_xyz $DATE
diff --git a/MY_RUN/INTEGRATION_CASES/HPC/PANAME/007_coupling2/PRE_REAL1.nam_mod2 b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/007_coupling2/PRE_REAL1.nam_mod2
new file mode 100644
index 0000000000000000000000000000000000000000..9167942fda8f0fdc4a2bc7483ffe3eb50b7c5c85
--- /dev/null
+++ b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/007_coupling2/PRE_REAL1.nam_mod2
@@ -0,0 +1,46 @@
+&NAM_CONFIO LCDF4=T,
+ LLFIREAD=F,
+ NIO_ABORT_LEVEL=1,
+ NGEN_ABORT_LEVEL=1/
+
+&NAM_FILE_NAMES HATMFILE='20230712_00H.spa01',
+ HATMFILETYPE='MESONH',
+ HPGDFILE='PGD_300m.neste1',
+ CINIFILE='20230712_00H_300m' /
+
+&NAM_REAL_CONF LCOUPLING=.FALSE.,
+ NVERB=5,
+ CEQNSYS='DUR',
+ CPRESOPT='ZRESI',
+ LRES=.TRUE.,
+ XRES=1.E07 /
+
+&NAM_PREP_SURF_ATM CFILE='INIT_SURF.20230712.00',
+ CFILETYPE='MESONH',
+ CFILEPGD='PGD_oper_41t1.01km30',
+ NHALO_PREP=50/
+
+!ECMWF
+!&NAM_PREP_SURF_ATM CFILE = 'EC.AN.SFC.20230712.0000.grib',
+! CFILETYPE='GRIB',
+! NHALO_PREP=50/
+
+
+&NAM_PREP_ISBA LISBA_CANOPY = F
+/
+
+&NAM_PREP_ISBA_SNOW CSNOW = '3-L' / ! pas le choix avec DIF
+
+&NAM_PREP_TEB LTEB_CANOPY = T / !Booleen qui ne sert plus par la suite, mais par soucis de clarte je mets a .T.
+
+&NAM_CONFZ MPI_BUFFER_SIZE=800/
+
+&NAM_VER_GRID NKMAX=68,
+ YZGRID_TYPE='FUNCTN',
+ ZDZGRD=4,
+ ZDZTOP=500,
+ ZZMAX_STRGRD=1000,
+ ZSTRGRD=17,
+ ZSTRTOP=4 /
+
+
diff --git a/MY_RUN/INTEGRATION_CASES/HPC/PANAME/007_coupling2/copie_automatique_namelist_COUPLING_MNH.sh b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/007_coupling2/copie_automatique_namelist_COUPLING_MNH.sh
new file mode 100755
index 0000000000000000000000000000000000000000..6349c3393f19f66d1af43ede81da77280fbc55e7
--- /dev/null
+++ b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/007_coupling2/copie_automatique_namelist_COUPLING_MNH.sh
@@ -0,0 +1,64 @@
+#!/bin/ksh
+#Programme qui copie la namelist PRE REAL MesoNH de façon automatique pour toutes les échéances voulues
+#adapte du programme de Marc Mandement Février 2020
+#dans /scratch/work/mandementm/SIMU_CAS_AUDE_2MODELS/004_arome2lfi_preal1
+# Clotilde Augros septembre 2020
+# adapte pour simu CEP
+
+AN=$(date +%Y)
+MM=$(date +%m)
+DD=$(date +%d)
+
+#DD=30
+DATE=$1
+HH=00
+
+cat <<EOT | tr '"' "'" > PRE_REAL1.nam_mod2
+&NAM_CONFIO LCDF4=T,
+ LLFIREAD=F,
+ NIO_ABORT_LEVEL=1,
+ NGEN_ABORT_LEVEL=1/
+
+&NAM_FILE_NAMES HATMFILE='${DATE}_00H.spa01',
+ HATMFILETYPE='MESONH',
+ HPGDFILE='PGD_300m.neste1',
+ CINIFILE='${DATE}_00H_300m' /
+
+&NAM_REAL_CONF LCOUPLING=.FALSE.,
+ NVERB=5,
+ CEQNSYS='DUR',
+ CPRESOPT='ZRESI',
+ LRES=.TRUE.,
+ XRES=1.E07 /
+
+&NAM_PREP_SURF_ATM CFILE='INIT_SURF.${DATE}.00',
+ CFILETYPE='MESONH',
+ CFILEPGD='PGD_oper_41t1.01km30',
+ NHALO_PREP=50/
+
+!ECMWF
+!&NAM_PREP_SURF_ATM CFILE = "EC.AN.SFC.${DATE}.${HH}00.grib",
+! CFILETYPE="GRIB",
+! NHALO_PREP=50/
+
+
+&NAM_PREP_ISBA LISBA_CANOPY = F
+/
+
+&NAM_PREP_ISBA_SNOW CSNOW = "3-L" / ! pas le choix avec DIF
+
+&NAM_PREP_TEB LTEB_CANOPY = T / !Booleen qui ne sert plus par la suite, mais par soucis de clarte je mets a .T.
+
+&NAM_CONFZ MPI_BUFFER_SIZE=800/
+
+&NAM_VER_GRID NKMAX=68,
+ YZGRID_TYPE='FUNCTN',
+ ZDZGRD=4,
+ ZDZTOP=500,
+ ZZMAX_STRGRD=1000,
+ ZSTRGRD=17,
+ ZSTRTOP=4 /
+
+
+EOT
+
diff --git a/MY_RUN/INTEGRATION_CASES/HPC/PANAME/007_coupling2/run_preal_xyz b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/007_coupling2/run_preal_xyz
new file mode 100755
index 0000000000000000000000000000000000000000..a56055853ffaef35c98032baa962222073dc7903
--- /dev/null
+++ b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/007_coupling2/run_preal_xyz
@@ -0,0 +1,46 @@
+#!/bin/sh
+#SBATCH -J PREP2
+#SBATCH -N 1 # nodes number
+#SBATCH -n 10 # CPUs number (on all nodes)
+#SBATCH -o real_2.eo%j
+#SBATCH -e real_2.eo%j
+#SBATCH -t 0:45:00 # time limit
+#SBATCH -p normal256
+module purge
+
+ulimit -s unlimited
+ulimit -c 0
+
+# Arret du job des la premiere erreur
+set -e
+
+#Date du jour
+DATE=$1
+
+#Creation namelist
+. copie_automatique_namelist_COUPLING_MNH.sh
+. ~/MNH-V5-6-2/conf/profile_mesonh-LXifort-R8I4-MNH-V5-6-2-ECRAD140-MPIAUTO-O2-PANAME_FOR571
+
+# Import fichiers PGD nestes
+ln -sf ../004_nest_pgd/PGD*.neste1.* .
+ln -sf ../005_coupling1/INIT_SURF.${DATE}.00* .
+ln -sf ../005_coupling1/PGD_oper_41t1.01km30* .
+
+# Import fichier initial
+ln -sf ../005_coupling1/${DATE}_00H* .
+
+# Import fichier spawning
+ln -sf ../006_spawning12/${DATE}_00H.spa01.* .
+
+cp PRE_REAL1.nam_mod2 PRE_REAL1.nam
+
+# Exécution
+export MPIRUN="Mpirun -np 10 "
+time ${MPIRUN} PREP_REAL_CASE${XYZ}
+
+# Nettoyage
+rm -f file_for_xtransfer PRE_REAL1.nam
+
+# Etape suivante
+cd ../008_spawning23/
+sbatch run_spawning_xyz_2_3 $DATE
diff --git a/MY_RUN/INTEGRATION_CASES/HPC/PANAME/008_spawning23/SPAWN1.nam_2_3 b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/008_spawning23/SPAWN1.nam_2_3
new file mode 100644
index 0000000000000000000000000000000000000000..753ecc0d63f5b27dc686c1e19a588568b6790612
--- /dev/null
+++ b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/008_spawning23/SPAWN1.nam_2_3
@@ -0,0 +1,8 @@
+
+&NAM_LUNIT2_SPA CINIFILE = '20230712_00H_300m',
+ CINIFILEPGD='PGD_300m.neste1',
+ YDOMAIN = 'PGD_100m.neste1',
+ YSPANBR = '01' /
+
+&NAM_SPAWN_SURF LSPAWN_SURF=F /
+
diff --git a/MY_RUN/INTEGRATION_CASES/HPC/PANAME/008_spawning23/creation_namelist_SPAWN.sh b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/008_spawning23/creation_namelist_SPAWN.sh
new file mode 100755
index 0000000000000000000000000000000000000000..d1c74e9dd9c8d9dd34df394eb81f0c954bd2bab6
--- /dev/null
+++ b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/008_spawning23/creation_namelist_SPAWN.sh
@@ -0,0 +1,24 @@
+#!/bin/ksh
+#Programme qui copie la namelist PRE REAL MesoNH de façon automatique pour toutes les échéances voulues
+#adapte du programme de Marc Mandement Février 2020
+#dans /scratch/work/mandementm/SIMU_CAS_AUDE_2MODELS/004_arome2lfi_preal1
+# Clotilde Augros septembre 2020
+# adapte pour simu CEP
+
+AN=$(date +%Y)
+MM=$(date +%m)
+DD=$(date +%d)
+
+DATE=$1
+
+cat <<EOT | tr '"' "'" > SPAWN1.nam_2_3
+
+&NAM_LUNIT2_SPA CINIFILE = '${DATE}_00H_300m',
+ CINIFILEPGD="PGD_300m.neste1",
+ YDOMAIN = 'PGD_100m.neste1',
+ YSPANBR = '01' /
+
+&NAM_SPAWN_SURF LSPAWN_SURF=F /
+
+EOT
+
diff --git a/MY_RUN/INTEGRATION_CASES/HPC/PANAME/008_spawning23/run_spawning_xyz_2_3 b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/008_spawning23/run_spawning_xyz_2_3
new file mode 100755
index 0000000000000000000000000000000000000000..9249a1ab55b1a08e1cfbda6920a49088944b8ccf
--- /dev/null
+++ b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/008_spawning23/run_spawning_xyz_2_3
@@ -0,0 +1,48 @@
+#!/bin/sh
+#
+#SBATCH -J SPAWN23
+#SBATCH -e spawn_1.eo%j
+#SBATCH -N 1 # nodes number
+#SBATCH -n 1 # CPUs number (on all nodes)
+#SBATCH -o spawn_1.eo%j
+#SBATCH -e spawn_1.eo%j
+#SBATCH -t 0:05:00 # time limit
+#SBATCH --no-requeue # afin d'éviter qu'un job coupable d'un dépassement mémoire soit remis en queue infiniment
+#SBATCH -p normal256
+
+module purge
+
+ulimit -s unlimited
+ulimit -c 0
+
+# Arret du job des la premiere erreur
+set -e
+
+# Echo des commandes
+set -x
+
+
+DATE=$1
+
+. creation_namelist_SPAWN.sh $DATE
+
+#. ~/MNH-V5-6-2/conf/profile_mesonh-LXifort-R8I4-MNH-V5-6-2-MPIAUTO-O2-PANAME_FOR571
+. ~/MNH-V5-6-2/conf/profile_mesonh-LXifort-R8I4-MNH-V5-6-2-ECRAD140-MPIAUTO-O2-PANAME_FOR571
+
+
+# Import fichiers PGD nestes
+ln -sf ../004_nest_pgd/PGD*.neste1.* .
+
+ln -sf ../007_coupling2/${DATE}_00H_300m* .
+
+export MPIRUN="Mpirun -np 1"
+
+cp SPAWN1.nam_2_3 SPAWN1.nam
+time ${MPIRUN} SPAWNING${XYZ}
+mv OUTPUT_LISTING0 OUTPUT_LISTING0_SPAWN_2_3
+rm SPAWN1.nam
+
+
+#job suivant
+cd ../009_coupling3/
+sbatch run_preal_xyz $DATE
diff --git a/MY_RUN/INTEGRATION_CASES/HPC/PANAME/009_coupling3/PRE_REAL1.nam_mod3 b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/009_coupling3/PRE_REAL1.nam_mod3
new file mode 100644
index 0000000000000000000000000000000000000000..4d209f6180a724c655e04db091d2472c45b5cf31
--- /dev/null
+++ b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/009_coupling3/PRE_REAL1.nam_mod3
@@ -0,0 +1,45 @@
+&NAM_CONFIO LCDF4=T,
+ LLFIREAD=F,
+ NIO_ABORT_LEVEL=1,
+ NGEN_ABORT_LEVEL=1/
+
+&NAM_FILE_NAMES HATMFILE='20230712_00H_300m.spa01',
+ HATMFILETYPE='MESONH',
+ HPGDFILE='PGD_100m.neste1',
+ CINIFILE='20230712_00H_100m' /
+
+&NAM_REAL_CONF LCOUPLING=.FALSE.,
+ NVERB=5,
+ CEQNSYS='DUR',
+ CPRESOPT='ZRESI',
+ LRES=.TRUE.,
+ XRES=1.E07 /
+
+!&NAM_PREP_SURF_ATM CFILE = 'EC.AN.SFC.20230712.0000.grib',
+! CFILETYPE='GRIB',
+! NHALO_PREP=50/
+
+
+&NAM_PREP_SURF_ATM CFILE='INIT_SURF.20230712.00',
+ CFILETYPE='MESONH',
+ CFILEPGD='PGD_oper_41t1.01km30',
+ NHALO_PREP=50/
+
+&NAM_PREP_ISBA LISBA_CANOPY = F
+/
+
+&NAM_PREP_ISBA_SNOW CSNOW = '3-L' / !Oblige avec DIF
+
+&NAM_PREP_TEB LTEB_CANOPY = T / !Booleen qui ne sert plus par la suite, mais par soucis de clarte je mets a .T.
+
+&NAM_CONFZ MPI_BUFFER_SIZE=800/
+
+&NAM_VER_GRID NKMAX=68,
+ YZGRID_TYPE='FUNCTN',
+ ZDZGRD=4,
+ ZDZTOP=500,
+ ZZMAX_STRGRD=1000,
+ ZSTRGRD=17,
+ ZSTRTOP=4 /
+
+
diff --git a/MY_RUN/INTEGRATION_CASES/HPC/PANAME/009_coupling3/copie_automatique_namelist_COUPLING_MNH.sh b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/009_coupling3/copie_automatique_namelist_COUPLING_MNH.sh
new file mode 100755
index 0000000000000000000000000000000000000000..9ca45a9676136dee5b396696d5cf0245b39fb7b4
--- /dev/null
+++ b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/009_coupling3/copie_automatique_namelist_COUPLING_MNH.sh
@@ -0,0 +1,65 @@
+#!/bin/ksh
+#Programme qui copie la namelist PRE REAL MesoNH de façon automatique pour toutes les échéances voulues
+#adapte du programme de Marc Mandement Février 2020
+#dans /scratch/work/mandementm/SIMU_CAS_AUDE_2MODELS/004_arome2lfi_preal1
+# Clotilde Augros septembre 2020
+# adapte pour simu CEP
+
+AN=$(date +%Y)
+MM=$(date +%m)
+DD=$(date +%d)
+
+#DD=30
+#MM=05
+HH=00
+
+DATE=$1
+
+cat <<EOT | tr '"' "'" > PRE_REAL1.nam_mod3
+&NAM_CONFIO LCDF4=T,
+ LLFIREAD=F,
+ NIO_ABORT_LEVEL=1,
+ NGEN_ABORT_LEVEL=1/
+
+&NAM_FILE_NAMES HATMFILE='${DATE}_00H_300m.spa01',
+ HATMFILETYPE='MESONH',
+ HPGDFILE='PGD_100m.neste1',
+ CINIFILE='${DATE}_00H_100m' /
+
+&NAM_REAL_CONF LCOUPLING=.FALSE.,
+ NVERB=5,
+ CEQNSYS='DUR',
+ CPRESOPT='ZRESI',
+ LRES=.TRUE.,
+ XRES=1.E07 /
+
+!&NAM_PREP_SURF_ATM CFILE = "EC.AN.SFC.${DATE}.${HH}00.grib",
+! CFILETYPE="GRIB",
+! NHALO_PREP=50/
+
+
+&NAM_PREP_SURF_ATM CFILE='INIT_SURF.${DATE}.00',
+ CFILETYPE='MESONH',
+ CFILEPGD='PGD_oper_41t1.01km30',
+ NHALO_PREP=50/
+
+&NAM_PREP_ISBA LISBA_CANOPY = F
+/
+
+&NAM_PREP_ISBA_SNOW CSNOW = "3-L" / !Oblige avec DIF
+
+&NAM_PREP_TEB LTEB_CANOPY = T / !Booleen qui ne sert plus par la suite, mais par soucis de clarte je mets a .T.
+
+&NAM_CONFZ MPI_BUFFER_SIZE=800/
+
+&NAM_VER_GRID NKMAX=68,
+ YZGRID_TYPE='FUNCTN',
+ ZDZGRD=4,
+ ZDZTOP=500,
+ ZZMAX_STRGRD=1000,
+ ZSTRGRD=17,
+ ZSTRTOP=4 /
+
+
+EOT
+
diff --git a/MY_RUN/INTEGRATION_CASES/HPC/PANAME/009_coupling3/run_preal_xyz b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/009_coupling3/run_preal_xyz
new file mode 100755
index 0000000000000000000000000000000000000000..8ad24102c9e9a5b766d278c7aa8da34b44e18417
--- /dev/null
+++ b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/009_coupling3/run_preal_xyz
@@ -0,0 +1,62 @@
+#!/bin/sh
+#SBATCH -J PREP3
+#SBATCH -N 1 # nodes number
+#SBATCH -n 20 # CPUs number (on all nodes)
+#SBATCH -o real_2_CEP.eo%j
+#SBATCH -e real_2_CEP.eo%j
+#SBATCH -t 0:45:00 # time limit
+#SBATCH -p normal256
+#SBATCH --no-requeue
+
+module purge
+
+ulimit -s unlimited
+ulimit -c 0
+
+# Arret du job des la premiere erreur
+set -e
+#Nom de la machine
+hostname
+
+
+
+HH=00
+
+
+DATE=$1
+
+#Creation namelist
+
+. copie_automatique_namelist_COUPLING_MNH.sh $DATE
+
+
+#. ~/MNH-V5-6-2/conf/profile_mesonh-LXifort-R8I4-MNH-V5-6-2-MPIAUTO-O2-PANAME_FOR571 #profile_mesonh-LXifort-R8I4-MNH-V5-6-2-MPIAUTO-O2-PANAME_FOR571
+. ~/MNH-V5-6-2/conf/profile_mesonh-LXifort-R8I4-MNH-V5-6-2-ECRAD140-MPIAUTO-O2-PANAME_FOR571
+
+
+# Import fichiers PGD nestes
+ln -sf ../004_nest_pgd/PGD*.neste1.* .
+
+ln -sf ../005_coupling1/INIT_SURF.${DATE}.00* .
+ln -sf ../005_coupling1/PGD_oper_41t1.01km30* .
+
+# Import fichier initial
+ln -sf ../007_coupling2/${DATE}_00H_300m* .
+ln -sf ../005_coupling1/${DATE}_00H* .
+
+# Import fichier spawning
+ln -sf ../008_spawning23/${DATE}_00H_300m.spa01.* .
+ln -sf ../006_spawning12/${DATE}_00H.spa01.* .
+
+cp PRE_REAL1.nam_mod3 PRE_REAL1.nam
+
+# Exécution
+export MPIRUN="Mpirun -np 20 "
+time ${MPIRUN} PREP_REAL_CASE${XYZ}
+
+# Nettoyage
+rm -f file_for_xtransfer PRE_REAL1.nam
+
+# Etape suivante
+cd ../010_run/
+sbatch run_mesonh_xyz_mod123 $DATE
diff --git a/MY_RUN/INTEGRATION_CASES/HPC/PANAME/010_run/EXSEG1.nam b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/010_run/EXSEG1.nam
new file mode 100644
index 0000000000000000000000000000000000000000..139377222a3eb90ca61d7ac22397e2f5790739c3
--- /dev/null
+++ b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/010_run/EXSEG1.nam
@@ -0,0 +1,140 @@
+&NAM_CONFIO LCDF4=T,
+ NIO_ABORT_LEVEL=1,
+ NGEN_ABORT_LEVEL=1 /
+
+&NAM_LUNITn CINIFILE = '20230712_00H',
+ CINIFILEPGD = 'PGD_1200m.neste1'
+ CCPLFILE(1) = '20230712_03H',
+ CCPLFILE(2) = '20230712_06H',
+ CCPLFILE(3) = '20230712_09H',
+ CCPLFILE(4) = '20230712_12H',
+ CCPLFILE(5) = '20230712_15H',
+ CCPLFILE(6) = '20230712_18H',
+ CCPLFILE(7) = '20230712_21H',
+ CCPLFILE(8) = '20230713_00H'
+/
+
+&NAM_DYN XSEGLEN = 86400.
+ XALKTOP = 0.01,
+ XALZBOT = 15000.,
+ LNUMDIFU = F, ! a desactiver quand on utilise WENO_K
+ LNUMDIFTH = F / ! si CMET_ADV_SCHEME = PPM_01 on peut desactiver la diffusion num pour varia meteo
+
+&NAM_BACKUP XBAK_TIME_FREQ(1) = 3600.,
+ XBAK_TIME_FREQ(2) = 3600.,
+ XBAK_TIME_FREQ(3) = 3600.,
+ XBAK_TIME_FREQ_FIRST(1) = 3600.,
+ XBAK_TIME_FREQ_FIRST(2) = 3600.,
+ XBAK_TIME_FREQ_FIRST(3) = 3600.,/
+
+&NAM_NESTING NDAD(2)=1,
+ NDAD(3)=2,
+ NDTRATIO(2)=4,
+ NDTRATIO(3)=3,
+ XWAY(2)=2 /
+ XWAY(3)=2 /
+
+&NAM_CONFZ MPI_BUFFER_SIZE=100,
+ NZ_VERB=1 /
+
+&NAM_CONF CEXP = 'PAN.00',
+ CCONF= 'START',
+ CSEG = 'SEG01',
+ NMODEL = 3,
+ NVERB = 1,
+ NHALO=3
+/
+
+&NAM_DYNn XTSTEP = 25.,
+ CPRESOPT = 'ZRESI',
+ NITR = 5,
+ LRES = T,
+ XRES = 1.E-07,
+ LITRADJ=T,
+ LHORELAX_UVWTH = F,
+ LHORELAX_RV = F,
+ LHORELAX_RC = F,
+ LHORELAX_RR = F,
+ LHORELAX_RS = F,
+ LHORELAX_RI = F,
+ LHORELAX_RG = F,
+ LVE_RELAX = T,
+ NRIMX = 5,
+ NRIMY = 5,
+ XRIMKMAX = 0.002 /
+
+NAM_COUPLING_LEVELSn NLEV_COUPLE=16 /
+
+&NAM_ADVn CUVW_ADV_SCHEME = 'WENO_K',
+ CMET_ADV_SCHEME = 'PPM_01',
+ CSV_ADV_SCHEME = 'PPM_01',
+ CTEMP_SCHEME = 'RK53',
+ NWENO_ORDER = 3,
+ LSPLIT_CFL = .TRUE.,
+ LSPLIT_WENO = .TRUE.,
+ XSPLIT_CFL = 0.7,
+ LCFL_WRIT = .FALSE. /
+
+&NAM_PARAMn CTURB = 'TKEL',
+ CRAD = 'ECMW',
+ CCLOUD = 'ICE3',
+ CDCONV = 'NONE',
+ CSCONV = 'EDKF'/
+
+&NAM_PARAM_RADn XDTRAD = 900.,
+ XDTRAD_CLONLY = 300.,
+ CAER = 'TEGE'
+/
+
+
+&NAM_LBCn CLBCX = 2*'OPEN',
+ CLBCY = 2*'OPEN' /
+
+&NAM_TURBn CTURBLEN = 'HM21',
+ CTURBDIM = '1DIM',
+
+
+ LLEONARD=.TRUE.,
+ XCOEFHGRADTHL=4.0,
+ XCOEFHGRADRM=4.0,
+ XCLDTHOLD=-1.0,
+ XALTHGRAD=2000.0 /
+&NAM_SURF_ATM
+XRIMAX=0
+LVERTSHIFT=F
+/
+
+&NAM_SEAFLUXn CSEA_FLUX = 'COARE3',
+ LPRECIP = T,
+ LPWEBB = T /
+
+&NAM_DIAG_SURFn N2M = 2,
+ LSURF_BUDGET = T,
+ LRAD_BUDGET = T,
+ LCOEF = F,
+ LSURF_VARS = F /
+
+&NAM_DIAG_SURF_ATMn LFRAC = F /
+
+
+&NAM_DIAG_ISBAn LPGD = F,
+ LSURF_EVAP_BUDGET = F,
+ LSURF_MISC_BUDGET = F,
+ LSURF_DIAG_ALBEDO = F,
+ LPATCH_BUDGET = F /
+
+&NAM_DIAG_TEBn
+ LPGD = F,
+ LSURF_DIAG_ALBEDO = F ,
+ LUTCI = F,
+ LSURF_MISC_BUDGET=F /
+
+
+
+&NAM_PARAM_ICEn CSUBG_AUCV_RC = 'CLFR', LRED=T,
+LSNOW_T=T
+/
+
+
+&NAM_NEBn LSUBG_COND = T /
+
diff --git a/MY_RUN/INTEGRATION_CASES/HPC/PANAME/010_run/EXSEG2.nam b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/010_run/EXSEG2.nam
new file mode 100644
index 0000000000000000000000000000000000000000..bf93ce5c6d532b96bf38d591c350bc8c742d7a20
--- /dev/null
+++ b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/010_run/EXSEG2.nam
@@ -0,0 +1,104 @@
+&NAM_LUNITn CINIFILE = '20230712_00H_300m',
+ CINIFILEPGD = 'PGD_300m.neste1' /
+
+&NAM_DYNn CPRESOPT = 'ZRESI',
+ NITR = 5,
+ LRES = T,
+ XRES = 1.E-07,
+ LITRADJ=T,
+ LHORELAX_UVWTH = F,
+ LHORELAX_RV = F,
+ LHORELAX_RC = F,
+ LHORELAX_RR = F,
+ LHORELAX_RS = F,
+ LHORELAX_RI = F,
+ LHORELAX_RG = F,
+ LVE_RELAX = T,
+ NRIMX = 5,
+ NRIMY = 5,
+ XRIMKMAX = 0.002 /
+&NAM_DYN
+ LNUMDIFTH=F/ ! conseille a false pour WENO
+
+&NAM_COUPLING_LEVELSn NLEV_COUPLE=16 / !20 nveaux pour Valery
+
+&NAM_ADVn CUVW_ADV_SCHEME = 'WENO_K',
+ CMET_ADV_SCHEME = 'PPM_01',
+ CSV_ADV_SCHEME = 'PPM_01',
+ CTEMP_SCHEME = 'RK53',
+ NWENO_ORDER = 3,
+ LSPLIT_CFL = .TRUE.,
+ LSPLIT_WENO = .TRUE.,
+ XSPLIT_CFL = 0.7,
+ LCFL_WRIT = .FALSE. /
+
+
+&NAM_PARAMn CTURB = 'TKEL',
+ CRAD = 'ECMW',
+ CCLOUD = 'ICE3',
+ CDCONV = 'NONE',
+ CSCONV = 'NONE'/
+
+&NAM_PARAM_RADn XDTRAD = 900.,
+ XDTRAD_CLONLY = 300.,
+ CAER = 'TEGE',
+ NRAD_AGG = 3
+/
+
+
+
+
+&NAM_LBCn CLBCX = 2*'OPEN',
+ CLBCY = 2*'OPEN' /
+
+&NAM_TURBn CTURBLEN = 'HM21',
+ CTURBDIM = '3DIM',
+
+
+ LLEONARD=.TRUE.,
+ XCOEFHGRADTHL=4.0,
+ XCOEFHGRADRM=4.0,
+ XCLDTHOLD=-1.0,
+ XALTHGRAD=2000.0 /
+
+&NAM_SURF_ATM
+XRIMAX=0
+LVERTSHIFT=F
+/
+
+&NAM_SEAFLUXn CSEA_FLUX = 'COARE3',
+ LPRECIP = T,
+ LPWEBB = T /
+
+&NAM_DIAG_SURFn N2M = 2,
+ LSURF_BUDGET = T,
+ LRAD_BUDGET = T,
+ LCOEF = F,
+ LSURF_VARS = F
+
+/
+
+
+&NAM_DIAG_SURF_ATMn LFRAC = F /
+
+&NAM_PARAM_ICEn CSUBG_AUCV_RC = 'CLFR', LRED=T,
+LSNOW_T=T
+ /
+
+&NAM_DIAG_ISBAn LPGD = F,
+ LSURF_EVAP_BUDGET = F,
+ LSURF_MISC_BUDGET = F,
+ LSURF_DIAG_ALBEDO = F,
+ LPATCH_BUDGET = F /
+
+&NAM_DIAG_TEBn
+ LPGD = F,
+ LSURF_DIAG_ALBEDO = F ,
+ LUTCI = F,
+ LSURF_MISC_BUDGET=T /
+
+
+
+
+&NAM_NEBn LSUBG_COND = T /
+
diff --git a/MY_RUN/INTEGRATION_CASES/HPC/PANAME/010_run/EXSEG3.nam b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/010_run/EXSEG3.nam
new file mode 100644
index 0000000000000000000000000000000000000000..11b1bd1307cd244078e2d97c5c39a128b68a1172
--- /dev/null
+++ b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/010_run/EXSEG3.nam
@@ -0,0 +1,96 @@
+&NAM_LUNITn CINIFILE = '20230712_00H_100m',
+ CINIFILEPGD = 'PGD_100m.neste1' /
+
+&NAM_DYNn CPRESOPT = 'ZRESI',
+ NITR = 5,
+ LRES = T,
+ XRES = 1.E-07, !-6
+ LITRADJ=T,
+ LHORELAX_UVWTH = F,
+ LHORELAX_RV = F,
+ LHORELAX_RC = F,
+ LHORELAX_RR = F,
+ LHORELAX_RS = F,
+ LHORELAX_RI = F,
+ LHORELAX_RG = F,
+ LVE_RELAX = T,
+ NRIMX = 5,
+ NRIMY = 5,
+ XRIMKMAX = 0.002 /
+&NAM_DYN
+ LNUMDIFU=T
+ LNUMDIFTH=F /
+
+&NAM_COUPLING_LEVELSn NLEV_COUPLE=19 / !20 nveaux pour Valery
+
+&NAM_ADVn CUVW_ADV_SCHEME = 'CEN4TH', !'WENO_K',
+ CMET_ADV_SCHEME = 'PPM_01',
+ CSV_ADV_SCHEME = 'PPM_01',
+ CTEMP_SCHEME = 'RKC4',
+ LSPLIT_CFL = .TRUE.,
+ XSPLIT_CFL = 0.9,
+ LCFL_WRIT = .FALSE. /
+
+&NAM_PARAMn CTURB = 'TKEL',
+ CRAD = 'ECMW',
+ CCLOUD = 'ICE3',
+ CDCONV = 'NONE',
+ CSCONV = 'NONE'/
+
+&NAM_PARAM_RADn XDTRAD = 900.,
+ XDTRAD_CLONLY = 300.,
+ CAER = 'TEGE',
+ NRAD_AGG = 3
+/
+
+
+&NAM_LBCn CLBCX = 2*'OPEN',
+ CLBCY = 2*'OPEN' /
+
+&NAM_TURBn CTURBLEN = 'HM21',
+ CTURBDIM = '3DIM',
+
+
+ LLEONARD=.TRUE.,
+ XCOEFHGRADTHL=4.0,
+ XCOEFHGRADRM=4.0,
+ XCLDTHOLD=-1.0,
+ XALTHGRAD=2000.0 /
+&NAM_SURF_ATM
+XRIMAX=0
+LVERTSHIFT=F
+/
+
+&NAM_SEAFLUXn CSEA_FLUX = 'COARE3',
+ LPRECIP = T,
+ LPWEBB = T /
+
+&NAM_DIAG_SURFn N2M = 2,
+ LSURF_BUDGET = T,
+ LRAD_BUDGET = T,
+ LCOEF = F,
+ LSURF_VARS = F /
+
+&NAM_DIAG_SURF_ATMn LFRAC = F /
+
+&NAM_PARAM_ICEn CSUBG_AUCV_RC = 'CLFR', LRED=T,
+LSNOW_T=T
+ /
+
+
+&NAM_DIAG_ISBAn LPGD = F,
+ LSURF_EVAP_BUDGET = F,
+ LSURF_MISC_BUDGET = F,
+ LSURF_DIAG_ALBEDO = F,
+ LPATCH_BUDGET = F /
+
+&NAM_DIAG_TEBn
+ LPGD = F,
+ LSURF_DIAG_ALBEDO = F ,
+ LUTCI = F,
+ LSURF_MISC_BUDGET=T/
+
+
+&NAM_NEBn LSUBG_COND = T /
+
+
diff --git a/MY_RUN/INTEGRATION_CASES/HPC/PANAME/010_run/run_mesonh_xyz_mod123 b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/010_run/run_mesonh_xyz_mod123
new file mode 100755
index 0000000000000000000000000000000000000000..fbbc08b89c7c074babf44ad0f9b56eb4993ae355
--- /dev/null
+++ b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/010_run/run_mesonh_xyz_mod123
@@ -0,0 +1,43 @@
+#!/bin/sh
+#SBATCH -J MESONH
+#SBATCH -N 11 # nodes number
+#SBATCH -n 1400 # CPUs number (on all nodes)
+#SBATCH -o RUN_PANAME.eo%j
+#SBATCH -e RUN_PANAME.eo%j
+#SBATCH -t 12:00:00
+
+# Arret du job des la premiere erreur (optionnel)
+set -e
+# Echo des commandes (optionnel)
+# set -x
+
+# Profil utilisateur MesoNH (obligatoire)
+. ~/MNH-V5-6-2/conf/profile_mesonh-LXifort-R8I4-MNH-V5-6-2-ECRAD140-MPIAUTO-O2-PANAME_FOR571
+
+# Import fichiers PGD nestes
+ln -sf ../004_nest_pgd/PGD*.neste1.* .
+
+# Import fichiers MESONH initiaux + couplage (obligatoire)
+DATE=20230712
+DATEP1=$(( DATE + 1 ))
+
+# ----- pour EXSEG1 (1200m)
+# initialisation et couplage
+ln -sf ../005_coupling1/${DATE}_??H.* .
+ln -sf ../005_coupling1/${DATEP1}_??H.* .
+
+#couplage
+# ----- pour EXSEG2 (300m)
+ln -sf ../007_coupling2/${DATE}_00H_300m.* .
+
+#couplage
+# ----- pour EXSEG2 (100m)
+ln -sf ../009_coupling3/${DATE}_00H_100m.* .
+
+#ECRAD
+ln -sf ${SRC_MESONH}/MNH-V5-6-2/src/LIB/RAD/ecrad-1.4.0/data/RAD* .
+ln -sf ${SRC_MESONH}/src/LIB/RAD/ecrad-1.4.0/data/*nc .
+
+# Exécution (obligatoire)
+export MPIRUN="Mpirun -np 1400 "
+time ${MPIRUN} MESONH${XYZ}
diff --git a/MY_RUN/INTEGRATION_CASES/HPC/PANAME/011_python/plot_PANAME_mod1.py b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/011_python/plot_PANAME_mod1.py
new file mode 100644
index 0000000000000000000000000000000000000000..25742d67e3b80e9acc04e81bce2695563a6e9926
--- /dev/null
+++ b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/011_python/plot_PANAME_mod1.py
@@ -0,0 +1,379 @@
+#!/usr/bin/env python3
+"""
+@author: Quentin Rodier
+Creation : 07/01/2021
+
+Last modifications
+"""
+
+import matplotlib as mpl
+mpl.use('Agg')
+from read_MNHfile import read_netcdf
+from Panel_Plot import PanelPlot
+import cartopy.crs as ccrs
+import numpy as np
+import os
+import cartopy.io.shapereader as shpreader
+
+os.system('rm -f tempgraph*')
+#
+# User's parameter / Namelist
+#
+path="../010_run/"
+LnameFiles = ['PAN.0.1.SEG01.004.nc', 'PAN.0.1.SEG01.008.nc', 'PAN.0.1.SEG01.012.nc', 'PAN.0.1.SEG01.016.nc', 'PAN.0.1.SEG01.020.nc', 'PAN.0.1.SEG01.024.nc']
+
+Dvar_input = {'f'+str(i):['T2M', 'TSRAD', 'HU2M','SWD','LWD','UT', 'CLDFR', 'TKET','WT', 'H','LE','level', 'longitude','latitude','time']
+ for i in range(1,7)}
+
+
+# Read the variables in the files
+Dvar = {}
+Dvar = read_netcdf(LnameFiles, Dvar_input, path=path, removeHALO=True)
+
+# Departements francais
+departements_shp='departements-20180101.shp'
+adm1_shapes = list(shpreader.Reader(departements_shp).geometries())
+
+titlepad=25
+minmaxpad=1.04
+timepad=-0.10
+colorbarpad=0.01
+labelcolorbarpad = 13
+colorbaraspect=20
+
+################################################################
+######### PANEL 1 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'H', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['H'][:,:] for i in range(1,7)]
+Lcbarlabel = ['W/m2']*len(Lplot)
+Lminval = [-100]*len(Lplot)
+Lmaxval = [1000]*len(Lplot)
+Lstep = [25]*len(Lplot)
+Lstepticks = [100]*len(Lplot)
+Lfacconv = [1]*len(Lplot)
+Lcolormap = ['rainbow']*len(Lplot)
+Ltime = [Dvar['f'+str(i)]['time']/3600.0 for i in range(1,7)]
+LaddWhite_cm = [True]*len(Lplot)
+LaxeX = [Dvar['f1']['longitude']]*len(Lplot)
+LaxeY = [Dvar['f1']['latitude']]*len(Lplot)
+Ltitle = ['']*len(Lplot)
+Lxlab = ['longitude']*len(Lplot)
+Lylab = ['latitude']*len(Lplot)
+Lpltype = ['cf']*len(Lplot)
+Lprojection = [ccrs.PlateCarree()]*len(Lplot)
+
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection, gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='gray', alpha=0.1)
+Panel1.save_graph(1, fig)
+
+################################################################
+######### PANEL 2 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'LE', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['LE'][:,:] for i in range(1,7)]
+Lcbarlabel = ['W/m2']*len(Lplot)
+Lminval = [0]*len(Lplot)
+Lmaxval = [600]*len(Lplot)
+Lstep = [25]*len(Lplot)
+Lstepticks = [100]*len(Lplot)
+Lfacconv = [1]*len(Lplot)
+Lcolormap = ['rainbow']*len(Lplot)
+LaddWhite_cm = [True]*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='gray', alpha=0.1)
+Panel1.save_graph(2, fig)
+
+
+################################################################
+######### PANEL 3 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'T2M', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['T2M'][:,:]-273.15 for i in range(1,7)]
+Lcbarlabel = ['°C']*len(Lplot)
+Lminval = [10]*len(Lplot)
+Lmaxval = [30]*len(Lplot)
+Lstep = [0.5]*len(Lplot)
+Lstepticks = [2.5]*len(Lplot)
+Lfacconv = [1]*len(Lplot)
+Lcolormap = ['rainbow']*len(Lplot)
+LaddWhite_cm = [False]*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='gray', alpha=0.1)
+Panel1.save_graph(3, fig)
+
+
+################################################################
+######### PANEL 4 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'HU2M', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['HU2M'][:,:]*100.0 for i in range(1,7)]
+Lcbarlabel = ['%']*len(Lplot)
+Lminval = [10]*len(Lplot)
+Lmaxval = [100]*len(Lplot)
+Lstep = [2]*len(Lplot)
+Lstepticks = [10]*len(Lplot)
+Lfacconv = [1]*len(Lplot)
+Lcolormap = ['rainbow_r']*len(Lplot)
+LaddWhite_cm = [False]*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='gray', alpha=0.1)
+Panel1.save_graph(4, fig)
+
+
+################################################################
+######### PANEL 5 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'TSRAD', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['TSRAD'][:,:]-273.15 for i in range(1,7)]
+Lcbarlabel = ['°C']*len(Lplot)
+Lminval = [10]*len(Lplot)
+Lmaxval = [50]*len(Lplot)
+Lstep = [0.5]*len(Lplot)
+Lstepticks = [5]*len(Lplot)
+Lfacconv = [1]*len(Lplot)
+Lcolormap = ['rainbow']*len(Lplot)
+LaddWhite_cm = [False]*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='gray', alpha=0.1)
+Panel1.save_graph(5, fig)
+
+
+################################################################
+######### PANEL 6 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'SWD', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['SWD'][:,:] for i in range(1,7)]
+Lcbarlabel = ['W/m2']*len(Lplot)
+Lminval = [0]*len(Lplot)
+Lmaxval = [800]*len(Lplot)
+Lstep = [25]*len(Lplot)
+Lstepticks = [100]*len(Lplot)
+Lfacconv = [1]*len(Lplot)
+Lcolormap = ['binary']*len(Lplot)
+LaddWhite_cm = [True]*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='grey', facecolor='white', alpha=0.1)
+Panel1.save_graph(6, fig)
+
+
+################################################################
+######### PANEL 7 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'LWD', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['LWD'][:,:] for i in range(1,7)]
+Lcbarlabel = ['W/m2']*len(Lplot)
+Lminval = [280]*len(Lplot)
+Lmaxval = [440]*len(Lplot)
+Lstep = [2.5]*len(Lplot)
+Lstepticks = [50]*len(Lplot)
+Lfacconv = [1]*len(Lplot)
+Lcolormap = ['binary_r']*len(Lplot)
+LaddWhite_cm = [False]*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='gray', alpha=0.1)
+Panel1.save_graph(7, fig)
+
+################################################################
+######### PANEL 8 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'W (k=20) 500m', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['WT'][20,:,:] for i in range(1,7)]
+Lcbarlabel = ['m/s']*len(Lplot)
+Lminval = [-5]*len(Lplot)
+Lmaxval = [5]*len(Lplot)
+Lstep = [0.1]*len(Lplot)
+Lstepticks = [1]*len(Lplot)
+Lcolormap = ['seismic']*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='gray', alpha=0.1)
+Panel1.save_graph(8, fig)
+
+################################################################
+######### PANEL 9 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'W (k=27) 1500m', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['WT'][27,:,:] for i in range(1,7)]
+Lcbarlabel = ['m/s']*len(Lplot)
+Lminval = [-5]*len(Lplot)
+Lmaxval = [5]*len(Lplot)
+Lstep = [0.1]*len(Lplot)
+Lstepticks = [1]*len(Lplot)
+Lcolormap = ['seismic']*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='gray', alpha=0.1)
+Panel1.save_graph(9, fig)
+
+
+################################################################
+######### PANEL 10 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'CLDFR 500m', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['CLDFR'][20,:,:] for i in range(1,7)]
+Lcbarlabel = ['-']*len(Lplot)
+Lminval = [0]*len(Lplot)
+Lmaxval = [1]*len(Lplot)
+Lstep = [0.025]*len(Lplot)
+Lstepticks = [0.2]*len(Lplot)
+Lcolormap = ['binary_r']*len(Lplot)
+LaddWhite_cm = [False]*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='white', facecolor='white', alpha=0.1)
+Panel1.save_graph(10, fig)
+
+################################################################
+######### PANEL 11 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'CLDFR 1500m', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['CLDFR'][27,:,:] for i in range(1,7)]
+Lcbarlabel = ['-']*len(Lplot)
+Lminval = [0]*len(Lplot)
+Lmaxval = [1]*len(Lplot)
+Lstep = [0.025]*len(Lplot)
+Lstepticks = [0.2]*len(Lplot)
+Lcolormap = ['binary_r']*len(Lplot)
+LaddWhite_cm = [False]*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='white', facecolor='white', alpha=0.1)
+Panel1.save_graph(11, fig)
+
+################################################################
+######### PANEL 12 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'CLDFR 10000m', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['CLDFR'][53,:,:] for i in range(1,7)]
+Lcbarlabel = ['-']*len(Lplot)
+Lminval = [0]*len(Lplot)
+Lmaxval = [1]*len(Lplot)
+Lstep = [0.025]*len(Lplot)
+Lstepticks = [0.2]*len(Lplot)
+Lcolormap = ['binary_r']*len(Lplot)
+LaddWhite_cm = [False]*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='white', facecolor='white', alpha=0.1)
+Panel1.save_graph(12, fig)
+
+
+################################################################
+######### PANEL 13 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'TKET k=20, 500m', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['TKET'][20,:,:] for i in range(1,7)]
+Lcbarlabel = ['m2/s2']*len(Lplot)
+Lminval = [0]*len(Lplot)
+Lmaxval = [2]*len(Lplot)
+Lstep = [0.1]*len(Lplot)
+Lstepticks = [1]*len(Lplot)
+Lcolormap = ['binary']*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='gray', alpha=0.1)
+Panel1.save_graph(13, fig)
+################################################################
+######### PANEL 14 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'TKET k=27, 1500m', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['TKET'][27,:,:] for i in range(1,7)]
+Lcbarlabel = ['m2/s2']*len(Lplot)
+Lminval = [0]*len(Lplot)
+Lmaxval = [2]*len(Lplot)
+Lstep = [0.1]*len(Lplot)
+Lstepticks = [1]*len(Lplot)
+Lcolormap = ['binary']*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='gray', alpha=0.1)
+Panel1.save_graph(14, fig)
+
+################################################################
+######### PANEL 15 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'U k=2, 6m', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['UT'][2,:,:] for i in range(1,7)]
+Lcbarlabel = ['m/s']*len(Lplot)
+Lminval = [-2]*len(Lplot)
+Lmaxval = [20]*len(Lplot)
+Lstep = [0.1]*len(Lplot)
+Lstepticks = [2.5]*len(Lplot)
+Lcolormap = ['twilight']*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='gray', alpha=0.1)
+Panel1.save_graph(15, fig)
+################################################################
+######### PANEL 16 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'U k=27, 1500m', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['UT'][27,:,:] for i in range(1,7)]
+Lcbarlabel = ['m/s']*len(Lplot)
+Lminval = [-2]*len(Lplot)
+Lmaxval = [20]*len(Lplot)
+Lstep = [0.1]*len(Lplot)
+Lstepticks = [2.5]*len(Lplot)
+Lcolormap = ['twilight']*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='gray', alpha=0.1)
+Panel1.save_graph(16, fig)
+
+################################################################
+######### PANEL 17 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'U k=53, 10000m', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['UT'][53,:,:] for i in range(1,7)]
+Lcbarlabel = ['m/s']*len(Lplot)
+Lminval = [10]*len(Lplot)
+Lmaxval = [40]*len(Lplot)
+Lstep = [0.1]*len(Lplot)
+Lstepticks = [5]*len(Lplot)
+Lcolormap = ['twilight']*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='gray', alpha=0.1)
+Panel1.save_graph(17, fig)
diff --git a/MY_RUN/INTEGRATION_CASES/HPC/PANAME/011_python/plot_PANAME_mod2.py b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/011_python/plot_PANAME_mod2.py
new file mode 100644
index 0000000000000000000000000000000000000000..ba3aeb486c3854307a6c89056607e6ba3ea21bf4
--- /dev/null
+++ b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/011_python/plot_PANAME_mod2.py
@@ -0,0 +1,379 @@
+#!/usr/bin/env python3
+"""
+@author: Quentin Rodier
+Creation : 07/01/2021
+
+Last modifications
+"""
+
+import matplotlib as mpl
+mpl.use('Agg')
+from read_MNHfile import read_netcdf
+from Panel_Plot import PanelPlot
+import cartopy.crs as ccrs
+import numpy as np
+import os
+import cartopy.io.shapereader as shpreader
+
+os.system('rm -f tempgraph*')
+#
+# User's parameter / Namelist
+#
+path="../010_run/"
+LnameFiles = ['PAN.0.2.SEG01.004.nc', 'PAN.0.2.SEG01.008.nc', 'PAN.0.2.SEG01.012.nc', 'PAN.0.2.SEG01.016.nc', 'PAN.0.2.SEG01.020.nc', 'PAN.0.2.SEG01.024.nc']
+
+Dvar_input = {'f'+str(i):['T2M', 'TSRAD', 'HU2M','SWD','LWD','UT', 'CLDFR', 'TKET','WT', 'H','LE','level', 'longitude','latitude','time']
+ for i in range(1,7)}
+
+
+# Read the variables in the files
+Dvar = {}
+Dvar = read_netcdf(LnameFiles, Dvar_input, path=path, removeHALO=True)
+
+# Departements francais
+departements_shp='departements-20180101.shp'
+adm1_shapes = list(shpreader.Reader(departements_shp).geometries())
+
+titlepad=25
+minmaxpad=1.04
+timepad=-0.10
+colorbarpad=0.01
+labelcolorbarpad = 13
+colorbaraspect=20
+
+################################################################
+######### PANEL 1 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'H', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['H'][:,:] for i in range(1,7)]
+Lcbarlabel = ['W/m2']*len(Lplot)
+Lminval = [-100]*len(Lplot)
+Lmaxval = [1000]*len(Lplot)
+Lstep = [25]*len(Lplot)
+Lstepticks = [100]*len(Lplot)
+Lfacconv = [1]*len(Lplot)
+Lcolormap = ['rainbow']*len(Lplot)
+Ltime = [Dvar['f'+str(i)]['time']/3600.0 for i in range(1,7)]
+LaddWhite_cm = [True]*len(Lplot)
+LaxeX = [Dvar['f1']['longitude']]*len(Lplot)
+LaxeY = [Dvar['f1']['latitude']]*len(Lplot)
+Ltitle = ['']*len(Lplot)
+Lxlab = ['longitude']*len(Lplot)
+Lylab = ['latitude']*len(Lplot)
+Lpltype = ['cf']*len(Lplot)
+Lprojection = [ccrs.PlateCarree()]*len(Lplot)
+
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection, gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='gray', alpha=0.1)
+Panel1.save_graph(1, fig)
+
+################################################################
+######### PANEL 2 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'LE', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['LE'][:,:] for i in range(1,7)]
+Lcbarlabel = ['W/m2']*len(Lplot)
+Lminval = [0]*len(Lplot)
+Lmaxval = [600]*len(Lplot)
+Lstep = [25]*len(Lplot)
+Lstepticks = [100]*len(Lplot)
+Lfacconv = [1]*len(Lplot)
+Lcolormap = ['rainbow']*len(Lplot)
+LaddWhite_cm = [True]*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='gray', alpha=0.1)
+Panel1.save_graph(2, fig)
+
+
+################################################################
+######### PANEL 3 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'T2M', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['T2M'][:,:]-273.15 for i in range(1,7)]
+Lcbarlabel = ['°C']*len(Lplot)
+Lminval = [10]*len(Lplot)
+Lmaxval = [30]*len(Lplot)
+Lstep = [0.5]*len(Lplot)
+Lstepticks = [2.5]*len(Lplot)
+Lfacconv = [1]*len(Lplot)
+Lcolormap = ['rainbow']*len(Lplot)
+LaddWhite_cm = [False]*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='gray', alpha=0.1)
+Panel1.save_graph(3, fig)
+
+
+################################################################
+######### PANEL 4 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'HU2M', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['HU2M'][:,:]*100.0 for i in range(1,7)]
+Lcbarlabel = ['%']*len(Lplot)
+Lminval = [10]*len(Lplot)
+Lmaxval = [100]*len(Lplot)
+Lstep = [2]*len(Lplot)
+Lstepticks = [10]*len(Lplot)
+Lfacconv = [1]*len(Lplot)
+Lcolormap = ['rainbow_r']*len(Lplot)
+LaddWhite_cm = [False]*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='gray', alpha=0.1)
+Panel1.save_graph(4, fig)
+
+
+################################################################
+######### PANEL 5 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'TSRAD', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['TSRAD'][:,:]-273.15 for i in range(1,7)]
+Lcbarlabel = ['°C']*len(Lplot)
+Lminval = [10]*len(Lplot)
+Lmaxval = [50]*len(Lplot)
+Lstep = [0.5]*len(Lplot)
+Lstepticks = [5]*len(Lplot)
+Lfacconv = [1]*len(Lplot)
+Lcolormap = ['rainbow']*len(Lplot)
+LaddWhite_cm = [False]*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='gray', alpha=0.1)
+Panel1.save_graph(5, fig)
+
+
+################################################################
+######### PANEL 6 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'SWD', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['SWD'][:,:] for i in range(1,7)]
+Lcbarlabel = ['W/m2']*len(Lplot)
+Lminval = [0]*len(Lplot)
+Lmaxval = [800]*len(Lplot)
+Lstep = [25]*len(Lplot)
+Lstepticks = [100]*len(Lplot)
+Lfacconv = [1]*len(Lplot)
+Lcolormap = ['binary']*len(Lplot)
+LaddWhite_cm = [True]*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='grey', facecolor='white', alpha=0.1)
+Panel1.save_graph(6, fig)
+
+
+################################################################
+######### PANEL 7 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'LWD', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['LWD'][:,:] for i in range(1,7)]
+Lcbarlabel = ['W/m2']*len(Lplot)
+Lminval = [280]*len(Lplot)
+Lmaxval = [440]*len(Lplot)
+Lstep = [2.5]*len(Lplot)
+Lstepticks = [50]*len(Lplot)
+Lfacconv = [1]*len(Lplot)
+Lcolormap = ['binary_r']*len(Lplot)
+LaddWhite_cm = [False]*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='gray', alpha=0.1)
+Panel1.save_graph(7, fig)
+
+################################################################
+######### PANEL 8 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'W (k=20) 500m', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['WT'][20,:,:] for i in range(1,7)]
+Lcbarlabel = ['m/s']*len(Lplot)
+Lminval = [-5]*len(Lplot)
+Lmaxval = [5]*len(Lplot)
+Lstep = [0.1]*len(Lplot)
+Lstepticks = [1]*len(Lplot)
+Lcolormap = ['seismic']*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='gray', alpha=0.1)
+Panel1.save_graph(8, fig)
+
+################################################################
+######### PANEL 9 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'W (k=27) 1500m', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['WT'][27,:,:] for i in range(1,7)]
+Lcbarlabel = ['m/s']*len(Lplot)
+Lminval = [-5]*len(Lplot)
+Lmaxval = [5]*len(Lplot)
+Lstep = [0.1]*len(Lplot)
+Lstepticks = [1]*len(Lplot)
+Lcolormap = ['seismic']*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='gray', alpha=0.1)
+Panel1.save_graph(9, fig)
+
+
+################################################################
+######### PANEL 10 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'CLDFR 500m', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['CLDFR'][20,:,:] for i in range(1,7)]
+Lcbarlabel = ['-']*len(Lplot)
+Lminval = [0]*len(Lplot)
+Lmaxval = [1]*len(Lplot)
+Lstep = [0.025]*len(Lplot)
+Lstepticks = [0.2]*len(Lplot)
+Lcolormap = ['binary_r']*len(Lplot)
+LaddWhite_cm = [False]*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='white', facecolor='white', alpha=0.1)
+Panel1.save_graph(10, fig)
+
+################################################################
+######### PANEL 11 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'CLDFR 1500m', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['CLDFR'][27,:,:] for i in range(1,7)]
+Lcbarlabel = ['-']*len(Lplot)
+Lminval = [0]*len(Lplot)
+Lmaxval = [1]*len(Lplot)
+Lstep = [0.025]*len(Lplot)
+Lstepticks = [0.2]*len(Lplot)
+Lcolormap = ['binary_r']*len(Lplot)
+LaddWhite_cm = [False]*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='white', facecolor='white', alpha=0.1)
+Panel1.save_graph(11, fig)
+
+################################################################
+######### PANEL 12 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'CLDFR 10000m', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['CLDFR'][53,:,:] for i in range(1,7)]
+Lcbarlabel = ['-']*len(Lplot)
+Lminval = [0]*len(Lplot)
+Lmaxval = [1]*len(Lplot)
+Lstep = [0.025]*len(Lplot)
+Lstepticks = [0.2]*len(Lplot)
+Lcolormap = ['binary_r']*len(Lplot)
+LaddWhite_cm = [False]*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='white', facecolor='white', alpha=0.1)
+Panel1.save_graph(12, fig)
+
+
+################################################################
+######### PANEL 13 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'TKET k=20, 500m', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['TKET'][20,:,:] for i in range(1,7)]
+Lcbarlabel = ['m2/s2']*len(Lplot)
+Lminval = [0]*len(Lplot)
+Lmaxval = [2]*len(Lplot)
+Lstep = [0.1]*len(Lplot)
+Lstepticks = [1]*len(Lplot)
+Lcolormap = ['binary']*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='gray', alpha=0.1)
+Panel1.save_graph(13, fig)
+################################################################
+######### PANEL 14 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'TKET k=27, 1500m', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['TKET'][27,:,:] for i in range(1,7)]
+Lcbarlabel = ['m2/s2']*len(Lplot)
+Lminval = [0]*len(Lplot)
+Lmaxval = [2]*len(Lplot)
+Lstep = [0.1]*len(Lplot)
+Lstepticks = [1]*len(Lplot)
+Lcolormap = ['binary']*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='gray', alpha=0.1)
+Panel1.save_graph(14, fig)
+
+################################################################
+######### PANEL 15 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'U k=2, 6m', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['UT'][2,:,:] for i in range(1,7)]
+Lcbarlabel = ['m/s']*len(Lplot)
+Lminval = [-2]*len(Lplot)
+Lmaxval = [20]*len(Lplot)
+Lstep = [0.1]*len(Lplot)
+Lstepticks = [2.5]*len(Lplot)
+Lcolormap = ['twilight']*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='gray', alpha=0.1)
+Panel1.save_graph(15, fig)
+################################################################
+######### PANEL 16 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'U k=27, 1500m', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['UT'][27,:,:] for i in range(1,7)]
+Lcbarlabel = ['m/s']*len(Lplot)
+Lminval = [-2]*len(Lplot)
+Lmaxval = [20]*len(Lplot)
+Lstep = [0.1]*len(Lplot)
+Lstepticks = [2.5]*len(Lplot)
+Lcolormap = ['twilight']*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='gray', alpha=0.1)
+Panel1.save_graph(16, fig)
+
+################################################################
+######### PANEL 17 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'U k=53, 10000m', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['UT'][53,:,:] for i in range(1,7)]
+Lcbarlabel = ['m/s']*len(Lplot)
+Lminval = [10]*len(Lplot)
+Lmaxval = [40]*len(Lplot)
+Lstep = [0.1]*len(Lplot)
+Lstepticks = [5]*len(Lplot)
+Lcolormap = ['twilight']*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='gray', alpha=0.1)
+Panel1.save_graph(17, fig)
diff --git a/MY_RUN/INTEGRATION_CASES/HPC/PANAME/011_python/plot_PANAME_mod3.py b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/011_python/plot_PANAME_mod3.py
new file mode 100644
index 0000000000000000000000000000000000000000..52822cee3642c120bc1c9b38de11a15377e7e01d
--- /dev/null
+++ b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/011_python/plot_PANAME_mod3.py
@@ -0,0 +1,379 @@
+#!/usr/bin/env python3
+"""
+@author: Quentin Rodier
+Creation : 07/01/2021
+
+Last modifications
+"""
+
+import matplotlib as mpl
+mpl.use('Agg')
+from read_MNHfile import read_netcdf
+from Panel_Plot import PanelPlot
+import cartopy.crs as ccrs
+import numpy as np
+import os
+import cartopy.io.shapereader as shpreader
+
+os.system('rm -f tempgraph*')
+#
+# User's parameter / Namelist
+#
+path="../010_run/"
+LnameFiles = ['PAN.0.3.SEG01.004.nc', 'PAN.0.3.SEG01.008.nc', 'PAN.0.3.SEG01.012.nc', 'PAN.0.3.SEG01.016.nc', 'PAN.0.3.SEG01.020.nc', 'PAN.0.3.SEG01.024.nc']
+
+Dvar_input = {'f'+str(i):['T2M', 'TSRAD', 'HU2M','SWD','LWD','UT', 'CLDFR', 'TKET','WT', 'H','LE','level', 'longitude','latitude','time']
+ for i in range(1,7)}
+
+
+# Read the variables in the files
+Dvar = {}
+Dvar = read_netcdf(LnameFiles, Dvar_input, path=path, removeHALO=True)
+
+# Departements francais
+departements_shp='departements-20180101.shp'
+adm1_shapes = list(shpreader.Reader(departements_shp).geometries())
+
+titlepad=25
+minmaxpad=1.04
+timepad=-0.10
+colorbarpad=0.01
+labelcolorbarpad = 13
+colorbaraspect=20
+
+################################################################
+######### PANEL 1 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'H', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['H'][:,:] for i in range(1,7)]
+Lcbarlabel = ['W/m2']*len(Lplot)
+Lminval = [-100]*len(Lplot)
+Lmaxval = [1000]*len(Lplot)
+Lstep = [25]*len(Lplot)
+Lstepticks = [100]*len(Lplot)
+Lfacconv = [1]*len(Lplot)
+Lcolormap = ['rainbow']*len(Lplot)
+Ltime = [Dvar['f'+str(i)]['time']/3600.0 for i in range(1,7)]
+LaddWhite_cm = [True]*len(Lplot)
+LaxeX = [Dvar['f1']['longitude']]*len(Lplot)
+LaxeY = [Dvar['f1']['latitude']]*len(Lplot)
+Ltitle = ['']*len(Lplot)
+Lxlab = ['longitude']*len(Lplot)
+Lylab = ['latitude']*len(Lplot)
+Lpltype = ['cf']*len(Lplot)
+Lprojection = [ccrs.PlateCarree()]*len(Lplot)
+
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection, gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='gray', alpha=0.1)
+Panel1.save_graph(1, fig)
+
+################################################################
+######### PANEL 2 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'LE', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['LE'][:,:] for i in range(1,7)]
+Lcbarlabel = ['W/m2']*len(Lplot)
+Lminval = [0]*len(Lplot)
+Lmaxval = [600]*len(Lplot)
+Lstep = [25]*len(Lplot)
+Lstepticks = [100]*len(Lplot)
+Lfacconv = [1]*len(Lplot)
+Lcolormap = ['rainbow']*len(Lplot)
+LaddWhite_cm = [True]*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='gray', alpha=0.1)
+Panel1.save_graph(2, fig)
+
+
+################################################################
+######### PANEL 3 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'T2M', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['T2M'][:,:]-273.15 for i in range(1,7)]
+Lcbarlabel = ['°C']*len(Lplot)
+Lminval = [10]*len(Lplot)
+Lmaxval = [30]*len(Lplot)
+Lstep = [0.5]*len(Lplot)
+Lstepticks = [2.5]*len(Lplot)
+Lfacconv = [1]*len(Lplot)
+Lcolormap = ['rainbow']*len(Lplot)
+LaddWhite_cm = [False]*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='gray', alpha=0.1)
+Panel1.save_graph(3, fig)
+
+
+################################################################
+######### PANEL 4 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'HU2M', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['HU2M'][:,:]*100.0 for i in range(1,7)]
+Lcbarlabel = ['%']*len(Lplot)
+Lminval = [10]*len(Lplot)
+Lmaxval = [100]*len(Lplot)
+Lstep = [2]*len(Lplot)
+Lstepticks = [10]*len(Lplot)
+Lfacconv = [1]*len(Lplot)
+Lcolormap = ['rainbow_r']*len(Lplot)
+LaddWhite_cm = [False]*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='gray', alpha=0.1)
+Panel1.save_graph(4, fig)
+
+
+################################################################
+######### PANEL 5 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'TSRAD', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['TSRAD'][:,:]-273.15 for i in range(1,7)]
+Lcbarlabel = ['°C']*len(Lplot)
+Lminval = [10]*len(Lplot)
+Lmaxval = [50]*len(Lplot)
+Lstep = [0.5]*len(Lplot)
+Lstepticks = [5]*len(Lplot)
+Lfacconv = [1]*len(Lplot)
+Lcolormap = ['rainbow']*len(Lplot)
+LaddWhite_cm = [False]*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='gray', alpha=0.1)
+Panel1.save_graph(5, fig)
+
+
+################################################################
+######### PANEL 6 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'SWD', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['SWD'][:,:] for i in range(1,7)]
+Lcbarlabel = ['W/m2']*len(Lplot)
+Lminval = [0]*len(Lplot)
+Lmaxval = [800]*len(Lplot)
+Lstep = [25]*len(Lplot)
+Lstepticks = [100]*len(Lplot)
+Lfacconv = [1]*len(Lplot)
+Lcolormap = ['binary']*len(Lplot)
+LaddWhite_cm = [True]*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='grey', facecolor='white', alpha=0.1)
+Panel1.save_graph(6, fig)
+
+
+################################################################
+######### PANEL 7 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'LWD', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['LWD'][:,:] for i in range(1,7)]
+Lcbarlabel = ['W/m2']*len(Lplot)
+Lminval = [280]*len(Lplot)
+Lmaxval = [440]*len(Lplot)
+Lstep = [2.5]*len(Lplot)
+Lstepticks = [50]*len(Lplot)
+Lfacconv = [1]*len(Lplot)
+Lcolormap = ['binary_r']*len(Lplot)
+LaddWhite_cm = [False]*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='gray', alpha=0.1)
+Panel1.save_graph(7, fig)
+
+################################################################
+######### PANEL 8 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'W (k=20) 500m', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['WT'][20,:,:] for i in range(1,7)]
+Lcbarlabel = ['m/s']*len(Lplot)
+Lminval = [-5]*len(Lplot)
+Lmaxval = [5]*len(Lplot)
+Lstep = [0.1]*len(Lplot)
+Lstepticks = [1]*len(Lplot)
+Lcolormap = ['seismic']*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='gray', alpha=0.1)
+Panel1.save_graph(8, fig)
+
+################################################################
+######### PANEL 9 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'W (k=27) 1500m', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['WT'][27,:,:] for i in range(1,7)]
+Lcbarlabel = ['m/s']*len(Lplot)
+Lminval = [-5]*len(Lplot)
+Lmaxval = [5]*len(Lplot)
+Lstep = [0.1]*len(Lplot)
+Lstepticks = [1]*len(Lplot)
+Lcolormap = ['seismic']*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='gray', alpha=0.1)
+Panel1.save_graph(9, fig)
+
+
+################################################################
+######### PANEL 10 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'CLDFR 500m', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['CLDFR'][20,:,:] for i in range(1,7)]
+Lcbarlabel = ['-']*len(Lplot)
+Lminval = [0]*len(Lplot)
+Lmaxval = [1]*len(Lplot)
+Lstep = [0.025]*len(Lplot)
+Lstepticks = [0.2]*len(Lplot)
+Lcolormap = ['binary_r']*len(Lplot)
+LaddWhite_cm = [False]*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='white', facecolor='white', alpha=0.1)
+Panel1.save_graph(10, fig)
+
+################################################################
+######### PANEL 11 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'CLDFR 1500m', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['CLDFR'][27,:,:] for i in range(1,7)]
+Lcbarlabel = ['-']*len(Lplot)
+Lminval = [0]*len(Lplot)
+Lmaxval = [1]*len(Lplot)
+Lstep = [0.025]*len(Lplot)
+Lstepticks = [0.2]*len(Lplot)
+Lcolormap = ['binary_r']*len(Lplot)
+LaddWhite_cm = [False]*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='white', facecolor='white', alpha=0.1)
+Panel1.save_graph(11, fig)
+
+################################################################
+######### PANEL 12 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'CLDFR 10000m', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['CLDFR'][53,:,:] for i in range(1,7)]
+Lcbarlabel = ['-']*len(Lplot)
+Lminval = [0]*len(Lplot)
+Lmaxval = [1]*len(Lplot)
+Lstep = [0.025]*len(Lplot)
+Lstepticks = [0.2]*len(Lplot)
+Lcolormap = ['binary_r']*len(Lplot)
+LaddWhite_cm = [False]*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='white', facecolor='white', alpha=0.1)
+Panel1.save_graph(12, fig)
+
+
+################################################################
+######### PANEL 13 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'TKET k=20, 500m', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['TKET'][20,:,:] for i in range(1,7)]
+Lcbarlabel = ['m2/s2']*len(Lplot)
+Lminval = [0]*len(Lplot)
+Lmaxval = [2]*len(Lplot)
+Lstep = [0.1]*len(Lplot)
+Lstepticks = [1]*len(Lplot)
+Lcolormap = ['binary']*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='gray', alpha=0.1)
+Panel1.save_graph(13, fig)
+################################################################
+######### PANEL 14 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'TKET k=27, 1500m', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['TKET'][27,:,:] for i in range(1,7)]
+Lcbarlabel = ['m2/s2']*len(Lplot)
+Lminval = [0]*len(Lplot)
+Lmaxval = [2]*len(Lplot)
+Lstep = [0.1]*len(Lplot)
+Lstepticks = [1]*len(Lplot)
+Lcolormap = ['binary']*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='gray', alpha=0.1)
+Panel1.save_graph(14, fig)
+
+################################################################
+######### PANEL 15 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'U k=2, 6m', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['UT'][2,:,:] for i in range(1,7)]
+Lcbarlabel = ['m/s']*len(Lplot)
+Lminval = [-2]*len(Lplot)
+Lmaxval = [20]*len(Lplot)
+Lstep = [0.1]*len(Lplot)
+Lstepticks = [2.5]*len(Lplot)
+Lcolormap = ['twilight']*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='gray', alpha=0.1)
+Panel1.save_graph(15, fig)
+################################################################
+######### PANEL 16 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'U k=27, 1500m', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['UT'][27,:,:] for i in range(1,7)]
+Lcbarlabel = ['m/s']*len(Lplot)
+Lminval = [-2]*len(Lplot)
+Lmaxval = [20]*len(Lplot)
+Lstep = [0.1]*len(Lplot)
+Lstepticks = [2.5]*len(Lplot)
+Lcolormap = ['twilight']*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='gray', alpha=0.1)
+Panel1.save_graph(16, fig)
+
+################################################################
+######### PANEL 17 #
+###############################################################
+Panel1 = PanelPlot(2,3, [25,14],'U k=53, 10000m', titlepad=titlepad, minmaxpad=minmaxpad, timepad=timepad, colorbarpad=colorbarpad, labelcolorbarpad = labelcolorbarpad , colorbaraspect=colorbaraspect)
+Lplot = [ Dvar['f'+str(i)]['UT'][53,:,:] for i in range(1,7)]
+Lcbarlabel = ['m/s']*len(Lplot)
+Lminval = [10]*len(Lplot)
+Lmaxval = [40]*len(Lplot)
+Lstep = [0.1]*len(Lplot)
+Lstepticks = [5]*len(Lplot)
+Lcolormap = ['twilight']*len(Lplot)
+fig = Panel1.psectionH(lon=LaxeX, lat=LaxeY, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval,
+ Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv,
+ colorbar=True, Ltime=Ltime, LaddWhite_cm=LaddWhite_cm, Lproj=Lprojection,gridLines=False)
+for i in range(6):
+ fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='gray', alpha=0.1)
+Panel1.save_graph(17, fig)
diff --git a/MY_RUN/INTEGRATION_CASES/HPC/PANAME/011_python/run_python b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/011_python/run_python
new file mode 100755
index 0000000000000000000000000000000000000000..2697cdded1a3b0a547c842187025592f38cccb25
--- /dev/null
+++ b/MY_RUN/INTEGRATION_CASES/HPC/PANAME/011_python/run_python
@@ -0,0 +1,30 @@
+#!/bin/sh
+#SBATCH -J py_PANAME
+#SBATCH -N 1 # nodes number
+#SBATCH -n 1 # CPUs number (on all nodes)
+#SBATCH -o tools.eo%j #
+#SBATCH -e tools.eo%j #
+#SBATCH -t 01:00:00 # time limit
+# Echo des commandes
+ulimit -c 0
+ulimit -s unlimited
+# Arrete du job des la premiere erreur
+set -e
+
+. ~rodierq/DEV_57/MNH-PHYEX070-b95d84d7/conf/profile_mesonh-LXifort-R8I4-MNH-V5-6-2-ECRAD140-MPIAUTO-O2
+ln -sf ${SRC_MESONH}/src/LIB/Python/* .
+ln -sf ~rodierq/SAVE/OUTILS/PYTHON/departements-20180101.* .
+
+module purge
+module load python/3.7.6
+
+python3 plot_PANAME_mod1.py
+convert *.png PANAMEmod1.pdf
+rm -f *.png
+
+python3 plot_PANAME_mod2.py
+convert *.png PANAMEmod2.pdf
+rm -f *.png
+
+python3 plot_PANAME_mod3.py
+convert *.png PANAMEmod3.pdf