adaptation WENO

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"

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
 ;=================================================;
 

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"
 

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