MY_RUN/INTEGRATION_CASES/HPC/16JAN/RUN/run_spa_real

+#!/bin/sh
+#SBATCH -J sparea_16jan
+#SBATCH -N 1            # nodes number (=NBP)
+#SBATCH -n 1            # CPUs number (on all nodes) (=NBP*TPN)
+#SBATCH -o sparea_16jan.eo%j   #
+#SBATCH -e sparea_16jan.eo%j   #
+#SBATCH -t 00:20: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"
+
+set -x
+set -e
+
+ls -lrt
+
+rm -f 16JA1.1.CEN4T.004.spa04.???
+time ${MPIRUN} SPAWNING${XYZ}
+mv OUTPUT_LISTING0  OUTPUT_LISTING0_spa
+mv OUTPUT_LISTING1  OUTPUT_LISTING1_spa
+mv OUTPUT_LISTING2  OUTPUT_LISTING2_spa
+
+ls -lrt
+
+
+rm -f 16JAN_00_MNH2.???
+cp PRE_REAL1.nam_dom2 PRE_REAL1.nam
+time ${MPIRUN} PREP_REAL_CASE${XYZ}
+mv OUTPUT_LISTING0  OUTPUT_LISTING0_prep
+mv OUTPUT_LISTING1  OUTPUT_LISTING1_prep
+
+ls -lrt
+
+rm -f file_for_xtransfer pipe_name
+rm -f PRE_REAL1.nam
+
+sbatch run_mesonh2
+
+ja

MY_RUN/INTEGRATION_CASES/HPC/ARMCU_LES/DEAR/EXSEG1.nam

+&NAM_CONFIO LCDF4=T LLFIOUT=F LLFIREAD=F /
+&NAM_LUNITn  CINIFILE = "ARM_3D_LF", CINIFILEPGD = "ARM_PGD" /
+&NAM_CONFn  LUSERV = .TRUE./
+&NAM_DYNn XTSTEP = 2.,CPRESOPT="ZRESI"  /
+&NAM_ADVn CUVW_ADV_SCHEME="CEN4TH" ,CTEMP_SCHEME='RKC4'/
+&NAM_PARAMn  CTURB = "TKEL", CRAD = "NONE",
+ CCLOUD = "ICE3", CDCONV="NONE", CSCONV="NONE" /
+&NAM_PARAM_RADn XDTRAD = 120., 
+                XDTRAD_CLONLY = 1., 
+                CLW='RRTM',
+                CAER = 'SURF' /
+&NAM_PARAM_GROUNDN /
+&NAM_SEAFLUXn  /
+&NAM_ISBAn /
+&NAM_DIAG_SURFn /
+&NAM_LBCn CLBCX = 2*"CYCL", CLBCY = 2*"CYCL" /
+&NAM_TURBn XIMPL = 1., CTURBLEN = "DEAR", CTURBDIM = "3DIM" , 
+           LTURB_DIAG = .TRUE., LRMC01=.TRUE.,
+           LTURB_FLX = .TRUE., 
+           LSIG_CONV=.FALSE. /
+&NAM_CH_MNHCn /
+&NAM_CONF  CCONF = "START", LFLAT = .TRUE., NMODEL = 1,NHALO=1
+                   NVERB = 4, CEQNSYS='DUR', lforcing=.true.,
+                   CEXP = "ARM__", CSEG = "CEN4T" /
+&NAM_DYN  XSEGLEN = 43200. , XASSELIN = 0.1, LCORIO = .true., 
+                    LNUMDIFU =T,
+                    XALKTOP = 0.01, XALZBOT = 3500. /
+&NAM_NESTING  NDAD(1) = 1 , NDAD(2) = 1 , NDTRATIO(1) = 1, NDTRATIO(2) = 2 ,
+              XWAY(1) = 1., XWAY(2) = 2. /
+&NAM_BACKUP  XBAK_TIME(1,1)= 43200. /
+&NAM_CONDSAMP LCONDSAMP=T, NCONDSAMP=3 /
+&NAM_LES    
+            LLES_MEAN=.TRUE., 
+            LLES_RESOLVED=.TRUE.,
+            LLES_SUBGRID=.TRUE.,
+            LLES_UPDRAFT=.FALSE.,
+            LLES_DOWNDRAFT=.FALSE.,CBL_HEIGHT_DEF='WTV',
+            CLES_NORM_TYPE='CONV',
+            XLES_TEMP_SAMPLING=300.,
+            LLES_NEB_MASK = .TRUE.,
+            LLES_CS_MASK = .TRUE.,
+            XLES_TEMP_MEAN_STEP=3600.,
+            XLES_TEMP_MEAN_START=3600.,
+            XLES_TEMP_MEAN_END=43200.  /
+&NAM_FRC LGEOST_UV_FRC=.TRUE., LTEND_THRV_FRC=.TRUE. /
+&NAM_CH_SOLVER /
+&NAM_PARAM_ICEn CSUBG_AUCV_RC='CLFR'   /
+&NAM_IDEAL_FLUX 
+NFORCT = 2,
+NFORCF = 31,
+XTIMET(1) = 0.,
+XTIMET(2) = 54000.,
+XTIMEF(1) = 0.,
+XTIMEF(2) = 1800.,
+XTIMEF(3) = 3600.,
+XTIMEF(4) = 5400.,
+XTIMEF(5) = 7200.,
+XTIMEF(6) = 9000.,
+XTIMEF(7) = 10800.,
+XTIMEF(8) = 12600.,
+XTIMEF(9) = 14400.,
+XTIMEF(10) = 16200.,
+XTIMEF(11) = 18000.,
+XTIMEF(12) = 19800.,
+XTIMEF(13) = 21600.,
+XTIMEF(14) = 23400.,
+XTIMEF(15) = 25200.,
+XTIMEF(16) = 27000.,
+XTIMEF(17) = 28800.,
+XTIMEF(18) = 30600.,
+XTIMEF(19) = 32400.,
+XTIMEF(20) = 34200.,
+XTIMEF(21) = 36000.,
+XTIMEF(22) = 37800.,
+XTIMEF(23) = 39600.,
+XTIMEF(24) = 41400.,
+XTIMEF(25) = 43200.,
+XTIMEF(26) = 45000.,
+XTIMEF(27) = 46800.,
+XTIMEF(28) = 48600.,
+XTIMEF(29) = 50400.,
+XTIMEF(30) = 52200.,
+XTIMEF(31) = 54000.,
+XSFTH(1) = -30.,
+XSFTH(2) = -15.0,
+XSFTH(3) = 0.,
+XSFTH(4) = 15.0,
+XSFTH(5) = 30.0,
+XSFTH(6) = 45.0,
+XSFTH(7) = 60.0,
+XSFTH(8) = 75.0,
+XSFTH(9) = 90.0,
+XSFTH(10) = 100.0,
+XSFTH(11) = 110.0,
+XSFTH(12) = 120.0,
+XSFTH(13) = 130.0,
+XSFTH(14) = 140.0,
+XSFTH(15) = 140.0,
+XSFTH(16) = 140.0,
+XSFTH(17) = 132.0,
+XSFTH(18) = 124.0,
+XSFTH(19) = 116.0,
+XSFTH(20) = 108.0,
+XSFTH(21) = 100.0,
+XSFTH(22) = 78.0,
+XSFTH(23) = 56.0,
+XSFTH(24) = 34.0,
+XSFTH(25) = 12.0,
+XSFTH(26) = -10.0,
+XSFTH(27) = -10.0,
+XSFTH(28) = -10.0,
+XSFTH(29) = -10.0,
+XSFTH(30) = -10.0,
+XSFTH(31) = -10.0,
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+CSFTQ='kg/m2/s'  
+XSFTQ(1)=1.99936020473448506E-006  
+XSFTQ(2)=1.42454414587332047E-005  
+XSFTQ(3)=2.64915227127319252E-005  
+XSFTQ(4)=3.87376039667306440E-005  
+XSFTQ(5)=5.09836852207293695E-005  
+XSFTQ(6)=6.32297664747280883E-005  
+XSFTQ(7)=7.54758477287268071E-005  
+XSFTQ(8)=8.77219289827255259E-005  
+XSFTQ(9)=9.99680102367242446E-005  
+XSFTQ(10)=1.15962891874600132E-004  
+XSFTQ(11)=1.31957773512476006E-004  
+XSFTQ(12)=1.47952655150351879E-004  
+XSFTQ(13)=1.63947536788227780E-004  
+XSFTQ(14)=1.79942418426103654E-004  
+XSFTQ(15)=1.89939219449776085E-004  
+XSFTQ(16)=1.99936020473448489E-004  
+XSFTQ(17)=1.93538067818298140E-004  
+XSFTQ(18)=1.87140115163147790E-004  
+XSFTQ(19)=1.80742162507997441E-004  
+XSFTQ(20)=1.74344209852847091E-004  
+XSFTQ(21)=1.67946257197696742E-004  
+XSFTQ(22)=1.48752399232245693E-004  
+XSFTQ(23)=1.29558541266794618E-004  
+XSFTQ(24)=1.10364683301343569E-004  
+XSFTQ(25)=9.11708253358925209E-005  
+XSFTQ(26)=7.19769673704414589E-005
+XSFTQ(27)=5.39827255278310908E-005  
+XSFTQ(28)=3.59884836852207294E-005  
+XSFTQ(29)=1.79942418426103647E-005   
+XSFTQ(30)=0.0000000000000000      
+XSFTQ(31)=-1.79942418426103647E-005 
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+XSFCO2(1) = 0.,
+XSFCO2(2) = 0.,
+XSFCO2(3) = 0.,
+XSFCO2(4) = 0.,
+XSFCO2(5) = 0.,
+XSFCO2(6) = 0.,
+XSFCO2(7) = 0.,
+XSFCO2(8) = 0.,
+XSFCO2(9) = 0.,
+XSFCO2(10) = 0.,
+XSFCO2(11) = 0.,
+XSFCO2(12) = 0.,
+XSFCO2(13) = 0.,
+XSFCO2(14) = 0.,
+XSFCO2(15) = 0.,
+XSFCO2(16) = 0.,
+XSFCO2(17) = 0.,
+XSFCO2(18) = 0.,
+XSFCO2(19) = 0.,
+XSFCO2(20) = 0.,
+XSFCO2(21) = 0.,
+XSFCO2(22) = 0.,
+XSFCO2(23) = 0.,
+XSFCO2(24) = 0.,
+XSFCO2(25) = 0.,
+XSFCO2(26) = 0.,
+XSFCO2(27) = 0.,
+XSFCO2(28) = 0.,
+XSFCO2(29) = 0.,
+XSFCO2(30) = 0.,
+XSFCO2(31) = 0.,
+
+CUSTARTYPE = 'Z0   ',
+
+XZ0=0.035,
+
+XALB   = 0.,
+XEMIS  = 1.,
+XTSRAD(1) = 273.16,
+XTSRAD(2) = 273.16
+/
+
+&NAM_NEBn LSIGMAS=.TRUE., LSUBG_COND =.FALSE.  /

MY_RUN/INTEGRATION_CASES/HPC/ARMCU_LES/DEAR/PRE_IDEA1.nam

+&NAM_CONFIO LCDF4=T LLFIOUT=F LLFIREAD=F /
+&NAM_REAL_PGD /
+&NAM_DIMn_PRE NIMAX=64, NJMAX=64 /
+&NAM_CONF_PRE LCARTESIAN=.TRUE., NVERB=10,
+              CIDEAL='RSOU',  CZS='FLAT', LFORCING=.TRUE., NVERB=10,
+              LBOUSS=.FALSE., CEQNSYS='DUR', LPERTURB=.TRUE.,JPHEXT=1 /
+&NAM_PERT_PRE CPERT_KIND='WH' /
+&NAM_CONFn LUSERV=.TRUE. /
+&NAM_GRID_PRE XLAT0=35.762 /
+&NAM_GRIDH_PRE XDELTAX=100., XDELTAY=100., /
+&NAM_VER_GRID  LTHINSHELL=.TRUE., NKMAX=100,
+	       ZDZGRD=40., ZDZTOP=40.,
+               ZZMAX_STRGRD=1000. , ZSTRGRD=0., ZSTRTOP=0. /
+&NAM_LUNITn CINIFILE='ARM_3D_LF',
+            CINIFILEPGD='ARM_PGD' /
+&NAM_POST_PRE /
+&NAM_DYNn_PRE  /
+&NAM_LBCn_PRE CLBCX=2*"CYCL", CLBCY=2*"CYCL" / 
+&NAM_VPROF_PRE /
+&NAM_GRn_PRE CSURF='EXTE'/
+&NAM_CH_MNHCn_PRE /
+&NAM_BLANK /
+&NAM_PGD_SCHEMES CSEA='FLUX  ' /
+&NAM_COVER XUNIF_COVER(1)=1. /
+RSOU  
+1997 6 21 41400.
+'ZUVTHDMR'
+0.
+97000.
+299.
+15.20E-03
+2
+0.    10. 0.  
+5500. 10. 0.
+8
+50.   301.50  15.17E-03
+350.  302.50  14.98E-03
+650.  303.53  14.80E-03
+700.  303.70  14.70E-03
+1300. 307.13  13.50E-03
+2500. 314.00   3.00E-03
+5500. 343.20   3.00E-03
+ZFRC
+6
+1997 6 21 41400.
+0.
+97000.
+299.
+1.51913E-02
+4
+0.     10. 0. 0. 0. 0. -3.47222E-05 2.22222E-08 0. 0.
+1000.  10. 0. 0. 0. 0. -3.47222E-05 2.22222E-08 0. 0.
+3000.  10. 0. 0. 0. 0. 0.           0. 0. 0.
+5000.  10. 0. 0. 0. 0. 0.           0. 0. 0.
+1997 6 21 52200.
+0.
+97000.
+299.
+1.51913E-02
+4
+0.     10. 0. 0. 0. 0. 0. 5.55555E-09 0. 0.
+1000.  10. 0. 0. 0. 0. 0. 5.55555E-09 0. 0.
+3000.  10. 0. 0. 0. 0. 0.           0. 0. 0.
+5000.  10. 0. 0. 0. 0. 0.           0. 0. 0.
+1997 6 21 63000.
+0.
+97000.
+299.
+1.51913E-02
+4
+0.     10. 0. 0. 0. 0. 0. -1.11111E-08 0. 0.
+1000.  10. 0. 0. 0. 0. 0. -1.11111E-08 0. 0.
+3000.  10. 0. 0. 0. 0. 0.           0. 0. 0.
+5000.  10. 0. 0. 0. 0. 0.           0. 0. 0.
+1997 6 21 73800.
+0.
+97000.
+299.
+1.51913E-02
+4
+0.     10. 0. 0. 0. 0. -2.22222E-05 -2.77778E-08 0. 0.
+1000.  10. 0. 0. 0. 0. -2.22222E-05 -2.77778E-08 0. 0.
+3000.  10. 0. 0. 0. 0. 0.           0. 0. 0.
+5000.  10. 0. 0. 0. 0. 0.           0. 0. 0.
+1997 6 21 84600.
+0.
+97000.
+299.
+1.51913E-02
+4
+0.     10. 0. 0. 0. 0. -4.44444E-05 -4.44444E-08 0. 0.
+1000.  10. 0. 0. 0. 0. -4.44444E-05 -4.44444E-08 0. 0.
+3000.  10. 0. 0. 0. 0. 0.           0. 0. 0.
+5000.  10. 0. 0. 0. 0. 0.           0. 0. 0.
+1997 6 22 9000.
+0.
+97000.
+299.
+1.51913E-02
+4
+0.     10. 0. 0. 0. 0. -7.77777E-05 -9.11111E-08 0. 0.
+1000.  10. 0. 0. 0. 0. -7.77777E-05 -9.11111E-08 0. 0.
+3000.  10. 0. 0. 0. 0. 0.           0. 0. 0.
+5000.  10. 0. 0. 0. 0. 0.           0. 0. 0.         
+
+
+
+
+

MY_RUN/INTEGRATION_CASES/HPC/ARMCU_LES/DEAR/run_mesonh

+#!/bin/sh
+#SBATCH -J run_armles
+#SBATCH -N 2            # nodes number (=NBP)
+#SBATCH -n 128            # CPUs number (on all nodes) (=NBP*TPN)
+#SBATCH -o run_armles.eo%j   #
+#SBATCH -e run_armles.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 
+
+
+. ~rodierq/DEV_57/MNH-PHYEX070-b95d84d7/conf/profile_mesonh-LXifort-R8I4-MNH-V5-6-2-ECRAD140-MPIAUTO-O2
+#. ~rodierq/DEV_57/MNH-PHYEX070-b95d84d7/conf/profile_mesonh-LXifort-R8I4-MNH-V5-6-2-ECRAD140-MPIAUTO-O2-ARMSV
+export MPIRUN="Mpirun -np 128"
+
+set -x
+set -e
+
+
+ls -lrt
+
+rm -f ARM__.1.CEN4T.*
+time ${MPIRUN} MESONH${XYZ}
+mv OUTPUT_LISTING0  OUTPUT_LISTING0_run
+mv OUTPUT_LISTING1  OUTPUT_LISTING1_run
+ 
+ls -lrt 
+
+rm -f file_for_xtransfer pipe_name PRESSURE REMAP*
+
+cd ../PYTHON
+sbatch run_python

MY_RUN/INTEGRATION_CASES/HPC/ARMCU_LES/DEAR/run_prep_ideal

+#!/bin/sh
+#SBATCH -J ideal_armles
+#SBATCH -N 1            # nodes number (=NBP)
+#SBATCH -n 32            # CPUs number (on all nodes) (=NBP*TPN)
+#SBATCH -o ideal_armles.eo%j   #
+#SBATCH -e ideal_armles.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
+# 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 32"
+
+set -x
+set -e
+
+
+ls -lrt
+
+rm -f ARM_3D_LF.??? ARM_PGD.???
+time ${MPIRUN} PREP_IDEAL_CASE${XYZ}
+mv OUTPUT_LISTING1  OUTPUT_LISTING1_ideal
+ 
+ls -lrt 
+
+rm -f file_for_xtransfer pipe_name
+
+ls -lrt 
+
+sbatch run_mesonh
+ja

MY_RUN/INTEGRATION_CASES/HPC/ARMCU_LES/HM21/EXSEG1.nam

+&NAM_CONFIO LCDF4=T LLFIOUT=F LLFIREAD=F /
+&NAM_LUNITn  CINIFILE = "ARM_3D_LF", CINIFILEPGD = "ARM_PGD" /
+&NAM_CONFn  LUSERV = .TRUE./
+&NAM_DYNn XTSTEP = 2.,CPRESOPT="ZRESI"  /
+&NAM_ADVn CUVW_ADV_SCHEME="CEN4TH" ,CTEMP_SCHEME='RKC4'/
+&NAM_PARAMn  CTURB = "TKEL", CRAD = "NONE",
+ CCLOUD = "ICE3", CDCONV="NONE", CSCONV="NONE" /
+&NAM_PARAM_RADn XDTRAD = 120., 
+                XDTRAD_CLONLY = 1., 
+                CLW='RRTM',
+                CAER = 'SURF' /
+&NAM_PARAM_GROUNDN /
+&NAM_SEAFLUXn  /
+&NAM_ISBAn /
+&NAM_DIAG_SURFn /
+&NAM_LBCn CLBCX = 2*"CYCL", CLBCY = 2*"CYCL" /
+&NAM_TURBn XIMPL = 1., CTURBLEN = "HM21", CTURBDIM = "3DIM" , 
+           LTURB_DIAG = .TRUE., LRMC01=.TRUE.,
+           LTURB_FLX = .TRUE., 
+           LSIG_CONV=.FALSE. /
+&NAM_CH_MNHCn /
+&NAM_CONF  CCONF = "START", LFLAT = .TRUE., NMODEL = 1,NHALO=1
+                   NVERB = 4, CEQNSYS='DUR', lforcing=.true.,
+                   CEXP = "ARM__", CSEG = "CEN4T" /
+&NAM_DYN  XSEGLEN = 43200. , XASSELIN = 0.1, LCORIO = .true., 
+                    LNUMDIFU =T,
+                    XALKTOP = 0.01, XALZBOT = 3500. /
+&NAM_NESTING  NDAD(1) = 1 , NDAD(2) = 1 , NDTRATIO(1) = 1, NDTRATIO(2) = 2 ,
+              XWAY(1) = 1., XWAY(2) = 2. /
+&NAM_BACKUP  XBAK_TIME(1,1)= 43200. /
+&NAM_CONDSAMP LCONDSAMP=T, NCONDSAMP=3 /
+&NAM_LES    
+            LLES_MEAN=.TRUE., 
+            LLES_RESOLVED=.TRUE.,
+            LLES_SUBGRID=.TRUE.,
+            LLES_UPDRAFT=.FALSE.,
+            LLES_DOWNDRAFT=.FALSE.,CBL_HEIGHT_DEF='WTV',
+            CLES_NORM_TYPE='CONV',
+            XLES_TEMP_SAMPLING=300.,
+            LLES_NEB_MASK = .TRUE.,
+            LLES_CS_MASK = .TRUE.,
+            XLES_TEMP_MEAN_STEP=3600.,
+            XLES_TEMP_MEAN_START=3600.,
+            XLES_TEMP_MEAN_END=43200.  /
+&NAM_FRC LGEOST_UV_FRC=.TRUE., LTEND_THRV_FRC=.TRUE. /
+&NAM_CH_SOLVER /
+&NAM_PARAM_ICEn CSUBG_AUCV_RC='CLFR'   /
+&NAM_IDEAL_FLUX 
+NFORCT = 2,
+NFORCF = 31,
+XTIMET(1) = 0.,
+XTIMET(2) = 54000.,
+XTIMEF(1) = 0.,
+XTIMEF(2) = 1800.,
+XTIMEF(3) = 3600.,
+XTIMEF(4) = 5400.,
+XTIMEF(5) = 7200.,
+XTIMEF(6) = 9000.,
+XTIMEF(7) = 10800.,
+XTIMEF(8) = 12600.,
+XTIMEF(9) = 14400.,
+XTIMEF(10) = 16200.,
+XTIMEF(11) = 18000.,
+XTIMEF(12) = 19800.,
+XTIMEF(13) = 21600.,
+XTIMEF(14) = 23400.,
+XTIMEF(15) = 25200.,
+XTIMEF(16) = 27000.,
+XTIMEF(17) = 28800.,
+XTIMEF(18) = 30600.,
+XTIMEF(19) = 32400.,
+XTIMEF(20) = 34200.,
+XTIMEF(21) = 36000.,
+XTIMEF(22) = 37800.,
+XTIMEF(23) = 39600.,
+XTIMEF(24) = 41400.,
+XTIMEF(25) = 43200.,
+XTIMEF(26) = 45000.,
+XTIMEF(27) = 46800.,
+XTIMEF(28) = 48600.,
+XTIMEF(29) = 50400.,
+XTIMEF(30) = 52200.,
+XTIMEF(31) = 54000.,
+XSFTH(1) = -30.,
+XSFTH(2) = -15.0,
+XSFTH(3) = 0.,
+XSFTH(4) = 15.0,
+XSFTH(5) = 30.0,
+XSFTH(6) = 45.0,
+XSFTH(7) = 60.0,
+XSFTH(8) = 75.0,
+XSFTH(9) = 90.0,
+XSFTH(10) = 100.0,
+XSFTH(11) = 110.0,
+XSFTH(12) = 120.0,
+XSFTH(13) = 130.0,
+XSFTH(14) = 140.0,
+XSFTH(15) = 140.0,
+XSFTH(16) = 140.0,
+XSFTH(17) = 132.0,
+XSFTH(18) = 124.0,
+XSFTH(19) = 116.0,
+XSFTH(20) = 108.0,
+XSFTH(21) = 100.0,
+XSFTH(22) = 78.0,
+XSFTH(23) = 56.0,
+XSFTH(24) = 34.0,
+XSFTH(25) = 12.0,
+XSFTH(26) = -10.0,
+XSFTH(27) = -10.0,
+XSFTH(28) = -10.0,
+XSFTH(29) = -10.0,
+XSFTH(30) = -10.0,
+XSFTH(31) = -10.0,
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+CSFTQ='kg/m2/s'  
+XSFTQ(1)=1.99936020473448506E-006  
+XSFTQ(2)=1.42454414587332047E-005  
+XSFTQ(3)=2.64915227127319252E-005  
+XSFTQ(4)=3.87376039667306440E-005  
+XSFTQ(5)=5.09836852207293695E-005  
+XSFTQ(6)=6.32297664747280883E-005  
+XSFTQ(7)=7.54758477287268071E-005  
+XSFTQ(8)=8.77219289827255259E-005  
+XSFTQ(9)=9.99680102367242446E-005  
+XSFTQ(10)=1.15962891874600132E-004  
+XSFTQ(11)=1.31957773512476006E-004  
+XSFTQ(12)=1.47952655150351879E-004  
+XSFTQ(13)=1.63947536788227780E-004  
+XSFTQ(14)=1.79942418426103654E-004  
+XSFTQ(15)=1.89939219449776085E-004  
+XSFTQ(16)=1.99936020473448489E-004  
+XSFTQ(17)=1.93538067818298140E-004  
+XSFTQ(18)=1.87140115163147790E-004  
+XSFTQ(19)=1.80742162507997441E-004  
+XSFTQ(20)=1.74344209852847091E-004  
+XSFTQ(21)=1.67946257197696742E-004  
+XSFTQ(22)=1.48752399232245693E-004  
+XSFTQ(23)=1.29558541266794618E-004  
+XSFTQ(24)=1.10364683301343569E-004  
+XSFTQ(25)=9.11708253358925209E-005  
+XSFTQ(26)=7.19769673704414589E-005
+XSFTQ(27)=5.39827255278310908E-005  
+XSFTQ(28)=3.59884836852207294E-005  
+XSFTQ(29)=1.79942418426103647E-005   
+XSFTQ(30)=0.0000000000000000      
+XSFTQ(31)=-1.79942418426103647E-005 
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+XSFCO2(1) = 0.,
+XSFCO2(2) = 0.,
+XSFCO2(3) = 0.,
+XSFCO2(4) = 0.,
+XSFCO2(5) = 0.,
+XSFCO2(6) = 0.,
+XSFCO2(7) = 0.,
+XSFCO2(8) = 0.,
+XSFCO2(9) = 0.,
+XSFCO2(10) = 0.,
+XSFCO2(11) = 0.,
+XSFCO2(12) = 0.,
+XSFCO2(13) = 0.,
+XSFCO2(14) = 0.,
+XSFCO2(15) = 0.,
+XSFCO2(16) = 0.,
+XSFCO2(17) = 0.,
+XSFCO2(18) = 0.,
+XSFCO2(19) = 0.,
+XSFCO2(20) = 0.,
+XSFCO2(21) = 0.,
+XSFCO2(22) = 0.,
+XSFCO2(23) = 0.,
+XSFCO2(24) = 0.,
+XSFCO2(25) = 0.,
+XSFCO2(26) = 0.,
+XSFCO2(27) = 0.,
+XSFCO2(28) = 0.,
+XSFCO2(29) = 0.,
+XSFCO2(30) = 0.,
+XSFCO2(31) = 0.,
+
+CUSTARTYPE = 'Z0   ',
+
+XZ0=0.035,
+
+XALB   = 0.,
+XEMIS  = 1.,
+XTSRAD(1) = 273.16,
+XTSRAD(2) = 273.16
+/
+
+&NAM_NEBn LSIGMAS=.TRUE., LSUBG_COND =.FALSE.  /

MY_RUN/INTEGRATION_CASES/HPC/ARMCU_LES/HM21/PRE_IDEA1.nam

+&NAM_CONFIO LCDF4=T LLFIOUT=F LLFIREAD=F /
+&NAM_REAL_PGD /
+&NAM_DIMn_PRE NIMAX=64, NJMAX=64 /
+&NAM_CONF_PRE LCARTESIAN=.TRUE., NVERB=10,
+              CIDEAL='RSOU',  CZS='FLAT', LFORCING=.TRUE., NVERB=10,
+              LBOUSS=.FALSE., CEQNSYS='DUR', LPERTURB=.TRUE.,JPHEXT=1 /
+&NAM_PERT_PRE CPERT_KIND='WH' /
+&NAM_CONFn LUSERV=.TRUE. /
+&NAM_GRID_PRE XLAT0=35.762 /
+&NAM_GRIDH_PRE XDELTAX=500., XDELTAY=500., /
+&NAM_VER_GRID  LTHINSHELL=.TRUE., NKMAX=100,
+	       ZDZGRD=40., ZDZTOP=40.,
+               ZZMAX_STRGRD=1000. , ZSTRGRD=0., ZSTRTOP=0. /
+&NAM_LUNITn CINIFILE='ARM_3D_LF',
+            CINIFILEPGD='ARM_PGD' /
+&NAM_POST_PRE /
+&NAM_DYNn_PRE  /
+&NAM_LBCn_PRE CLBCX=2*"CYCL", CLBCY=2*"CYCL" / 
+&NAM_VPROF_PRE /
+&NAM_GRn_PRE CSURF='EXTE'/
+&NAM_CH_MNHCn_PRE /
+&NAM_BLANK /
+&NAM_PGD_SCHEMES CSEA='FLUX  ' /
+&NAM_COVER XUNIF_COVER(1)=1. /
+RSOU  
+1997 6 21 41400.
+'ZUVTHDMR'
+0.
+97000.
+299.
+15.20E-03
+2
+0.    10. 0.  
+5500. 10. 0.
+8
+50.   301.50  15.17E-03
+350.  302.50  14.98E-03
+650.  303.53  14.80E-03
+700.  303.70  14.70E-03
+1300. 307.13  13.50E-03
+2500. 314.00   3.00E-03
+5500. 343.20   3.00E-03
+ZFRC
+6
+1997 6 21 41400.
+0.
+97000.
+299.
+1.51913E-02
+4
+0.     10. 0. 0. 0. 0. -3.47222E-05 2.22222E-08 0. 0.
+1000.  10. 0. 0. 0. 0. -3.47222E-05 2.22222E-08 0. 0.
+3000.  10. 0. 0. 0. 0. 0.           0. 0. 0.
+5000.  10. 0. 0. 0. 0. 0.           0. 0. 0.
+1997 6 21 52200.
+0.
+97000.
+299.
+1.51913E-02
+4
+0.     10. 0. 0. 0. 0. 0. 5.55555E-09 0. 0.
+1000.  10. 0. 0. 0. 0. 0. 5.55555E-09 0. 0.
+3000.  10. 0. 0. 0. 0. 0.           0. 0. 0.
+5000.  10. 0. 0. 0. 0. 0.           0. 0. 0.
+1997 6 21 63000.
+0.
+97000.
+299.
+1.51913E-02
+4
+0.     10. 0. 0. 0. 0. 0. -1.11111E-08 0. 0.
+1000.  10. 0. 0. 0. 0. 0. -1.11111E-08 0. 0.
+3000.  10. 0. 0. 0. 0. 0.           0. 0. 0.
+5000.  10. 0. 0. 0. 0. 0.           0. 0. 0.
+1997 6 21 73800.
+0.
+97000.
+299.
+1.51913E-02
+4
+0.     10. 0. 0. 0. 0. -2.22222E-05 -2.77778E-08 0. 0.
+1000.  10. 0. 0. 0. 0. -2.22222E-05 -2.77778E-08 0. 0.
+3000.  10. 0. 0. 0. 0. 0.           0. 0. 0.
+5000.  10. 0. 0. 0. 0. 0.           0. 0. 0.
+1997 6 21 84600.
+0.
+97000.
+299.
+1.51913E-02
+4
+0.     10. 0. 0. 0. 0. -4.44444E-05 -4.44444E-08 0. 0.
+1000.  10. 0. 0. 0. 0. -4.44444E-05 -4.44444E-08 0. 0.
+3000.  10. 0. 0. 0. 0. 0.           0. 0. 0.
+5000.  10. 0. 0. 0. 0. 0.           0. 0. 0.
+1997 6 22 9000.
+0.
+97000.
+299.
+1.51913E-02
+4
+0.     10. 0. 0. 0. 0. -7.77777E-05 -9.11111E-08 0. 0.
+1000.  10. 0. 0. 0. 0. -7.77777E-05 -9.11111E-08 0. 0.
+3000.  10. 0. 0. 0. 0. 0.           0. 0. 0.
+5000.  10. 0. 0. 0. 0. 0.           0. 0. 0.         
+
+
+
+
+

MY_RUN/INTEGRATION_CASES/HPC/ARMCU_LES/HM21/run_mesonh

+#!/bin/sh
+#SBATCH -J run_armles
+#SBATCH -N 2            # nodes number (=NBP)
+#SBATCH -n 128            # CPUs number (on all nodes) (=NBP*TPN)
+#SBATCH -o run_armles.eo%j   #
+#SBATCH -e run_armles.eo%j   #
+#SBATCH -t 06: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 
+
+
+. ~rodierq/DEV_57/MNH-PHYEX070-b95d84d7/conf/profile_mesonh-LXifort-R8I4-MNH-V5-6-2-ECRAD140-MPIAUTO-O2
+
+export MPIRUN="Mpirun -np 128"
+
+set -x
+set -e
+
+
+ls -lrt
+
+rm -f ARM__.1.CEN4T.*
+time ${MPIRUN} MESONH${XYZ}
+mv OUTPUT_LISTING0  OUTPUT_LISTING0_run
+mv OUTPUT_LISTING1  OUTPUT_LISTING1_run
+ 
+ls -lrt 
+
+rm -f file_for_xtransfer pipe_name PRESSURE REMAP*
+
+cd ../PYTHON
+sbatch run_python

MY_RUN/INTEGRATION_CASES/HPC/ARMCU_LES/HM21/run_prep_ideal

+#!/bin/sh
+#SBATCH -J ideal_armles
+#SBATCH -N 1            # nodes number (=NBP)
+#SBATCH -n 32            # CPUs number (on all nodes) (=NBP*TPN)
+#SBATCH -o ideal_armles.eo%j   #
+#SBATCH -e ideal_armles.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
+# 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 32"
+
+set -x
+set -e
+
+
+ls -lrt
+
+rm -f ARM_3D_LF.??? ARM_PGD.???
+time ${MPIRUN} PREP_IDEAL_CASE${XYZ}
+mv OUTPUT_LISTING1  OUTPUT_LISTING1_ideal
+ 
+ls -lrt 
+
+rm -f file_for_xtransfer pipe_name
+
+ls -lrt 
+
+sbatch run_mesonh
+ja

MY_RUN/INTEGRATION_CASES/HPC/ARMCU_LES/PYTHON/plot_ARMCU_LES_DEAR.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
+import os
+
+os.system('rm -f tempgraph*')
+#
+#  User's parameter / Namelist
+#
+path="../DEAR/"
+
+LnameFiles = ['ARM__.1.CEN4T.000.nc' ]
+
+LG_CART = '/LES_budgets/Mean/Cartesian/Not_time_averaged/Not_normalized/cart/'
+LG_NEB = '/LES_budgets/Mean/Cartesian/Not_time_averaged/Not_normalized/neb/'
+LG_CLEAR = '/LES_budgets/Mean/Cartesian/Not_time_averaged/Not_normalized/clear/'
+LG_CS1 = '/LES_budgets/Mean/Cartesian/Not_time_averaged/Not_normalized/cs1/'
+LG_CS2 = '/LES_budgets/Mean/Cartesian/Not_time_averaged/Not_normalized/cs2/'
+LG_CS3 = '/LES_budgets/Mean/Cartesian/Not_time_averaged/Not_normalized/cs3/'
+LG_SBGCART = '/LES_budgets/Subgrid/Cartesian/Not_time_averaged/Not_normalized/cart/'
+LG_SBGNEB = '/LES_budgets/Subgrid/Cartesian/Not_time_averaged/Not_normalized/neb/'
+LG_SBGCLEAR = '/LES_budgets/Subgrid/Cartesian/Not_time_averaged/Not_normalized/clear/'
+LG_SBGCS1 = '/LES_budgets/Subgrid/Cartesian/Not_time_averaged/Not_normalized/cs1/'
+LG_SBGCS2 = '/LES_budgets/Subgrid/Cartesian/Not_time_averaged/Not_normalized/cs2/'
+LG_SBGCS3 = '/LES_budgets/Subgrid/Cartesian/Not_time_averaged/Not_normalized/cs3/'
+
+Dvar_input = {'f1':['time_les','level_les',
+(LG_CART, 'MEAN_RC'),  (LG_NEB, 'MEAN_RC'),  (LG_CLEAR, 'MEAN_RC'),  (LG_CS1, 'MEAN_RC'),  (LG_CS2, 'MEAN_RC'),  (LG_CS3, 'MEAN_RC'),
+(LG_CART, 'MEAN_RR'),  (LG_NEB, 'MEAN_RR'),  (LG_CLEAR, 'MEAN_RR'),  (LG_CS1, 'MEAN_RR'),  (LG_CS2, 'MEAN_RR'),  (LG_CS3, 'MEAN_RR'),
+(LG_CART, 'MEAN_U'),   (LG_NEB, 'MEAN_U'),   (LG_CLEAR, 'MEAN_U'),   (LG_CS1, 'MEAN_U'),   (LG_CS2, 'MEAN_U'),   (LG_CS3, 'MEAN_U'),
+(LG_CART, 'MEAN_V'),   (LG_NEB, 'MEAN_V'),   (LG_CLEAR, 'MEAN_V'),   (LG_CS1, 'MEAN_V'),   (LG_CS2, 'MEAN_V'),   (LG_CS3, 'MEAN_V'),
+(LG_CART, 'MEAN_W'),   (LG_NEB, 'MEAN_W'),   (LG_CLEAR, 'MEAN_W'),   (LG_CS1, 'MEAN_W'),   (LG_CS2, 'MEAN_W'),   (LG_CS3, 'MEAN_W'),
+(LG_CART, 'MEAN_THL'), (LG_NEB, 'MEAN_THL'), (LG_CLEAR, 'MEAN_THL'), (LG_CS1, 'MEAN_THL'), (LG_CS2, 'MEAN_THL'), (LG_CS3, 'MEAN_THL'),
+(LG_CART, 'MEAN_RT'),  (LG_NEB, 'MEAN_RT'),  (LG_CLEAR, 'MEAN_RT'),  (LG_CS1, 'MEAN_RT'),  (LG_CS2, 'MEAN_RT'),  (LG_CS3, 'MEAN_RT'),
+(LG_CART, 'MEAN_MF'),  (LG_NEB, 'MEAN_MF'),  (LG_CLEAR, 'MEAN_MF'),  (LG_CS1, 'MEAN_MF'),  (LG_CS2, 'MEAN_MF'),  (LG_CS3, 'MEAN_MF'),
+(LG_CART, 'MEAN_SV'),  (LG_NEB, 'MEAN_SV'),  (LG_CLEAR, 'MEAN_SV'),  (LG_CS1, 'MEAN_SV'),  (LG_CS2, 'MEAN_SV'),  (LG_CS3, 'MEAN_SV'),
+(LG_SBGCART, 'SBG_WTHL'), (LG_SBGNEB, 'SBG_WTHL'), (LG_SBGCLEAR, 'SBG_WTHL'), (LG_SBGCS1, 'SBG_WTHL'), (LG_SBGCS2, 'SBG_WTHL'), (LG_SBGCS3, 'SBG_WTHL'),
+(LG_SBGCART, 'SBG_WRT'),  (LG_SBGNEB, 'SBG_WRT'),  (LG_SBGCLEAR, 'SBG_WRT'),  (LG_SBGCS1, 'SBG_WRT'),  (LG_SBGCS2, 'SBG_WRT'),  (LG_SBGCS3, 'SBG_WRT')]
+}
+
+#  Read the variables in the files
+Dvar = {}
+Dvar = read_netcdf(LnameFiles, Dvar_input, path=path, removeHALO=False)
+
+################################################################
+#########          PANEL  1
+###############################################################
+Panel = PanelPlot(2,3, [25,14],'MEAN_RC', titlepad=25, minmaxpad=1.04, timepad=-0.07, colorbarpad=0.03, labelcolorbarpad = 13, colorbaraspect=40)
+var='MEAN_RC'
+Lplot = [Dvar['f1'][(LG_CART,var)],Dvar['f1'][(LG_NEB,var)],Dvar['f1'][(LG_CLEAR,var)],
+         Dvar['f1'][(LG_CS1,var)],Dvar['f1'][(LG_CS2,var)],Dvar['f1'][(LG_CS3,var)]]
+
+LaxeX = [Dvar['f1']['time_les']/3600.]*len(Lplot)
+LaxeZ = [Dvar['f1']['level_les']]*len(Lplot)
+Ltitle = ['cart','neb','clear','cs1','cs2','cs3']
+Lcbarlabel = ['g/kg']*len(Lplot)
+Lxlab = ['time (h)']*len(Lplot)
+Lylab = ['altitude (m)']*len(Lplot)
+Lylim = [(0,3000)]*len(Lplot)
+Lminval = [-0.009]*len(Lplot)
+Lmaxval = [0.2]*len(Lplot)
+Lstep = [0.01]*len(Lplot)
+Lstepticks = Lstep
+Lfacconv = [1000]*len(Lplot)
+Lcolormap = ['gist_rainbow_r']*len(Lplot)
+LaddWhite = [True]*len(Lplot)
+
+fig = Panel.psectionV(Lxx=LaxeX, Lzz=LaxeZ, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, 
+                                Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv, 
+                                LaddWhite_cm=LaddWhite, Lylim=Lylim)
+Panel.save_graph(1,fig)
+################################################################
+#########          PANEL 2
+###############################################################
+Panel = PanelPlot(2,3, [25,14],'MEAN_RR', titlepad=25, minmaxpad=1.04, timepad=-0.07, colorbarpad=0.03, labelcolorbarpad = 13, colorbaraspect=40)
+var='MEAN_RR'
+Lplot = [Dvar['f1'][(LG_CART,var)],Dvar['f1'][(LG_NEB,var)],Dvar['f1'][(LG_CLEAR,var)],
+         Dvar['f1'][(LG_CS1,var)],Dvar['f1'][(LG_CS2,var)],Dvar['f1'][(LG_CS3,var)]]
+
+Lcbarlabel = ['g/kg']*len(Lplot)
+Lminval = [-0.0009]*len(Lplot)
+Lmaxval = [0.02]*len(Lplot)
+Lstep = [0.001]*len(Lplot)
+Lstepticks = Lstep
+Lfacconv = [1000]*len(Lplot)
+LaddWhite = [True]*len(Lplot)
+
+fig = Panel.psectionV(Lxx=LaxeX, Lzz=LaxeZ, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, 
+                                Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv, 
+                                LaddWhite_cm=LaddWhite, Lylim=Lylim)
+Panel.save_graph(2,fig)
+################################################################
+#########          PANEL 3
+###############################################################
+Panel = PanelPlot(2,3, [25,14],'MEAN_U', titlepad=25, minmaxpad=1.04, timepad=-0.07, colorbarpad=0.03, labelcolorbarpad = 13, colorbaraspect=40)
+var='MEAN_U'
+Lplot = [Dvar['f1'][(LG_CART,var)],Dvar['f1'][(LG_NEB,var)],Dvar['f1'][(LG_CLEAR,var)],
+         Dvar['f1'][(LG_CS1,var)],Dvar['f1'][(LG_CS2,var)],Dvar['f1'][(LG_CS3,var)]]
+
+Lcbarlabel = ['m/s']*len(Lplot)
+Lminval = [3]*len(Lplot)
+Lmaxval = [11]*len(Lplot)
+Lstep = [0.5]*len(Lplot)
+Lstepticks = Lstep
+Lfacconv = [1]*len(Lplot)
+LaddWhite = [False]*len(Lplot)
+
+fig = Panel.psectionV(Lxx=LaxeX, Lzz=LaxeZ, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, 
+                                Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv, 
+                                LaddWhite_cm=LaddWhite, Lylim=Lylim)
+Panel.save_graph(3,fig)
+################################################################
+#########          PANEL 4
+###############################################################
+Panel = PanelPlot(2,3, [25,14],'MEAN_V', titlepad=25, minmaxpad=1.04, timepad=-0.07, colorbarpad=0.03, labelcolorbarpad = 13, colorbaraspect=40)
+var='MEAN_V'
+Lplot = [Dvar['f1'][(LG_CART,var)],Dvar['f1'][(LG_NEB,var)],Dvar['f1'][(LG_CLEAR,var)],
+         Dvar['f1'][(LG_CS1,var)],Dvar['f1'][(LG_CS2,var)],Dvar['f1'][(LG_CS3,var)]]
+
+Lcbarlabel = ['m/s']*len(Lplot)
+Lminval = [0]*len(Lplot)
+Lmaxval = [3.2]*len(Lplot)
+Lstep = [0.2]*len(Lplot)
+Lstepticks = Lstep
+Lfacconv = [1]*len(Lplot)
+LaddWhite = [False]*len(Lplot)
+
+fig = Panel.psectionV(Lxx=LaxeX, Lzz=LaxeZ, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, 
+                                Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv, 
+                                LaddWhite_cm=LaddWhite, Lylim=Lylim)
+Panel.save_graph(4,fig)
+################################################################
+#########          PANEL 5
+###############################################################
+Panel = PanelPlot(2,3, [25,14],'MEAN_W', titlepad=25, minmaxpad=1.04, timepad=-0.07, colorbarpad=0.03, labelcolorbarpad = 13, colorbaraspect=40)
+var='MEAN_W'
+Lplot = [Dvar['f1'][(LG_CART,var)],Dvar['f1'][(LG_NEB,var)],Dvar['f1'][(LG_CLEAR,var)],
+         Dvar['f1'][(LG_CS1,var)],Dvar['f1'][(LG_CS2,var)],Dvar['f1'][(LG_CS3,var)]]
+
+Lcbarlabel = ['m/s']*len(Lplot)
+Lminval = [0., 0., -0.18, 0., -1.25, -1.2]
+Lmaxval = [4.25, 4.25, 0., 5.25, 2.75, 0.8]
+Lstep = [0.25, 0.25, 0.01, 0.25, 0.25, 0.1]
+Lstepticks = Lstep
+Lfacconv = [1]*len(Lplot)
+LaddWhite = [False]*len(Lplot)
+
+fig = Panel.psectionV(Lxx=LaxeX, Lzz=LaxeZ, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, 
+                                Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv, 
+                                LaddWhite_cm=LaddWhite, Lylim=Lylim)
+Panel.save_graph(5,fig)
+################################################################
+#########          PANEL 6
+###############################################################
+Panel = PanelPlot(2,3, [25,14],'MEAN_THL', titlepad=25, minmaxpad=1.04, timepad=-0.07, colorbarpad=0.03, labelcolorbarpad = 13, colorbaraspect=40)
+var='MEAN_THL'
+Lplot = [Dvar['f1'][(LG_CART,var)],Dvar['f1'][(LG_NEB,var)],Dvar['f1'][(LG_CLEAR,var)],
+         Dvar['f1'][(LG_CS1,var)],Dvar['f1'][(LG_CS2,var)],Dvar['f1'][(LG_CS3,var)]]
+
+Lcbarlabel = ['K']*len(Lplot)
+Lminval = [300]*len(Lplot)
+Lmaxval = [320]*len(Lplot)
+Lstep = [0.5]*len(Lplot)
+Lstepticks = Lstep
+Lfacconv = [1]*len(Lplot)
+LaddWhite = [False]*len(Lplot)
+
+fig = Panel.psectionV(Lxx=LaxeX, Lzz=LaxeZ, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, 
+                                Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv, 
+                                LaddWhite_cm=LaddWhite, Lylim=Lylim)
+Panel.save_graph(6,fig)
+################################################################
+#########          PANEL 7
+###############################################################
+Panel = PanelPlot(2,3, [25,14],'MEAN_SV 1', titlepad=25, minmaxpad=1.04, timepad=-0.07, colorbarpad=0.03, labelcolorbarpad = 13, colorbaraspect=40)
+var='MEAN_SV'
+Lplot = [Dvar['f1'][(LG_CART,var)][:,:,0],Dvar['f1'][(LG_NEB,var)][:,:,0],Dvar['f1'][(LG_CLEAR,var)][:,:,0],
+         Dvar['f1'][(LG_CS1,var)][:,:,0] ,Dvar['f1'][(LG_CS2,var)][:,:,0],Dvar['f1'][(LG_CS3,var)][:,:,0]]
+
+Lcbarlabel = ['kg/kg']*len(Lplot)
+Lminval = [2.5,9, 2.5, 2.5, 2, 2 ]
+Lmaxval = [50, 35, 50, 50, 34, 34]
+Lstep = [2.5, 1, 2.5, 2.5, 2, 2]
+Lstepticks = Lstep
+Lfacconv = [1]*len(Lplot)
+LaddWhite = [False]*len(Lplot)
+
+fig = Panel.psectionV(Lxx=LaxeX, Lzz=LaxeZ, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, 
+                                Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv, 
+                                LaddWhite_cm=LaddWhite, Lylim=Lylim)
+Panel.save_graph(7,fig)
+################################################################
+#########          PANEL 8
+###############################################################
+Panel = PanelPlot(2,3, [25,14],'MEAN_SV 2', titlepad=25, minmaxpad=1.04, timepad=-0.07, colorbarpad=0.03, labelcolorbarpad = 13, colorbaraspect=40)
+var='MEAN_SV'
+Lplot = [Dvar['f1'][(LG_CART,var)][:,:,1],Dvar['f1'][(LG_NEB,var)][:,:,1],Dvar['f1'][(LG_CLEAR,var)][:,:,1],
+         Dvar['f1'][(LG_CS1,var)][:,:,1] ,Dvar['f1'][(LG_CS2,var)][:,:,1],Dvar['f1'][(LG_CS3,var)][:,:,1]]
+
+Lminval = [0.1, 0.3, 0.1, 0.1, 0.1, 0.1 ]
+Lmaxval = [2.1, 1.6, 2.1, 1.7, 2.1, 1.75]
+Lstep = [0.1, 0.05, 0.1, 0.1, 0.1, 0.05]
+Lstepticks = Lstep
+Lcolormap = ['gist_rainbow_r']*len(Lplot)
+fig = Panel.psectionV(Lxx=LaxeX, Lzz=LaxeZ, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, 
+                                Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv, 
+                                LaddWhite_cm=LaddWhite, Lylim=Lylim)
+Panel.save_graph(8,fig)
+################################################################
+#########          PANEL 9
+###############################################################
+Panel = PanelPlot(2,3, [25,14],'MEAN_SV 3', titlepad=25, minmaxpad=1.04, timepad=-0.07, colorbarpad=0.03, labelcolorbarpad = 13, colorbaraspect=40)
+var='MEAN_SV'
+Lplot = [Dvar['f1'][(LG_CART,var)][:,:,2],Dvar['f1'][(LG_NEB,var)][:,:,2],Dvar['f1'][(LG_CLEAR,var)][:,:,2],
+         Dvar['f1'][(LG_CS1,var)][:,:,2] ,Dvar['f1'][(LG_CS2,var)][:,:,2],Dvar['f1'][(LG_CS3,var)][:,:,2]]
+
+Lminval = [0.1,0.05, 0.1, 0.05, 0.1, 0.1 ]
+Lmaxval = [2.3, 1.05, 2.3, 1.05, 3.5, 3.1]
+Lstep = [0.1, 0.05, 0.1, 0.05, 0.1, 0.1]
+Lstepticks = Lstep
+fig = Panel.psectionV(Lxx=LaxeX, Lzz=LaxeZ, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, 
+                                Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv, 
+                                LaddWhite_cm=LaddWhite, Lylim=Lylim)
+Panel.save_graph(9,fig)
+################################################################
+#########          PANEL 10
+###############################################################
+Panel = PanelPlot(2,3, [25,14],'MEAN_RT', titlepad=25, minmaxpad=1.04, timepad=-0.07, colorbarpad=0.03, labelcolorbarpad = 13, colorbaraspect=40)
+var='MEAN_RT'
+Lplot = [Dvar['f1'][(LG_CART,var)],Dvar['f1'][(LG_NEB,var)],Dvar['f1'][(LG_CLEAR,var)],
+         Dvar['f1'][(LG_CS1,var)],Dvar['f1'][(LG_CS2,var)],Dvar['f1'][(LG_CS3,var)]]
+
+Lcbarlabel = ['kg/kg']*len(Lplot)
+Lminval = [0.01]*len(Lplot)
+Lmaxval = [0.019]*len(Lplot)
+Lstep = [0.001]*len(Lplot)
+Lstepticks = Lstep
+Lfacconv = [1]*len(Lplot)
+LaddWhite = [False]*len(Lplot)
+
+fig = Panel.psectionV(Lxx=LaxeX, Lzz=LaxeZ, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, 
+                                Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv, 
+                                LaddWhite_cm=LaddWhite, Lylim=Lylim)
+Panel.save_graph(10,fig)
+################################################################
+#########          PANEL 11
+###############################################################
+Panel = PanelPlot(2,3, [25,14],'MEAN_MF', titlepad=25, minmaxpad=1.04, timepad=-0.07, colorbarpad=0.03, labelcolorbarpad = 13, colorbaraspect=40)
+var='MEAN_MF'
+Lplot = [Dvar['f1'][(LG_CART,var)],Dvar['f1'][(LG_NEB,var)],Dvar['f1'][(LG_CLEAR,var)],
+         Dvar['f1'][(LG_CS1,var)],Dvar['f1'][(LG_CS2,var)],Dvar['f1'][(LG_CS3,var)]]
+
+Lcbarlabel = ['kg/s']*len(Lplot)
+Lminval = [-0.5E-3, 0.2, -0.18, 0.2, -1.2, -1.2]
+Lmaxval = [0.2E-3, 4.2, 0.,4.2, 2.5, 1.1 ]
+Lstep = [0.05E-3, 0.2, 0.01, 0.2, 0.1, 0.1]
+Lstepticks = Lstep
+Lfacconv = [1]*len(Lplot)
+LaddWhite = [False]*len(Lplot)
+
+fig = Panel.psectionV(Lxx=LaxeX, Lzz=LaxeZ, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, 
+                                Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv, 
+                                LaddWhite_cm=LaddWhite, Lylim=Lylim)
+Panel.save_graph(11,fig)
+################################################################
+#########          PANEL 12
+###############################################################
+Panel = PanelPlot(2,3, [25,14],'SBG_WTHL', titlepad=25, minmaxpad=1.04, timepad=-0.07, colorbarpad=0.03, labelcolorbarpad = 13, colorbaraspect=40)
+var='SBG_WTHL'
+Lplot = [Dvar['f1'][(LG_SBGCART,var)],Dvar['f1'][(LG_SBGNEB,var)],Dvar['f1'][(LG_SBGCLEAR,var)],
+         Dvar['f1'][(LG_SBGCS1,var)],Dvar['f1'][(LG_SBGCS2,var)],Dvar['f1'][(LG_SBGCS3,var)]]
+
+Lcbarlabel = ['g/kg']*len(Lplot)
+Lminval = [-0.035,-0.8, -0.035, -0.8, -0.4, -0.04 ]
+Lmaxval = [0.105, 0., 0.105, 0., 0.05, 0.085]
+Lstep = [0.005, 0.05, 0.005, 0.05, 0.005, 0.005]
+Lstepticks = Lstep
+Lfacconv = [1]*len(Lplot)
+LaddWhite = [False]*len(Lplot)
+
+fig = Panel.psectionV(Lxx=LaxeX, Lzz=LaxeZ, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, 
+                                Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv, 
+                                LaddWhite_cm=LaddWhite, Lylim=Lylim)
+Panel.save_graph(12,fig)
+################################################################
+#########          PANEL 13
+###############################################################
+Panel = PanelPlot(2,3, [25,14],'SBG_WRT', titlepad=25, minmaxpad=1.04, timepad=-0.07, colorbarpad=0.03, labelcolorbarpad = 13, colorbaraspect=40)
+var='SBG_WRT'
+Lplot = [Dvar['f1'][(LG_SBGCART,var)],Dvar['f1'][(LG_SBGNEB,var)],Dvar['f1'][(LG_SBGCLEAR,var)],
+         Dvar['f1'][(LG_SBGCS1,var)],Dvar['f1'][(LG_SBGCS2,var)],Dvar['f1'][(LG_SBGCS3,var)]]
+
+Lcbarlabel = ['m/s kg/kg']*len(Lplot)
+Lminval = [0.]*len(Lplot)
+Lmaxval = [0.15E-3, 0.16E-2, 0.15E-3, 0.16E-2, 0.85E-3, 0.105E-3]
+Lstep = [0.1E-4, 0.05E-3, 0.1E-4,0.05E-3, 0.05E-3, 0.05E-4 ]
+Lstepticks = Lstep
+Lfacconv = [1]*len(Lplot)
+LaddWhite = [True]*len(Lplot)
+
+fig = Panel.psectionV(Lxx=LaxeX, Lzz=LaxeZ, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, 
+                                Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv, 
+                                LaddWhite_cm=LaddWhite, Lylim=Lylim)
+Panel.save_graph(13,fig)

MY_RUN/INTEGRATION_CASES/HPC/ARMCU_LES/PYTHON/plot_ARMCU_LES_HM21.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
+import os
+
+os.system('rm -f tempgraph*')
+#
+#  User's parameter / Namelist
+#
+path="../HM21/"
+
+LnameFiles = ['ARM__.1.CEN4T.000.nc' ]
+
+LG_CART = '/LES_budgets/Mean/Cartesian/Not_time_averaged/Not_normalized/cart/'
+LG_NEB = '/LES_budgets/Mean/Cartesian/Not_time_averaged/Not_normalized/neb/'
+LG_CLEAR = '/LES_budgets/Mean/Cartesian/Not_time_averaged/Not_normalized/clear/'
+LG_CS1 = '/LES_budgets/Mean/Cartesian/Not_time_averaged/Not_normalized/cs1/'
+LG_CS2 = '/LES_budgets/Mean/Cartesian/Not_time_averaged/Not_normalized/cs2/'
+LG_CS3 = '/LES_budgets/Mean/Cartesian/Not_time_averaged/Not_normalized/cs3/'
+LG_SBGCART = '/LES_budgets/Subgrid/Cartesian/Not_time_averaged/Not_normalized/cart/'
+LG_SBGNEB = '/LES_budgets/Subgrid/Cartesian/Not_time_averaged/Not_normalized/neb/'
+LG_SBGCLEAR = '/LES_budgets/Subgrid/Cartesian/Not_time_averaged/Not_normalized/clear/'
+LG_SBGCS1 = '/LES_budgets/Subgrid/Cartesian/Not_time_averaged/Not_normalized/cs1/'
+LG_SBGCS2 = '/LES_budgets/Subgrid/Cartesian/Not_time_averaged/Not_normalized/cs2/'
+LG_SBGCS3 = '/LES_budgets/Subgrid/Cartesian/Not_time_averaged/Not_normalized/cs3/'
+
+Dvar_input = {'f1':['time_les','level_les',
+(LG_CART, 'MEAN_RC'),  (LG_NEB, 'MEAN_RC'),  (LG_CLEAR, 'MEAN_RC'),  (LG_CS1, 'MEAN_RC'),  (LG_CS2, 'MEAN_RC'),  (LG_CS3, 'MEAN_RC'),
+(LG_CART, 'MEAN_RR'),  (LG_NEB, 'MEAN_RR'),  (LG_CLEAR, 'MEAN_RR'),  (LG_CS1, 'MEAN_RR'),  (LG_CS2, 'MEAN_RR'),  (LG_CS3, 'MEAN_RR'),
+(LG_CART, 'MEAN_U'),   (LG_NEB, 'MEAN_U'),   (LG_CLEAR, 'MEAN_U'),   (LG_CS1, 'MEAN_U'),   (LG_CS2, 'MEAN_U'),   (LG_CS3, 'MEAN_U'),
+(LG_CART, 'MEAN_V'),   (LG_NEB, 'MEAN_V'),   (LG_CLEAR, 'MEAN_V'),   (LG_CS1, 'MEAN_V'),   (LG_CS2, 'MEAN_V'),   (LG_CS3, 'MEAN_V'),
+(LG_CART, 'MEAN_W'),   (LG_NEB, 'MEAN_W'),   (LG_CLEAR, 'MEAN_W'),   (LG_CS1, 'MEAN_W'),   (LG_CS2, 'MEAN_W'),   (LG_CS3, 'MEAN_W'),
+(LG_CART, 'MEAN_THL'), (LG_NEB, 'MEAN_THL'), (LG_CLEAR, 'MEAN_THL'), (LG_CS1, 'MEAN_THL'), (LG_CS2, 'MEAN_THL'), (LG_CS3, 'MEAN_THL'),
+(LG_CART, 'MEAN_RT'),  (LG_NEB, 'MEAN_RT'),  (LG_CLEAR, 'MEAN_RT'),  (LG_CS1, 'MEAN_RT'),  (LG_CS2, 'MEAN_RT'),  (LG_CS3, 'MEAN_RT'),
+(LG_CART, 'MEAN_MF'),  (LG_NEB, 'MEAN_MF'),  (LG_CLEAR, 'MEAN_MF'),  (LG_CS1, 'MEAN_MF'),  (LG_CS2, 'MEAN_MF'),  (LG_CS3, 'MEAN_MF'),
+(LG_CART, 'MEAN_SV'),  (LG_NEB, 'MEAN_SV'),  (LG_CLEAR, 'MEAN_SV'),  (LG_CS1, 'MEAN_SV'),  (LG_CS2, 'MEAN_SV'),  (LG_CS3, 'MEAN_SV'),
+(LG_SBGCART, 'SBG_WTHL'), (LG_SBGNEB, 'SBG_WTHL'), (LG_SBGCLEAR, 'SBG_WTHL'), (LG_SBGCS1, 'SBG_WTHL'), (LG_SBGCS2, 'SBG_WTHL'), (LG_SBGCS3, 'SBG_WTHL'),
+(LG_SBGCART, 'SBG_WRT'),  (LG_SBGNEB, 'SBG_WRT'),  (LG_SBGCLEAR, 'SBG_WRT'),  (LG_SBGCS1, 'SBG_WRT'),  (LG_SBGCS2, 'SBG_WRT'),  (LG_SBGCS3, 'SBG_WRT')]
+}
+
+#  Read the variables in the files
+Dvar = {}
+Dvar = read_netcdf(LnameFiles, Dvar_input, path=path, removeHALO=False)
+
+################################################################
+#########          PANEL  1
+###############################################################
+Panel = PanelPlot(2,3, [25,14],'MEAN_RC', titlepad=25, minmaxpad=1.04, timepad=-0.07, colorbarpad=0.03, labelcolorbarpad = 13, colorbaraspect=40)
+var='MEAN_RC'
+Lplot = [Dvar['f1'][(LG_CART,var)],Dvar['f1'][(LG_NEB,var)],Dvar['f1'][(LG_CLEAR,var)],
+         Dvar['f1'][(LG_CS1,var)],Dvar['f1'][(LG_CS2,var)],Dvar['f1'][(LG_CS3,var)]]
+
+LaxeX = [Dvar['f1']['time_les']/3600.]*len(Lplot)
+LaxeZ = [Dvar['f1']['level_les']]*len(Lplot)
+Ltitle = ['cart','neb','clear','cs1','cs2','cs3']
+Lcbarlabel = ['g/kg']*len(Lplot)
+Lxlab = ['time (h)']*len(Lplot)
+Lylab = ['altitude (m)']*len(Lplot)
+Lylim = [(0,3000)]*len(Lplot)
+Lminval = [-0.009]*len(Lplot)
+Lmaxval = [0.2]*len(Lplot)
+Lstep = [0.01]*len(Lplot)
+Lstepticks = Lstep
+Lfacconv = [1000]*len(Lplot)
+Lcolormap = ['gist_rainbow_r']*len(Lplot)
+LaddWhite = [True]*len(Lplot)
+
+fig = Panel.psectionV(Lxx=LaxeX, Lzz=LaxeZ, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, 
+                                Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv, 
+                                LaddWhite_cm=LaddWhite, Lylim=Lylim)
+Panel.save_graph(1,fig)
+################################################################
+#########          PANEL 2
+###############################################################
+Panel = PanelPlot(2,3, [25,14],'MEAN_RR', titlepad=25, minmaxpad=1.04, timepad=-0.07, colorbarpad=0.03, labelcolorbarpad = 13, colorbaraspect=40)
+var='MEAN_RR'
+Lplot = [Dvar['f1'][(LG_CART,var)],Dvar['f1'][(LG_NEB,var)],Dvar['f1'][(LG_CLEAR,var)],
+         Dvar['f1'][(LG_CS1,var)],Dvar['f1'][(LG_CS2,var)],Dvar['f1'][(LG_CS3,var)]]
+
+Lcbarlabel = ['g/kg']*len(Lplot)
+Lminval = [-0.0009]*len(Lplot)
+Lmaxval = [0.02]*len(Lplot)
+Lstep = [0.001]*len(Lplot)
+Lstepticks = Lstep
+Lfacconv = [1000]*len(Lplot)
+LaddWhite = [True]*len(Lplot)
+
+fig = Panel.psectionV(Lxx=LaxeX, Lzz=LaxeZ, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, 
+                                Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv, 
+                                LaddWhite_cm=LaddWhite, Lylim=Lylim)
+Panel.save_graph(2,fig)
+################################################################
+#########          PANEL 3
+###############################################################
+Panel = PanelPlot(2,3, [25,14],'MEAN_U', titlepad=25, minmaxpad=1.04, timepad=-0.07, colorbarpad=0.03, labelcolorbarpad = 13, colorbaraspect=40)
+var='MEAN_U'
+Lplot = [Dvar['f1'][(LG_CART,var)],Dvar['f1'][(LG_NEB,var)],Dvar['f1'][(LG_CLEAR,var)],
+         Dvar['f1'][(LG_CS1,var)],Dvar['f1'][(LG_CS2,var)],Dvar['f1'][(LG_CS3,var)]]
+
+Lcbarlabel = ['m/s']*len(Lplot)
+Lminval = [3]*len(Lplot)
+Lmaxval = [11]*len(Lplot)
+Lstep = [0.5]*len(Lplot)
+Lstepticks = Lstep
+Lfacconv = [1]*len(Lplot)
+LaddWhite = [False]*len(Lplot)
+
+fig = Panel.psectionV(Lxx=LaxeX, Lzz=LaxeZ, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, 
+                                Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv, 
+                                LaddWhite_cm=LaddWhite, Lylim=Lylim)
+Panel.save_graph(3,fig)
+################################################################
+#########          PANEL 4
+###############################################################
+Panel = PanelPlot(2,3, [25,14],'MEAN_V', titlepad=25, minmaxpad=1.04, timepad=-0.07, colorbarpad=0.03, labelcolorbarpad = 13, colorbaraspect=40)
+var='MEAN_V'
+Lplot = [Dvar['f1'][(LG_CART,var)],Dvar['f1'][(LG_NEB,var)],Dvar['f1'][(LG_CLEAR,var)],
+         Dvar['f1'][(LG_CS1,var)],Dvar['f1'][(LG_CS2,var)],Dvar['f1'][(LG_CS3,var)]]
+
+Lcbarlabel = ['m/s']*len(Lplot)
+Lminval = [0]*len(Lplot)
+Lmaxval = [3.2]*len(Lplot)
+Lstep = [0.2]*len(Lplot)
+Lstepticks = Lstep
+Lfacconv = [1]*len(Lplot)
+LaddWhite = [False]*len(Lplot)
+
+fig = Panel.psectionV(Lxx=LaxeX, Lzz=LaxeZ, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, 
+                                Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv, 
+                                LaddWhite_cm=LaddWhite, Lylim=Lylim)
+Panel.save_graph(4,fig)
+################################################################
+#########          PANEL 5
+###############################################################
+Panel = PanelPlot(2,3, [25,14],'MEAN_W', titlepad=25, minmaxpad=1.04, timepad=-0.07, colorbarpad=0.03, labelcolorbarpad = 13, colorbaraspect=40)
+var='MEAN_W'
+Lplot = [Dvar['f1'][(LG_CART,var)],Dvar['f1'][(LG_NEB,var)],Dvar['f1'][(LG_CLEAR,var)],
+         Dvar['f1'][(LG_CS1,var)],Dvar['f1'][(LG_CS2,var)],Dvar['f1'][(LG_CS3,var)]]
+
+Lcbarlabel = ['m/s']*len(Lplot)
+Lminval = [0., 0., -0.18, 0., -1.25, -1.2]
+Lmaxval = [4.25, 4.25, 0., 5.25, 2.75, 0.8]
+Lstep = [0.25, 0.25, 0.01, 0.25, 0.25, 0.1]
+Lstepticks = Lstep
+Lfacconv = [1]*len(Lplot)
+LaddWhite = [False]*len(Lplot)
+
+fig = Panel.psectionV(Lxx=LaxeX, Lzz=LaxeZ, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, 
+                                Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv, 
+                                LaddWhite_cm=LaddWhite, Lylim=Lylim)
+Panel.save_graph(5,fig)
+################################################################
+#########          PANEL 6
+###############################################################
+Panel = PanelPlot(2,3, [25,14],'MEAN_THL', titlepad=25, minmaxpad=1.04, timepad=-0.07, colorbarpad=0.03, labelcolorbarpad = 13, colorbaraspect=40)
+var='MEAN_THL'
+Lplot = [Dvar['f1'][(LG_CART,var)],Dvar['f1'][(LG_NEB,var)],Dvar['f1'][(LG_CLEAR,var)],
+         Dvar['f1'][(LG_CS1,var)],Dvar['f1'][(LG_CS2,var)],Dvar['f1'][(LG_CS3,var)]]
+
+Lcbarlabel = ['K']*len(Lplot)
+Lminval = [300]*len(Lplot)
+Lmaxval = [320]*len(Lplot)
+Lstep = [0.5]*len(Lplot)
+Lstepticks = Lstep
+Lfacconv = [1]*len(Lplot)
+LaddWhite = [False]*len(Lplot)
+
+fig = Panel.psectionV(Lxx=LaxeX, Lzz=LaxeZ, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, 
+                                Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv, 
+                                LaddWhite_cm=LaddWhite, Lylim=Lylim)
+Panel.save_graph(6,fig)
+################################################################
+#########          PANEL 7
+###############################################################
+Panel = PanelPlot(2,3, [25,14],'MEAN_SV 1', titlepad=25, minmaxpad=1.04, timepad=-0.07, colorbarpad=0.03, labelcolorbarpad = 13, colorbaraspect=40)
+var='MEAN_SV'
+Lplot = [Dvar['f1'][(LG_CART,var)][:,:,0],Dvar['f1'][(LG_NEB,var)][:,:,0],Dvar['f1'][(LG_CLEAR,var)][:,:,0],
+         Dvar['f1'][(LG_CS1,var)][:,:,0] ,Dvar['f1'][(LG_CS2,var)][:,:,0],Dvar['f1'][(LG_CS3,var)][:,:,0]]
+
+Lcbarlabel = ['kg/kg']*len(Lplot)
+Lminval = [2.5,9, 2.5, 2.5, 2, 2 ]
+Lmaxval = [50, 35, 50, 50, 34, 34]
+Lstep = [2.5, 1, 2.5, 2.5, 2, 2]
+Lstepticks = Lstep
+Lfacconv = [1]*len(Lplot)
+LaddWhite = [False]*len(Lplot)
+
+fig = Panel.psectionV(Lxx=LaxeX, Lzz=LaxeZ, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, 
+                                Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv, 
+                                LaddWhite_cm=LaddWhite, Lylim=Lylim)
+Panel.save_graph(7,fig)
+################################################################
+#########          PANEL 8
+###############################################################
+Panel = PanelPlot(2,3, [25,14],'MEAN_SV 2', titlepad=25, minmaxpad=1.04, timepad=-0.07, colorbarpad=0.03, labelcolorbarpad = 13, colorbaraspect=40)
+var='MEAN_SV'
+Lplot = [Dvar['f1'][(LG_CART,var)][:,:,1],Dvar['f1'][(LG_NEB,var)][:,:,1],Dvar['f1'][(LG_CLEAR,var)][:,:,1],
+         Dvar['f1'][(LG_CS1,var)][:,:,1] ,Dvar['f1'][(LG_CS2,var)][:,:,1],Dvar['f1'][(LG_CS3,var)][:,:,1]]
+
+Lminval = [0.1, 0.3, 0.1, 0.1, 0.1, 0.1 ]
+Lmaxval = [2.1, 1.6, 2.1, 1.7, 2.1, 1.75]
+Lstep = [0.1, 0.05, 0.1, 0.1, 0.1, 0.05]
+Lstepticks = Lstep
+Lcolormap = ['gist_rainbow_r']*len(Lplot)
+fig = Panel.psectionV(Lxx=LaxeX, Lzz=LaxeZ, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, 
+                                Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv, 
+                                LaddWhite_cm=LaddWhite, Lylim=Lylim)
+Panel.save_graph(8,fig)
+################################################################
+#########          PANEL 9
+###############################################################
+Panel = PanelPlot(2,3, [25,14],'MEAN_SV 3', titlepad=25, minmaxpad=1.04, timepad=-0.07, colorbarpad=0.03, labelcolorbarpad = 13, colorbaraspect=40)
+var='MEAN_SV'
+Lplot = [Dvar['f1'][(LG_CART,var)][:,:,2],Dvar['f1'][(LG_NEB,var)][:,:,2],Dvar['f1'][(LG_CLEAR,var)][:,:,2],
+         Dvar['f1'][(LG_CS1,var)][:,:,2] ,Dvar['f1'][(LG_CS2,var)][:,:,2],Dvar['f1'][(LG_CS3,var)][:,:,2]]
+
+Lminval = [0.1,0.05, 0.1, 0.05, 0.1, 0.1 ]
+Lmaxval = [2.3, 1.05, 2.3, 1.05, 3.5, 3.1]
+Lstep = [0.1, 0.05, 0.1, 0.05, 0.1, 0.1]
+Lstepticks = Lstep
+fig = Panel.psectionV(Lxx=LaxeX, Lzz=LaxeZ, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, 
+                                Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv, 
+                                LaddWhite_cm=LaddWhite, Lylim=Lylim)
+Panel.save_graph(9,fig)
+################################################################
+#########          PANEL 10
+###############################################################
+Panel = PanelPlot(2,3, [25,14],'MEAN_RT', titlepad=25, minmaxpad=1.04, timepad=-0.07, colorbarpad=0.03, labelcolorbarpad = 13, colorbaraspect=40)
+var='MEAN_RT'
+Lplot = [Dvar['f1'][(LG_CART,var)],Dvar['f1'][(LG_NEB,var)],Dvar['f1'][(LG_CLEAR,var)],
+         Dvar['f1'][(LG_CS1,var)],Dvar['f1'][(LG_CS2,var)],Dvar['f1'][(LG_CS3,var)]]
+
+Lcbarlabel = ['kg/kg']*len(Lplot)
+Lminval = [0.01]*len(Lplot)
+Lmaxval = [0.019]*len(Lplot)
+Lstep = [0.001]*len(Lplot)
+Lstepticks = Lstep
+Lfacconv = [1]*len(Lplot)
+LaddWhite = [False]*len(Lplot)
+
+fig = Panel.psectionV(Lxx=LaxeX, Lzz=LaxeZ, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, 
+                                Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv, 
+                                LaddWhite_cm=LaddWhite, Lylim=Lylim)
+Panel.save_graph(10,fig)
+################################################################
+#########          PANEL 11
+###############################################################
+Panel = PanelPlot(2,3, [25,14],'MEAN_MF', titlepad=25, minmaxpad=1.04, timepad=-0.07, colorbarpad=0.03, labelcolorbarpad = 13, colorbaraspect=40)
+var='MEAN_MF'
+Lplot = [Dvar['f1'][(LG_CART,var)],Dvar['f1'][(LG_NEB,var)],Dvar['f1'][(LG_CLEAR,var)],
+         Dvar['f1'][(LG_CS1,var)],Dvar['f1'][(LG_CS2,var)],Dvar['f1'][(LG_CS3,var)]]
+
+Lcbarlabel = ['kg/s']*len(Lplot)
+Lminval = [-0.5E-3, 0.2, -0.18, 0.2, -1.2, -1.2]
+Lmaxval = [0.2E-3, 4.2, 0.,4.2, 2.5, 1.1 ]
+Lstep = [0.05E-3, 0.2, 0.01, 0.2, 0.1, 0.1]
+Lstepticks = Lstep
+Lfacconv = [1]*len(Lplot)
+LaddWhite = [False]*len(Lplot)
+
+fig = Panel.psectionV(Lxx=LaxeX, Lzz=LaxeZ, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, 
+                                Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv, 
+                                LaddWhite_cm=LaddWhite, Lylim=Lylim)
+Panel.save_graph(11,fig)
+################################################################
+#########          PANEL 12
+###############################################################
+Panel = PanelPlot(2,3, [25,14],'SBG_WTHL', titlepad=25, minmaxpad=1.04, timepad=-0.07, colorbarpad=0.03, labelcolorbarpad = 13, colorbaraspect=40)
+var='SBG_WTHL'
+Lplot = [Dvar['f1'][(LG_SBGCART,var)],Dvar['f1'][(LG_SBGNEB,var)],Dvar['f1'][(LG_SBGCLEAR,var)],
+         Dvar['f1'][(LG_SBGCS1,var)],Dvar['f1'][(LG_SBGCS2,var)],Dvar['f1'][(LG_SBGCS3,var)]]
+
+Lcbarlabel = ['g/kg']*len(Lplot)
+Lminval = [-0.035,-0.8, -0.035, -0.8, -0.4, -0.04 ]
+Lmaxval = [0.105, 0., 0.105, 0., 0.05, 0.085]
+Lstep = [0.005, 0.05, 0.005, 0.05, 0.005, 0.005]
+Lstepticks = Lstep
+Lfacconv = [1]*len(Lplot)
+LaddWhite = [False]*len(Lplot)
+
+fig = Panel.psectionV(Lxx=LaxeX, Lzz=LaxeZ, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, 
+                                Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv, 
+                                LaddWhite_cm=LaddWhite, Lylim=Lylim)
+Panel.save_graph(12,fig)
+################################################################
+#########          PANEL 13
+###############################################################
+Panel = PanelPlot(2,3, [25,14],'SBG_WRT', titlepad=25, minmaxpad=1.04, timepad=-0.07, colorbarpad=0.03, labelcolorbarpad = 13, colorbaraspect=40)
+var='SBG_WRT'
+Lplot = [Dvar['f1'][(LG_SBGCART,var)],Dvar['f1'][(LG_SBGNEB,var)],Dvar['f1'][(LG_SBGCLEAR,var)],
+         Dvar['f1'][(LG_SBGCS1,var)],Dvar['f1'][(LG_SBGCS2,var)],Dvar['f1'][(LG_SBGCS3,var)]]
+
+Lcbarlabel = ['m/s kg/kg']*len(Lplot)
+Lminval = [0.]*len(Lplot)
+Lmaxval = [0.15E-3, 0.16E-2, 0.15E-3, 0.16E-2, 0.85E-3, 0.105E-3]
+Lstep = [0.1E-4, 0.05E-3, 0.1E-4,0.05E-3, 0.05E-3, 0.05E-4 ]
+Lstepticks = Lstep
+Lfacconv = [1]*len(Lplot)
+LaddWhite = [True]*len(Lplot)
+
+fig = Panel.psectionV(Lxx=LaxeX, Lzz=LaxeZ, Lvar=Lplot, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, 
+                                Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lfacconv=Lfacconv, 
+                                LaddWhite_cm=LaddWhite, Lylim=Lylim)
+Panel.save_graph(13,fig)

MY_RUN/INTEGRATION_CASES/HPC/ARMCU_LES/PYTHON/run_python

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

MY_RUN/INTEGRATION_CASES/HPC/AZF_2M/001_pgd1/PRE_PGD1.nam

+&NAM_CONFIO  LCDF4=T, LLFIOUT=F, LLFIREAD=F /
+&NAM_PGDFILE CPGDFILE='PGD_AZF_2.5km' /       
+&NAM_CONF_PGD JPHEXT=1,NHALO_MNH=1 /     
+&NAM_CONF_PROJ
+ XLAT0=43.567, XLON0=1.439, XRPK=0.7, XBETA=0. /
+&NAM_CONF_PROJ_GRID
+ XLATCEN=43.567,  XLONCEN=1.439,
+ NIMAX=50, NJMAX=50,
+ XDX=2500., XDY=2500. /
+&NAM_COVER     YCOVER='ECOCLIMAP_v2.0', YCOVERFILETYPE='DIRECT' /
+&NAM_ZS        YZS='gtopo30',    YZSFILETYPE='DIRECT' /
+&NAM_ISBA      YCLAY='clay_fao', YCLAYFILETYPE='DIRECT' ,
+               YSAND='sand_fao', YSANDFILETYPE='DIRECT' /

MY_RUN/INTEGRATION_CASES/HPC/AZF_2M/001_pgd1/clean_prep_pgd_xyz

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

MY_RUN/INTEGRATION_CASES/HPC/AZF_2M/001_pgd1/run_prep_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"
+#
+
+ln -sf $MESONH/PGD/ECOCLIMAP_v2.0.* .
+ln -sf $MESONH/PGD/gtopo* .
+ln -sf $MESONH/PGD/clay* .
+ln -sf $MESONH/PGD/sand* .
+
+set -x
+set -e
+
+time ${MPIRUN} PREP_PGD${XYZ}
+
+rm -f *.dir *.hdr
+cd ../002_pgd2
+sbatch run_prep_pgd_xyz

MY_RUN/INTEGRATION_CASES/HPC/AZF_2M/002_pgd2/PRE_PGD1.nam

+&NAM_CONFIO  LCDF4=T, LLFIOUT=F, LLFIREAD=F /
+&NAM_CONF_PGD JPHEXT=1,NHALO_MNH=1 /    
+&NAM_PGDFILE CPGDFILE='PGD_AZF_500m' NHALO=30  /            
+&NAM_PGD_GRID YINIFILE='PGD_AZF_2.5km', YINIFILETYPE='MESONH'/
+&NAM_INIFILE_CONF_PROJ IXOR=9,  
+                       IYOR=22,   
+                       IXSIZE=24 ,
+                       IYSIZE=18 ,
+                       IDXRATIO=5 ,
+                       IDYRATIO=5 /
+&NAM_FRAC              LECOCLIMAP = T
+/
+&NAM_PGD_SCHEMES CNATURE='ISBA',
+                 CTOWN='TEB' ,
+                 CWATER='WATFLX' /
+&NAM_COVER     YCOVER='ECOCLIMAP_v2.0', YCOVERFILETYPE='DIRECT' /
+&NAM_ZS        YZS='france250',    YZSFILETYPE='DIRECT' /
+&NAM_ISBA      YCLAY='clay_fao', YCLAYFILETYPE='DIRECT' ,
+               YSAND='sand_fao', YSANDFILETYPE='DIRECT' /

MY_RUN/INTEGRATION_CASES/HPC/AZF_2M/002_pgd2/clean_prep_pgd_xyz

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

MY_RUN/INTEGRATION_CASES/HPC/AZF_2M/002_pgd2/run_prep_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"
+#
+
+ln -sf $MESONH/PGD/ECOCLIMAP_v2.0.* .
+ln -sf $MESONH/PGD/france250* .
+ln -sf $MESONH/PGD/clay* .
+ln -sf $MESONH/PGD/sand* .
+
+ln -sf ../001_pgd1/PGD_AZF* .
+
+set -x
+set -e
+
+time ${MPIRUN} PREP_PGD${XYZ}
+
+rm -f *.hdr *.dir
+cd ../003_nest
+sbatch run_prep_nest_pgd_xyz

MY_RUN/INTEGRATION_CASES/HPC/AZF_2M/003_nest/PRE_NEST_PGD1.nam

+&NAM_CONFIO  LCDF4=T, LLFIOUT=F, LLFIREAD=F /
+&NAM_PGD1 YPGD1= 'PGD_AZF_2.5km' /       
+&NAM_PGD2 YPGD2= 'PGD_AZF_500m', IDAD = 1 /
+&NAM_PGD3 /
+&NAM_PGD4 /
+&NAM_PGD5 /
+&NAM_PGD6 /
+&NAM_PGD7 /
+&NAM_PGD8 /
+&NAM_NEST_PGD YNEST= 'ne' /
+&NAM_CONF_NEST JPHEXT=1,NHALO_MNH=1 / 

MY_RUN/INTEGRATION_CASES/HPC/AZF_2M/003_nest/clean_prep_nest_pgd_xyz

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