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
Hide whitespace changes
Inline Side-by-side
Showing
with 1654 additions and 0 deletions
MY_RUN/INTEGRATION_CASES/HPC/16JAN/RUN/run_spa_real 0 → 100755
View file @ 6bb35ca8
#!/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 0 → 100644
View file @ 6bb35ca8
&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 0 → 100644
View file @ 6bb35ca8
&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 0 → 100755
View file @ 6bb35ca8
#!/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 0 → 100755
View file @ 6bb35ca8
#!/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 0 → 100644
View file @ 6bb35ca8
&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 0 → 100644
View file @ 6bb35ca8
&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 0 → 100755
View file @ 6bb35ca8
#!/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 0 → 100755
View file @ 6bb35ca8
#!/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 0 → 100644
View file @ 6bb35ca8
#!/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 0 → 100644
View file @ 6bb35ca8
#!/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 0 → 100755
View file @ 6bb35ca8
#!/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 0 → 100644
View file @ 6bb35ca8
&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 0 → 100644
View file @ 6bb35ca8
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 0 → 100755
View file @ 6bb35ca8
#!/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 0 → 100644
View file @ 6bb35ca8
&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 0 → 100644
View file @ 6bb35ca8
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 0 → 100755
View file @ 6bb35ca8
#!/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 0 → 100644
View file @ 6bb35ca8
&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 0 → 100644
View file @ 6bb35ca8
set -x
rm -f PGD_AZF* OUTPUT_LISTING* OUTPUT_TRANSFER pipe* *.tex

Documentation | Mentions légales | CGU | Accessibilité : non conforme