From 828c39a4d146dcc6a2f7a4557ae2d13caced0afa Mon Sep 17 00:00:00 2001
From: Christine Lac <christine.lac@meteo.fr>
Date: Fri, 13 Dec 2013 16:12:06 +0000
Subject: [PATCH] adaptation WENO
---
.../005_ncl_lfi2cdf/mnh_CrossSection2.ncl | 6 +++---
.../004_Reunion/005_ncl_lfi2cdf/mnh_ModelLevels2.ncl | 10 +++++-----
.../005_ncl_lfi2cdf/mnh_PressureLevel2.ncl | 12 ++++++------
3 files changed, 14 insertions(+), 14 deletions(-)
diff --git a/MY_RUN/KTEST/004_Reunion/005_ncl_lfi2cdf/mnh_CrossSection2.ncl b/MY_RUN/KTEST/004_Reunion/005_ncl_lfi2cdf/mnh_CrossSection2.ncl
index da775b04b..47b058931 100644
--- a/MY_RUN/KTEST/004_Reunion/005_ncl_lfi2cdf/mnh_CrossSection2.ncl
+++ b/MY_RUN/KTEST/004_Reunion/005_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)
diff --git a/MY_RUN/KTEST/004_Reunion/005_ncl_lfi2cdf/mnh_ModelLevels2.ncl b/MY_RUN/KTEST/004_Reunion/005_ncl_lfi2cdf/mnh_ModelLevels2.ncl
index 4d37194d4..cf09ccd0b 100644
--- a/MY_RUN/KTEST/004_Reunion/005_ncl_lfi2cdf/mnh_ModelLevels2.ncl
+++ b/MY_RUN/KTEST/004_Reunion/005_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/004_Reunion/005_ncl_lfi2cdf/mnh_PressureLevel2.ncl b/MY_RUN/KTEST/004_Reunion/005_ncl_lfi2cdf/mnh_PressureLevel2.ncl
index 033d55d55..ba6f30ad6 100644
--- a/MY_RUN/KTEST/004_Reunion/005_ncl_lfi2cdf/mnh_PressureLevel2.ncl
+++ b/MY_RUN/KTEST/004_Reunion/005_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
--
GitLab