Skip to content
Snippets Groups Projects

Compare revisions

Changes are shown as if the source revision was being merged into the target revision. Learn more about comparing revisions.

Source

Select target project
No results found

Target

Select target project
  • mesonh/mesonh-code
  • quentin.rodier/mesonh-code-fork
  • benoit.vie/mesonh-code
  • joris.pianezze/mesonh-code
  • 8qq4g5s7/mesonh-code
  • jean.baptiste.filippi/meso-nh-fire-code
  • fdl68d9p/mesonh-code-sophia
7 results
Show changes
Showing
with 1798 additions and 0 deletions
#!/bin/sh
#SBATCH -J run1_16jan
#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 04: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 16"
set -x
set -e
ls -lrt
rm -f 16JA1.1.CEN4T.*
cp EXSEG1.nam_run1 EXSEG1.nam
time ${MPIRUN} MESONH${XYZ}
mv OUTPUT_LISTING0 OUTPUT_LISTING0_run1
mv OUTPUT_LISTING1 OUTPUT_LISTING1_run1
ls -lrt
rm -f file_for_xtransfer pipe_name PRESSURE REMAP*
rm -f EXSEG1.nam
ls -lrt
sbatch run_spa_real
ja
#!/bin/sh
#SBATCH -J run2_16jan
#SBATCH -N 1 # nodes number (=NBP)
#SBATCH -n 16 # CPUs number (on all nodes) (=NBP*TPN)
#SBATCH -o run2_16jan.eo%j #
#SBATCH -e run2_16jan.eo%j #
#SBATCH -t 04: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 16"
set -x
set -e
ls -lrt
rm -f 16JA2.?.CEN4T.*
cp EXSEG1.nam_run2 EXSEG1.nam
cp EXSEG2.nam_run2 EXSEG2.nam
time ${MPIRUN} MESONH${XYZ}
mv OUTPUT_LISTING0 OUTPUT_LISTING0_run2
mv OUTPUT_LISTING1 OUTPUT_LISTING1_run2
mv OUTPUT_LISTING2 OUTPUT_LISTING2_run2
ls -lrt
rm -f EXSEG1.nam EXSEG2.nam
rm -f file_for_xtransfer pipe_name PRESSURE REMAP*
ls -lrt
sbatch run_diag
ja
#!/bin/sh
#SBATCH -J pgd_16jan
#SBATCH -N 1 # nodes number
#SBATCH -n 4 # CPUs number (on all nodes)
#SBATCH -o pgd_16jan.eo%j #
#SBATCH -e pgd_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
. ~rodierq/DEV_57/MNH-PHYEX070-b95d84d7/conf/profile_mesonh-LXifort-R8I4-MNH-V5-6-2-ECRAD140-MPIAUTO-O2
export MPIRUN="Mpirun -np 4"
set -x
set -e
ln -sf $MESONH/PGD/ECOCLIMAP_v2.0.* .
ln -sf $MESONH/PGD/gtopo30.* .
ln -sf $MESONH/PGD/CLAY_HWSD_MOY.* .
ln -sf $MESONH/PGD/SAND_HWSD_MOY.* .
ls -lrt
rm -f 16JAN98_36km.???
cp PRE_PGD1.nam_1 PRE_PGD1.nam
time ${MPIRUN} PREP_PGD${XYZ}
mv OUTPUT_LISTING0 OUTPUT_LISTING0_pgd36km
ls -lrt
rm -f 16JAN98_9km.???
cp PRE_PGD1.nam_2 PRE_PGD1.nam
time ${MPIRUN} PREP_PGD${XYZ}
mv OUTPUT_LISTING0 OUTPUT_LISTING0_pgd9km
ls -lrt
rm -f *.neste1.*
time ${MPIRUN} PREP_NEST_PGD${XYZ}
mv OUTPUT_LISTING0 OUTPUT_LISTING0_nest
ls -lrt
rm -f *.dir *.hdr
rm -f PRE_PGD1.nam
rm -f file_for_xtransfer pipe_name
sbatch run_prep_real_case
ja
#!/bin/sh
#SBATCH -J prep_16jan
#SBATCH -N 1 # nodes number (=NBP)
#SBATCH -n 20 # CPUs number (on all nodes) (=NBP*TPN)
#SBATCH -o prep_16jan.eo%j #
#SBATCH -e prep_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
. ~rodierq/DEV_57/MNH-PHYEX070-b95d84d7/conf/profile_mesonh-LXifort-R8I4-MNH-V5-6-2-ECRAD140-MPIAUTO-O2
export MPIRUN="Mpirun -np 20"
set -x
set -e
ln -sf $HOME/SAVE/GRIB_KTEST/ecmwf.EI.1998011* .
ls -lrt
rm -f 15JAN_12_MNH.???
cp PRE_REAL1.nam_1 PRE_REAL1.nam
time ${MPIRUN} PREP_REAL_CASE${XYZ}
mv OUTPUT_LISTING0 OUTPUT_LISTING0_prep1
ls -lrt
rm -f 15JAN_18_MNH.???
cp PRE_REAL1.nam_2 PRE_REAL1.nam
time ${MPIRUN} PREP_REAL_CASE${XYZ}
mv OUTPUT_LISTING0 OUTPUT_LISTING0_prep2
ls -lrt
rm -f 16JAN_00_MNH.???
cp PRE_REAL1.nam_3 PRE_REAL1.nam
time ${MPIRUN} PREP_REAL_CASE${XYZ}
mv OUTPUT_LISTING0 OUTPUT_LISTING0_prep3
ls -lrt
rm -f 16JAN_06_MNH.???
cp PRE_REAL1.nam_4 PRE_REAL1.nam
time ${MPIRUN} PREP_REAL_CASE${XYZ}
mv OUTPUT_LISTING0 OUTPUT_LISTING0_prep4
ls -lrt
rm -f PRE_REAL1.nam
rm -f file_for_xtransfer pipe_name
rm -f ecmwf.*
sbatch run_mesonh1
ja
#!/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
&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. /
&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.
#!/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
#!/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
&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. /
&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.
#!/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
#!/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
#!/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)
#!/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)
#!/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
&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' /
set -x
rm -f *.dir *.hdr *.dat *.bin
rm -f PGD* OUTPUT_LISTING* OUTPUT_TRANSFER pipe* *.tex
#!/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
&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' /