From ce1a28f6bc33f94cbdd21ba439fe3ae98938af78 Mon Sep 17 00:00:00 2001
From: Gaelle DELAUTIER <gaelle.delautier@meteo.fr>
Date: Thu, 17 May 2018 13:07:44 +0200
Subject: [PATCH] Gaelle 17/5/2018 : new ktest 005_ARM
---
.../005_ARM/001_prep_ideal/PRE_IDEA1.nam | 104 ++++
.../001_prep_ideal/clean_prep_ideal_case_xyz | 1 +
.../001_prep_ideal/run_prep_ideal_case_xyz | 6 +
MY_RUN/KTEST/005_ARM/002_mesonh/EXSEG1.nam | 198 +++++++
.../KTEST/005_ARM/002_mesonh/clean_mesonh_xyz | 2 +
.../KTEST/005_ARM/002_mesonh/run_mesonh_xyz | 6 +
.../KTEST/005_ARM/003_diaprog/clean_diaprog | 3 +
MY_RUN/KTEST/005_ARM/003_diaprog/dir_arm | 77 +++
MY_RUN/KTEST/005_ARM/003_diaprog/run_diaprog | 7 +
MY_RUN/KTEST/005_ARM/004_ncl/clean_ncl | 2 +
MY_RUN/KTEST/005_ARM/004_ncl/plot_arm.ncl | 513 ++++++++++++++++++
MY_RUN/KTEST/005_ARM/004_ncl/run_ncl | 34 ++
MY_RUN/KTEST/005_ARM/Makefile | 11 +
MY_RUN/KTEST/005_ARM/run_005_ARM | 3 +
src/Makefile | 2 +-
15 files changed, 968 insertions(+), 1 deletion(-)
create mode 100755 MY_RUN/KTEST/005_ARM/001_prep_ideal/PRE_IDEA1.nam
create mode 100755 MY_RUN/KTEST/005_ARM/001_prep_ideal/clean_prep_ideal_case_xyz
create mode 100755 MY_RUN/KTEST/005_ARM/001_prep_ideal/run_prep_ideal_case_xyz
create mode 100755 MY_RUN/KTEST/005_ARM/002_mesonh/EXSEG1.nam
create mode 100755 MY_RUN/KTEST/005_ARM/002_mesonh/clean_mesonh_xyz
create mode 100755 MY_RUN/KTEST/005_ARM/002_mesonh/run_mesonh_xyz
create mode 100755 MY_RUN/KTEST/005_ARM/003_diaprog/clean_diaprog
create mode 100755 MY_RUN/KTEST/005_ARM/003_diaprog/dir_arm
create mode 100755 MY_RUN/KTEST/005_ARM/003_diaprog/run_diaprog
create mode 100755 MY_RUN/KTEST/005_ARM/004_ncl/clean_ncl
create mode 100644 MY_RUN/KTEST/005_ARM/004_ncl/plot_arm.ncl
create mode 100755 MY_RUN/KTEST/005_ARM/004_ncl/run_ncl
create mode 100755 MY_RUN/KTEST/005_ARM/Makefile
create mode 100755 MY_RUN/KTEST/005_ARM/run_005_ARM
diff --git a/MY_RUN/KTEST/005_ARM/001_prep_ideal/PRE_IDEA1.nam b/MY_RUN/KTEST/005_ARM/001_prep_ideal/PRE_IDEA1.nam
new file mode 100755
index 000000000..ae229db2b
--- /dev/null
+++ b/MY_RUN/KTEST/005_ARM/001_prep_ideal/PRE_IDEA1.nam
@@ -0,0 +1,104 @@
+&NAM_CONFIO LCDF4=T, LLFIOUT=T, LLFIREAD=F/
+&NAM_REAL_PGD /
+&NAM_DIMn_PRE NIMAX=1, NJMAX=1 /
+&NAM_CONF_PRE LCARTESIAN=.TRUE., NVERB=5,
+ CIDEAL='RSOU', CZS='FLAT', LFORCING=.TRUE., LPACK=.FALSE.,
+ LBOUSS=.FALSE., CEQNSYS='DUR', LPERTURB=.FALSE.,
+ JPHEXT=1,NHALO=1 /
+&NAM_PERT_PRE /
+&NAM_CONFn LUSERV=.TRUE. /
+&NAM_GRID_PRE XLAT0=35.762 XLON0=-96. /
+&NAM_GRIDH_PRE XDELTAX=40000., XDELTAY=40000. /
+&NAM_VER_GRID LTHINSHELL=.TRUE., NKMAX=100,
+ ZDZGRD=40., ZDZTOP=40.,
+ ZZMAX_STRGRD=1000. , ZSTRGRD=0., ZSTRTOP=0. /
+&NAM_LUNITn CINIFILE='eurocs_noground', CINIFILEPGD='eurocs_noground_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.
diff --git a/MY_RUN/KTEST/005_ARM/001_prep_ideal/clean_prep_ideal_case_xyz b/MY_RUN/KTEST/005_ARM/001_prep_ideal/clean_prep_ideal_case_xyz
new file mode 100755
index 000000000..3b366f16a
--- /dev/null
+++ b/MY_RUN/KTEST/005_ARM/001_prep_ideal/clean_prep_ideal_case_xyz
@@ -0,0 +1 @@
+rm -f eurocs* OUTPUT_LISTING* pipe* *.tex
diff --git a/MY_RUN/KTEST/005_ARM/001_prep_ideal/run_prep_ideal_case_xyz b/MY_RUN/KTEST/005_ARM/001_prep_ideal/run_prep_ideal_case_xyz
new file mode 100755
index 000000000..19645e4d7
--- /dev/null
+++ b/MY_RUN/KTEST/005_ARM/001_prep_ideal/run_prep_ideal_case_xyz
@@ -0,0 +1,6 @@
+#!/bin/sh
+set -x
+set -e
+rm -f eurocs* OUTPUT_LISTING* pipe* *.tex
+time ${MONORUN} PREP_IDEAL_CASE${XYZ}
+#ddd --directory=~/DEV/MNH.V4.6.2/src/dir_obj_bug2 PREP_PGD${XYZ}
diff --git a/MY_RUN/KTEST/005_ARM/002_mesonh/EXSEG1.nam b/MY_RUN/KTEST/005_ARM/002_mesonh/EXSEG1.nam
new file mode 100755
index 000000000..5633d913d
--- /dev/null
+++ b/MY_RUN/KTEST/005_ARM/002_mesonh/EXSEG1.nam
@@ -0,0 +1,198 @@
+&NAM_CONFIO LCDF4=T, LLFIOUT=T, LLFIREAD=F /
+&NAM_LUNITn CINIFILE = "eurocs_noground", CINIFILEPGD = "eurocs_noground_pgd" /
+&NAM_CONFn LUSERV = .TRUE./
+&NAM_DYNn XTSTEP = 60. /
+&NAM_ADVn CUVW_ADV_SCHEME = "CEN4TH" ,CTEMP_SCHEME='LEFR'/
+&NAM_PARAMn CTURB = "TKEL", CRAD = "NONE",
+ CCLOUD = "ICE3", CDCONV="NONE", CSCONV="EDKF" /
+&NAM_PARAM_RADn /
+&NAM_PARAM_KAFRn XDTCONV= 60., LDIAGCONV=.TRUE./
+&NAM_PARAM_GROUNDN /
+&NAM_SEAFLUXn /
+&NAM_ISBAn /
+&NAM_DIAG_SURFn /
+&NAM_LBCn /
+&NAM_TURBn XIMPL = 1., CTURBLEN = "BL89", CTURBDIM = "1DIM" ,
+ LTURB_DIAG = .TRUE., LRMC01=.TRUE.,
+ LTURB_FLX = .TRUE., LSUBG_COND =.TRUE., LSIGMAS=.TRUE.,
+ LSIG_CONV=.FALSE., CSUBG_AUCV='CLFR' /
+&NAM_PARAM_MFSHALLn CMF_UPDRAFT='EDKF',CMF_CLOUD='DIRE', LMIXUV=T, LMF_FLX=T /
+&NAM_CH_MNHCn /
+&NAM_CONF CCONF = "START", LFLAT = .TRUE., NMODEL = 1,
+ NVERB = 4, CEQNSYS='DUR', lforcing=.true.,
+ CEXP = "ARM__", CSEG = "CEN4T" /
+&NAM_DYN XSEGLEN = 54000. , XASSELIN = 0.1, LCORIO = .true.,
+ LNUMDIFU = .FALSE.,
+ 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)= 36000. /
+&NAM_BUDGET CBUTYPE = 'NONE' /
+&NAM_BU_RU LBU_RU = .FALSE. /
+&NAM_BU_RV LBU_RV = .FALSE. /
+&NAM_BU_RW LBU_RW = .FALSE. /
+&NAM_BU_RTH LBU_RTH = .FALSE. /
+&NAM_BU_RTKE LBU_RTKE = .FALSE. /
+&NAM_BU_RRV LBU_RRV = .FALSE. /
+&NAM_BU_RRC LBU_RRC = .FALSE. /
+&NAM_BU_RRR LBU_RRR = .FALSE. /
+&NAM_BU_RRI LBU_RRI = .FALSE. /
+&NAM_BU_RRS LBU_RRS = .FALSE. /
+&NAM_BU_RRG LBU_RRG = .FALSE. /
+&NAM_BU_RRH LBU_RRH = .FALSE. /
+&NAM_BU_RSV LBU_RSV = .FALSE. /
+&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.,
+ XLES_TEMP_MEAN_STEP=3600.,
+ XLES_TEMP_MEAN_START=3600., XLES_TEMP_MEAN_END=32400. /
+&NAM_BLANK /
+&NAM_FRC LGEOST_UV_FRC=.TRUE., LTEND_THRV_FRC=.TRUE. /
+&NAM_CH_SOLVER /
+&NAM_PARAM_ICE /
+&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
+/
+
diff --git a/MY_RUN/KTEST/005_ARM/002_mesonh/clean_mesonh_xyz b/MY_RUN/KTEST/005_ARM/002_mesonh/clean_mesonh_xyz
new file mode 100755
index 000000000..424f6d31b
--- /dev/null
+++ b/MY_RUN/KTEST/005_ARM/002_mesonh/clean_mesonh_xyz
@@ -0,0 +1,2 @@
+rm -f ARM* pipe_name
+rm -f eurocs* OUT*
diff --git a/MY_RUN/KTEST/005_ARM/002_mesonh/run_mesonh_xyz b/MY_RUN/KTEST/005_ARM/002_mesonh/run_mesonh_xyz
new file mode 100755
index 000000000..c8aeaf943
--- /dev/null
+++ b/MY_RUN/KTEST/005_ARM/002_mesonh/run_mesonh_xyz
@@ -0,0 +1,6 @@
+set -x
+set -e
+ln -fs ../001_prep_ideal/eurocs_noground.* .
+ln -fs ../001_prep_ideal/eurocs_noground_pgd.* .
+rm -f ARM* OUT*
+time ${MONORUN} MESONH${XYZ}
diff --git a/MY_RUN/KTEST/005_ARM/003_diaprog/clean_diaprog b/MY_RUN/KTEST/005_ARM/003_diaprog/clean_diaprog
new file mode 100755
index 000000000..f02e52272
--- /dev/null
+++ b/MY_RUN/KTEST/005_ARM/003_diaprog/clean_diaprog
@@ -0,0 +1,3 @@
+rm -f ARM**
+rm -f gmeta OUT*
+rm -f rm dir.*:*
diff --git a/MY_RUN/KTEST/005_ARM/003_diaprog/dir_arm b/MY_RUN/KTEST/005_ARM/003_diaprog/dir_arm
new file mode 100755
index 000000000..56cfecac4
--- /dev/null
+++ b/MY_RUN/KTEST/005_ARM/003_diaprog/dir_arm
@@ -0,0 +1,77 @@
+visu
+_file1_'ARM__.1.CEN4T.000'
+print groups
+lcolarea=t
+lminmax=t
+liso=f
+xhmax=3000
+LCOLZERO=t
+NCOLZERO=1
+LMINMAX=t
+LVPTUSER=T
+XVPTVL=0.20
+XVPTR=0.80
+NIMNMX=1
+LTITDEF=F
+LVPTUSER=T
+XVPTVL=0.20
+XVPTR=0.80
+CTITB1='white'
+CTITB2='white'
+CTITB3='white'
+CTITB4='white'
+CTITT1='white'
+CTITT2='white'
+CTITVAR2='white'
+LDATFILE=.FALSE.
+XISOMIN_SBG_TKE=0
+XISOMAX_SBG_TKE=2.5
+XDIAINT_SBG_TKE=0.25
+XISOMIN_MEAN_TH=298
+XISOMAX_MEAN_TH=316
+XDIAINT_MEAN_TH=1
+XISOMIN_MEAN_U=3
+XISOMAX_MEAN_U=11.
+XDIAINT_MEAN_U=0.5
+XISOMIN_MEAN_V=0
+XISOMAX_MEAN_V=3.2
+XDIAINT_MEAN_V=0.2
+XISOMIN_MEAN_RC=0
+XISOMAX_MEAN_RC=0.2E-03
+XDIAINT_MEAN_RC=0.1E-04
+XISOMIN_MEAN_RR=0
+XISOMAX_MEAN_RR=0.2E-04
+XDIAINT_MEAN_RR=0.1E-05
+XISOMIN_SBG_WTHL=-0.3
+XISOMAX_SBG_WTHL=0.12
+XDIAINT_SBG_WTHL=0.02
+XISOMIN_SBG_WRT=-0.1E-03
+XISOMAX_SBG_WRT=0.4E-03
+XDIAINT_SBG_WRT=0.05E-03
+SBG_TKE_pvt_
+MEAN_TH_pvt_
+MEAN_U_pvt_
+MEAN_V_pvt_
+MEAN_RC_pvt_
+MEAN_RR_pvt_
+SBG_WTHL_pvt_
+SBG_WRT_pvt_
+THLUP_MF_pvt_
+RTUP_MF_pvt_
+RVUP_MF_pvt_
+RCUP_MF_pvt_
+RIUP_MF_pvt_
+WUP_MF_pvt_
+MAFLX_MF_pvt_
+DETR_MF_pvt_
+ENTR_MF_pvt_
+FRCUP_MF_pvt_
+THVUP_MF_pvt_
+WTHL_MF_pvt_
+WRT_MF_pvt_
+WTHV_MF_pvt_
+WU_MF_pvt_
+WV_MF_pvt_
+quit
+
+
diff --git a/MY_RUN/KTEST/005_ARM/003_diaprog/run_diaprog b/MY_RUN/KTEST/005_ARM/003_diaprog/run_diaprog
new file mode 100755
index 000000000..e5a25c963
--- /dev/null
+++ b/MY_RUN/KTEST/005_ARM/003_diaprog/run_diaprog
@@ -0,0 +1,7 @@
+set -x
+set -e
+ln -sf ../002_mesonh/ARM__.1.CEN4T.000.* .
+rm -f dir.0*
+${POSTRUN} diaprog < dir_arm
+
+
diff --git a/MY_RUN/KTEST/005_ARM/004_ncl/clean_ncl b/MY_RUN/KTEST/005_ARM/004_ncl/clean_ncl
new file mode 100755
index 000000000..ed94e798a
--- /dev/null
+++ b/MY_RUN/KTEST/005_ARM/004_ncl/clean_ncl
@@ -0,0 +1,2 @@
+
+rm -f ARM__.1.CEN4T.000* *.png
diff --git a/MY_RUN/KTEST/005_ARM/004_ncl/plot_arm.ncl b/MY_RUN/KTEST/005_ARM/004_ncl/plot_arm.ncl
new file mode 100644
index 000000000..5d315db0f
--- /dev/null
+++ b/MY_RUN/KTEST/005_ARM/004_ncl/plot_arm.ncl
@@ -0,0 +1,513 @@
+;================================================;
+load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
+load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
+load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
+load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
+; ================================================;
+;================================================;
+begin
+;=================================================;
+; open file and read in data
+;=================================================;
+ fichier1 = addfile("ARM__.1.CEN4T.000.nc", "r")
+;==================================================;
+; Open the workstation
+;==================================================;
+ type = "png"
+ wks = gsn_open_wks(type,"visu_ARM")
+;=================================================;
+; lecture des différents champs
+;=================================================;
+kmax=100
+zhat = fichier1->ZHAT(:) ; ZHAT
+; Unstagger zhat (from grid 4 to 1)
+ nzh=new(kmax,double)
+ do k=0,kmax-2
+ nzh(k)=(zhat(k)+zhat(k+1))/2.
+ end do
+ nzh(kmax-1)=2*zhat(kmax-1)-zhat(kmax-2)
+TIME=new((/180/),double)
+TIME(0)=300
+do j=1,179
+ TIME(j)=TIME(j-1)+300
+end do
+
+
+
+vname="SBG_TKE___PROC1"
+SBGTKE = fichier1->$vname$(0,:,:,0,0)
+SBGTKE@long_name="SBG_TKE"
+SBGTKE!1="nzh"
+SBGTKE!0="TIME"
+SBG_TKE=transpose(SBGTKE)
+
+vname="MEAN_TH___PROC1"
+MEANTH = fichier1->$vname$(0,:,:,0,0) ; MEAN_TH
+MEANTH@long_name="MEAN_TH"
+MEANTH@units="K"
+MEANTH!1="nzh"
+MEANTH!0="TIME"
+MEAN_TH=transpose(MEANTH)
+
+
+vname="MEAN_U___PROC1"
+MEANU = fichier1->$vname$(0,:,:,0,0)
+MEANU!1="nzh"
+MEANU!0="TIME"
+MEAN_U=transpose(MEANU)
+MEAN_U@long_name="MEAN_U"
+
+vname="MEAN_V___PROC1"
+MEANV = fichier1->$vname$(0,:,:,0,0)
+MEANV!1="nzh"
+MEANV!0="TIME"
+MEAN_V=transpose(MEANV)
+MEAN_V@long_name="MEAN_V"
+
+vname="MEAN_RC___PROC1"
+MEANRC = fichier1->$vname$(0,:,:,0,0)
+MEANRC!1="nzh"
+MEANRC!0="TIME"
+MEAN_RC=transpose(MEANRC)
+MEAN_RC@long_name="MEAN_RC"
+
+vname="MEAN_RR___PROC1"
+MEANRR = fichier1->$vname$(0,:,:,0,0)
+MEANRR!1="nzh"
+MEANRR!0="TIME"
+MEAN_RR=transpose(MEANRR)
+MEAN_RR@long_name="MEAN_RR"
+
+vname="SBG_WTHL___PROC1"
+SBGWTHL = fichier1->$vname$(0,:,:,0,0)
+SBGWTHL!1="nzh"
+SBGWTHL!0="TIME"
+SBG_WTHL=transpose(SBGWTHL)
+SBG_WTHL@long_name="SBG_WTHL"
+
+vname="SBG_WRT___PROC1"
+SBGWRT = fichier1->$vname$(0,:,:,0,0)
+SBGWRT!1="nzh"
+SBGWRT!0="TIME"
+SBG_WRT=transpose(SBGWRT)
+SBG_WRT@long_name="SBG_WRT"
+
+
+vname="THLUP_MF___PROC1"
+ THLUPMF= fichier1->$vname$(0,:,:,0,0)
+THLUPMF!1="nzh"
+THLUPMF!0="TIME"
+THLUP_MF=transpose(THLUPMF)
+THLUP_MF@long_name="THLUP_MF"
+
+vname="RTUP_MF___PROC1"
+RTUPMF = fichier1->$vname$(0,:,:,0,0)
+RTUPMF!1="nzh"
+RTUPMF!0="TIME"
+RTUP_MF=transpose(RTUPMF)
+RTUP_MF@long_name="RTUP_MF"
+
+vname="RVUP_MF___PROC1"
+RVUPMF = fichier1->$vname$(0,:,:,0,0)
+RVUPMF!1="nzh"
+RVUPMF!0="TIME"
+RVUP_MF=transpose(RVUPMF)
+RVUP_MF@long_name="RVUP_MF"
+
+vname="RCUP_MF___PROC1"
+RCUPMF = fichier1->$vname$(0,:,:,0,0)
+RCUPMF!1="nzh"
+RCUPMF!0="TIME"
+RCUP_MF=transpose(RCUPMF)
+RCUP_MF@long_name="RCUP_MF"
+
+vname="RIUP_MF___PROC1"
+RIUPMF = fichier1->$vname$(0,:,:,0,0)
+RIUPMF!1="nzh"
+RIUPMF!0="TIME"
+RIUP_MF=transpose(RIUPMF)
+RIUP_MF@long_name="RIUP_MF"
+
+vname="WUP_MF___PROC1"
+WUPMF = fichier1->$vname$(0,:,:,0,0)
+WUPMF!1="nzh"
+WUPMF!0="TIME"
+WUP_MF=transpose(WUPMF)
+WUP_MF@long_name="WUP_MF"
+
+vname="MAFLX_MF___PROC1"
+MAFLXMF = fichier1->$vname$(0,:,:,0,0)
+MAFLXMF!1="nzh"
+MAFLXMF!0="TIME"
+MAFLX_MF=transpose(MAFLXMF)
+MAFLX_MF@long_name="MAFLX_MF"
+
+vname="DETR_MF___PROC1"
+DETRMF = fichier1->$vname$(0,:,:,0,0)
+DETRMF!1="nzh"
+DETRMF!0="TIME"
+DETR_MF=transpose(DETRMF)
+DETR_MF@long_name="DETR_MF"
+
+vname="ENTR_MF___PROC1"
+ENTRMF = fichier1->$vname$(0,:,:,0,0)
+ENTRMF!1="nzh"
+ENTRMF!0="TIME"
+ENTR_MF=transpose(ENTRMF)
+ENTR_MF@long_name="ENTR_MF"
+
+vname="FRCUP_MF___PROC1"
+FRCUPMF = fichier1->$vname$(0,:,:,0,0)
+FRCUPMF!1="nzh"
+FRCUPMF!0="TIME"
+FRCUP_MF=transpose(FRCUPMF)
+FRCUP_MF@long_name="FRCUP_MF"
+
+vname="THVUP_MF___PROC1"
+THVUPMF = fichier1->$vname$(0,:,:,0,0)
+THVUPMF!1="nzh"
+THVUPMF!0="TIME"
+THVUP_MF=transpose(THVUPMF)
+THVUP_MF@long_name="THVUP_MF"
+
+vname="WTHL_MF___PROC1"
+WTHLMF = fichier1->$vname$(0,:,:,0,0)
+WTHLMF!1="nzh"
+WTHLMF!0="TIME"
+WTHL_MF=transpose(WTHLMF)
+WTHL_MF@long_name="WTHL_MF"
+
+vname="WRT_MF___PROC1"
+WRTMF = fichier1->$vname$(0,:,:,0,0)
+WRTMF!1="nzh"
+WRTMF!0="TIME"
+WRT_MF=transpose(WRTMF)
+WRT_MF@long_name="WRT_MF"
+
+
+vname="WTHV_MF___PROC1"
+WTHVMF = fichier1->$vname$(0,:,:,0,0)
+WTHVMF!1="nzh"
+WTHVMF!0="TIME"
+WTHV_MF=transpose(WTHVMF)
+WTHV_MF@long_name="WTHV_MF"
+
+
+vname="WU_MF___PROC1"
+WUMF = fichier1->$vname$(0,:,:,0,0)
+WUMF!1="nzh"
+WUMF!0="TIME"
+WU_MF=transpose(WUMF)
+WU_MF@long_name="WU_MF"
+
+
+vname="WV_MF___PROC1"
+WVMF = fichier1->$vname$(0,:,:,0,0)
+WVMF!1="nzh"
+WVMF!0="TIME"
+WV_MF=transpose(WVMF)
+WV_MF@long_name="WV_MF"
+
+;=================================================;
+; Altitude des niveaux modèles
+;=================================================;
+; Unstagger zhat (from grid 4 to 1)
+ nzh=new(kmax,double)
+ do k=0,kmax-2
+ nzh(k)=(zhat(k)+zhat(k+1))/2.
+ end do
+ nzh(kmax-1)=2*zhat(kmax-1)-zhat(kmax-2)
+;=================================================;
+; Set some other basic resources
+;=================================================;
+ resmap = True
+ resmap@gsnFrame = False
+ resmap@gsnDraw = False
+ resmap@gsnMaximize = True
+ resmap@gsnPaperOrientation = "portrait"
+ resmap@gsnSpreadColors = True ; use full range of colormap
+ resmap@tiYAxisString =" "
+ resmap@cnFillOn = True ; turn on color fill
+ resmap@cnLinesOn = False ; turn off contour lines
+ ;resmap@tmXBLabelStride = 2 ; to reduce the number of labels on xaxis
+ ;resmap@lbLabelStride = 2. ; to reduce the number of labels on labelbar
+ resmap@sfYArray = nzh(0:76)
+ resmap@sfXArray = TIME
+ resmap@tmXBTickSpacingF = 10800.
+ resmap@tmXBMode = "Explicit"
+ resmap@tmXBValues = (/10800,21600,32400,43200,54000/)
+ resmap@tmXBLabels =(/3,6,9,12,15/)
+resmap@tiXAxisPosition="Left"
+resmap@tiXAxisFontHeightF=0.015
+;=================================================;
+; TRACE
+;=================================================;
+ gsn_define_colormap(wks,"rainbow") ; Choose colormap
+ colors = gsn_retrieve_colormap(wks) ; retrieve color map for editing
+ colors(2,:) = (/ 1, 1, 1 /) ; replace the first color with white color (les deux premières sont background /foreground donc c'est bien la 2 qu'il faut changer)
+ gsn_define_colormap(wks,colors) ; redefine colormap to workstation, color map now includes a gray
+
+
+
+opts=resmap
+opts@cnLevelSelectionMode = "ManualLevels"
+opts@cnLevelSpacingF=0.25
+opts@cnMinLevelValF=0
+opts@cnMaxLevelValF=2.5
+
+plot = gsn_csm_contour(wks,SBG_TKE(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+opts=resmap
+opts@cnLevelSelectionMode = "ManualLevels"
+opts@cnLevelSpacingF=1
+opts@cnMinLevelValF=298
+opts@cnMaxLevelValF=316
+
+plot = gsn_csm_contour(wks,MEAN_TH(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+
+opts=resmap
+opts@cnLevelSelectionMode = "ManualLevels"
+opts@cnLevelSpacingF=0.5
+opts@cnMinLevelValF=3
+opts@cnMaxLevelValF=11
+
+plot = gsn_csm_contour(wks,MEAN_U(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+opts=resmap
+opts@cnLevelSelectionMode = "ManualLevels"
+opts@cnLevelSpacingF=0.2
+opts@cnMinLevelValF=0.
+opts@cnMaxLevelValF=3.2
+
+plot = gsn_csm_contour(wks,MEAN_V(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+opts=resmap
+opts@cnLevelSelectionMode = "ManualLevels"
+opts@cnLevelSpacingF=0.00001
+opts@cnMinLevelValF=0
+opts@cnMaxLevelValF=0.0002
+
+plot = gsn_csm_contour(wks,MEAN_RC(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+
+opts=resmap
+opts@cnLevelSelectionMode = "ManualLevels"
+opts@cnLevelSpacingF=0.000001
+opts@cnMinLevelValF=0
+opts@cnMaxLevelValF=0.00002
+
+plot = gsn_csm_contour(wks,MEAN_RR(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+
+opts=resmap
+opts@cnLevelSelectionMode = "ManualLevels"
+opts@cnLevelSpacingF=0.02
+opts@cnMinLevelValF=-0.3
+opts@cnMaxLevelValF=0.12
+
+plot = gsn_csm_contour(wks,SBG_WTHL(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+opts=resmap
+opts@cnLevelSelectionMode = "ManualLevels"
+opts@cnLevelSpacingF=0.00005
+opts@cnMinLevelValF=-0.0001
+opts@cnMaxLevelValF=0.0004
+
+plot = gsn_csm_contour(wks,SBG_WRT(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+
+opts=resmap
+opts@cnLevelSelectionMode = "ManualLevels"
+opts@cnLevelSpacingF=1
+opts@cnMinLevelValF=301
+opts@cnMaxLevelValF=318
+
+plot = gsn_csm_contour(wks,THLUP_MF(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+opts=resmap
+opts@cnLevelSelectionMode = "ManualLevels"
+opts@cnLevelSpacingF=0.001
+opts@cnMinLevelValF=0.0035
+opts@cnMaxLevelValF=0.017
+plot = gsn_csm_contour(wks,RTUP_MF(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+opts=resmap
+opts@cnLevelSelectionMode = "ManualLevels"
+opts@cnLevelSpacingF=0.001
+opts@cnMinLevelValF=0.001
+opts@cnMaxLevelValF=0.017
+
+plot = gsn_csm_contour(wks,RVUP_MF(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+opts=resmap
+opts@cnLevelSelectionMode = "ManualLevels"
+opts@cnLevelSpacingF=0.00025
+opts@cnMinLevelValF=0.00025
+opts@cnMaxLevelValF=0.00425
+
+plot = gsn_csm_contour(wks,RCUP_MF(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+opts=resmap
+opts@cnLevelSelectionMode = "AutomaticLevels"
+plot = gsn_csm_contour(wks,RIUP_MF(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+opts=resmap
+opts@cnLevelSelectionMode = "ManualLevels"
+opts@cnLevelSpacingF=0.25
+opts@cnMinLevelValF=0.250
+opts@cnMaxLevelValF=5.5
+
+plot = gsn_csm_contour(wks,WUP_MF(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+opts=resmap
+opts@cnLevelSelectionMode = "ManualLevels"
+opts@cnLevelSpacingF=0.02
+opts@cnMinLevelValF=0.02
+opts@cnMaxLevelValF=0.32
+
+plot = gsn_csm_contour(wks,MAFLX_MF(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+opts=resmap
+opts@cnLevelSelectionMode = "ManualLevels"
+opts@cnLevelSpacingF=0.025
+opts@cnMinLevelValF=0.025
+opts@cnMaxLevelValF=0.45
+
+plot = gsn_csm_contour(wks,DETR_MF(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+opts=resmap
+opts@cnLevelSelectionMode = "ManualLevels"
+opts@cnLevelSpacingF=0.00025
+opts@cnMinLevelValF=0.00025
+opts@cnMaxLevelValF=0.005
+
+plot = gsn_csm_contour(wks,ENTR_MF(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+opts=resmap
+opts@cnLevelSelectionMode = "ManualLevels"
+opts@cnLevelSpacingF=0.01
+opts@cnMinLevelValF=0.01
+opts@cnMaxLevelValF=0.16
+
+plot = gsn_csm_contour(wks,FRCUP_MF(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+opts=resmap
+opts@cnLevelSelectionMode = "ManualLevels"
+opts@cnLevelSpacingF=1
+opts@cnMinLevelValF=303
+opts@cnMaxLevelValF=319
+
+plot = gsn_csm_contour(wks,THVUP_MF(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+opts=resmap
+opts@cnLevelSelectionMode = "ManualLevels"
+opts@cnLevelSpacingF=0.005
+opts@cnMinLevelValF=-0.075
+opts@cnMaxLevelValF=0.055
+
+plot = gsn_csm_contour(wks,WTHL_MF(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+opts=resmap
+opts@cnLevelSelectionMode = "ManualLevels"
+opts@cnLevelSpacingF=0.00001
+opts@cnMinLevelValF=0.00001
+opts@cnMaxLevelValF=0.00019
+
+plot = gsn_csm_contour(wks,WRT_MF(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+opts=resmap
+opts@cnLevelSelectionMode = "ManualLevels"
+opts@cnLevelSpacingF=0.005
+opts@cnMinLevelValF=-0.016
+opts@cnMaxLevelValF=0.075
+plot = gsn_csm_contour(wks,WTHV_MF(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+opts=resmap
+opts@cnLevelSelectionMode = "ManualLevels"
+opts@cnLevelSpacingF=0.005
+opts@cnMinLevelValF=-0.120
+opts@cnMaxLevelValF=-0.005
+
+plot = gsn_csm_contour(wks,WU_MF(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+opts=resmap
+opts@cnLevelSelectionMode = "ManualLevels"
+opts@cnLevelSpacingF=0.004
+opts@cnMinLevelValF=-0.024
+opts@cnMaxLevelValF=0.048
+plot = gsn_csm_contour(wks,WV_MF(0:76,:),opts)
+draw(plot)
+frame(wks)
+delete(opts)
+
+end
diff --git a/MY_RUN/KTEST/005_ARM/004_ncl/run_ncl b/MY_RUN/KTEST/005_ARM/004_ncl/run_ncl
new file mode 100755
index 000000000..5fd5dd23c
--- /dev/null
+++ b/MY_RUN/KTEST/005_ARM/004_ncl/run_ncl
@@ -0,0 +1,34 @@
+#MNH_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
+#MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
+#MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
+#MNH_LIC for details. version 1.
+set -x
+
+ln -sf ../002_mesonh/ARM__.1.CEN4T.000.* .
+
+ncl plot_arm.ncl
+display visu_ARM.000001.png
+display visu_ARM.000002.png
+display visu_ARM.000003.png
+display visu_ARM.000004.png
+display visu_ARM.000005.png
+display visu_ARM.000006.png
+display visu_ARM.000007.png
+display visu_ARM.000008.png
+display visu_ARM.000009.png
+display visu_ARM.000010.png
+display visu_ARM.000011.png
+display visu_ARM.000012.png
+display visu_ARM.000013.png
+display visu_ARM.000014.png
+display visu_ARM.000015.png
+display visu_ARM.000016.png
+display visu_ARM.000017.png
+display visu_ARM.000018.png
+display visu_ARM.000019.png
+display visu_ARM.000020.png
+display visu_ARM.000021.png
+display visu_ARM.000022.png
+display visu_ARM.000023.png
+display visu_ARM.000024.png
+exit 0
diff --git a/MY_RUN/KTEST/005_ARM/Makefile b/MY_RUN/KTEST/005_ARM/Makefile
new file mode 100755
index 000000000..15aadf169
--- /dev/null
+++ b/MY_RUN/KTEST/005_ARM/Makefile
@@ -0,0 +1,11 @@
+all:
+ cd 001_prep_ideal && run_prep_ideal_case_xyz
+ cd 002_mesonh && run_mesonh_xyz
+ cd 003_diaprog && run_diaprog
+ cd 004_ncl && run_ncl
+
+clean:
+ cd 001_prep_ideal && clean_prep_ideal_case_xyz
+ cd 002_mesonh && clean_mesonh_xyz
+ cd 003_diaprog && clean_diaprog
+ cd 004_ncl && clean_ncl
diff --git a/MY_RUN/KTEST/005_ARM/run_005_ARM b/MY_RUN/KTEST/005_ARM/run_005_ARM
new file mode 100755
index 000000000..9917191a7
--- /dev/null
+++ b/MY_RUN/KTEST/005_ARM/run_005_ARM
@@ -0,0 +1,3 @@
+set -e
+make clean
+make
\ No newline at end of file
diff --git a/src/Makefile b/src/Makefile
index 317f4015c..c85c2c023 100644
--- a/src/Makefile
+++ b/src/Makefile
@@ -420,7 +420,7 @@ endif
PWD_KTEST=$(SRC_MESONH)/MY_RUN/KTEST
-EXAMPLES=001_2Drelief 002_3Drelief 003_KW78 004_Reunion 007_16janvier 014_LIMA
+EXAMPLES=001_2Drelief 002_3Drelief 003_KW78 004_Reunion 005_ARM 007_16janvier 014_LIMA
examples : $(EXAMPLES)
EXAMPLES_CHIMIE=009_ICARTT 011_KW78CHEM
--
GitLab