MY_RUN/INTEGRATION_CASES/HPC/AZF_2M/003_nest/run_prep_nest_pgd_xyz

+#!/bin/sh
+#SBATCH -J azf_pgd
+#SBATCH -N 1            # nodes number (=NBP)
+#SBATCH -n 4            # CPUs number (on all nodes) (=NBP*TPN)
+#SBATCH -o run1_16jan.eo%j   #
+#SBATCH -e run1_16jan.eo%j   #
+#SBATCH -t 00:30:00    # time limit
+#
+## Echo des commandes
+ulimit -c 0
+ulimit -s unlimited
+## Arrete du job des la premiere erreur
+set -e
+## Nom de la machine
+hostname
+#
+#
+. ~rodierq/DEV_57/MNH-PHYEX070-b95d84d7/conf/profile_mesonh-LXifort-R8I4-MNH-V5-6-2-ECRAD140-MPIAUTO-O2
+
+#
+export MPIRUN="Mpirun -np 4"
+#
+rm -f *ne.??? OUTPUT_LISTING* pipe* *.tex
+
+ln -sf ../002_pgd2/PGD_AZF_500m.* .
+ln -sf ../001_pgd1/PGD_AZF_2.5km.* .
+
+set -x
+set -e
+
+time ${MPIRUN} PREP_NEST_PGD${XYZ}
+cd ../004_arp2lfi/
+sbatch run_arp2lfi_xyz

MY_RUN/INTEGRATION_CASES/HPC/AZF_2M/004_arp2lfi/PRE_REAL1.nam_09

+&NAM_CONFIO  LCDF4=T, LLFIOUT=F, LLFIREAD=F /
+&NAM_FILE_NAMES HATMFILE='ecmwf.OD.20010921.00.09',HATMFILETYPE='GRIBEX',
+                HPGDFILE='PGD_AZF_2.5km.nestne',
+                CINIFILE='CPL_J1_20010921.09' /
+&NAM_REAL_CONF NVERB=5, CEQNSYS='DUR',JPHEXT=1 , NHALO=1  /
+&NAM_VER_GRID YZGRID_TYPE='MANUAL', NKMAX=40  /
+ZHAT
+.0000000E+00
+.7213016E+02
+.1504758E+03
+.2378252E+03
+.3367744E+03
+.4497268E+03
+.5788933E+03
+.7262926E+03
+.8937506E+03
+.1082901E+04
+.1295185E+04
+.1531851E+04
+.1793956E+04
+.2082363E+04
+.2397744E+04
+.2740577E+04
+.3111149E+04
+.3509554E+04
+.3935693E+04
+.4389276E+04
+.4869819E+04
+.5376646E+04
+.5908888E+04
+.6465486E+04
+.7045184E+04
+.7646539E+04
+.8269539E+04
+.8914539E+04
+.9582539E+04
+.10275039E+05
+.10992639E+05
+.11736239E+05
+.12506739E+05
+.13305139E+05
+.14132439E+05
+.14989699E+05
+.15877999E+05
+.16798459E+05
+.17752259E+05
+.18740559E+05
+.19764669E+05
+

MY_RUN/INTEGRATION_CASES/HPC/AZF_2M/004_arp2lfi/PRE_REAL1.nam_12

+&NAM_CONFIO  LCDF4=T, LLFIOUT=F, LLFIREAD=F /
+&NAM_FILE_NAMES HATMFILE='ecmwf.OD.20010921.00.12',HATMFILETYPE='GRIBEX',
+                HPGDFILE='PGD_AZF_2.5km.nestne',
+                CINIFILE='CPL_J1_20010921.12' /
+&NAM_REAL_CONF NVERB=5, CEQNSYS='DUR' , JPHEXT=1 , NHALO=1/
+&NAM_VER_GRID YZGRID_TYPE='MANUAL', NKMAX=40  /
+ZHAT
+.0000000E+00
+.7213016E+02
+.1504758E+03
+.2378252E+03
+.3367744E+03
+.4497268E+03
+.5788933E+03
+.7262926E+03
+.8937506E+03
+.1082901E+04
+.1295185E+04
+.1531851E+04
+.1793956E+04
+.2082363E+04
+.2397744E+04
+.2740577E+04
+.3111149E+04
+.3509554E+04
+.3935693E+04
+.4389276E+04
+.4869819E+04
+.5376646E+04
+.5908888E+04
+.6465486E+04
+.7045184E+04
+.7646539E+04
+.8269539E+04
+.8914539E+04
+.9582539E+04
+.10275039E+05
+.10992639E+05
+.11736239E+05
+.12506739E+05
+.13305139E+05
+.14132439E+05
+.14989699E+05
+.15877999E+05
+.16798459E+05
+.17752259E+05
+.18740559E+05
+.19764669E+05
+

MY_RUN/INTEGRATION_CASES/HPC/AZF_2M/004_arp2lfi/clean_arp2lfi_xyz

+set -x
+
+rm -f OUTPUT_LISTING* OUTPUT_TRANSFER pipe* *.tex *.dat *.bin PGD* CPL*
+

MY_RUN/INTEGRATION_CASES/HPC/AZF_2M/004_arp2lfi/run_arp2lfi_xyz

+#!/bin/sh
+#SBATCH -J azf_prep
+#SBATCH -N 1            # nodes number (=NBP)
+#SBATCH -n 1            # CPUs number (on all nodes) (=NBP*TPN)
+#SBATCH -o prep_azf.eo%j   #
+#SBATCH -e prep_azf.eo%j   #
+#SBATCH -t 01:30:00    # time limit
+
+# Echo des commandes
+ulimit -c 0
+ulimit -s unlimited
+# Arrete du job des la premiere erreur
+set -e
+# Nom de la machine
+hostname
+
+
+. ~rodierq/DEV_57/MNH-PHYEX070-b95d84d7/conf/profile_mesonh-LXifort-R8I4-MNH-V5-6-2-ECRAD140-MPIAUTO-O2
+
+ln -sf ../003_nest/PGD_AZF* .
+ln -sf $HOME/SAVE/GRIB_KTEST/ecmwf.OD.20010921.00.09
+ln -sf $HOME/SAVE/GRIB_KTEST/ecmwf.OD.20010921.00.12
+rm -f  OUTPUT_LISTING* pipe* *.tex
+
+export MPIRUN="Mpirun -np 1"
+
+set -x
+set -e
+
+cp -f PRE_REAL1.nam_09 PRE_REAL1.nam
+time ${MPIRUN} PREP_REAL_CASE${XYZ}
+mv OUTPUT_LISTING0 OUTPUT_LISTING0_09
+cp -f PRE_REAL1.nam_12 PRE_REAL1.nam
+time ${MPIRUN} PREP_REAL_CASE${XYZ}
+mv OUTPUT_LISTING0 OUTPUT_LISTING0_12
+
+rm -f PRE_REAL1.nam
+cd ../005_spa_mod1_mod2/
+sbatch run_spawning_xyz

MY_RUN/INTEGRATION_CASES/HPC/AZF_2M/005_spa_mod1_mod2/SPAWN1.nam

+&NAM_CONFIO  LCDF4=T, LLFIOUT=F, LLFIREAD=F /
+&NAM_GRID2_SPA GBAL_ONLY=.FALSE./
+&NAM_LUNIT2_SPA CINIFILE='CPL_J1_20010921.09',
+                CINIFILEPGD='PGD_AZF_2.5km.nestne'
+                YDOMAIN='PGD_AZF_500m.nestne',           
+                YSPAFILE='CPL_J1_20010921.09.spa2' /
+&NAM_SPAWN_SURF LSPAWN_SURF=T /
+&NAM_BLANK /

MY_RUN/INTEGRATION_CASES/HPC/AZF_2M/005_spa_mod1_mod2/clean_spawning_xyz

+set -x
+rm -f PGD* CPL* OUTPUT_LISTING* OUTPUT_TRANSFER pipe* *.tex
+

MY_RUN/INTEGRATION_CASES/HPC/AZF_2M/005_spa_mod1_mod2/run_spawning_xyz

+#!/bin/sh
+#SBATCH -J azf
+#SBATCH -N 1            # nodes number (=NBP)
+#SBATCH -n 1            # CPUs number (on all nodes) (=NBP*TPN)
+#SBATCH -o run1_16jan.eo%j   #
+#SBATCH -e run1_16jan.eo%j   #
+#SBATCH -t 00:30:00    # time limit
+
+# Echo des commandes
+ulimit -c 0
+ulimit -s unlimited
+# Arrete du job des la premiere erreur
+set -e
+# Nom de la machine
+hostname
+
+
+. ~rodierq/DEV_57/MNH-PHYEX070-b95d84d7/conf/profile_mesonh-LXifort-R8I4-MNH-V5-6-2-ECRAD140-MPIAUTO-O2
+
+export MPIRUN="Mpirun -np 1"
+
+rm -f *spa04.* OUTPUT_LISTING* pipe* *.tex
+ln -sf  ../003_nest/PGD_AZF_2.5km.nestne.* .
+ln -sf  ../003_nest/PGD_AZF_500m.nestne.* .
+ln -sf  ../004_arp2lfi/CPL_J1_20010921.09.* .
+
+set -x
+set -e
+
+time ${MPIRUN} SPAWNING${XYZ}
+cd ../006_preal/
+sbatch run_preal_xyz

MY_RUN/INTEGRATION_CASES/HPC/AZF_2M/006_preal/PRE_REAL1.nam

+&NAM_CONFIO  LCDF4=T, LLFIOUT=F, LLFIREAD=F /
+&NAM_FILE_NAMES HATMFILE='CPL_J1_20010921.09.spa2', 
+                HATMFILETYPE='MESONH',
+                HPGDFILE='PGD_AZF_500m.nestne',             
+                CINIFILE='CPL_J1_20010921.09.2' /
+&NAM_VER_GRID   YZGRID_TYPE='SAMEGR' /
+&NAM_REAL_CONF  CEQNSYS='DUR' ,
+                CPRESOPT= 'RICHA',
+                NVERB=5 ,JPHEXT=1,NHALO=1/
+&NAM_BLANK /
+&NAM_PREP_SURF_ATM CFILE='CPL_J1_20010921.09',
+                   CFILEPGD="PGD_AZF_2.5km.nestne",
+                   CFILETYPE='MESONH' /
+&NAM_DIAG_ISBAn LPGD=.TRUE. /
+&NAM_DIAG_SURF_ATMn LFRAC=.TRUE. /
+
+            

MY_RUN/INTEGRATION_CASES/HPC/AZF_2M/006_preal/clean_preal_xyz

+set -x
+rm -f CPL* PGD* OUTPUT_LISTING* OUTPUT_TRANSFER pipe* *.tex
+

MY_RUN/INTEGRATION_CASES/HPC/AZF_2M/006_preal/run_preal_xyz

+#!/bin/sh
+#SBATCH -J azf_prep
+#SBATCH -N 1            # nodes number (=NBP)
+#SBATCH -n 8            # CPUs number (on all nodes) (=NBP*TPN)
+#SBATCH -o prep_azf.eo%j   #
+#SBATCH -e prep_azf.eo%j   #
+#SBATCH -t 01:30:00    # time limit
+
+# Echo des commandes
+ulimit -c 0
+ulimit -s unlimited
+# Arrete du job des la premiere erreur
+set -e
+# Nom de la machine
+hostname
+
+
+. ~rodierq/DEV_57/MNH-PHYEX070-b95d84d7/conf/profile_mesonh-LXifort-R8I4-MNH-V5-6-2-ECRAD140-MPIAUTO-O2
+
+ln -sf ../003_nest/PGD_AZF* .
+ln -sf ../004_arp2lfi/CPL_J1_20010921.09.* .
+ln -sf ../004_arp2lfi/CPL_J3_20010921.09.* .
+ln -sf ../005_spa_mod1_mod2/CPL_J1_20010921.09.spa2.* .
+ln -sf ../005_spa_mod1_mod2/CPL_J3_20010921.09.spa2.* .
+
+export MPIRUN="Mpirun -np 8"
+
+time ${MPIRUN} PREP_REAL_CASE${XYZ}
+cd ../007_run
+sbatch run_mesonh

MY_RUN/INTEGRATION_CASES/HPC/AZF_2M/007_run/EXSEG1.nam

+&NAM_CONFIO  LCDF4=T, LLFIOUT=F, LLFIREAD=F, NIO_VERB=5 /
+&NAM_LUNITn CINIFILE = "CPL_J1_20010921.09",  
+            CINIFILEPGD='PGD_AZF_2.5km.nestne',       
+            CCPLFILE(1) = "CPL_J1_20010921.12" /
+&NAM_CONFZ MPI_BUFFER_SIZE = 80 /
+&NAM_CONFn LUSERV = T, LUSERC = T, LUSERR = T ,
+           LUSERI = T, LUSERS = T, LUSERG = T, LUSERH = F, LUSECI= T,
+           NSV_USER =0 /
+&NAM_DYNn XTSTEP=10., 
+      NITR = 12, CPRESOPT="RICHA", XRELAX=1.0,  LITRADJ=.TRUE.,
+      LVE_RELAX = .TRUE., 
+      LHORELAX_UVWTH=.TRUE., LHORELAX_RV=.TRUE., 
+      LHORELAX_RC=.FALSE., LHORELAX_RR=.FALSE., 
+      LHORELAX_RS=.FALSE., LHORELAX_RI=.FALSE., LHORELAX_RG=.FALSE.,
+      LHORELAX_SVPP=.F., 
+      NRIMX = 5, NRIMY = 5, 
+      XRIMKMAX = .002, XT4DIFU=1800. /
+&NAM_ADVn CUVW_ADV_SCHEME="CEN4TH", CMET_ADV_SCHEME="PPM_01", CSV_ADV_SCHEME="PPM_01" ,CTEMP_SCHEME='RKC4'/
+&NAM_PARAMn  CTURB   = "TKEL", CRAD   = "ECMW",   
+             CCLOUD = "ICE3" ,  CSCONV = "EDKF" /
+&NAM_PARAM_RADn XDTRAD = 900., XDTRAD_CLONLY = 300., 
+                LCLEAR_SKY=.FALSE. /
+&NAM_PARAM_MFSHALLn CMF_UPDRAFT='EDKF',CMF_CLOUD='DIRE', LMIXUV=T, LMF_FLX=F /
+&NAM_LBCn  CLBCX = 2*"OPEN", CLBCY = 2*"OPEN", XCPHASE = 5. /
+&NAM_TURBn CTURBDIM="1DIM", CTURBLEN="BL89",
+           LTURB_DIAG=.TRUE.,  LTURB_FLX=.TRUE.  /
+&NAM_CH_MNHCn LUSECHEM = .FALSE. /
+&NAM_CONF CCONF = "START", LFLAT = F, CEQNSYS = "DUR",
+          NMODEL = 2, NVERB = 10, CEXP = "AZF02", CSEG = "CEN4T" /
+&NAM_DYN XSEGLEN = 3600, XASSELIN = 0.2, LCORIO = T, LNUMDIFU = T,
+         XALKTOP = 0.003, XALZBOT = 15000. /
+&NAM_NESTING NDAD(2) = 1,NDTRATIO(2) =5 ,XWAY(2) = 2.  /
+&NAM_BACKUP XBAK_TIME(1,1)=10., XBAK_TIME(1,2)=60.,XBAK_TIME(1,3)=600.,XBAK_TIME(1,4)=900.,XBAK_TIME(1,5)=1800.,XBAK_TIME(1,6)=2700.,XBAK_TIME(1,7)=3600. /
+&NAM_PASPOL
+LPASPOL   =  T  ,
+NRELEASE     =   2        ,
+CPPINIT(1)   = "1PT"      ,
+XPPLAT(1)    =  43.567 ,
+XPPLON(1)    =   1.439 ,
+XPPBOT(1)    =   10.0     ,
+XPPTOP(1)    =   500.0     ,
+XPPMASS(1)   = 10000000.  ,
+CPPT1(1)     = "20010921090000",
+CPPT2(1)     = "20010921090000",
+CPPT3(1)     = "20010921091500",
+CPPT4(1)     = "20010921091500", 
+CPPINIT(2)   = "1PT"      ,
+XPPLAT(2)    =  43.9   ,
+XPPLON(2)    =   1.5   ,
+XPPBOT(2)    =   10.0     ,
+XPPTOP(2)    =   1000.0     ,
+XPPMASS(2)   = 10000000.  ,
+CPPT1(2)     = "20010921090000",
+CPPT2(2)     = "20010921090000",
+CPPT3(2)     = "20010921090100",
+CPPT4(2)     = "20010921090100" /
+&NAM_ISBA 
+       CALBEDO = "DRY"
+/
+&NAM_ISBAn
+       CSCOND  = "NP89",
+       CC1DRY     = 'DEF',
+       CSOILFRZ   = 'DEF',
+       CDIFSFCOND = 'DEF',
+       CSNOWRES   = 'DEF'
+         /
+&NAM_SGH_ISBAn
+       CRUNOFF = "WSAT" /
+&NAM_SEAFLUXn   CSEA_FLUX="ECUME", CSEA_ALB="UNIF" /
+&NAM_DIAG_SURFn N2M=1 LSURF_BUDGET=T /
+&NAM_BUDGET CBUTYPE='NONE' /
+&NAM_AERO_CONF LDUST=F /
+
+&NAM_FLYERS NAIRCRAFTS=1 /
+
+&NAM_AIRCRAFTS
+  NMODEL(1) = 1,
+  CMODEL(1) = 'MOB',
+  CTITLE(1) = 'AVION',
+  TLAUNCH(1)%nyear  = 2001,
+  TLAUNCH(1)%nmonth = 09,
+  TLAUNCH(1)%nday   = 21,
+  TLAUNCH(1)%xtime  = 32400.,
+  XTSTEP(1)     = 10.,
+  NPOS(1) = 9,
+  LALTDEF(1) = .TRUE.,
+  CFILE(1) = 'aircraft_data.csv' /
+&NAM_NEBn LSUBG_COND=.T.  /

MY_RUN/INTEGRATION_CASES/HPC/AZF_2M/007_run/EXSEG2.nam

+&NAM_LUNITn CINIFILE = "CPL_J1_20010921.09.2",CINIFILEPGD='PGD_AZF_500m.nestne' /      
+&NAM_CONFn LUSERV = T, LUSERC = T, LUSERR = T ,
+           LUSERI = T, LUSERS = T, LUSERG = T, LUSERH = F, LUSECI= T,
+           NSV_USER =0 /
+&NAM_DYNn 
+      NITR = 8, LITRADJ=.TRUE.,
+      LVE_RELAX = .TRUE., 
+      LHORELAX_UVWTH=.FALSE.,  LHORELAX_RV=.FALSE., 
+      LHORELAX_RC=.FALSE., LHORELAX_RR=.FALSE., 
+      LHORELAX_RS=.FALSE., LHORELAX_RI=.FALSE., LHORELAX_RG=.FALSE., 
+      NRIMX = 5, NRIMY = 5, 
+      XRIMKMAX = .002, XT4DIFU=1800. /
+&NAM_ADVn CUVW_ADV_SCHEME="CEN4TH4", CMET_ADV_SCHEME="PPM_01", CSV_ADV_SCHEME="PPM_01" ,CTEMP_SCHEME='RKC4'/
+&NAM_PARAMn  CTURB   = "TKEL", CRAD   = "ECMW",   
+             CCLOUD = "ICE3" ,  CSCONV = "NONE" /
+&NAM_PARAM_RADn XDTRAD = 900., XDTRAD_CLONLY = 300., 
+                LCLEAR_SKY=.FALSE. /
+&NAM_LBCn  CLBCX = 2*"OPEN", CLBCY = 2*"OPEN", XCPHASE = 5. /
+&NAM_TURBn CTURBDIM="1DIM", CTURBLEN="BL89",
+           LTURB_DIAG=.FALSE.,  LTURB_FLX=.FALSE.  /
+&NAM_CH_MNHCn LUSECHEM = .FALSE. /
+&NAM_ISBA 
+       CALBEDO = "DRY"
+/
+&NAM_ISBAn
+       CSCOND  = "NP89",
+       CC1DRY     = 'DEF',
+       CSOILFRZ   = 'DEF',
+       CDIFSFCOND = 'DEF',
+       CSNOWRES   = 'DEF'
+         /
+&NAM_SGH_ISBAn
+       CRUNOFF = "WSAT" /
+&NAM_SEAFLUXn   CSEA_FLUX="ECUME", CSEA_ALB="UNIF" /
+&NAM_DIAG_SURFn N2M=1 LSURF_BUDGET=T /
+&NAM_NEBn LSUBG_COND=.T.  /

MY_RUN/INTEGRATION_CASES/HPC/AZF_2M/007_run/aircraft_data.csv

+Time Lat Lon Alt(hPa)
+0.   43.2 1. 900.
+60.  43.4 1. 850.
+120. 43.4 1.2 700.
+180. 43.6 1.2 700.
+240. 43.6 1.4 700.
+300. 43.8 1.4 700.
+360. 43.8 1.6 700.
+420. 44.0 1.6 850.
+480. 44.0 1.8 900.

MY_RUN/INTEGRATION_CASES/HPC/AZF_2M/007_run/clean_mesonh_xyz

+set -x
+rm -f PGD* CPL* AZF* OUTPUT_LISTING* OUTPUT_TRANSFER pipe* *.tex *.dat *.bin
+

MY_RUN/INTEGRATION_CASES/HPC/AZF_2M/007_run/run_mesonh

+#!/bin/sh
+#SBATCH -J azf
+#SBATCH -N 1            # nodes number (=NBP)
+#SBATCH -n 16            # CPUs number (on all nodes) (=NBP*TPN)
+#SBATCH -o run1_16jan.eo%j   #
+#SBATCH -e run1_16jan.eo%j   #
+#SBATCH -t 02:00:00    # time limit
+
+# Echo des commandes
+ulimit -c 0
+ulimit -s unlimited
+# Arrete du job des la premiere erreur
+set -e
+# Nom de la machine
+hostname 
+
+ln -sf ../006_preal/PGD* .
+ln -sf ../006_preal/CPL* .
+ln -sf ../004_arp2lfi/CPL* .
+. ~rodierq/DEV_57/MNH-PHYEX070-b95d84d7/conf/profile_mesonh-LXifort-R8I4-MNH-V5-6-2-ECRAD140-MPIAUTO-O2
+
+export MPIRUN="Mpirun -np 8"
+
+set -x
+set -e
+ls -lrt
+
+time ${MPIRUN} MESONH${XYZ}
+cd ../009_python
+sbatch run_python

MY_RUN/INTEGRATION_CASES/HPC/AZF_2M/009_python/plot_AZF2M.py

+#!/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
+from misc_functions import mean_operator, convert_date
+import cartopy.crs as ccrs
+import numpy as np
+import os
+import cartopy.io.shapereader as shpreader
+import matplotlib.patches as mpatches
+os.system('rm -f tempgraph*')
+#
+#  User's parameter / Namelist
+#
+path=""
+
+LnameFiles = ['AZF02.1.CEN4T.001.nc', 'AZF02.1.CEN4T.002.nc', 'AZF02.1.CEN4T.003.nc',
+              'AZF02.1.CEN4T.004.nc', 'AZF02.1.CEN4T.005.nc', 'AZF02.1.CEN4T.007.nc',
+              'AZF02.2.CEN4T.001.nc', 'AZF02.2.CEN4T.002.nc', 'AZF02.2.CEN4T.003.nc',
+              'AZF02.2.CEN4T.004.nc', 'AZF02.2.CEN4T.005.nc', 'AZF02.2.CEN4T.007.nc',
+              'AZF02.1.CEN4T.000.nc']
+LG_AVION='/Flyers/Aircrafts/AVION/'
+LG_AVIONT='/Flyers/Aircrafts/AVION/Point/'
+LG_AVIONZT='/Flyers/Aircrafts/AVION/Vertical_profile/'
+
+Dvar_input = {
+'f1':['SVPP001','SVPP002','ATC001','ATC002','UT','VT','latitude','longitude','level'],
+'f2':['SVPP001','SVPP002','ATC001','ATC002','UT','VT','latitude','longitude','level'],
+'f3':['SVPP001','SVPP002','ATC001','ATC002','UT','VT','latitude','longitude','level'],
+'f4':['SVPP001','SVPP002','ATC001','ATC002','UT','VT','latitude','longitude','level'],
+'f5':['SVPP001','SVPP002','ATC001','ATC002','UT','VT','latitude','longitude','level'],
+'f6':['SVPP001','SVPP002','ATC001','ATC002','UT','VT','latitude','longitude','level'],
+'f7':['SVPP001','SVPP002','ATC001','ATC002','UT','VT','latitude','longitude','level','LONOR','LATOR','LAT','LON'],
+'f8':['SVPP001','SVPP002','ATC001','ATC002','UT','VT','latitude','longitude','level'],
+'f9':['SVPP001','SVPP002','ATC001','ATC002','UT','VT','latitude','longitude','level'],
+'f10':['SVPP001','SVPP002','ATC001','ATC002','UT','VT','latitude','longitude','level'],
+'f11':['SVPP001','SVPP002','ATC001','ATC002','UT','VT','latitude','longitude','level'],
+'f12':['SVPP001','SVPP002','ATC001','ATC002','UT','VT','latitude','longitude','level'],
+'f13':[(LG_AVION,'time_flyer'),(LG_AVIONT,'ZS'), (LG_AVIONT,'P'), (LG_AVIONT,'LON'),(LG_AVIONT,'MER_WIND'),
+       (LG_AVIONT,'ZON_WIND'),(LG_AVIONT,'W'), (LG_AVIONT,'Th'), (LG_AVIONT,'Rv'),(LG_AVIONT,'Tke'),
+       (LG_AVIONT,'H_FLUX'),(LG_AVIONT,'LE_FLUX'), (LG_AVIONT,'Tke_Diss'), (LG_AVIONT,'Tsrad')]
+}
+
+#  Read the variables in the files
+Dvar = {}
+Dvar = read_netcdf(LnameFiles, Dvar_input, path=path, removeHALO=True)
+
+################################################################
+#########          PANEL 1
+###############################################################
+Panel = PanelPlot(2,3, [25,14],'Domaine 1 SV 001', titlepad=25, minmaxpad=1.04, timepad=-0.07, colorbarpad=0.03, labelcolorbarpad = 13, colorbaraspect=22)
+Lplot = [Dvar['f1']['SVPP001'], Dvar['f2']['SVPP001'], Dvar['f3']['SVPP001'],
+         Dvar['f4']['SVPP001'], Dvar['f5']['SVPP001'], Dvar['f6']['SVPP001']]
+
+lon = [Dvar['f1']['longitude']]*len(Lplot)
+lat = [Dvar['f1']['latitude']]*len(Lplot)
+Ltitle = ['SVPP001']*len(Lplot)
+Lcbarlabel = ['kg/kg']*len(Lplot)
+Lxlab = ['longitude']*len(Lplot)
+Lylab = ['latitude']*len(Lplot)
+Lminval = [0]*len(Lplot)
+Lmaxval = [0.15E-5]*len(Lplot)
+Lstep = [0.05E-6]*len(Lplot)
+Lstepticks = [0.2E-6]*len(Lplot) 
+Lcolormap = ['gist_rainbow_r']*len(Lplot)
+Lprojection = [ccrs.PlateCarree()]*len(Lplot)
+LaddWhite = [True]*len(Lplot)
+Llevel = [0]*len(Lplot)
+Ltime = [Dvar['f1']['date'], Dvar['f2']['date'], Dvar['f3']['date'], Dvar['f4']['date'], Dvar['f5']['date'], Dvar['f6']['date']]
+Lcbformatlabel=[True]*len(Lplot)
+fig = Panel.psectionH(lon=lon, lat=lat, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, 
+                                Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, 
+                                Ltime=Ltime, LaddWhite_cm=LaddWhite, Lproj=Lprojection, Llevel=Llevel, Lcbformatlabel=Lcbformatlabel)
+
+Lplot1 = [ Dvar['f1']['UT'], Dvar['f2']['UT'], Dvar['f3']['UT'], Dvar['f4']['UT'], Dvar['f5']['UT'], Dvar['f6']['UT']]
+Lplot2 = [ Dvar['f1']['VT'], Dvar['f2']['VT'], Dvar['f3']['VT'], Dvar['f4']['VT'], Dvar['f5']['VT'], Dvar['f6']['VT']]
+Ltitle = ['wind vectors at K=2']*len(Lplot)
+Llegendval = [7.5]*len(Lplot)
+Lcbarlabel = ['(m/s)']*len(Lplot1)
+Larrowstep = [2]*len(Lplot1)
+Lwidth = [0.002]*len(Lplot1)
+Lcolor = ['black']*len(Lplot1)
+Lscale = [100]*len(Lplot1)
+fig = Panel.pvector(Lxx=lon, Lyy=lat, Llevel=Llevel, Lvar1=Lplot1, Lvar2=Lplot2, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lwidth=Lwidth, Larrowstep=Larrowstep,
+                      Llegendval=Llegendval, Llegendlabel=Lcbarlabel, Lproj=Lprojection, Lid_overlap=[0,2,4,6,8,10], ax=fig.axes, Lscale=Lscale)
+
+# Departements francais
+departements_shp='departements-20180101.shp'
+adm1_shapes = list(shpreader.Reader(departements_shp).geometries())
+
+# Add departements to each axes + scatter point of emission source
+loncar, latcar = [1.439,1.5], [43.567, 43.9]
+label=['AZF1','AZF2']
+for i in range(len(Lplot)):
+  fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='white', alpha=0.2)
+  fig.axes[i*2].scatter(loncar,latcar)
+  for lab, txt in enumerate(label):
+    fig.axes[i*2].annotate(label[lab], (loncar[lab], latcar[lab]), color='black',size=10, weight="bold")
+
+# Add a Rectangle displaying the domain of the model 2
+for i in range(len(Lplot)):
+  fig.axes[i*2].add_patch(mpatches.Rectangle(xy=[Dvar['f7']['LONOR'], Dvar['f7']['LATOR']], width=Dvar['f7']['LON'][-1,-1]-Dvar['f7']['LONOR'], height=Dvar['f7']['LAT'][-1,-1]-Dvar['f7']['LATOR'],
+   facecolor='blue', alpha=0.15, transform=ccrs.PlateCarree()))
+
+fig.tight_layout()
+Panel.save_graph(1,fig)
+
+################################################################
+#########          PANEL 2
+###############################################################
+Panel = PanelPlot(2,3, [25,14],'Domaine 1 SV 002', titlepad=25, minmaxpad=1.04, timepad=-0.07, colorbarpad=0.03, labelcolorbarpad = 13, colorbaraspect=22)
+Lplot = [Dvar['f1']['SVPP002'], Dvar['f2']['SVPP002'], Dvar['f3']['SVPP002'],
+         Dvar['f4']['SVPP002'], Dvar['f5']['SVPP002'], Dvar['f6']['SVPP002']]
+
+lon = [Dvar['f1']['longitude']]*len(Lplot)
+lat = [Dvar['f1']['latitude']]*len(Lplot)
+Ltitle = ['SVPP001']*len(Lplot)
+Lcbarlabel = ['kg/kg']*len(Lplot)
+Lminval = [0]*len(Lplot)
+Lmaxval = [0.15E-5]*len(Lplot)
+Lstep = [0.05E-6]*len(Lplot)
+Lstepticks = [0.2E-6]*len(Lplot) 
+Lcbformatlabel=[True]*len(Lplot)
+fig = Panel.psectionH(lon=lon, lat=lat, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, 
+                                Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, 
+                                Ltime=Ltime, LaddWhite_cm=LaddWhite, Lproj=Lprojection, Llevel=Llevel, Lcbformatlabel=Lcbformatlabel)
+
+fig = Panel.pvector(Lxx=lon, Lyy=lat, Llevel=Llevel, Lvar1=Lplot1, Lvar2=Lplot2, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lwidth=Lwidth, Larrowstep=Larrowstep,
+                      Llegendval=Llegendval, Llegendlabel=Lcbarlabel, Lproj=Lprojection, Lid_overlap=[0,2,4,6,8,10], ax=fig.axes, Lscale=Lscale)
+
+# Add departements to each axes + scatter point of emission source
+for i in range(len(Lplot)):
+  fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='white', alpha=0.2)
+  fig.axes[i*2].scatter(loncar,latcar)
+  for lab, txt in enumerate(label):
+    fig.axes[i*2].annotate(label[lab], (loncar[lab], latcar[lab]), color='black',size=10, weight="bold")
+
+# Add a Rectangle displaying the domain of the model 2
+for i in range(len(Lplot)):
+  fig.axes[i*2].add_patch(mpatches.Rectangle(xy=[Dvar['f7']['LONOR'], Dvar['f7']['LATOR']], width=Dvar['f7']['LON'][-1,-1]-Dvar['f7']['LONOR'], height=Dvar['f7']['LAT'][-1,-1]-Dvar['f7']['LATOR'],
+   facecolor='blue', alpha=0.15, transform=ccrs.PlateCarree()))
+
+fig.tight_layout()
+Panel.save_graph(2,fig)
+
+################################################################
+#########          PANEL 3
+###############################################################
+Panel = PanelPlot(2,3, [25,14],'Domaine 1 ATC 001', titlepad=25, minmaxpad=1.04, timepad=-0.07, colorbarpad=0.03, labelcolorbarpad = 13, colorbaraspect=22)
+Lplot = [Dvar['f1']['ATC001'], Dvar['f2']['ATC001'], Dvar['f3']['ATC001'],
+         Dvar['f4']['ATC001'], Dvar['f5']['ATC001'], Dvar['f6']['ATC001']]
+
+Ltitle = ['ATC001']*len(Lplot)
+Lcbarlabel = ['$m^{-3}$']*len(Lplot)
+Lminval = [0]*len(Lplot)
+Lmaxval = [0.15E-6]*len(Lplot)
+Lstep = [0.05E-7]*len(Lplot)
+Lstepticks = [0.2E-7]*len(Lplot) 
+fig = Panel.psectionH(lon=lon, lat=lat, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, 
+                                Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, 
+                                Ltime=Ltime, LaddWhite_cm=LaddWhite, Lproj=Lprojection, Llevel=Llevel, Lcbformatlabel=Lcbformatlabel)
+
+# Add departements to each axes + scatter point of emission source
+for i in range(len(Lplot)):
+  fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='white', alpha=0.2)
+  fig.axes[i*2].scatter(loncar,latcar)
+  for lab, txt in enumerate(label):
+    fig.axes[i*2].annotate(label[lab], (loncar[lab], latcar[lab]), color='black',size=10, weight="bold")
+
+# Add a Rectangle displaying the domain of the model 2
+for i in range(len(Lplot)):
+  fig.axes[i*2].add_patch(mpatches.Rectangle(xy=[Dvar['f7']['LONOR'], Dvar['f7']['LATOR']], width=Dvar['f7']['LON'][-1,-1]-Dvar['f7']['LONOR'], height=Dvar['f7']['LAT'][-1,-1]-Dvar['f7']['LATOR'],
+   facecolor='blue', alpha=0.15, transform=ccrs.PlateCarree()))
+
+fig.tight_layout()
+Panel.save_graph(3,fig)
+
+################################################################
+#########          PANEL 4
+###############################################################
+Panel = PanelPlot(2,3, [25,14],'Domaine 1 ATC 002', titlepad=25, minmaxpad=1.04, timepad=-0.07, colorbarpad=0.03, labelcolorbarpad = 13, colorbaraspect=22)
+Lplot = [Dvar['f1']['ATC002'], Dvar['f2']['ATC002'], Dvar['f3']['ATC002'],
+         Dvar['f4']['ATC002'], Dvar['f5']['ATC002'], Dvar['f6']['ATC002']]
+
+Ltitle = ['ATC002']*len(Lplot)
+Lcbarlabel = ['$m^{-3}$']*len(Lplot)
+Lminval = [0]*len(Lplot)
+Lmaxval = [0.15E-6]*len(Lplot)
+Lstep = [0.05E-7]*len(Lplot)
+Lstepticks = [0.2E-7]*len(Lplot) 
+fig = Panel.psectionH(lon=lon, lat=lat, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, 
+                                Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, 
+                                Ltime=Ltime, LaddWhite_cm=LaddWhite, Lproj=Lprojection, Llevel=Llevel, Lcbformatlabel=Lcbformatlabel)
+
+# Add departements to each axes + scatter point of emission source
+for i in range(len(Lplot)):
+  fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='white', alpha=0.2)
+  fig.axes[i*2].scatter(loncar,latcar)
+  for lab, txt in enumerate(label):
+    fig.axes[i*2].annotate(label[lab], (loncar[lab], latcar[lab]), color='black',size=10, weight="bold")
+
+# Add a Rectangle displaying the domain of the model 2
+for i in range(len(Lplot)):
+  fig.axes[i*2].add_patch(mpatches.Rectangle(xy=[Dvar['f7']['LONOR'], Dvar['f7']['LATOR']], width=Dvar['f7']['LON'][-1,-1]-Dvar['f7']['LONOR'], height=Dvar['f7']['LAT'][-1,-1]-Dvar['f7']['LATOR'],
+   facecolor='blue', alpha=0.15, transform=ccrs.PlateCarree()))
+
+fig.tight_layout()
+Panel.save_graph(4,fig)
+################################################################
+#########          PANEL 5 : Domaine fils
+###############################################################
+Panel = PanelPlot(2,3, [25,14],'Domaine 2 SV 001', titlepad=25, minmaxpad=1.04, timepad=-0.07, colorbarpad=0.03, labelcolorbarpad = 13, colorbaraspect=18)
+Lplot = [Dvar['f7']['SVPP001'], Dvar['f8']['SVPP001'], Dvar['f9']['SVPP001'],
+         Dvar['f10']['SVPP001'], Dvar['f11']['SVPP001'], Dvar['f12']['SVPP001']]
+
+lon = [Dvar['f7']['longitude']]*len(Lplot)
+lat = [Dvar['f7']['latitude']]*len(Lplot)
+Ltitle = ['SVPP001']*len(Lplot)
+Lcbarlabel = ['kg/kg']*len(Lplot)
+Lxlab = ['longitude']*len(Lplot)
+Lylab = ['latitude']*len(Lplot)
+Lminval = [0]*len(Lplot)
+Lmaxval = [0.15E-5]*len(Lplot)
+Lstep = [0.05E-6]*len(Lplot)
+Lstepticks = [0.2E-6]*len(Lplot) 
+Lcolormap = ['gist_rainbow_r']*len(Lplot)
+Lprojection = [ccrs.PlateCarree()]*len(Lplot)
+LaddWhite = [True]*len(Lplot)
+Llevel = [0]*len(Lplot)
+Ltime = [Dvar['f7']['date'], Dvar['f8']['date'], Dvar['f9']['date'], Dvar['f10']['date'], Dvar['f11']['date'], Dvar['f12']['date']]
+Lcbformatlabel=[True]*len(Lplot)
+fig = Panel.psectionH(lon=lon, lat=lat, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, 
+                                Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, 
+                                Ltime=Ltime, LaddWhite_cm=LaddWhite, Lproj=Lprojection, Llevel=Llevel, Lcbformatlabel=Lcbformatlabel)
+
+Lplot1 = [ Dvar['f7']['UT'], Dvar['f8']['UT'], Dvar['f9']['UT'], Dvar['f10']['UT'], Dvar['f11']['UT'], Dvar['f12']['UT']]
+Lplot2 = [ Dvar['f7']['VT'], Dvar['f8']['VT'], Dvar['f9']['VT'], Dvar['f10']['VT'], Dvar['f11']['VT'], Dvar['f12']['VT']]
+Ltitle = ['wind vectors at K=2']*len(Lplot)
+Llegendval = [7.5]*len(Lplot)
+Lcbarlabel = ['(m/s)']*len(Lplot1)
+Larrowstep = [4]*len(Lplot1)
+Lwidth = [0.002]*len(Lplot1)
+Lcolor = ['black']*len(Lplot1)
+Lscale = [75]*len(Lplot1)
+fig = Panel.pvector(Lxx=lon, Lyy=lat, Llevel=Llevel, Lvar1=Lplot1, Lvar2=Lplot2, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lwidth=Lwidth, Larrowstep=Larrowstep,
+                      Llegendval=Llegendval, Llegendlabel=Lcbarlabel, Lproj=Lprojection, Lid_overlap=[0,2,4,6,8,10], ax=fig.axes, Lscale=Lscale)
+
+# Departements francais
+departements_shp='departements-20180101.shp'
+adm1_shapes = list(shpreader.Reader(departements_shp).geometries())
+
+# Add departements to each axes + scatter point of emission source
+loncar, latcar = [1.439,1.5], [43.567, 43.9]
+label=['AZF1','AZF2']
+for i in range(len(Lplot)):
+  fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='white', alpha=0.2)
+  fig.axes[i*2].scatter(loncar,latcar)
+  for lab, txt in enumerate(label):
+    fig.axes[i*2].annotate(label[lab], (loncar[lab], latcar[lab]), color='black',size=10, weight="bold")
+
+fig.tight_layout()
+Panel.save_graph(5,fig)
+
+################################################################
+#########          PANEL 6
+###############################################################
+Panel = PanelPlot(2,3, [25,14],'Domaine 2 ATC 001', titlepad=25, minmaxpad=1.04, timepad=-0.07, colorbarpad=0.03, labelcolorbarpad = 13, colorbaraspect=18)
+Lplot = [Dvar['f7']['ATC001'], Dvar['f8']['ATC001'], Dvar['f9']['ATC001'],
+         Dvar['f10']['ATC001'], Dvar['f11']['ATC001'], Dvar['f12']['ATC001']]
+Ltitle = ['ATC001']*len(Lplot)
+Lcbarlabel = ['$m^{-3}$']*len(Lplot)
+Lminval = [0]*len(Lplot)
+Lmaxval = [0.6E-6]*len(Lplot)
+Lstep = [0.01E-6]*len(Lplot)
+Lstepticks = [0.1E-6]*len(Lplot) 
+fig = Panel.psectionH(lon=lon, lat=lat, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, 
+                                Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, 
+                                Ltime=Ltime, LaddWhite_cm=LaddWhite, Lproj=Lprojection, Llevel=Llevel, Lcbformatlabel=Lcbformatlabel)
+
+# Add departements to each axes + scatter point of emission source
+for i in range(len(Lplot)):
+  fig.axes[i*2].add_geometries(adm1_shapes, ccrs.PlateCarree(),edgecolor='black', facecolor='white', alpha=0.2)
+  fig.axes[i*2].scatter(loncar,latcar)
+  for lab, txt in enumerate(label):
+    fig.axes[i*2].annotate(label[lab], (loncar[lab], latcar[lab]), color='black',size=10, weight="bold")
+
+fig.tight_layout()
+Panel.save_graph(6,fig)
+
+################################################################
+#########          PANEL 7
+###############################################################
+Panel = PanelPlot(8,2, [14,20],'Time series from Aircraft', titlepad=25, minmaxpad=1.04, timepad=-0.07, colorbarpad=0.03, labelcolorbarpad = 13, colorbaraspect=18)
+
+Lplot = [ Dvar['f13'][(LG_AVIONT,'ZS')]]
+Ltime = [Dvar['f13'][(LG_AVION,'time_flyer')]/3600.0]
+Ltitle = ['Orography']
+Lxlab = ['Time (h)']
+Lylab = ['ZS (m)']
+Lylim = [(0, 350)]
+Lxlim = [(9.0, 9.2)]
+fig = Panel.pXY_lines(Lyy=Lplot, Lxx=Ltime, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lylim=Lylim, Lxlim=Lxlim)
+
+Lplot = [ Dvar['f13'][(LG_AVIONT,'P')]]
+Ltitle = ['Pressure']
+Lylab = ['P (Pa)']
+Lylim = [(0, 95000)]
+fig = Panel.pXY_lines(Lyy=Lplot, Lxx=Ltime, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lylim=Lylim, Lxlim=Lxlim, ax=fig.axes) 
+
+Lplot = [ Dvar['f13'][(LG_AVIONT,'LON')]]
+Ltitle = ['Longitude']
+Lylab = ['Longitude']
+Lylim = [(0, 2.5)]
+fig = Panel.pXY_lines(Lyy=Lplot, Lxx=Ltime, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lylim=Lylim, Lxlim=Lxlim, ax=fig.axes)
+
+Lplot = [ Dvar['f13'][(LG_AVIONT,'ZON_WIND')]]
+Ltitle = ['Zonal wind']
+Lylab = ['u (m/s)']
+Lylim = [(-1, 11)]
+fig = Panel.pXY_lines(Lyy=Lplot, Lxx=Ltime, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lylim=Lylim, Lxlim=Lxlim, ax=fig.axes)
+
+Lplot = [ Dvar['f13'][(LG_AVIONT,'MER_WIND')]]
+Ltitle = ['Meridional wind']
+Lylab = ['v (m/s)']
+Lylim = [(-3, 3)]
+fig = Panel.pXY_lines(Lyy=Lplot, Lxx=Ltime, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lylim=Lylim, Lxlim=Lxlim, ax=fig.axes)
+
+Lplot = [ Dvar['f13'][(LG_AVIONT,'W')]]
+Ltitle = ['Vertical velocity']
+Lylab = ['w (m/s)']
+Lylim = [(-0.1, 0.1)]
+fig = Panel.pXY_lines(Lyy=Lplot, Lxx=Ltime, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lylim=Lylim, Lxlim=Lxlim, ax=fig.axes)
+
+Lplot = [ Dvar['f13'][(LG_AVIONT,'Th')]]
+Ltitle = ['Potential Temperature']
+Lylab = [r'$\theta$ (K)']
+Lylim = [(290, 305)]
+fig = Panel.pXY_lines(Lyy=Lplot, Lxx=Ltime, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lylim=Lylim, Lxlim=Lxlim, ax=fig.axes)
+
+Lplot = [ Dvar['f13'][(LG_AVIONT,'Rv')]]
+Ltitle = ['Water vapor mixing ratio']
+Lylab = ['Rv (kg/kg))']
+Lylim = [(0, 0.01)]
+fig = Panel.pXY_lines(Lyy=Lplot, Lxx=Ltime, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lylim=Lylim, Lxlim=Lxlim, ax=fig.axes)
+
+Lplot = [ Dvar['f13'][(LG_AVIONT,'Tke')]]
+Ltitle = ['Turbulent Kinetic Energy']
+Lylab = ['TKE ($m^2s^{-2}$)']
+Lylim = [(0, 0.1)]
+fig = Panel.pXY_lines(Lyy=Lplot, Lxx=Ltime, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lylim=Lylim, Lxlim=Lxlim, ax=fig.axes)
+
+Lplot = [ Dvar['f13'][(LG_AVIONT,'Tke_Diss')]]
+Ltitle = ['Turbulent Kinetic Energy Dissipation']
+Lylab = ['TKE Diss ($m^2s^{-2}$']
+Lylim = [(0, 1000)]
+fig = Panel.pXY_lines(Lyy=Lplot, Lxx=Ltime, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lylim=Lylim, Lxlim=Lxlim, ax=fig.axes)
+
+Lplot = [ Dvar['f13'][(LG_AVIONT,'H_FLUX')]]
+Ltitle = ['Sensible Heat Flux H']
+Lylab = ['H ($W/m^2$)']
+Lylim = [(-0.7, 0.)]
+fig = Panel.pXY_lines(Lyy=Lplot, Lxx=Ltime, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lylim=Lylim, Lxlim=Lxlim, ax=fig.axes)
+
+Lplot = [ Dvar['f13'][(LG_AVIONT,'LE_FLUX')]]
+Ltitle = ['Latent Heat Flux LE']
+Lylab = ['LE ($W/m^2$)']
+Lylim = [(0, 2.0)]
+fig = Panel.pXY_lines(Lyy=Lplot, Lxx=Ltime, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lylim=Lylim, Lxlim=Lxlim, ax=fig.axes)
+
+Lplot = [ Dvar['f13'][(LG_AVIONT,'Tsrad')]]
+Ltitle = ['Radiative surface temperature']
+Lylab = ['Tsrad (K))']
+Lylim = [(250, 1000)]
+fig = Panel.pXY_lines(Lyy=Lplot, Lxx=Ltime, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lylim=Lylim, Lxlim=Lxlim, ax=fig.axes)
+
+fig.tight_layout()
+Panel.save_graph(7,fig)

MY_RUN/INTEGRATION_CASES/HPC/AZF_2M/009_python/run_python

+#!/bin/sh
+#SBATCH -J tools
+#SBATCH -N 1            # nodes number
+#SBATCH -n 1            # CPUs number (on all nodes) 
+#SBATCH -o tools.eo%j   #
+#SBATCH -e tools.eo%j   #
+#SBATCH -t 02: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.* .
+ln -sf ../007_run/AZF02.*.CEN4T.*.nc .
+module purge
+module load python/3.7.6
+
+python3 plot_AZF2M.py
+convert *.png AZF_2M.pdf

MY_RUN/INTEGRATION_CASES/HPC/BIOMAIDO_DMS/001_prep_pgd1/PRE_PGD1.nam

+&NAM_CONFIO LCDF4=.TRUE.,
+            LLFIOUT=.FALSE.,
+            LLFIREAD=.FALSE. /
+
+&NAM_PGDFILE CPGDFILE='BIOMAIDO_CH_8km', NHALO=0 /
+
+&NAM_COVER YCOVER='ECOCLIMAP_v2.0',
+           YCOVERFILETYPE='DIRECT' /
+
+&NAM_ZS YZS='gtopo30',
+        YZSFILETYPE='DIRECT' /
+
+&NAM_ZSFILTER NZSFILTER=0 /
+
+&NAM_PGD_GRID CGRID='CONF PROJ' /
+
+&NAM_CONF_PROJ XLAT0=-20.93,
+               XLON0=55.370,
+               XRPK=0.,
+               XBETA=0. /
+
+&NAM_CONF_PROJ_GRID XLATCEN=-20.93,
+                    XLONCEN=55.370,
+                    NIMAX=60,
+                    NJMAX=60,
+                    XDX=8000.,
+                    XDY=8000. /
+
+&NAM_ISBA CPHOTO='NIT',
+          YCLAY='CLAY_HWSD_MOY_v2',
+          YCLAYFILETYPE='DIRECT',
+          YSAND='SAND_HWSD_MOY_v2',
+          YSANDFILETYPE='DIRECT',
+          CISBA ="3-L",
+          NPATCH=12 /
+
+&NAM_CH_EMISSIONS CCH_EMIS = 'AGGR', CCH_BIOEMIS = 'MEGA', CCH_DMSEMIS= 'DMSD' /
+
+&NAM_DMS_PGD    NDMS_NBR = 12,
+                CDMS_NAME(1) = 'DMS_JAN',
+                CDMS_FILE(1) = 'data_DMS_january.dat',
+                CDMS_FILETYPE(1) = 'ASCLLV',
+                CDMS_AREA(1) = 'SEA',
+                CDMS_NAME(2) = 'DMS_FEV',
+                CDMS_FILE(2) = 'data_DMS_february.dat',
+                CDMS_FILETYPE(2) = 'ASCLLV',
+                CDMS_AREA(2) = 'SEA',
+                CDMS_NAME(3) = 'DMS_MAR',
+                CDMS_FILE(3) = 'data_DMS_march.dat',
+                CDMS_FILETYPE(3) = 'ASCLLV',
+                CDMS_AREA(3) = 'SEA',
+                CDMS_NAME(4) = 'DMS_APR',
+                CDMS_FILE(4) = 'data_DMS_april.dat',
+                CDMS_FILETYPE(4) = 'ASCLLV',
+                CDMS_AREA(4) = 'SEA',
+                CDMS_NAME(5) = 'DMS_MAY',
+                CDMS_FILE(5) = 'data_DMS_may.dat',
+                CDMS_FILETYPE(5) = 'ASCLLV',
+                CDMS_AREA(5) = 'SEA',
+                CDMS_NAME(6) = 'DMS_JUN',
+                CDMS_FILE(6) = 'data_DMS_june.dat',
+                CDMS_FILETYPE(6) = 'ASCLLV',
+                CDMS_AREA(6) = 'SEA',
+                CDMS_NAME(7) = 'DMS_JUL',
+                CDMS_FILE(7) = 'data_DMS_july.dat',
+                CDMS_FILETYPE(7) = 'ASCLLV',
+                CDMS_AREA(7) = 'SEA',
+                CDMS_NAME(8) = 'DMS_AUG',
+                CDMS_FILE(8) = 'data_DMS_august.dat',
+                CDMS_FILETYPE(8) = 'ASCLLV',
+                CDMS_AREA(9) = 'SEA',
+                CDMS_NAME(9) = 'DMS_SEP',
+                CDMS_FILE(9) = 'data_DMS_september.dat',
+                CDMS_FILETYPE(9) = 'ASCLLV',
+                CDMS_AREA(10) = 'SEA',
+                CDMS_NAME(10) = 'DMS_OCT',
+                CDMS_FILE(10) = 'data_DMS_october.dat',
+                CDMS_FILETYPE(10) = 'ASCLLV',
+                CDMS_AREA(11) = 'SEA',
+                CDMS_NAME(11) = 'DMS_NOV',
+                CDMS_FILE(11) = 'data_DMS_november.dat',
+                CDMS_FILETYPE(11) = 'ASCLLV',
+                CDMS_AREA(11) = 'SEA',
+                CDMS_NAME(12) = 'DMS_DEC',
+                CDMS_FILE(12) = 'data_DMS_december.dat',
+                CDMS_FILETYPE(12) = 'ASCLLV',
+                CDMS_AREA(12) = 'SEA' /
+
+
+&NAM_CH_EMIS_PGD    NEMIS_PGD_NBR = 38, 
+                    CEMIS_PGD_NAME(1) = 'COAN',
+                    NEMIS_PGD_TIME(1) = 0,
+                    CEMIS_PGD_FILE(1) = 'MACCity_CO_juil2014.asc',
+                    CEMIS_PGD_FILETYPE(1) = 'ASCLLV',
+                    CEMIS_PGD_AREA(1) = 'LAN',
+                    CEMIS_PGD_NAME(2) = 'NOXAN',
+                    NEMIS_PGD_TIME(2) = 0,
+                    CEMIS_PGD_FILE(2) = 'MACCity_NOX_juil2014.asc',
+                    CEMIS_PGD_FILETYPE(2) = 'ASCLLV',
+                    CEMIS_PGD_AREA(2) = 'LAN',
+                    CEMIS_PGD_NAME(3) = 'C2H6AN',
+                    NEMIS_PGD_TIME(3) = 0,
+                    CEMIS_PGD_FILE(3) = 'MACCity_C2H6_juil2014.asc',
+                    CEMIS_PGD_FILETYPE(3) = 'ASCLLV',
+                    CEMIS_PGD_AREA(3) = 'LAN',
+                    CEMIS_PGD_NAME(4) = 'C2H4AN',
+                    NEMIS_PGD_TIME(4) = 0,
+                    CEMIS_PGD_FILE(4) = 'MACCity_C2H4_juil2014.asc',
+                    CEMIS_PGD_FILETYPE(4) = 'ASCLLV',
+                    CEMIS_PGD_AREA(4) = 'LAN',
+                    CEMIS_PGD_NAME(5) = 'C3H8AN',
+                    NEMIS_PGD_TIME(5) = 0,
+                    CEMIS_PGD_FILE(5) = 'MACCity_C3H8_juil2014.asc',
+                    CEMIS_PGD_FILETYPE(5) = 'ASCLLV',
+                    CEMIS_PGD_AREA(5) = 'LAN',
+                    CEMIS_PGD_NAME(6) = 'C3H6AN',
+                    NEMIS_PGD_TIME(6) = 0,
+                    CEMIS_PGD_FILE(6) = 'MACCity_C3H6_juil2014.asc',
+                    CEMIS_PGD_FILETYPE(6) = 'ASCLLV',
+                    CEMIS_PGD_AREA(6) = 'LAN',
+                    CEMIS_PGD_NAME(7) = 'BALKAN',
+                    NEMIS_PGD_TIME(7) = 0,
+                    CEMIS_PGD_FILE(7) = 'MACCity_BIGALK_juil2014.asc',
+                    CEMIS_PGD_FILETYPE(7) = 'ASCLLV',
+                    CEMIS_PGD_AREA(7) = 'LAN',
+                    CEMIS_PGD_NAME(8) = 'BENEAN',
+                    NEMIS_PGD_TIME(8) = 0,
+                    CEMIS_PGD_FILE(8) = 'MACCity_BIGENE_juil2014.asc',
+                    CEMIS_PGD_FILETYPE(8) = 'ASCLLV',
+                    CEMIS_PGD_AREA(8) = 'LAN',
+                    CEMIS_PGD_NAME(9) = 'CH2OAN',
+                    NEMIS_PGD_TIME(9) = 0,
+                    CEMIS_PGD_FILE(9) = 'MACCity_CH2O_juil2014.asc',
+                    CEMIS_PGD_FILETYPE(9) = 'ASCLLV',
+                    CEMIS_PGD_AREA(9) = 'LAN',
+                    CEMIS_PGD_NAME(10) = 'BALDAN',
+                    NEMIS_PGD_TIME(10) = 0,
+                    CEMIS_PGD_FILE(10) = 'MACCity_BIGALD_juil2014.asc',
+                    CEMIS_PGD_FILETYPE(10) = 'ASCLLV',
+                    CEMIS_PGD_AREA(10) = 'LAN',
+                    CEMIS_PGD_NAME(11) = 'MTOLAN',
+                    NEMIS_PGD_TIME(11) = 0,
+                    CEMIS_PGD_FILE(11) = 'MACCity_CH3OH_juil2014.asc',
+                    CEMIS_PGD_FILETYPE(11) = 'ASCLLV',
+                    CEMIS_PGD_AREA(11) = 'LAN',
+                    CEMIS_PGD_NAME(12) = 'BOLAN',
+                    NEMIS_PGD_TIME(12) = 0,
+                    CEMIS_PGD_FILE(12) = 'MACCity_BIGOL_juil2014.asc',
+                    CEMIS_PGD_FILETYPE(12) = 'ASCLLV',
+                    CEMIS_PGD_AREA(12) = 'LAN',
+                    CEMIS_PGD_NAME(13) = 'ACETAN',
+                    NEMIS_PGD_TIME(13) = 0,
+                    CEMIS_PGD_FILE(13) = 'MACCity_ACETONE_juil2014.asc',
+                    CEMIS_PGD_FILETYPE(13) = 'ASCLLV',
+                    CEMIS_PGD_AREA(13) = 'LAN',
+                    CEMIS_PGD_NAME(14) = 'BONEAN',
+                    NEMIS_PGD_TIME(14) = 0,
+                    CEMIS_PGD_FILE(14) = 'MACCity_BIGONE_juil2014.asc',
+                    CEMIS_PGD_FILETYPE(14) = 'ASCLLV',
+                    CEMIS_PGD_AREA(14) = 'LAN',
+                    CEMIS_PGD_NAME(15) = 'AROAN',
+                    NEMIS_PGD_TIME(15) = 0,
+                    CEMIS_PGD_FILE(15) = 'MACCity_ARO_juil2014.asc',
+                    CEMIS_PGD_FILETYPE(15) = 'ASCLLV',
+                    CEMIS_PGD_AREA(15) = 'LAN',
+                    CEMIS_PGD_NAME(16) = 'SO2AN',
+                    NEMIS_PGD_TIME(16) = 0,
+                    CEMIS_PGD_FILE(16) = 'MACCity_SO2_juil2014.asc',
+                    CEMIS_PGD_FILETYPE(16) = 'ASCLLV',
+                    CEMIS_PGD_AREA(16) = 'LAN',
+                    NEMIS_PGD_TIME(17) = 0,
+                    CEMIS_PGD_NAME(17) = 'NH3AN',
+                    CEMIS_PGD_FILE(17) = 'MACCity_NH3_juil2014.asc',
+                    CEMIS_PGD_FILETYPE(17) = 'ASCLLV',
+                    CEMIS_PGD_AREA(17) = 'LAN',
+                    CEMIS_PGD_NAME(18) = 'BCAN',
+                    NEMIS_PGD_TIME(18) = 0,
+                    CEMIS_PGD_FILE(18) = 'MACCity_BC_juil2014.asc',
+                    CEMIS_PGD_FILETYPE(18) = 'ASCLLV',
+                    CEMIS_PGD_AREA(18) = 'LAN',
+                    CEMIS_PGD_NAME(19) = 'OCAN',
+                    NEMIS_PGD_TIME(19) = 0,
+                    CEMIS_PGD_FILE(19) = 'MACCity_OC_juil2014.asc',
+                    CEMIS_PGD_FILETYPE(19) = 'ASCLLV',
+                    CEMIS_PGD_AREA(19) = 'LAN',
+                    CEMIS_PGD_NAME(20) = 'COANSEA',
+                    NEMIS_PGD_TIME(20) = 0,
+                    CEMIS_PGD_FILE(20) = 'MACCity_CO_juil2014.asc',
+                    CEMIS_PGD_FILETYPE(20) = 'ASCLLV',
+                    CEMIS_PGD_AREA(20) = 'SEA',
+                    CEMIS_PGD_NAME(21) = 'NOXANSEA',
+                    NEMIS_PGD_TIME(21) = 0,
+                    CEMIS_PGD_FILE(21) = 'MACCity_NOX_juil2014.asc',
+                    CEMIS_PGD_FILETYPE(21) = 'ASCLLV',
+                    CEMIS_PGD_AREA(21) = 'SEA',
+                    CEMIS_PGD_NAME(22) = 'C2H6ANSEA',
+                    NEMIS_PGD_TIME(22) = 0,
+                    CEMIS_PGD_FILE(22) = 'MACCity_C2H6_juil2014.asc',
+                    CEMIS_PGD_FILETYPE(22) = 'ASCLLV',
+                    CEMIS_PGD_AREA(22) = 'SEA',
+                    CEMIS_PGD_NAME(23) = 'C2H4ANSEA',
+                    NEMIS_PGD_TIME(23) = 0,
+                    CEMIS_PGD_FILE(23) = 'MACCity_C2H4_juil2014.asc',
+                    CEMIS_PGD_FILETYPE(23) = 'ASCLLV',
+                    CEMIS_PGD_AREA(23) = 'SEA',
+                    CEMIS_PGD_NAME(24) = 'C3H8ANSEA',
+                    NEMIS_PGD_TIME(24) = 0,
+                    CEMIS_PGD_FILE(24) = 'MACCity_C3H8_juil2014.asc',
+                    CEMIS_PGD_FILETYPE(24) = 'ASCLLV',
+                    CEMIS_PGD_AREA(24) = 'SEA',
+                    CEMIS_PGD_NAME(25) = 'C3H6ANSEA',
+                    NEMIS_PGD_TIME(25) = 0,
+                    CEMIS_PGD_FILE(25) = 'MACCity_C3H6_juil2014.asc',
+                    CEMIS_PGD_FILETYPE(25) = 'ASCLLV',
+                    CEMIS_PGD_AREA(25) = 'SEA',
+                    CEMIS_PGD_NAME(26) = 'BALKANSEA',
+                    NEMIS_PGD_TIME(26) = 0,
+                    CEMIS_PGD_FILE(26) = 'MACCity_BIGALK_juil2014.asc',
+                    CEMIS_PGD_FILETYPE(26) = 'ASCLLV',
+                    CEMIS_PGD_AREA(26) = 'SEA',
+                    CEMIS_PGD_NAME(27) = 'BENEANSEA',
+                    NEMIS_PGD_TIME(27) = 0,
+                    CEMIS_PGD_FILE(27) = 'MACCity_BIGENE_juil2014.asc',
+                    CEMIS_PGD_FILETYPE(27) = 'ASCLLV',
+                    CEMIS_PGD_AREA(27) = 'SEA',
+                    CEMIS_PGD_NAME(28) = 'CH2OANSEA',
+                    NEMIS_PGD_TIME(28) = 0,
+                    CEMIS_PGD_FILE(28) = 'MACCity_CH2O_juil2014.asc',
+                    CEMIS_PGD_FILETYPE(28) = 'ASCLLV',
+                    CEMIS_PGD_AREA(28) = 'SEA',
+                    CEMIS_PGD_NAME(29) = 'BALDANSEA',
+                    NEMIS_PGD_TIME(29) = 0,
+                    CEMIS_PGD_FILE(29) = 'MACCity_BIGALD_juil2014.asc',
+                    CEMIS_PGD_FILETYPE(29) = 'ASCLLV',
+                    CEMIS_PGD_AREA(29) = 'SEA',
+                    CEMIS_PGD_NAME(30) = 'MTOLANSEA',
+                    NEMIS_PGD_TIME(30) = 0,
+                    CEMIS_PGD_FILE(30) = 'MACCity_CH3OH_juil2014.asc',
+                    CEMIS_PGD_FILETYPE(30) = 'ASCLLV',
+                    CEMIS_PGD_AREA(30) = 'SEA',
+                    CEMIS_PGD_NAME(31) = 'BOLANSEA',
+                    NEMIS_PGD_TIME(31) = 0,
+                    CEMIS_PGD_FILE(31) = 'MACCity_BIGOL_juil2014.asc',
+                    CEMIS_PGD_FILETYPE(31) = 'ASCLLV',
+                    CEMIS_PGD_AREA(31) = 'SEA',
+                    CEMIS_PGD_NAME(32) = 'ACETANSEA',
+                    NEMIS_PGD_TIME(32) = 0,
+                    CEMIS_PGD_FILE(32) = 'MACCity_ACETONE_juil2014.asc',
+                    CEMIS_PGD_FILETYPE(32) = 'ASCLLV',
+                    CEMIS_PGD_AREA(32) = 'SEA',
+                    CEMIS_PGD_NAME(33) = 'BONEANSEA',
+                    NEMIS_PGD_TIME(33) = 0,
+                    CEMIS_PGD_FILE(33) = 'MACCity_BIGONE_juil2014.asc',
+                    CEMIS_PGD_FILETYPE(33) = 'ASCLLV',
+                    CEMIS_PGD_AREA(33) = 'SEA',
+                    CEMIS_PGD_NAME(34) = 'AROANSEA',
+                    NEMIS_PGD_TIME(34) = 0,
+                    CEMIS_PGD_FILE(34) = 'MACCity_ARO_juil2014.asc',
+                    CEMIS_PGD_FILETYPE(34) = 'ASCLLV',
+                    CEMIS_PGD_AREA(34) = 'SEA',
+                    CEMIS_PGD_NAME(35) = 'SO2ANSEA',
+                    NEMIS_PGD_TIME(35) = 0,
+                    CEMIS_PGD_FILE(35) = 'MACCity_SO2_juil2014.asc',
+                    CEMIS_PGD_FILETYPE(35) = 'ASCLLV',
+                    CEMIS_PGD_AREA(35) = 'SEA',
+                    CEMIS_PGD_NAME(36) = 'NH3ANSEA',
+                    NEMIS_PGD_TIME(36) = 0,
+                    CEMIS_PGD_FILE(36) = 'MACCity_NH3_juil2014.asc',
+                    CEMIS_PGD_FILETYPE(36) = 'ASCLLV',
+                    CEMIS_PGD_AREA(36) = 'SEA',
+                    CEMIS_PGD_NAME(37) = 'BCANSEA',
+                    NEMIS_PGD_TIME(37) = 0,
+                    CEMIS_PGD_FILE(37) = 'MACCity_BC_juil2014.asc',
+                    CEMIS_PGD_FILETYPE(37) = 'ASCLLV',
+                    CEMIS_PGD_AREA(37) = 'SEA',
+                    CEMIS_PGD_NAME(38) = 'OCANSEA',
+                    NEMIS_PGD_TIME(38) = 0,
+                    CEMIS_PGD_FILE(38) = 'MACCity_OC_juil2014.asc',
+                    CEMIS_PGD_FILETYPE(38) = 'ASCLLV',
+                    CEMIS_PGD_AREA(38) = 'SEA' /
+
+&NAM_MEGAN_PGD    NMEGAN_NBR = 15,
+                  CMEGAN_NAME(1) = 'EFISOP',
+                  CMEGAN_FILE(1) = 'isop100.asc',
+                  CMEGAN_FILETYPE(1) = 'ASCLLV',
+                  CMEGAN_AREA(1) = 'LAN',
+                  CMEGAN_NAME(2) = 'EFAPIN',
+                  CMEGAN_FILE(2) = 'apin100.asc',
+                  CMEGAN_FILETYPE(2) = 'ASCLLV',
+                  CMEGAN_AREA(2) = 'LAN',
+                  CMEGAN_NAME(3) = 'EFBPIN',
+                  CMEGAN_FILE(3) = 'bpin100.asc',
+                  CMEGAN_FILETYPE(3) = 'ASCLLV',
+                  CMEGAN_AREA(3) = 'LAN',
+                  CMEGAN_NAME(4) = 'EFCARE',
+                  CMEGAN_FILE(4) = 'bcar100.asc',
+                  CMEGAN_FILETYPE(4) = 'ASCLLV',
+                  CMEGAN_AREA(4) = 'LAN',
+                  CMEGAN_NAME(5) = 'EFLIMO',
+                  CMEGAN_FILE(5) = 'limo100.asc',
+                  CMEGAN_FILETYPE(5) = 'ASCLLV',
+                  CMEGAN_AREA(5) = 'LAN',
+                  CMEGAN_NAME(6) = 'EFMYRC',
+                  CMEGAN_FILE(6) = 'myrc100.asc',
+                  CMEGAN_FILETYPE(6) = 'ASCLLV',
+                  CMEGAN_AREA(6) = 'LAN',
+                  CMEGAN_NAME(7) = 'EFOCIM',
+                  CMEGAN_FILE(7) = 'ocim100.asc',
+                  CMEGAN_FILETYPE(7) = 'ASCLLV',
+                  CMEGAN_AREA(7) = 'LAN',
+                  CMEGAN_NAME(8) = 'EFSABI',
+                  CMEGAN_FILE(8) = 'sabi100.asc',
+                  CMEGAN_FILETYPE(8) = 'ASCLLV',
+                  CMEGAN_AREA(8) = 'LAN',
+                  CMEGAN_NAME(10) = 'LAI',
+                  CMEGAN_FILE(10) = 'lai100.asc',
+                  CMEGAN_FILETYPE(10) = 'ASCLLV',
+                  CMEGAN_AREA(10) = 'LAN'
+                  CMEGAN_NAME(11) = 'PFT1',
+                  CMEGAN_FILE(11) = 'pft1.asc',
+                  CMEGAN_FILETYPE(11) = 'ASCLLV',
+                  CMEGAN_AREA(11) = 'LAN'
+                  CMEGAN_NAME(12) = 'PFT4',
+                  CMEGAN_FILE(12) = 'pft4.asc',
+                  CMEGAN_FILETYPE(12) = 'ASCLLV',
+                  CMEGAN_AREA(12) = 'LAN'
+                  CMEGAN_NAME(13) = 'PFT9',
+                  CMEGAN_FILE(13) = 'pft9.asc',
+                  CMEGAN_FILETYPE(13) = 'ASCLLV',
+                  CMEGAN_AREA(13) = 'LAN'
+                  CMEGAN_NAME(14) = 'PFT13',
+                  CMEGAN_FILE(14) = 'pft13.asc',
+                  CMEGAN_FILETYPE(14) = 'ASCLLV',
+                  CMEGAN_AREA(14) = 'LAN'
+                  CMEGAN_NAME(15) = 'PFT14',
+                  CMEGAN_FILE(15) = 'pft14.asc',
+                  CMEGAN_FILETYPE(15) = 'ASCLLV',
+                  CMEGAN_AREA(15) = 'LAN'
+                  CMEGAN_NAME(9) = 'PFTru',
+                  CMEGAN_FILE(9) = 'pftru.asc',
+                  CMEGAN_FILETYPE(9) = 'ASCLLV',
+                  CMEGAN_AREA(9) = 'LAN' /
+

MY_RUN/INTEGRATION_CASES/HPC/BIOMAIDO_DMS/001_prep_pgd1/run_pgd.sh

+#!/bin/sh
+#SBATCH -J pgd_megan1
+#SBATCH -N 1            # nodes number (=NBP)
+#SBATCH -n 8            # CPUs number (on all nodes) (=NBP*TPN)
+#SBATCH -o run_megan.eo%j   #
+#SBATCH -e run_megan.eo%j   #
+#SBATCH -t 00:30:00    # time limit
+
+# Echo des commandes
+ulimit -c 0
+ulimit -s unlimited
+# Arrete du job des la premiere erreur
+set -e
+# Nom de la machine
+hostname
+
+. ~/DEV_57/MNH-PHYEX070-b95d84d7/conf/profile_mesonh-LXifort-R8I4-MNH-V5-6-2-ECRAD140-MPIAUTO-O2
+ln -sf ~/SAVE/mesonh/PGD/* .
+ln -sf ~/SAVE/EMISSIONS_MEGAN/DATA/DATA_EMIS/* .
+ln -sf ~/SAVE/CHIMIE_FILES/EMISSIONS/* .
+ln -sf ~/SAVE/EMISSIONS_MEGAN/* .
+ln -sf ~/SAVE/EMISSIONS_DMS/* .
+
+export MPIRUN="Mpirun -np 8"
+
+set -x
+set -e
+
+
+time ${MPIRUN} PREP_PGD${XYZ}
+cd ../002_prep_pgd2
+sbatch run_pgd.sh