From 585aab8340f2bef94aa9486173cefeec8911c00b Mon Sep 17 00:00:00 2001
From: Christine Lac <christine.lac@meteo.fr>
Date: Fri, 13 Dec 2013 16:21:16 +0000
Subject: [PATCH] adaptation WENO
---
.../KTEST/004_Reunion/003_mesonh/EXSEG1.nam | 6 +--
.../004_Reunion/003_mesonh/EXSEG1.nam_CEN4TH | 38 +++++++++++++++++++
.../004_Reunion/003_mesonh/EXSEG1.nam_WENO | 38 +++++++++++++++++++
.../KTEST/004_Reunion/005_diaprog/dir_Reunion | 14 +++----
.../005_ncl_extractdia/dir_exnetcdf | 18 ++++-----
.../005_ncl_extractdia/mnh_Cloud1.ncl | 8 ++--
.../005_ncl_extractdia/mnh_PressureLevel1.ncl | 10 ++---
7 files changed, 104 insertions(+), 28 deletions(-)
create mode 100644 MY_RUN/KTEST/004_Reunion/003_mesonh/EXSEG1.nam_CEN4TH
create mode 100644 MY_RUN/KTEST/004_Reunion/003_mesonh/EXSEG1.nam_WENO
diff --git a/MY_RUN/KTEST/004_Reunion/003_mesonh/EXSEG1.nam b/MY_RUN/KTEST/004_Reunion/003_mesonh/EXSEG1.nam
index fd46b178f..f588e4075 100644
--- a/MY_RUN/KTEST/004_Reunion/003_mesonh/EXSEG1.nam
+++ b/MY_RUN/KTEST/004_Reunion/003_mesonh/EXSEG1.nam
@@ -3,12 +3,12 @@
/
&NAM_LUNITn CINIFILE = "REUNION_IDEA_410",
CINIFILEPGD="REUNION_PGD_1km5_410" /
-&NAM_DYNn XTSTEP = 5.,
+&NAM_DYNn XTSTEP = 20.,
!CPRESOPT = "ZRESI",
LITRADJ = T,
LHORELAX_UVWTH = T, LHORELAX_RV = T, LVE_RELAX = T,
NRIMX = 5, NRIMY = 5, XRIMKMAX = 0.01, XT4DIFU = 500. /
-&NAM_ADVn CUVW_ADV_SCHEME = "CEN4TH" , CMET_ADV_SCHEME = "FCT2ND" /
+&NAM_ADVn CUVW_ADV_SCHEME = "WENO_K" , CMET_ADV_SCHEME = "PPM_01" /
&NAM_PARAMn CTURB = "TKEL", CRAD = "FIXE",
CCLOUD = "NONE", CDCONV = "NONE" /
&NAM_LBCn CLBCX = 2*"OPEN", CLBCY = 2*"OPEN" /
@@ -20,7 +20,7 @@
NVERB = 5, CEXP = "REUNI", CSEG = "00A20" ,
CSPLIT="BSPLITTING" /
&NAM_DYN XSEGLEN = 40.,
- LCORIO = F, LNUMDIFU = T,
+ LCORIO = F, LNUMDIFU = F,
XALKTOP = 0.01, XALZBOT = 14000. /
&NAM_FMOUT XFMOUT(1,1) = 10., XFMOUT(1,2) = 20.,
XFMOUT(1,3) = 30., XFMOUT(1,4) = 40. /
diff --git a/MY_RUN/KTEST/004_Reunion/003_mesonh/EXSEG1.nam_CEN4TH b/MY_RUN/KTEST/004_Reunion/003_mesonh/EXSEG1.nam_CEN4TH
new file mode 100644
index 000000000..cf5d50f09
--- /dev/null
+++ b/MY_RUN/KTEST/004_Reunion/003_mesonh/EXSEG1.nam_CEN4TH
@@ -0,0 +1,38 @@
+&NAM_CONFZ
+ ! NZ_VERB=5 , NB_PROCIO_R=8 , NB_PROCIO_W=1
+/
+&NAM_LUNITn CINIFILE = "REUNION_IDEA_410",
+ CINIFILEPGD="REUNION_PGD_1km5_410" /
+&NAM_DYNn XTSTEP = 5.,
+ !CPRESOPT = "ZRESI",
+ LITRADJ = T,
+ LHORELAX_UVWTH = T, LHORELAX_RV = T, LVE_RELAX = T,
+ NRIMX = 5, NRIMY = 5, XRIMKMAX = 0.01, XT4DIFU = 500. /
+&NAM_ADVn CUVW_ADV_SCHEME = "CEN4TH" , CMET_ADV_SCHEME = "PPM_01" /
+&NAM_PARAMn CTURB = "TKEL", CRAD = "FIXE",
+ CCLOUD = "NONE", CDCONV = "NONE" /
+&NAM_LBCn CLBCX = 2*"OPEN", CLBCY = 2*"OPEN" /
+&NAM_TURBn CTURBLEN = "BL89", CTURBDIM = "1DIM", LSUBG_COND = F /
+&NAM_SERIESn NIBOXL=25 NIBOXH=50 NJBOXL=20 NJBOXH=40
+ NKCLS=3 NKLOW=12 NKUP=18 NKMID=14
+ NBJSLICE=1 NJSLICEL(1)=30 NJSLICEH(1)=35 /
+&NAM_CONF CCONF = "START", NMODEL = 1,
+ NVERB = 5, CEXP = "REUNI", CSEG = "00A20" ,
+ CSPLIT="BSPLITTING" /
+&NAM_DYN XSEGLEN = 40.,
+ LCORIO = F, LNUMDIFU = T,
+ XALKTOP = 0.01, XALZBOT = 14000. /
+&NAM_FMOUT XFMOUT(1,1) = 10., XFMOUT(1,2) = 20.,
+ XFMOUT(1,3) = 30., XFMOUT(1,4) = 40. /
+&NAM_BUDGET CBUTYPE='CART', XBULEN=200.,
+ NBUKL=2, NBUKH=51, NBUIL=15, NBUIH=64, NBUJL=15, NBUJH=68,
+ LBU_KCP=.TRUE., LBU_ICP=.FALSE., LBU_JCP=.FALSE. /
+&NAM_BU_RTH LBU_RTH=.FALSE. /
+&NAM_BU_RTKE LBU_RTKE=.FALSE.,
+ NDPTKE=1, NTPTKE=1, NDISSTKE=0, NTRTKE=0 /
+&NAM_SERIES LSERIES=.TRUE. LMASKLANDSEA=.TRUE. LWMINMAX=.FALSE. /
+&NAM_SEAFLUXn CSEA_ALB = 'UNIF' /
+&NAM_ISBAn /
+&NAM_DIAG_SURFn LSURF_BUDGET=F /
+&NAM_DIAG_ISBAn /
+
diff --git a/MY_RUN/KTEST/004_Reunion/003_mesonh/EXSEG1.nam_WENO b/MY_RUN/KTEST/004_Reunion/003_mesonh/EXSEG1.nam_WENO
new file mode 100644
index 000000000..f588e4075
--- /dev/null
+++ b/MY_RUN/KTEST/004_Reunion/003_mesonh/EXSEG1.nam_WENO
@@ -0,0 +1,38 @@
+&NAM_CONFZ
+ ! NZ_VERB=5 , NB_PROCIO_R=8 , NB_PROCIO_W=1
+/
+&NAM_LUNITn CINIFILE = "REUNION_IDEA_410",
+ CINIFILEPGD="REUNION_PGD_1km5_410" /
+&NAM_DYNn XTSTEP = 20.,
+ !CPRESOPT = "ZRESI",
+ LITRADJ = T,
+ LHORELAX_UVWTH = T, LHORELAX_RV = T, LVE_RELAX = T,
+ NRIMX = 5, NRIMY = 5, XRIMKMAX = 0.01, XT4DIFU = 500. /
+&NAM_ADVn CUVW_ADV_SCHEME = "WENO_K" , CMET_ADV_SCHEME = "PPM_01" /
+&NAM_PARAMn CTURB = "TKEL", CRAD = "FIXE",
+ CCLOUD = "NONE", CDCONV = "NONE" /
+&NAM_LBCn CLBCX = 2*"OPEN", CLBCY = 2*"OPEN" /
+&NAM_TURBn CTURBLEN = "BL89", CTURBDIM = "1DIM", LSUBG_COND = F /
+&NAM_SERIESn NIBOXL=25 NIBOXH=50 NJBOXL=20 NJBOXH=40
+ NKCLS=3 NKLOW=12 NKUP=18 NKMID=14
+ NBJSLICE=1 NJSLICEL(1)=30 NJSLICEH(1)=35 /
+&NAM_CONF CCONF = "START", NMODEL = 1,
+ NVERB = 5, CEXP = "REUNI", CSEG = "00A20" ,
+ CSPLIT="BSPLITTING" /
+&NAM_DYN XSEGLEN = 40.,
+ LCORIO = F, LNUMDIFU = F,
+ XALKTOP = 0.01, XALZBOT = 14000. /
+&NAM_FMOUT XFMOUT(1,1) = 10., XFMOUT(1,2) = 20.,
+ XFMOUT(1,3) = 30., XFMOUT(1,4) = 40. /
+&NAM_BUDGET CBUTYPE='CART', XBULEN=200.,
+ NBUKL=2, NBUKH=51, NBUIL=15, NBUIH=64, NBUJL=15, NBUJH=68,
+ LBU_KCP=.TRUE., LBU_ICP=.FALSE., LBU_JCP=.FALSE. /
+&NAM_BU_RTH LBU_RTH=.FALSE. /
+&NAM_BU_RTKE LBU_RTKE=.FALSE.,
+ NDPTKE=1, NTPTKE=1, NDISSTKE=0, NTRTKE=0 /
+&NAM_SERIES LSERIES=.TRUE. LMASKLANDSEA=.TRUE. LWMINMAX=.FALSE. /
+&NAM_SEAFLUXn CSEA_ALB = 'UNIF' /
+&NAM_ISBAn /
+&NAM_DIAG_SURFn LSURF_BUDGET=F /
+&NAM_DIAG_ISBAn /
+
diff --git a/MY_RUN/KTEST/004_Reunion/005_diaprog/dir_Reunion b/MY_RUN/KTEST/004_Reunion/005_diaprog/dir_Reunion
index 67da61d2a..fa9ef813d 100644
--- a/MY_RUN/KTEST/004_Reunion/005_diaprog/dir_Reunion
+++ b/MY_RUN/KTEST/004_Reunion/005_diaprog/dir_Reunion
@@ -7,25 +7,25 @@ _file_'REUNI.1.00A20.004dg'
! ouverture de la fenetre d affichage
visu
! module du vent et vecteurs du niveau K=2
-MUMVM_K_2_ON_UMVM_K_2
+MUTVT_K_2_ON_UTVT_K_2
! temperature potentielle a 1500m
-THM_Z_1500
+THT_Z_1500
pression et vent en vecteurs a 1500m
-PABSM_Z_1500_ON_UMVM_Z_1500
+PABST_Z_1500_ON_UTVT_Z_1500
! definition d une coupe verticale
NIDEBCOU=35 NJDEBCOU=2 NLANGLE=90 NLMAX=80
! trace de la coupe verticale dans le plan horizontal
LTRACECV=T
! temperature potentielle
-THM_CV_
+THT_CV_
! plus de trace de la coupe verticale
LTRACECV=F
! difference pour la temperature potentielle avec le champ initial
-THM_CV__MINUS_LSTHM_CV_
+THT_CV__MINUS_LSTHM_CV_
! difference pour la composante du vent selon Y avec le champ initial
-VM_CV__MINUS_LSVM_CV_
+VT_CV__MINUS_LSVM_CV_
! difference pour la vitesse verticale avec le champ initial
-WM_CV__ON_ULMWM_CV_
+WT_CV__ON_ULMWM_CV_
!
! a vous d indiquer vos directives
quit
diff --git a/MY_RUN/KTEST/004_Reunion/005_ncl_extractdia/dir_exnetcdf b/MY_RUN/KTEST/004_Reunion/005_ncl_extractdia/dir_exnetcdf
index 020304e70..f39dba49c 100644
--- a/MY_RUN/KTEST/004_Reunion/005_ncl_extractdia/dir_exnetcdf
+++ b/MY_RUN/KTEST/004_Reunion/005_ncl_extractdia/dir_exnetcdf
@@ -14,16 +14,16 @@ KCDL
0
0
CONF
- UM
- VM
- WM
+ UT
+ VT
+ WT
ALT
-THM
-PABSM
-RVM
-RCM
-RRM
-RIM
+THT
+PABST
+RVT
+RCT
+RRT
+RIT
ACPRR
ZSBIS
END
diff --git a/MY_RUN/KTEST/004_Reunion/005_ncl_extractdia/mnh_Cloud1.ncl b/MY_RUN/KTEST/004_Reunion/005_ncl_extractdia/mnh_Cloud1.ncl
index 0517e0dfd..38e12ff80 100644
--- a/MY_RUN/KTEST/004_Reunion/005_ncl_extractdia/mnh_Cloud1.ncl
+++ b/MY_RUN/KTEST/004_Reunion/005_ncl_extractdia/mnh_Cloud1.ncl
@@ -43,8 +43,8 @@ begin
;=================================================;
; Read water wapor mixing ratio if it exists
; -----------------------------------------
- if(isfilevar(a,"RVM"))
- qv = a->RVM(0,:,:,:)
+ if(isfilevar(a,"RVT"))
+ qv = a->RVT(0,:,:,:)
qv = qv*1000.
qv@units = "g/kg"
qv@description="Water vapor mixing ratio"
@@ -70,8 +70,8 @@ begin
; Read ice mixing ratio if it exists
; -----------------------------------------
- if(isfilevar(a,"RIM"))
- qi = a->RIM(0,:,:,:)
+ if(isfilevar(a,"RIT"))
+ qi = a->RIT(0,:,:,:)
qi = qi*1000.
qi@units = "g/kg"
qi@description="Ice mixing ratio"
diff --git a/MY_RUN/KTEST/004_Reunion/005_ncl_extractdia/mnh_PressureLevel1.ncl b/MY_RUN/KTEST/004_Reunion/005_ncl_extractdia/mnh_PressureLevel1.ncl
index 90f709bdb..2667ac88b 100644
--- a/MY_RUN/KTEST/004_Reunion/005_ncl_extractdia/mnh_PressureLevel1.ncl
+++ b/MY_RUN/KTEST/004_Reunion/005_ncl_extractdia/mnh_PressureLevel1.ncl
@@ -34,17 +34,17 @@ begin
;=================================================;
; Read the variables we need
;=================================================;
- tk = a->THM(0,:,:,:) ; potential temperature
+ tk = a->THT(0,:,:,:) ; potential temperature
tk@description="Potential Temperature"
- u = a->UM(0,:,:,:) ; u
+ u = a->UT(0,:,:,:) ; u
- v = a->VM(0,:,:,:) ; v
+ v = a->VT(0,:,:,:) ; v
- p = a->PABSM(0,:,:,:) ; pressure
+ p = a->PABST(0,:,:,:) ; pressure
p@description="Pressure"
- rv = a->RVM(0,:,:,:)*1000 ; water vapor mixing ratio
+ rv = a->RVT(0,:,:,:)*1000 ; water vapor mixing ratio
rv@description="Water vapor mixing ratio"
ZS = a->ZSBIS(0,:,:) ; terrain
--
GitLab