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 3a5857c57e46ebb1ed25e66f549d48c757a701f7..41eb6d0eabd1d0c4bba10d705dac87a3bd68e58b 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 25e4449e367d7920fb68fd664210b85f57198531..6706552180c7485a15750f37f1e56294925f794e 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 149b786d2b9effb019d4f7f67b3d86ad6bb4b912..d721d9aae79491ac8197f610e61909a7f1a7404c 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 bb885331e5c759a3af72f0d4a0b54cfbc5ac45a3..26bc2d37386a86c4e79f245f1b755d66b8653bff 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 731f3e8d5d16b0b465770aa0c9fabab24169993b..73483697141487eaa1724dda02f50e1d52cecedf 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 4f9f88b34d84b74809865c5cfef064550a4ecef1..ec3585d714f6f3da190ab8572ddc1a0754f0222a 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