diff --git a/MY_RUN/KTEST/007_16janvier/011_ncl_lfi2cdf/mnh_Cloud2.ncl b/MY_RUN/KTEST/007_16janvier/011_ncl_lfi2cdf/mnh_Cloud2.ncl
index a26f5164aca6277d75e11ffca1ad68b674a775c8..a46a6512923ecfdd69d73f1bcdd7a2898c6f1cb8 100644
--- a/MY_RUN/KTEST/007_16janvier/011_ncl_lfi2cdf/mnh_Cloud2.ncl
+++ b/MY_RUN/KTEST/007_16janvier/011_ncl_lfi2cdf/mnh_Cloud2.ncl
@@ -28,7 +28,7 @@ begin
; Get informations on variable sizes
;=================================================;
- mdims = getfilevardimsizes(a,"RVM") ; get some dimension sizes for the file
+ mdims = getfilevardimsizes(a,"RVT") ; get some dimension sizes for the file
nd = dimsizes(mdims)
imax=mdims(nd-1)-2
jmax=mdims(nd-2)-2
@@ -41,8 +41,8 @@ begin
; Read water wapor mixing ratio if it exists
; -----------------------------------------
- if(isfilevar(a,"RVM"))
- qv = a->RVM(1:kmax,1:jmax,1:imax)
+ if(isfilevar(a,"RVT"))
+ qv = a->RVT(1:kmax,1:jmax,1:imax)
qv = qv*1000.
qv@units = "g/kg"
qv@description="Water vapor mixing ratio"
@@ -50,8 +50,8 @@ begin
; Read cloud mixing ratio if it exists
; -----------------------------------------
- if(isfilevar(a,"RCM"))
- qc = a->RCM(1:kmax,1:jmax,1:imax)
+ if(isfilevar(a,"RCT"))
+ qc = a->RCT(1:kmax,1:jmax,1:imax)
qc = qc*1000.
qc@units = "g/kg"
qc@description="Cloud mixing ratio"
@@ -59,8 +59,8 @@ begin
; Read rain mixing ratio if it exists
; -----------------------------------------
- if(isfilevar(a,"RRM"))
- qr = a->RRM(1:kmax,1:jmax,1:imax)
+ if(isfilevar(a,"RRT"))
+ qr = a->RRT(1:kmax,1:jmax,1:imax)
qr = qr*1000.
qr@units = "g/kg"
qr@description="Rain mixing ratio"
@@ -68,8 +68,8 @@ begin
; Read ice mixing ratio if it exists
; -----------------------------------------
- if(isfilevar(a,"RIM"))
- qi = a->RIM(1:kmax,1:jmax,1:imax)
+ if(isfilevar(a,"RIT"))
+ qi = a->RIT(1:kmax,1:jmax,1:imax)
qi = qi*1000.
qi@units = "g/kg"
qi@description="Ice mixing ratio"
diff --git a/MY_RUN/KTEST/007_16janvier/011_ncl_lfi2cdf/mnh_CrossSection2.ncl b/MY_RUN/KTEST/007_16janvier/011_ncl_lfi2cdf/mnh_CrossSection2.ncl
index b0198978d6d7e68014b17d0a1defce67f4558089..86f3062fd91e737d8844d41692377d28c826e72c 100644
--- a/MY_RUN/KTEST/007_16janvier/011_ncl_lfi2cdf/mnh_CrossSection2.ncl
+++ b/MY_RUN/KTEST/007_16janvier/011_ncl_lfi2cdf/mnh_CrossSection2.ncl
@@ -28,7 +28,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)
imax=mdims(nd-1)-2
jmax=mdims(nd-2)-2
@@ -39,11 +39,11 @@ begin
;=================================================;
ZS =a->ZS(1:jmax,1:imax)
- th = a->THM(1:kmax,1:jmax,1:imax)
+ th = a->THT(1:kmax,1:jmax,1:imax)
th@long_name="Potential Temperature"
th@description="TH"
- um = a->UM(1:kmax,1:jmax,1:imax+1)
+ um = a->UT(1:kmax,1:jmax,1:imax+1)
um@description="U"
zh= a->ZHAT(1:kmax+1)
@@ -86,7 +86,7 @@ begin
nz = floattoint(zmax + 1)
;=================================================;
-; Use MESONHtools procedure to set map projection
+; Use TESONHtools procedure to set map projection
; and return latitude and longitude values
;=================================================;
diff --git a/MY_RUN/KTEST/007_16janvier/011_ncl_lfi2cdf/mnh_ModelLevels2.ncl b/MY_RUN/KTEST/007_16janvier/011_ncl_lfi2cdf/mnh_ModelLevels2.ncl
index 5cc8d6dd0dd077f0f9992b5745f4343b8188fec2..110fd69a806d8203cc082b6a466e55ad45758411 100644
--- a/MY_RUN/KTEST/007_16janvier/011_ncl_lfi2cdf/mnh_ModelLevels2.ncl
+++ b/MY_RUN/KTEST/007_16janvier/011_ncl_lfi2cdf/mnh_ModelLevels2.ncl
@@ -29,7 +29,7 @@ begin
; Get informations on variable sizes
; dims are dims-2 to remove non-physical values
;=================================================;
- mdims = getfilevardimsizes(a,"THM") ; get dimension sizes
+ mdims = getfilevardimsizes(a,"THT") ; get dimension sizes
nd = dimsizes(mdims)
imax=mdims(nd-1)-2
jmax=mdims(nd-2)-2
@@ -38,20 +38,20 @@ begin
;=================================================;
; Read the variables we need
;=================================================;
- th = a->THM(1:kmax,1:jmax,1:imax) ; theta
+ th = a->THT(1:kmax,1:jmax,1:imax) ; theta
th@long_name="Potential Temperature"
th@units = "K"
- qv = a->RVM(1:kmax,1:jmax,1:imax) ; Qv
+ qv = a->RVT(1:kmax,1:jmax,1:imax) ; Qv
qv = qv*1000.
qv@long_name="Water vapor mixing ratio"
qv@units = "g/kg"
- um = a->UM(1:kmax,1:jmax,1:imax+1) ; u
+ um = a->UT(1:kmax,1:jmax,1:imax+1) ; u
um@long_name="U"
um@units="m/s"
- vm = a->VM(1:kmax,1:jmax+1,1:imax) ; v
+ vm = a->VT(1:kmax,1:jmax+1,1:imax) ; v
vm@long_name="V"
vm@units="m/s"
diff --git a/MY_RUN/KTEST/007_16janvier/011_ncl_lfi2cdf/mnh_PressureLevel2.ncl b/MY_RUN/KTEST/007_16janvier/011_ncl_lfi2cdf/mnh_PressureLevel2.ncl
index 65688047564656f43a5c19fdc801c17b4e1b0af7..f4eb2129ac56fc8e97e3f84ba246f1156b7d99e7 100644
--- a/MY_RUN/KTEST/007_16janvier/011_ncl_lfi2cdf/mnh_PressureLevel2.ncl
+++ b/MY_RUN/KTEST/007_16janvier/011_ncl_lfi2cdf/mnh_PressureLevel2.ncl
@@ -37,7 +37,7 @@ begin
; Get informations on variable sizes
; dims are dims-2 to remove non-physical values
;=================================================;
- mdims = getfilevardimsizes(a,"UM") ; get dimension sizes
+ mdims = getfilevardimsizes(a,"UT") ; get dimension sizes
nd = dimsizes(mdims)
imax=mdims(nd-1)-2
jmax=mdims(nd-2)-2
@@ -47,11 +47,11 @@ begin
; Read the variables we need
;=================================================;
- tk = a->THM(1:kmax,1:jmax,1:imax) ; potential temperature
+ tk = a->THT(1:kmax,1:jmax,1:imax) ; potential temperature
tk@description="Potential Temperature"
- u = a->UM(1:kmax,1:jmax,1:imax+1) ; u
- v = a->VM(1:kmax,1:jmax+1,1:imax) ; v
+ u = a->UT(1:kmax,1:jmax,1:imax+1) ; u
+ v = a->VT(1:kmax,1:jmax+1,1:imax) ; v
; Unstagger u,v
u!0="Z"
@@ -64,10 +64,10 @@ begin
v!2="X"
v_unst = wrf_user_unstagger(v,"Y")
- p = a->PABSM(1:kmax,1:jmax,1:imax) ; pressure
+ p = a->PABST(1:kmax,1:jmax,1:imax) ; pressure
p@description="Pressure"
- rv = a->RVM(1:kmax,1:jmax,1:imax)*1000 ; water vapor mixing ratio
+ rv = a->RVT(1:kmax,1:jmax,1:imax)*1000 ; water vapor mixing ratio
rv@description="Water vapor mixing ratio"
zs = a->ZS(1:jmax,1:imax) ; terrain