adaptation WENO

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)

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"
 

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