From 19088a1e1ff97f3ba3d8d2515e37cef9c15767df Mon Sep 17 00:00:00 2001
From: Gaelle DELAUTIER <gaelle.delautier@meteo.fr>
Date: Mon, 23 Apr 2018 14:36:32 +0200
Subject: [PATCH] Gaelle 23/04/2018 : bug ncl script corrections
---
.../004_Reunion/006_ncl/plot_Reunion.ncl | 12 ++--
.../KTEST/007_16janvier/011_ncl/plot_16j.ncl | 12 ++--
.../007_16janvier/011_ncl/plot_16j_2.ncl | 12 ++--
.../KTEST/009_ICARTT/006_ncl/plot_ICARTT.ncl | 58 +++++++++----------
.../009_ICARTT/006_ncl/plot_ICARTT_budget.ncl | 19 +++---
.../011_KW78CHEM/005_ncl/plot_KW78_chem.ncl | 58 +++++++++----------
MY_RUN/KTEST/012_dust/007_ncl/plot_dust.ncl | 20 +++----
MY_RUN/KTEST/014_LIMA/003_ncl/plot_LIMA.ncl | 12 ++--
8 files changed, 103 insertions(+), 100 deletions(-)
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 747d0c3c4..3898b26b9 100644
--- a/MY_RUN/KTEST/004_Reunion/006_ncl/plot_Reunion.ncl
+++ b/MY_RUN/KTEST/004_Reunion/006_ncl/plot_Reunion.ncl
@@ -102,28 +102,28 @@ load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
zs@lat2d = lat2d
zs@lon2d = lon2d
- 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
tht= fichier1->THT(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; THT
tht@long_name="Potential Temperature"
tht@units="K"
tht@lat2d = lat2d
tht@lon2d = lon2d
- lsthm= fichier1->LSTHM(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1) ;LSTHM
+ lsthm= fichier1->LSTHM(0,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1) ;LSTHM
lsthm@long_name="Large SCale Potential Temperature"
lsthm@units="K"
lsthm@lat2d = lat2d
lsthm@lon2d = lon2d
- lsvm= fichier1->LSVM(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; LSVM
+ lsvm= fichier1->LSVM(0,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; LSVM
lsvm@long_name="Large SCale Merdian Wind"
lsvm@units="m/s"
lsvm@lat2d = lat2d
lsvm@lon2d = lon2d
- pabst= fichier1->PABST(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; PABST
+ pabst= fichier1->PABST(0,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; PABST
pabst@long_name="Pressure"
pabst@units="Pa"
pabst@lat2d = lat2d
diff --git a/MY_RUN/KTEST/007_16janvier/011_ncl/plot_16j.ncl b/MY_RUN/KTEST/007_16janvier/011_ncl/plot_16j.ncl
index cd16dc4f3..961b4a278 100644
--- a/MY_RUN/KTEST/007_16janvier/011_ncl/plot_16j.ncl
+++ b/MY_RUN/KTEST/007_16janvier/011_ncl/plot_16j.ncl
@@ -110,11 +110,11 @@ zs@lat2d = lat2d
zs@lon2d = lon2d
-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
-tht1= fichier1->THT(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; THT
+tht1= fichier1->THT(0,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; THT
tht1@long_name="Potential Temperature"
tht1@units="K"
tht1@lat2d = lat2d
@@ -122,13 +122,13 @@ tht1@lon2d = lon2d
-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"
pabst1@lat2d = lat2d
pabst1@lon2d = lon2d
-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=rvt1*1000
rvt1@long_name="Vapor mixing ratio"
rvt1@units="kg/kg"
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 5de8a4a62..ffef5341c 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
@@ -109,11 +109,11 @@ zs@lat2d = lat2d
zs@lon2d = lon2d
-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
-tht1= fichier1->THT(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; THT
+tht1= fichier1->THT(0,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; THT
tht1@long_name="Potential Temperature"
tht1@units="K"
tht1@lat2d = lat2d
@@ -121,13 +121,13 @@ tht1@lon2d = lon2d
-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"
pabst1@lat2d = lat2d
pabst1@lon2d = lon2d
-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=rvt1*1000
rvt1@long_name="Vapor mixing ratio"
rvt1@units="kg/kg"
diff --git a/MY_RUN/KTEST/009_ICARTT/006_ncl/plot_ICARTT.ncl b/MY_RUN/KTEST/009_ICARTT/006_ncl/plot_ICARTT.ncl
index 40e2405c6..0ae28c72e 100644
--- a/MY_RUN/KTEST/009_ICARTT/006_ncl/plot_ICARTT.ncl
+++ b/MY_RUN/KTEST/009_ICARTT/006_ncl/plot_ICARTT.ncl
@@ -38,7 +38,7 @@ zs@lon2d = lon2d
printMinMax(zs,0)
- rc_t1 = a->MRC(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
+ rc_t1 = a->MRC(:,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
rc_t1@long_name="Cloud mixing ratio"
rc_t1@units="g/kg"
rc_t1@lat2d=lat2d
@@ -46,21 +46,21 @@ printMinMax(zs,0)
printMinMax(rc_t1,0)
;
- o3_t1 = a->O3T(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
+ o3_t1 = a->O3T(:,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
o3_t1@long_name="Ozone"
o3_t1@units="ppbv"
o3_t1@lat2d=lat2d
o3_t1@lon2d=lon2d
;
- co_t1 = a->COT(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
+ co_t1 = a->COT(:,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
co_t1@long_name="carbon monoxide"
co_t1@units="ppbv"
co_t1@lat2d=lat2d
co_t1@lon2d=lon2d
;
- o3p_t1=a->O3_PROD(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
+ o3p_t1=a->O3_PROD(:,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
o3p_t1 = o3p_t1*1e9*3600
o3p_t1@long_name="ozone production"
o3p_t1@units="ppbv/h"
@@ -68,7 +68,7 @@ printMinMax(rc_t1,0)
o3p_t1@lon2d=lon2d
;
- o3l_t1=a->O3_LOSS(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
+ o3l_t1=a->O3_LOSS(:,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
o3l_t1 = o3l_t1*1e9*3600
o3l_t1@long_name="ozone loss"
o3l_t1@units="ppbv/h"
@@ -76,7 +76,7 @@ printMinMax(rc_t1,0)
o3l_t1@lon2d=lon2d
;
- cop_t1=a->CO_PROD(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
+ cop_t1=a->CO_PROD(:,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
cop_t1 = cop_t1*1e9*3600
cop_t1@long_name="carbon monoxide production"
cop_t1@units="ppbv/h"
@@ -84,7 +84,7 @@ printMinMax(rc_t1,0)
cop_t1@lon2d=lon2d
;
- col_t1=a->CO_LOSS(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
+ col_t1=a->CO_LOSS(:,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
col_t1 = col_t1*1e9*3600
col_t1@long_name="carbon monoxide loss"
col_t1@units="ppbv/h"
@@ -94,28 +94,28 @@ printMinMax(rc_t1,0)
;
;
- rc_t2 = a2->MRC(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
+ rc_t2 = a2->MRC(:,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
rc_t2@long_name="Cloud mixing ratio"
rc_t2@units="g/kg"
rc_t2@lat2d=lat2d
rc_t2@lon2d=lon2d
;
- o3_t2 = a2->O3T(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
+ o3_t2 = a2->O3T(:,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
o3_t2@long_name="Ozone"
o3_t2@units="ppbv"
o3_t2@lat2d=lat2d
o3_t2@lon2d=lon2d
;
- co_t2 = a2->COT(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
+ co_t2 = a2->COT(:,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
co_t2@long_name="carbon monoxide"
co_t2@units="ppbv"
co_t2@lat2d=lat2d
co_t2@lon2d=lon2d
;
- o3p_t2=a2->O3_PROD(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
+ o3p_t2=a2->O3_PROD(:,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
o3p_t2=o3p_t2*1e9*3600
o3p_t2@long_name="ozone production"
o3p_t2@units="pptv/h"
@@ -123,7 +123,7 @@ printMinMax(rc_t1,0)
o3p_t2@lon2d=lon2d
;
- o3l_t2=a2->O3_LOSS(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
+ o3l_t2=a2->O3_LOSS(:,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
o3l_t2=o3l_t2*1e9*3600
o3l_t2@long_name="ozone loss"
o3l_t2@units="pptv/h"
@@ -131,7 +131,7 @@ printMinMax(rc_t1,0)
o3l_t2@lon2d=lon2d
;
- cop_t2=a2->CO_PROD(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
+ cop_t2=a2->CO_PROD(:,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
cop_t2=cop_t2*1e9*3600
cop_t2@long_name="carbon monoxide production"
cop_t2@units="pptv/h"
@@ -139,7 +139,7 @@ printMinMax(rc_t1,0)
cop_t2@lon2d=lon2d
;
- col_t2=a2->CO_LOSS(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
+ col_t2=a2->CO_LOSS(:,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
col_t2=col_t2*1e9*3600
col_t2@long_name="carbon monoxide loss"
col_t2@units="ppbv/h"
@@ -167,7 +167,7 @@ zhat= a2->ZHAT(1:kmax+1)
do i=0,imax-1
do j=0,jmax-1
- alt(:,j,i) = nzhat*zcoef(j,i)+zs(j,i)
+ alt(0,:,j,i) = nzhat*zcoef(j,i)+zs(j,i)
end do
end do
@@ -250,23 +250,23 @@ alt@lon2d = lon2d
; PLOT
;=================================================;
; interpolation des champs a 1250 m
-rc_t1_plane = wrf_user_intrp3d(rc_t1,alt,"h",1250,0.,False)
+rc_t1_plane = wrf_user_intrp3d(rc_t1(0,:,:,:),alt(0,:,:,:),"h",1250,0.,False)
printMinMax(rc_t1_plane,0)
printMinMax(alt,0)
-rc_t2_plane = wrf_user_intrp3d(rc_t2,alt,"h",1250,0.,False)
-co_t1_plane = wrf_user_intrp3d(co_t1,alt,"h",1250,0.,False)
-co_t2_plane = wrf_user_intrp3d(co_t2,alt,"h",1250,0.,False)
-cop_t1_plane= wrf_user_intrp3d(cop_t1,alt,"h",1250,0.,False)
-cop_t2_plane= wrf_user_intrp3d(cop_t2,alt,"h",1250,0.,False)
-col_t1_plane= wrf_user_intrp3d(col_t1,alt,"h",1250,0.,False)
-col_t2_plane= wrf_user_intrp3d(col_t2,alt,"h",1250,0.,False)
-o3_t1_plane = wrf_user_intrp3d(o3_t1,alt,"h",1250,0.,False)
-o3_t2_plane = wrf_user_intrp3d(o3_t2,alt,"h",1250,0.,False)
-o3l_t1_plane= wrf_user_intrp3d(o3l_t1,alt,"h",1250,0.,False)
-o3l_t2_plane= wrf_user_intrp3d(o3l_t2,alt,"h",1250,0.,False)
-o3p_t1_plane= wrf_user_intrp3d(o3p_t1,alt,"h",1250,0.,False)
-o3p_t2_plane= wrf_user_intrp3d(o3p_t2,alt,"h",1250,0.,False)
+rc_t2_plane = wrf_user_intrp3d(rc_t2(0,:,:,:),alt(0,:,:,:),"h",1250,0.,False)
+co_t1_plane = wrf_user_intrp3d(co_t1(0,:,:,:),alt(0,:,:,:),"h",1250,0.,False)
+co_t2_plane = wrf_user_intrp3d(co_t2(0,:,:,:),alt(0,:,:,:),"h",1250,0.,False)
+cop_t1_plane= wrf_user_intrp3d(cop_t1(0,:,:,:),alt(0,:,:,:),"h",1250,0.,False)
+cop_t2_plane= wrf_user_intrp3d(cop_t2(0,:,:,:),alt(0,:,:,:),"h",1250,0.,False)
+col_t1_plane= wrf_user_intrp3d(col_t1(0,:,:,:),alt(0,:,:,:),"h",1250,0.,False)
+col_t2_plane= wrf_user_intrp3d(col_t2(0,:,:,:),alt(0,:,:,:),"h",1250,0.,False)
+o3_t1_plane = wrf_user_intrp3d(o3_t1(0,:,:,:),alt(0,:,:,:),"h",1250,0.,False)
+o3_t2_plane = wrf_user_intrp3d(o3_t2(0,:,:,:),alt(0,:,:,:),"h",1250,0.,False)
+o3l_t1_plane= wrf_user_intrp3d(o3l_t1(0,:,:,:),alt(0,:,:,:),"h",1250,0.,False)
+o3l_t2_plane= wrf_user_intrp3d(o3l_t2(0,:,:,:),alt(0,:,:,:),"h",1250,0.,False)
+o3p_t1_plane= wrf_user_intrp3d(o3p_t1(0,:,:,:),alt(0,:,:,:),"h",1250,0.,False)
+o3p_t2_plane= wrf_user_intrp3d(o3p_t2(0,:,:,:),alt(0,:,:,:),"h",1250,0.,False)
figname ="zsection_1250"
diff --git a/MY_RUN/KTEST/009_ICARTT/006_ncl/plot_ICARTT_budget.ncl b/MY_RUN/KTEST/009_ICARTT/006_ncl/plot_ICARTT_budget.ncl
index 4350a7ae8..b3c9aac40 100644
--- a/MY_RUN/KTEST/009_ICARTT/006_ncl/plot_ICARTT_budget.ncl
+++ b/MY_RUN/KTEST/009_ICARTT/006_ncl/plot_ICARTT_budget.ncl
@@ -27,41 +27,44 @@ begin
;-------------------------------------------------;
; Liste de toutes les réactions impliquant O3
- o3_list_t1=a->O3_CHREACLIST
+ o3_list_t1=a->O3_CHREACLIST(0,:)
; Tableau 4D (reac,Z,Y,X) regroupant les termes de prod. et destr. de O3
; Niveau k=13 approx. 1250m
- o3_budget_t1=a->O3_BUDGET(:,13,jphext:jmax+jphext-1,jphext:imax+jphext-1)
+ o3_budget_t1=a->O3_BUDGET(0,:,13,jphext:jmax+jphext-1,jphext:imax+jphext-1)
+ printVarSummary(o3_budget_t1)
o3_bud_t1=dim_avg_n(o3_budget_t1,(/1,2/))
+ printVarSummary(o3_bud_t1)
+
o3_bud_t1=o3_bud_t1*1e9*3600
o3_bud_t1@units="ppbv/h"
; Liste de toutes les réactions impliquant CO
- co_list_t1=a->CO_CHREACLIST
+ co_list_t1=a->CO_CHREACLIST(0,:)
; Tableau 4D (reac,Z,Y,X) regroupant les termes de prod. et destr. de CO
; Niveau k=13 approx. 1250m
- co_budget_t1=a->CO_BUDGET(:,13,jphext:jmax+jphext-1,jphext:imax+jphext-1)
+ co_budget_t1=a->CO_BUDGET(0,:,13,jphext:jmax+jphext-1,jphext:imax+jphext-1)
co_bud_t1=dim_avg_n(co_budget_t1,(/1,2/))
co_bud_t1=co_bud_t1*1e9*3600
co_bud_t1@units="ppbv/h"
; Liste de toutes les réactions impliquant O3
- o3_list_t2=a2->O3_CHREACLIST
+ o3_list_t2=a2->O3_CHREACLIST(0,:)
; Tableau 4D (reac,Z,Y,X) regroupant les termes de prod. et destr. de O3
; Niveau k=13 approx. 1250m
- o3_budget_t2=a2->O3_BUDGET(:,13,jphext:jmax+jphext-1,jphext:imax+jphext-1)
+ o3_budget_t2=a2->O3_BUDGET(0,:,13,jphext:jmax+jphext-1,jphext:imax+jphext-1)
o3_bud_t2=dim_avg_n(o3_budget_t2,(/1,2/))
o3_bud_t2=o3_bud_t2*1e9*3600
o3_bud_t2@units="ppbv/h"
; Liste de toutes les réactions impliquant CO
- co_list_t2=a2->CO_CHREACLIST
+ co_list_t2=a2->CO_CHREACLIST(0,:)
; Tableau 4D (reac,Z,Y,X) regroupant les termes de prod. et destr. de CO
; Niveau k=13 approx. 1250m
- co_budget_t2=a2->CO_BUDGET(:,13,jphext:jmax+jphext-1,jphext:imax+jphext-1)
+ co_budget_t2=a2->CO_BUDGET(0,:,13,jphext:jmax+jphext-1,jphext:imax+jphext-1)
co_bud_t2=dim_avg_n(co_budget_t2,(/1,2/))
co_bud_t2=co_bud_t2*1e9*3600
co_bud_t2@units="ppbv/h"
diff --git a/MY_RUN/KTEST/011_KW78CHEM/005_ncl/plot_KW78_chem.ncl b/MY_RUN/KTEST/011_KW78CHEM/005_ncl/plot_KW78_chem.ncl
index 0b46e4c93..9c7bdef5b 100644
--- a/MY_RUN/KTEST/011_KW78CHEM/005_ncl/plot_KW78_chem.ncl
+++ b/MY_RUN/KTEST/011_KW78CHEM/005_ncl/plot_KW78_chem.ncl
@@ -33,94 +33,94 @@ zhat = fichier1->ZHAT(1:kmax+1) ; ZHAT
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"
-rst1= fichier1->RST(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; RRT
+rst1= fichier1->RST(0,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; RRT
rst1@long_name="Snow mixing ratio"
rst1@units="kg/kg"
-rgt1= fichier1->RGT(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; RRT
+rgt1= fichier1->RGT(0,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; RRT
rgt1@long_name="Graupel mixing ratio"
rgt1@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"
-cot1= fichier1->COT(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; COT
+cot1= fichier1->COT(0,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; COT
cot1@long_name="carbon monoxide"
cot1@units="ppp"
-o3t1= fichier1->O3T(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; O3T
+o3t1= fichier1->O3T(0,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; O3T
o3t1@long_name="Ozone"
o3t1@units="ppp"
-hchot1= fichier1->HCHOT(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; HCHOT
+hchot1= fichier1->HCHOT(0,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; HCHOT
hchot1@long_name="HCHO"
hchot1@units="ppp"
-wc_hchot1= fichier1->WC_HCHOT(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
+wc_hchot1= fichier1->WC_HCHOT(0,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
wc_hchot1@long_name="WC_HCHO"
wc_hchot1@units="ppp"
-wr_hchot1= fichier1->WR_HCHOT(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
+wr_hchot1= fichier1->WR_HCHOT(0,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
wr_hchot1@long_name="WR_HCHO"
wr_hchot1@units="ppp"
-ic_hchot1= fichier1->IC_HCHOT(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
+ic_hchot1= fichier1->IC_HCHOT(0,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
ic_hchot1@long_name="IC_HCHO"
ic_hchot1@units="ppp"
-so2t1= fichier1->SO2T(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; SO2T
+so2t1= fichier1->SO2T(0,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; SO2T
so2t1@long_name="SO2"
so2t1@units="ppp"
-wc_so2t1= fichier1->WC_SO2T(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
+wc_so2t1= fichier1->WC_SO2T(0,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
wc_so2t1@long_name="WC_SO2"
wc_so2t1@units="ppp"
-wr_so2t1= fichier1->WR_SO2T(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
+wr_so2t1= fichier1->WR_SO2T(0,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
wr_so2t1@long_name="WR_SO2"
wr_so2t1@units="ppp"
-ic_so2t1= fichier1->IC_SO2T(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
+ic_so2t1= fichier1->IC_SO2T(0,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
ic_so2t1@long_name="IC_SO2"
ic_so2t1@units="ppp"
-hno3t1= fichier1->HNO3T(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; HNO3T
+hno3t1= fichier1->HNO3T(0,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; HNO3T
hno3t1@long_name="HNO3"
hno3t1@units="ppp"
-wc_hno3t1= fichier1->WC_HNO3T(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
+wc_hno3t1= fichier1->WC_HNO3T(0,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
wc_hno3t1@long_name="WC_HNO3"
wc_hno3t1@units="ppp"
-wr_hno3t1= fichier1->WR_HNO3T(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
+wr_hno3t1= fichier1->WR_HNO3T(0,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
wr_hno3t1@long_name="WR_HNO3"
wr_hno3t1@units="ppp"
-ic_hno3t1= fichier1->IC_HNO3T(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
+ic_hno3t1= fichier1->IC_HNO3T(0,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
ic_hno3t1@long_name="IC_HNO3"
ic_hno3t1@units="ppp"
-phc1= fichier1->PHC(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
+phc1= fichier1->PHC(0,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
phc1@long_name="PHC"
phc1@units=" "
-phr1= fichier1->PHR(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
+phr1= fichier1->PHR(0,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
phr1@long_name="PHR"
phr1@units=" "
diff --git a/MY_RUN/KTEST/012_dust/007_ncl/plot_dust.ncl b/MY_RUN/KTEST/012_dust/007_ncl/plot_dust.ncl
index c3553f1dd..48570d40c 100644
--- a/MY_RUN/KTEST/012_dust/007_ncl/plot_dust.ncl
+++ b/MY_RUN/KTEST/012_dust/007_ncl/plot_dust.ncl
@@ -111,11 +111,11 @@ BETA = fichier1->BETA
;=================================================;
; lecture des différents champs
;=================================================;
-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
-tht= fichier1->THT(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; THT
+tht= fichier1->THT(0,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; THT
tht@long_name="Potential Temperature"
tht@units="K"
tht@lat2d = lat2d
@@ -123,37 +123,37 @@ tht@lon2d = lon2d
-dstm03t= fichier1->DSTM03T(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
+dstm03t= fichier1->DSTM03T(0,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
dstm03t@long_name="DSTM03T"
dstm03t@units="ppp"
dstm03t@lat2d = lat2d
dstm03t@lon2d = lon2d
-dstm33t= fichier1->DSTM33T(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
+dstm33t= fichier1->DSTM33T(0,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
dstm33t@long_name="DSTM33T"
dstm33t@units="ppp"
dstm33t@lat2d = lat2d
dstm33t@lon2d = lon2d
-dstm02t= fichier1->DSTM02T(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
+dstm02t= fichier1->DSTM02T(0,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
dstm02t@long_name="DSTM02T"
dstm02t@units="ppp"
dstm02t@lat2d = lat2d
dstm02t@lon2d = lon2d
-dstm32t= fichier1->DSTM32T(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
+dstm32t= fichier1->DSTM32T(0,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
dstm32t@long_name="DSTM32T"
dstm32t@units="ppp"
dstm32t@lat2d = lat2d
dstm32t@lon2d = lon2d
-dstm01t= fichier1->DSTM01T(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
+dstm01t= fichier1->DSTM01T(0,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
dstm01t@long_name="DSTM01T"
dstm01t@units="ppp"
dstm01t@lat2d = lat2d
dstm01t@lon2d = lon2d
-dstm31t= fichier1->DSTM31T(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
+dstm31t= fichier1->DSTM31T(0,1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1)
dstm31t@long_name="DSTM31T"
dstm31t@units="ppp"
dstm31t@lat2d = lat2d
diff --git a/MY_RUN/KTEST/014_LIMA/003_ncl/plot_LIMA.ncl b/MY_RUN/KTEST/014_LIMA/003_ncl/plot_LIMA.ncl
index cf5400f7e..aca671e9a 100644
--- a/MY_RUN/KTEST/014_LIMA/003_ncl/plot_LIMA.ncl
+++ b/MY_RUN/KTEST/014_LIMA/003_ncl/plot_LIMA.ncl
@@ -16,32 +16,32 @@ begin
;=================================================;
; lecture des différents champs
;=================================================;
-vname = "RI___0001--PROC3"
+vname = "RI___0001___AVEF_BU_RRI"
RI___0001_P3= fichier1->$vname$(0,0,:,0,:)
RI___0001_P3@long_name="ice water content"
RI___0001_P3@units="g/kg"
-vname2 = "SV0050001--PROC3"
+vname2 = "SV0050001___AVEF_BU_RSV"
NI___0001_P3= fichier1->$vname2$(0,0,:,0,:)
NI___0001_P3@long_name="ice concentration"
NI___0001_P3@units="/kg"
-vname3 = "RS___0001--PROC3"
+vname3 = "RS___0001___AVEF_BU_RRS"
RS___0001_P3= fichier1->$vname3$(0,0,:,0,:)
RS___0001_P3@long_name="snow water content"
RS___0001_P3@units="g/kg"
-vname4 = "RG___0001--PROC3"
+vname4 = "RG___0001___AVEF_BU_RRG"
RG___0001_P3= fichier1->$vname4$(0,0,:,0,:)
RG___0001_P3@long_name="graupel water content"
RG___0001_P3@units="g/kg"
-vname5 = "SV0060001--PROC3"
+vname5 = "SV0060001___AVEF_BU_RSV"
N_IFN_FREE_P3= fichier1->$vname5$(0,0,:,0,:)
N_IFN_FREE_P3@long_name="concentration of free IFN"
N_IFN_FREE_P3@units="/kg"
-vname6 = "SV0070001--PROC3"
+vname6 = "SV0070001___AVEF_BU_RSV"
N_IFN_NUCL_P3= fichier1->$vname6$(0,0,:,0,:)
N_IFN_NUCL_P3@long_name="concentration of nucleated IFN"
N_IFN_NUCL_P3@units="/kg"
--
GitLab