diff --git a/MY_RUN/KTEST/003_KW78/005_ncl/plot_KW78.ncl b/MY_RUN/KTEST/003_KW78/005_ncl/plot_KW78.ncl
index d721d9aae79491ac8197f610e61909a7f1a7404c..a25a8a407774a675ff16edb746d49d0c8c28ae7b 100644
--- a/MY_RUN/KTEST/003_KW78/005_ncl/plot_KW78.ncl
+++ b/MY_RUN/KTEST/003_KW78/005_ncl/plot_KW78.ncl
@@ -35,6 +35,8 @@
   yhat  = fichier1->YHAT(1:jmax+jphext-1+1) ; YHAT
   xhat  = fichier1->XHAT(1:imax+jphext-1+1) ; XHAT
 
+  vtype = typeof(zs)
+
   wt= fichier1->WT(0,1:kmax+1,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; WT
   vt= fichier1->VT(0,1:kmax,jphext:jmax+jphext-1+1,jphext:imax+jphext-1) ; VT
   ut= fichier1->UT(0,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1+1) ; UT
@@ -106,14 +108,14 @@
 ; Altitude des niveaux modèles
 ;=================================================;
 ; Unstagger zhat (from grid 4 to 1)
-  nzh=new(kmax,double)
+  nzh=new(kmax,vtype)
   do k=0,kmax-1
     nzh(k)=(zhat(k)+zhat(k+1))/2.
   end do
 
 ; Create alt (altitude des niveaux modèle)
-  alt=new(dimsizes(tht1),double)
-  zcoef=new((/jmax,imax/),double)
+  alt=new(dimsizes(tht1),vtype)
+  zcoef=new((/jmax,imax/),vtype)
   zcoef=1.-zs/nzh(kmax-1)
   do i=0,imax-1
     do j=0,jmax-1
diff --git a/MY_RUN/KTEST/004_Reunion/006_ncl/plot_Reunion.ncl b/MY_RUN/KTEST/004_Reunion/006_ncl/plot_Reunion.ncl
index 2e7e9813a3efc76dc8f6fffe7afb56b4b7793a68..f13d0b901dfa169c109a999f0d0beb14e2cc1e49 100644
--- a/MY_RUN/KTEST/004_Reunion/006_ncl/plot_Reunion.ncl
+++ b/MY_RUN/KTEST/004_Reunion/006_ncl/plot_Reunion.ncl
@@ -102,6 +102,8 @@ load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
   zs@lat2d = lat2d
   zs@lon2d = lon2d
 
+  vtype = typeof(zs)
+
   wt= fichier1->WT(0,1:kmax+1,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; WT
   vt= fichier1->VT(0,1:kmax,jphext:jmax+jphext-1+1,jphext:imax+jphext-1); VT
   ut= fichier1->UT(0,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1+1) ; UT
@@ -156,13 +158,13 @@ load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
   zhat= fichier1->ZHAT(1:kmax+1)
 
 ; Unstagger zhat (from grid 4 to 1)
-  nzhat=new(kmax,double)
+  nzhat=new(kmax,vtype)
   do k=0,kmax-1
     nzhat(k)=(zhat(k)+zhat(k+1))/2.
   end do
 
 ; Create Z3D == ALT
-  alt=new(dimsizes(tht),double)
+  alt=new(dimsizes(tht),vtype)
   zcoef=1.-zs/nzhat(kmax-1)
 
   do i=0,imax-1
diff --git a/MY_RUN/KTEST/007_16janvier/011_ncl/plot_16j_2.ncl b/MY_RUN/KTEST/007_16janvier/011_ncl/plot_16j_2.ncl
index ffef5341c8c8b704b61e6b69daa7fceb3f501b47..0fedb8530ca39ec1886b44d8ea159efc98bca4c6 100644
--- a/MY_RUN/KTEST/007_16janvier/011_ncl/plot_16j_2.ncl
+++ b/MY_RUN/KTEST/007_16janvier/011_ncl/plot_16j_2.ncl
@@ -108,6 +108,7 @@ zs@units="m"
 zs@lat2d = lat2d
 zs@lon2d = lon2d
 
+vtype = typeof(zs)
 
 wt= fichier1->WT(0,1:kmax+1,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; WT
 vt= fichier1->VT(0,1:kmax,jphext:jmax+jphext-1+1,jphext:imax+jphext-1) ; VT
@@ -160,13 +161,13 @@ wt1@lon2d = lon2d
 zhat= fichier1->ZHAT(1:kmax+1)
 
 ; Unstagger zhat (from grid 4 to 1)
-    nzhat=new(kmax,double)
+    nzhat=new(kmax,vtype)
     do k=0,kmax-1
      nzhat(k)=(zhat(k)+zhat(k+1))/2.
     end do
 
 ; Create Z3D == ALT
-    alt=new(dimsizes(tht1),double)
+    alt=new(dimsizes(tht1),vtype)
     zcoef=1.-zs/nzhat(kmax-1)
 
     do i=0,imax-1