adaptation WENO

MY_RUN/KTEST/007_16janvier/008_run2/EXSEG1.nam.src

@@ -8,7 +8,7 @@
            NITR = 8,
            LHORELAX_UVWTH = T, LHORELAX_RV = T, LVE_RELAX = T,
            NRIMX = 5, NRIMY = 5, XRIMKMAX = 0.0083, XT4DIFU = 5000. /
-&NAM_ADVn  CMET_ADV_SCHEME = "FCT2ND", CSV_ADV_SCHEME = "FCT2ND" /
+&NAM_ADVn  CUVW_ADV_SCHEME="CEN4TH",CMET_ADV_SCHEME = "PPM_01", CSV_ADV_SCHEME = "PPM_01" /
 &NAM_PARAMn  CCLOUD = "KESS", CTURB = "TKEL", CRAD = "ECMW",
              CDCONV = "KAFR",CSCONV="KAFR" /
 &NAM_PARAM_RADn  XDTRAD = 3600., XDTRAD_CLONLY = 3600., NRAD_COLNBR = 400 /

MY_RUN/KTEST/007_16janvier/008_run2/EXSEG1.nam.src_CEN4TH

+&NAM_CONFZ
+ NZ_VERB=5
+ ! NZ_VERB=5 , NZ_PROC=0 , NB_PROCIO_R=1 , NB_PROCIO_W=8 
+/
+&NAM_LUNITn  CINIFILE = "16JAN_06_MNH",CINIFILEPGD="16JAN98_36km.neste1"  /
+&NAM_DYNn  XTSTEP = 60.,  
+           ! CPRESOPT = "ZRESI", 
+           NITR = 8,
+           LHORELAX_UVWTH = T, LHORELAX_RV = T, LVE_RELAX = T,
+           NRIMX = 5, NRIMY = 5, XRIMKMAX = 0.0083, XT4DIFU = 5000. /
+&NAM_ADVn  CUVW_ADV_SCHEME="CEN4TH",CMET_ADV_SCHEME = "PPM_01", CSV_ADV_SCHEME = "PPM_01" /
+&NAM_PARAMn  CCLOUD = "KESS", CTURB = "TKEL", CRAD = "ECMW",
+             CDCONV = "KAFR",CSCONV="KAFR" /
+&NAM_PARAM_RADn  XDTRAD = 3600., XDTRAD_CLONLY = 3600., NRAD_COLNBR = 400 /
+&NAM_PARAM_KAFRn XDTCONV = 300., NICE = 1 LREFRESH_ALL = T,LDOWN = T /
+&NAM_LBCn  CLBCX = 2*"OPEN", CLBCY = 2*"OPEN" /
+&NAM_TURBn  CTURBLEN = "BL89", CTURBDIM = "1DIM", LSUBG_COND = F /
+&NAM_CH_MNHCn LUSECHEM = F /
+&NAM_CONF  CCONF = "START", NVERB=5, NMODEL = 2,
+            CEXP = "16JAN", CSEG = "12B18" ,
+            CSPLIT='BSPLITTING' /
+&NAM_DYN  XSEGLEN = 300., LCORIO = T, LNUMDIFU = T,
+           XALKTOP = 0.001, XALZBOT = 14500. /
+&NAM_NESTING NDAD(2) = 1, NDTRATIO(2) = 4, XWAY(2) = 2. /
+&NAM_FMOUT XFMOUT(1,1) = 300. /
+&NAM_DIAG_SURFn  /
+&NAM_ISBAn CSCOND="NP89", CALBEDO="DRY",CROUGH='NONE', 
+       CC1DRY='DEF', CSOILFRZ='DEF', CDIFSFCOND='DEF', CSNOWRES='DEF' /
+&NAM_SSOn CROUGH='Z01D' /
+&NAM_DIAG_ISBAn /
+&NAM_SEAFLUXn CSEA_ALB="UNIF" /
+&NAM_BLANK /
+

MY_RUN/KTEST/007_16janvier/008_run2/EXSEG1.nam.src_WENO

+&NAM_CONFZ
+ NZ_VERB=5
+ ! NZ_VERB=5 , NZ_PROC=0 , NB_PROCIO_R=1 , NB_PROCIO_W=8 
+/
+&NAM_LUNITn  CINIFILE = "16JAN_06_MNH",CINIFILEPGD="16JAN98_36km.neste1"  /
+&NAM_DYNn  XTSTEP = 150.,  
+           ! CPRESOPT = "ZRESI", 
+           NITR = 8,
+           LHORELAX_UVWTH = T, LHORELAX_RV = T, LVE_RELAX = T,
+           NRIMX = 5, NRIMY = 5, XRIMKMAX = 0.0083, XT4DIFU = 5000. /
+&NAM_ADVn  CUVW_ADV_SCHEME="WENO_K",CMET_ADV_SCHEME = "PPM_01", CSV_ADV_SCHEME = "PPM_01" /
+&NAM_PARAMn  CCLOUD = "KESS", CTURB = "TKEL", CRAD = "ECMW",
+             CDCONV = "KAFR",CSCONV="KAFR" /
+&NAM_PARAM_RADn  XDTRAD = 3600., XDTRAD_CLONLY = 3600., NRAD_COLNBR = 400 /
+&NAM_PARAM_KAFRn XDTCONV = 300., NICE = 1 LREFRESH_ALL = T,LDOWN = T /
+&NAM_LBCn  CLBCX = 2*"OPEN", CLBCY = 2*"OPEN" /
+&NAM_TURBn  CTURBLEN = "BL89", CTURBDIM = "1DIM", LSUBG_COND = F /
+&NAM_CH_MNHCn LUSECHEM = F /
+&NAM_CONF  CCONF = "START", NVERB=5, NMODEL = 2,
+            CEXP = "16JAN", CSEG = "12B18" ,
+            CSPLIT='BSPLITTING' /
+&NAM_DYN  XSEGLEN = 300., LCORIO = T, LNUMDIFU = F,
+           XALKTOP = 0.001, XALZBOT = 14500. /
+&NAM_NESTING NDAD(2) = 1, NDTRATIO(2) = 3, XWAY(2) = 2. /
+&NAM_FMOUT XFMOUT(1,1) = 300. /
+&NAM_DIAG_SURFn  /
+&NAM_ISBAn CSCOND="NP89", CALBEDO="DRY",CROUGH='NONE', 
+       CC1DRY='DEF', CSOILFRZ='DEF', CDIFSFCOND='DEF', CSNOWRES='DEF' /
+&NAM_SSOn CROUGH='Z01D' /
+&NAM_DIAG_ISBAn /
+&NAM_SEAFLUXn CSEA_ALB="UNIF" /
+&NAM_BLANK /
+

MY_RUN/KTEST/007_16janvier/008_run2/EXSEG2.nam.src

@@ -3,7 +3,7 @@
            LHORELAX_UVWTH = F, LHORELAX_RV = F, LHORELAX_RC= F, LHORELAX_RR= F,
            LHORELAX_RS= F, LHORELAX_RI= F, LHORELAX_RG= F, LHORELAX_TKE= F,
            NRIMX = 0, NRIMY = 0, LVE_RELAX = T, XT4DIFU = 1500. /
-&NAM_ADVn CMET_ADV_SCHEME = "FCT2ND", CSV_ADV_SCHEME = "FCT2ND" /
+&NAM_ADVn CUVW_ADV_SCHEME="CEN4TH",CMET_ADV_SCHEME = "PPM_01", CSV_ADV_SCHEME = "PPM_01" /
 &NAM_PARAMn  CCLOUD = "KESS", CTURB = "TKEL", CRAD = "ECMW", CDCONV = "KAFR" /
 &NAM_PARAM_RADn  XDTRAD = 1800., XDTRAD_CLONLY = 1800., LCLEAR_SKY = F,
              NRAD_COLNBR = 400 /

MY_RUN/KTEST/007_16janvier/008_run2/EXSEG2.nam.src_CEN4TH

+&NAM_LUNITn  CINIFILE = "16JAN_06_MNH2" ,CINIFILEPGD="16JAN98_9km.neste1" /
+&NAM_DYNn CPRESOPT = "RICHA", NITR = 8, XRELAX = 1.,
+           LHORELAX_UVWTH = F, LHORELAX_RV = F, LHORELAX_RC= F, LHORELAX_RR= F,
+           LHORELAX_RS= F, LHORELAX_RI= F, LHORELAX_RG= F, LHORELAX_TKE= F,
+           NRIMX = 0, NRIMY = 0, LVE_RELAX = T, XT4DIFU = 1500. /
+&NAM_ADVn CUVW_ADV_SCHEME="CEN4TH",CMET_ADV_SCHEME = "PPM_01", CSV_ADV_SCHEME = "PPM_01" /
+&NAM_PARAMn  CCLOUD = "KESS", CTURB = "TKEL", CRAD = "ECMW", CDCONV = "KAFR" /
+&NAM_PARAM_RADn  XDTRAD = 1800., XDTRAD_CLONLY = 1800., LCLEAR_SKY = F,
+             NRAD_COLNBR = 400 /
+&NAM_PARAM_KAFRn XDTCONV = 300., NICE = 1, LREFRESH_ALL = T,
+              LDOWN = T /
+&NAM_LBCn  CLBCX = 2*"OPEN", CLBCY = 2*"OPEN", XCPHASE = 20. /
+&NAM_TURBn  XIMPL = 1., CTURBLEN = "BL89", CTURBDIM = "1DIM" ,
+            LTURB_FLX = F, LTURB_DIAG = F, LSUBG_COND = F /
+&NAM_CH_MNHCn LUSECHEM = F /
+&NAM_DIAG_SURFn  /
+&NAM_SSOn CROUGH='Z01D' /
+&NAM_ISBAn CSCOND  = "NP89", CALBEDO = "DRY",CROUGH='NONE',
+       CC1DRY = 'DEF', CSOILFRZ = 'DEF', CDIFSFCOND = 'DEF', CSNOWRES = 'DEF'  /
+&NAM_DIAG_ISBAn /
+&NAM_SEAFLUXn CSEA_ALB="UNIF" /
+&NAM_BLANK /
+

MY_RUN/KTEST/007_16janvier/008_run2/EXSEG2.nam.src_WENO

+&NAM_LUNITn  CINIFILE = "16JAN_06_MNH2" ,CINIFILEPGD="16JAN98_9km.neste1" /
+&NAM_DYNn CPRESOPT = "RICHA", NITR = 8, XRELAX = 1.,
+           LHORELAX_UVWTH = F, LHORELAX_RV = F, LHORELAX_RC= F, LHORELAX_RR= F,
+           LHORELAX_RS= F, LHORELAX_RI= F, LHORELAX_RG= F, LHORELAX_TKE= F,
+           NRIMX = 0, NRIMY = 0, LVE_RELAX = T, XT4DIFU = 1500. /
+&NAM_ADVn CUVW_ADV_SCHEME="WENO_K",CMET_ADV_SCHEME = "PPM_01", CSV_ADV_SCHEME = "PPM_01" /
+&NAM_PARAMn  CCLOUD = "KESS", CTURB = "TKEL", CRAD = "ECMW", CDCONV = "KAFR" /
+&NAM_PARAM_RADn  XDTRAD = 1800., XDTRAD_CLONLY = 1800., LCLEAR_SKY = F,
+             NRAD_COLNBR = 400 /
+&NAM_PARAM_KAFRn XDTCONV = 300., NICE = 1, LREFRESH_ALL = T,
+              LDOWN = T /
+&NAM_LBCn  CLBCX = 2*"OPEN", CLBCY = 2*"OPEN", XCPHASE = 20. /
+&NAM_TURBn  XIMPL = 1., CTURBLEN = "BL89", CTURBDIM = "1DIM" ,
+            LTURB_FLX = F, LTURB_DIAG = F, LSUBG_COND = F /
+&NAM_CH_MNHCn LUSECHEM = F /
+&NAM_DIAG_SURFn  /
+&NAM_SSOn CROUGH='Z01D' /
+&NAM_ISBAn CSCOND  = "NP89", CALBEDO = "DRY",CROUGH='NONE',
+       CC1DRY = 'DEF', CSOILFRZ = 'DEF', CDIFSFCOND = 'DEF', CSNOWRES = 'DEF'  /
+&NAM_DIAG_ISBAn /
+&NAM_SEAFLUXn CSEA_ALB="UNIF" /
+&NAM_BLANK /
+

MY_RUN/KTEST/007_16janvier/011_diaprog/dir_16janv

@@ -10,15 +10,15 @@ ZS
 LCOLAREA=T LISO=F                                                               
 LCOLZERO=T NCOLZERO=1                                                           
 LCOLZERO=F
-NIMNMX=0 XDIAINT_THM=1.                                                         
-THM_PR_850_ON_UMVM_PR_850                                                       
+NIMNMX=0 XDIAINT_THT=1.                                                         
+THT_PR_850_ON_UTVT_PR_850                                                       
 NIMNMX=-1                                                                       
-RVM(*1000)_PR_700                                                               
-PABSM_Z_9000_ON_UMVM_Z_9000                                                     
+RVT(*1000)_PR_700                                                               
+PABST_Z_9000_ON_UTVT_Z_9000                                                     
 _file2_
-NIMNMX=0 XDIAINT_THM=1.                                                         
-THM_PR_850_ON_UMVM_PR_850                                                       
+NIMNMX=0 XDIAINT_THT=1.                                                         
+THT_PR_850_ON_UTVT_PR_850                                                       
 NIMNMX=-1                                                                       
-RVM(*1000)_PR_700                                                               
-PABSM_Z_9000_ON_UMVM_Z_9000                                                     
+RVT(*1000)_PR_700                                                               
+PABST_Z_9000_ON_UTVT_Z_9000                                                     
 quit                           

MY_RUN/KTEST/007_16janvier/011_ncl_extractdia/dir_exnetcdf

@@ -14,16 +14,16 @@ KCDL
   0
   0
 CONF
- UM
- VM
- WM
+ UT
+ VT
+ WT
 ALT
-THM
-PABSM
-RVM
+THT
+PABST
+RVT
 MRC
 MRR
-RIM
+RIT
 ACPRR
 ZSBIS
 END

MY_RUN/KTEST/007_16janvier/011_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"

MY_RUN/KTEST/007_16janvier/011_ncl_extractdia/mnh_CrossSection1.ncl

@@ -36,7 +36,7 @@ begin
 ;=================================================;
 ; Get informations on variable sizes
 ;=================================================;
-  mdims = getfilevardimsizes(a,"UM") ; get some dimension sizes for the file
+  mdims = getfilevardimsizes(a,"UT") ; get some dimension sizes for the file
   nd = dimsizes(mdims)
 
 ;=================================================;
@@ -44,11 +44,11 @@ begin
 ;=================================================;
   ZS =a->ZSBIS(0,:,:)
 
-  th = a->THM(0,:,:,:)
+  th = a->THT(0,:,:,:)
   th@long_name="Potential Temperature"
   th@description="TH"
 
-  um = a->UM(0,:,:,:)
+  um = a->UT(0,:,:,:)
   um@description="U"
 
   z= a->ALT(0,:,:,:)

MY_RUN/KTEST/007_16janvier/011_ncl_extractdia/mnh_ModelLevels1.ncl

@@ -24,18 +24,18 @@ begin
 ;=================================================;
 ; Read the variables we need
 ;=================================================;
-  th  = a->THM(0,:,:,:)   ; theta
+  th  = a->THT(0,:,:,:)   ; theta
   th@description="Potential Temperature"
 
-  qv  = a->RVM(0,:,:,:)   ; Qv
+  qv  = a->RVT(0,:,:,:)   ; Qv
   qv = qv*1000.
   qv@units = "g/kg"
   qv@description="Water vapor mixing ratio"
 
-  u   = a->UM(0,:,:,:)        ; u 
+  u   = a->UT(0,:,:,:)        ; u 
   u@description="U"
 
-  v   = a->VM(0,:,:,:)        ; v 
+  v   = a->VT(0,:,:,:)        ; v 
   v@description="V"
 
 ; Calculate wind speed

MY_RUN/KTEST/007_16janvier/011_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