Juan 4/02/2021, For real4 execution , correction of some ncl script using "double" variables

(cherry picked from commit 7b64bb25)

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

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

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