From 7b64bb25a3525bc9474368b436c5c328115afe63 Mon Sep 17 00:00:00 2001
From: Juan Escobar <escj@aero.obs-mip.fr>
Date: Thu, 4 Feb 2021 11:06:52 +0100
Subject: [PATCH] Juan 4/02/2021, For real4 execution , correction of some ncl
 script using "double" variables

---
 MY_RUN/KTEST/003_KW78/005_ncl/plot_KW78.ncl       | 8 +++++---
 MY_RUN/KTEST/004_Reunion/006_ncl/plot_Reunion.ncl | 6 ++++--
 MY_RUN/KTEST/007_16janvier/011_ncl/plot_16j_2.ncl | 5 +++--
 3 files changed, 12 insertions(+), 7 deletions(-)

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 d721d9aae..a25a8a407 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 2e7e9813a..f13d0b901 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 ffef5341c..0fedb8530 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
-- 
GitLab