Juan 20/04/2018: ncl script corrections

MY_RUN/KTEST/001_2Drelief/005_ncl/plot_2Drelief.ncl

@@ -31,27 +31,27 @@ zs  = fichier1->ZS(0,jphext:imax+jphext-1) ; ZS
 zhat  = fichier1->ZHAT(1:kmax+1) ; ZHAT
 xhat  = fichier1->XHAT(jphext:imax+jphext-1+1) ; XHAT
 
-lsthm1 = fichier1->LSTHM(1:kmax,0,jphext:imax+jphext-1) ; LSTHM
+lsthm1 = fichier1->LSTHM(0,1:kmax,0,jphext:imax+jphext-1) ; LSTHM
 lsthm1@long_name="LSTHM"
 lsthm1@units="K"
 
-lsum1_old = fichier1->LSUM(1:kmax,0,jphext:imax+jphext-1+1) ; LSUM
-wt1_old= fichier1->WT(1:kmax+1,0,jphext:imax+jphext-1) ; WT
-ut1_old= fichier1->UT(1:kmax,0,jphext:imax+jphext-1+1) ; UT
+lsum1_old = fichier1->LSUM(0,1:kmax,0,jphext:imax+jphext-1+1) ; LSUM
+wt1_old= fichier1->WT(0,1:kmax+1,0,jphext:imax+jphext-1) ; WT
+ut1_old= fichier1->UT(0,1:kmax,0,jphext:imax+jphext-1+1) ; UT
 
-tht1= fichier1->THT(1:kmax,0,jphext:imax+jphext-1) ; THT
+tht1= fichier1->THT(0,1:kmax,0,jphext:imax+jphext-1) ; THT
 tht1@long_name="Potential Temperature"
 tht1@units="K"
 
-rvt1= fichier1->RVT(1:kmax,0,jphext:imax+jphext-1) ; RVT
+rvt1= fichier1->RVT(0,1:kmax,0,jphext:imax+jphext-1) ; RVT
 rvt1@long_name="Vapor mixing ratio"
 rvt1@units="g/kg"
 
-cflu1= fichier1->CFLU(1:kmax,0,jphext:imax+jphext-1) ; CFLU
+cflu1= fichier1->CFLU(0,1:kmax,0,jphext:imax+jphext-1) ; CFLU
 cflu1@long_name="CFLU"
 cflu1@units="-"
 
-cflw1= fichier1->CFLW(1:kmax,0,jphext:imax+jphext-1) ; CFLW
+cflw1= fichier1->CFLW(0,1:kmax,0,jphext:imax+jphext-1) ; CFLW
 cflw1@long_name="CFLW"
 cflw1@units=" "
 
@@ -60,27 +60,27 @@ SS1=fichier1->DTCUR
 ;YYYYDDMM2= fichier2->DTCUR__TDATE
 SS2=fichier2->DTCUR
 
-lsthm2 = fichier2->LSTHM(1:kmax,0,jphext:imax+jphext-1) ; LSTHM
+lsthm2 = fichier2->LSTHM(0,1:kmax,0,jphext:imax+jphext-1) ; LSTHM
 lsthm2@long_name="LSTHM"
 lsthm2@units="K"
 
-lsum2_old = fichier2->LSUM(1:kmax,0,jphext:imax+jphext-1+1) ; LSUM
-wt2_old= fichier2->WT(1:kmax+1,0,jphext:imax+jphext-1) ; WT
-ut2_old= fichier2->UT(1:kmax,0,jphext:imax+jphext-1+1) ; UT
+lsum2_old = fichier2->LSUM(0,1:kmax,0,jphext:imax+jphext-1+1) ; LSUM
+wt2_old= fichier2->WT(0,1:kmax+1,0,jphext:imax+jphext-1) ; WT
+ut2_old= fichier2->UT(0,1:kmax,0,jphext:imax+jphext-1+1) ; UT
 
-tht2= fichier2->THT(1:kmax,0,jphext:imax+jphext-1) ; THT
+tht2= fichier2->THT(0,1:kmax,0,jphext:imax+jphext-1) ; THT
 tht2@long_name="Potential Temperature"
 tht2@units="K"
 
-rvt2= fichier2->RVT(1:kmax,0,jphext:imax+jphext-1) ; RVT
+rvt2= fichier2->RVT(0,1:kmax,0,jphext:imax+jphext-1) ; RVT
 rvt2@long_name="Vapor mixing ratio"
 rvt2@units="g/kg"
 
-cflu2= fichier2->CFLU(1:kmax,0,jphext:imax+jphext-1) ; CFLU
+cflu2= fichier2->CFLU(0,1:kmax,0,jphext:imax+jphext-1) ; CFLU
 cflu2@long_name="CFLU"
 cflu2@units="-"
 
-cflw2= fichier2->CFLW(1:kmax,0,jphext:imax+jphext-1) ; CFLW
+cflw2= fichier2->CFLW(0,1:kmax,0,jphext:imax+jphext-1) ; CFLW
 cflw2@long_name="CFLW"
 cflw2@units=" "
 

MY_RUN/KTEST/002_3Drelief/005_ncl/plot_3Drelief.ncl

@@ -27,7 +27,7 @@
 
 
 
-  ut= fichier1->UT(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1+1) ; UT
+  ut= fichier1->UT(0,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; UT
   ut@long_name="Zonal wind"
   ut@units="m/s"
 

MY_RUN/KTEST/003_KW78/005_ncl/plot_KW78.ncl

@@ -35,33 +35,33 @@
   yhat  = fichier1->YHAT(1:jmax+jphext-1+1) ; YHAT
   xhat  = fichier1->XHAT(1:imax+jphext-1+1) ; XHAT
 
-  wt= fichier1->WT(1:kmax+1,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; WT
-  vt= fichier1->VT(1:kmax,jphext:jmax+jphext-1+1,jphext:imax+jphext-1) ; VT
-  ut= fichier1->UT(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1+1) ; UT
-  tht1= fichier1->THT(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; THT
+  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
+  tht1= fichier1->THT(0,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; THT
   tht1@long_name="Potential Temperature"
   tht1@units="K"
-  lsthm1= fichier1->LSTHM(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1) 
+  lsthm1= fichier1->LSTHM(0,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1) 
   lsthm1@long_name="Large SCale Potential Temperature"
   lsthm1@units="K"
 
-  pabst1= fichier1->PABST(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; PABST
+  pabst1= fichier1->PABST(0,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; PABST
   pabst1@long_name="Pressure"
   pabst1@units="Pa"
-  rvt1= fichier1->RVT(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; RVT
+  rvt1= fichier1->RVT(0,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; RVT
   rvt1@long_name="Vapor mixing ratio"
   rvt1@units="kg/kg"
-  rct1= fichier1->RCT(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; RCT
+  rct1= fichier1->RCT(0,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; RCT
   rct1@long_name="Cloud mixing ratio"
   rct1@units="kg/kg"
-  rrt1= fichier1->RRT(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; RRT
+  rrt1= fichier1->RRT(0,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; RRT
   rrt1@long_name="Rain mixing ratio"
   rrt1@units="kg/kg"
 
-  inprr= fichier1->INPRR(jphext:jmax+jphext-1,jphext:imax+jphext-1) 
+  inprr= fichier1->INPRR(0,jphext:jmax+jphext-1,jphext:imax+jphext-1) 
   inprr@long_name="Instantaneous Precipitation Rain Rate"
   inprr@units="mm/h"
-  acprr= fichier1->ACPRR(jphext:jmax+jphext-1,jphext:imax+jphext-1) 
+  acprr= fichier1->ACPRR(0,jphext:jmax+jphext-1,jphext:imax+jphext-1) 
   acprr@long_name="Accumulated Precipitation Rain Rate"
   acprr@units="mm"
 

MY_RUN/KTEST/003_KW78/005_ncl/run_ncl

@@ -6,6 +6,8 @@ set -x
 
 ln -sf ../002_mesonh/KWRAI.1.SEG01.004.nc .
 
+rm -f visu_KW78.*.png
+
 ncl  plot_KW78.ncl
 display visu_KW78.000001.png
 display visu_KW78.000002.png

MY_RUN/KTEST/003_KW78/007_ncl_diag/plot_KW78_diag.ncl

@@ -36,23 +36,23 @@
   yhat  = fichier1->YHAT(1:jmax+jphext-1+1) ; YHAT
   xhat  = fichier1->XHAT(1:imax+jphext-1+1) ; XHAT
 
-  wt= fichier1->WT(1:kmax+1,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; WT
-  vt= fichier1->VT(1:kmax,jphext:jmax+jphext-1+1,jphext:imax+jphext-1) ; VT
-  ut= fichier1->UT(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1+1) ; UT
+  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
   
-  idt1= fichier2->CLUSTERID(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; CLUSTERID
+  idt1= fichier2->CLUSTERID(0,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; CLUSTERID
   idt1@long_name="Cluster Identity Number"
   idt1@units=""
 
-  YYYYDDMM1= fichier1->DTCUR__TDATE
-  SS1=fichier1->DTCUR__TIME
+  YYYYDDMM1= fichier1->DTCUR
+  SS1=fichier1->DTCUR
 
 ;=================================================;
 ; Récupération de la date 
 ;=================================================;
-  year1=YYYYDDMM1(0)
-  day1=YYYYDDMM1(1)
-  min1=YYYYDDMM1(2)
+  ;year1=YYYYDDMM1(0)
+  ;day1=YYYYDDMM1(1)
+  ;min1=YYYYDDMM1(2)
   time1=SS1
 ;=================================================;
 ; On mets toutes les variables sur la grille 1 

MY_RUN/KTEST/003_KW78/007_ncl_diag/run_ncl

@@ -7,6 +7,8 @@ set -x
 ln -sf ../002_mesonh/KWRAI.1.SEG01.004.nc .
 ln -sf ../006_diag/KWRAI.1.SEG01.004dia.nc .
 
+rm -f visu_KW78_diag.png
+
 ncl  plot_KW78_diag.ncl
 display visu_KW78_diag.png