From 821bd27a2d50fdd3be55784d7e5680110347f5e9 Mon Sep 17 00:00:00 2001 From: Philippe WAUTELET <philippe.wautelet@aero.obs-mip.fr> Date: Fri, 20 Apr 2018 15:01:59 +0200 Subject: [PATCH] Juan 20/04/2018: ncl script corrections --- .../001_2Drelief/005_ncl/plot_2Drelief.ncl | 32 +++++++++---------- .../002_3Drelief/005_ncl/plot_3Drelief.ncl | 2 +- MY_RUN/KTEST/003_KW78/005_ncl/plot_KW78.ncl | 22 ++++++------- MY_RUN/KTEST/003_KW78/005_ncl/run_ncl | 2 ++ .../003_KW78/007_ncl_diag/plot_KW78_diag.ncl | 18 +++++------ MY_RUN/KTEST/003_KW78/007_ncl_diag/run_ncl | 2 ++ 6 files changed, 41 insertions(+), 37 deletions(-) diff --git a/MY_RUN/KTEST/001_2Drelief/005_ncl/plot_2Drelief.ncl b/MY_RUN/KTEST/001_2Drelief/005_ncl/plot_2Drelief.ncl index 3a5857c57..41eb6d0ea 100644 --- a/MY_RUN/KTEST/001_2Drelief/005_ncl/plot_2Drelief.ncl +++ b/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=" " diff --git a/MY_RUN/KTEST/002_3Drelief/005_ncl/plot_3Drelief.ncl b/MY_RUN/KTEST/002_3Drelief/005_ncl/plot_3Drelief.ncl index 25e4449e3..670655218 100644 --- a/MY_RUN/KTEST/002_3Drelief/005_ncl/plot_3Drelief.ncl +++ b/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" 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 149b786d2..d721d9aae 100644 --- a/MY_RUN/KTEST/003_KW78/005_ncl/plot_KW78.ncl +++ b/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" diff --git a/MY_RUN/KTEST/003_KW78/005_ncl/run_ncl b/MY_RUN/KTEST/003_KW78/005_ncl/run_ncl index bb885331e..26bc2d373 100755 --- a/MY_RUN/KTEST/003_KW78/005_ncl/run_ncl +++ b/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 diff --git a/MY_RUN/KTEST/003_KW78/007_ncl_diag/plot_KW78_diag.ncl b/MY_RUN/KTEST/003_KW78/007_ncl_diag/plot_KW78_diag.ncl index 731f3e8d5..734836971 100755 --- a/MY_RUN/KTEST/003_KW78/007_ncl_diag/plot_KW78_diag.ncl +++ b/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 diff --git a/MY_RUN/KTEST/003_KW78/007_ncl_diag/run_ncl b/MY_RUN/KTEST/003_KW78/007_ncl_diag/run_ncl index 4f9f88b34..ec3585d71 100755 --- a/MY_RUN/KTEST/003_KW78/007_ncl_diag/run_ncl +++ b/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 -- GitLab