Skip to content
Snippets Groups Projects
Commit 19088a1e authored by Gaelle DELAUTIER's avatar Gaelle DELAUTIER
Browse files

Gaelle 23/04/2018 : bug ncl script corrections

parent 14972ad8
No related branches found
Tags PACK-MNH-V5-4--0-09bis_IO
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Showing
with 103 additions and 100 deletions
MY_RUN/KTEST/004_Reunion/006_ncl/plot_Reunion.ncl
+ 6
− 6
View file @ 19088a1e
...@@ -102,28 +102,28 @@ load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl" ...@@ -102,28 +102,28 @@ load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
zs@lat2d = lat2d zs@lat2d = lat2d
zs@lon2d = lon2d zs@lon2d = lon2d
wt= fichier1->WT(1:kmax+1,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; WT wt= fichier1->WT(0,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 vt= fichier1->VT(0,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 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(1:kmax,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; THT
tht@long_name="Potential Temperature" tht@long_name="Potential Temperature"
tht@units="K" tht@units="K"
tht@lat2d = lat2d tht@lat2d = lat2d
tht@lon2d = lon2d 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@long_name="Large SCale Potential Temperature"
lsthm@units="K" lsthm@units="K"
lsthm@lat2d = lat2d lsthm@lat2d = lat2d
lsthm@lon2d = lon2d 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@long_name="Large SCale Merdian Wind"
lsvm@units="m/s" lsvm@units="m/s"
lsvm@lat2d = lat2d lsvm@lat2d = lat2d
lsvm@lon2d = lon2d 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@long_name="Pressure"
pabst@units="Pa" pabst@units="Pa"
pabst@lat2d = lat2d pabst@lat2d = lat2d
......
MY_RUN/KTEST/007_16janvier/011_ncl/plot_16j.ncl
+ 6
− 6
View file @ 19088a1e
...@@ -110,11 +110,11 @@ zs@lat2d = lat2d ...@@ -110,11 +110,11 @@ zs@lat2d = lat2d
zs@lon2d = lon2d zs@lon2d = lon2d
wt= fichier1->WT(1:kmax+1,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; WT wt= fichier1->WT(0,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 vt= fichier1->VT(0,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 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@long_name="Potential Temperature"
tht1@units="K" tht1@units="K"
tht1@lat2d = lat2d tht1@lat2d = lat2d
...@@ -122,13 +122,13 @@ tht1@lon2d = lon2d ...@@ -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@long_name="Pressure"
pabst1@units="Pa" pabst1@units="Pa"
pabst1@lat2d = lat2d pabst1@lat2d = lat2d
pabst1@lon2d = lon2d 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=rvt1*1000
rvt1@long_name="Vapor mixing ratio" rvt1@long_name="Vapor mixing ratio"
rvt1@units="kg/kg" rvt1@units="kg/kg"
......
MY_RUN/KTEST/007_16janvier/011_ncl/plot_16j_2.ncl
+ 6
− 6
View file @ 19088a1e
...@@ -109,11 +109,11 @@ zs@lat2d = lat2d ...@@ -109,11 +109,11 @@ zs@lat2d = lat2d
zs@lon2d = lon2d zs@lon2d = lon2d
wt= fichier1->WT(1:kmax+1,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; WT wt= fichier1->WT(0,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 vt= fichier1->VT(0,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 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@long_name="Potential Temperature"
tht1@units="K" tht1@units="K"
tht1@lat2d = lat2d tht1@lat2d = lat2d
...@@ -121,13 +121,13 @@ tht1@lon2d = lon2d ...@@ -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@long_name="Pressure"
pabst1@units="Pa" pabst1@units="Pa"
pabst1@lat2d = lat2d pabst1@lat2d = lat2d
pabst1@lon2d = lon2d 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=rvt1*1000
rvt1@long_name="Vapor mixing ratio" rvt1@long_name="Vapor mixing ratio"
rvt1@units="kg/kg" rvt1@units="kg/kg"
......
MY_RUN/KTEST/009_ICARTT/006_ncl/plot_ICARTT.ncl
+ 29
− 29
View file @ 19088a1e
...@@ -38,7 +38,7 @@ zs@lon2d = lon2d ...@@ -38,7 +38,7 @@ zs@lon2d = lon2d
printMinMax(zs,0) 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@long_name="Cloud mixing ratio"
rc_t1@units="g/kg" rc_t1@units="g/kg"
rc_t1@lat2d=lat2d rc_t1@lat2d=lat2d
...@@ -46,21 +46,21 @@ printMinMax(zs,0) ...@@ -46,21 +46,21 @@ printMinMax(zs,0)
printMinMax(rc_t1,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@long_name="Ozone"
o3_t1@units="ppbv" o3_t1@units="ppbv"
o3_t1@lat2d=lat2d o3_t1@lat2d=lat2d
o3_t1@lon2d=lon2d 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@long_name="carbon monoxide"
co_t1@units="ppbv" co_t1@units="ppbv"
co_t1@lat2d=lat2d co_t1@lat2d=lat2d
co_t1@lon2d=lon2d 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 = o3p_t1*1e9*3600
o3p_t1@long_name="ozone production" o3p_t1@long_name="ozone production"
o3p_t1@units="ppbv/h" o3p_t1@units="ppbv/h"
...@@ -68,7 +68,7 @@ printMinMax(rc_t1,0) ...@@ -68,7 +68,7 @@ printMinMax(rc_t1,0)
o3p_t1@lon2d=lon2d 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 = o3l_t1*1e9*3600
o3l_t1@long_name="ozone loss" o3l_t1@long_name="ozone loss"
o3l_t1@units="ppbv/h" o3l_t1@units="ppbv/h"
...@@ -76,7 +76,7 @@ printMinMax(rc_t1,0) ...@@ -76,7 +76,7 @@ printMinMax(rc_t1,0)
o3l_t1@lon2d=lon2d 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 = cop_t1*1e9*3600
cop_t1@long_name="carbon monoxide production" cop_t1@long_name="carbon monoxide production"
cop_t1@units="ppbv/h" cop_t1@units="ppbv/h"
...@@ -84,7 +84,7 @@ printMinMax(rc_t1,0) ...@@ -84,7 +84,7 @@ printMinMax(rc_t1,0)
cop_t1@lon2d=lon2d 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 = col_t1*1e9*3600
col_t1@long_name="carbon monoxide loss" col_t1@long_name="carbon monoxide loss"
col_t1@units="ppbv/h" col_t1@units="ppbv/h"
...@@ -94,28 +94,28 @@ printMinMax(rc_t1,0) ...@@ -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@long_name="Cloud mixing ratio"
rc_t2@units="g/kg" rc_t2@units="g/kg"
rc_t2@lat2d=lat2d rc_t2@lat2d=lat2d
rc_t2@lon2d=lon2d 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@long_name="Ozone"
o3_t2@units="ppbv" o3_t2@units="ppbv"
o3_t2@lat2d=lat2d o3_t2@lat2d=lat2d
o3_t2@lon2d=lon2d 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@long_name="carbon monoxide"
co_t2@units="ppbv" co_t2@units="ppbv"
co_t2@lat2d=lat2d co_t2@lat2d=lat2d
co_t2@lon2d=lon2d 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=o3p_t2*1e9*3600
o3p_t2@long_name="ozone production" o3p_t2@long_name="ozone production"
o3p_t2@units="pptv/h" o3p_t2@units="pptv/h"
...@@ -123,7 +123,7 @@ printMinMax(rc_t1,0) ...@@ -123,7 +123,7 @@ printMinMax(rc_t1,0)
o3p_t2@lon2d=lon2d 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=o3l_t2*1e9*3600
o3l_t2@long_name="ozone loss" o3l_t2@long_name="ozone loss"
o3l_t2@units="pptv/h" o3l_t2@units="pptv/h"
...@@ -131,7 +131,7 @@ printMinMax(rc_t1,0) ...@@ -131,7 +131,7 @@ printMinMax(rc_t1,0)
o3l_t2@lon2d=lon2d 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=cop_t2*1e9*3600
cop_t2@long_name="carbon monoxide production" cop_t2@long_name="carbon monoxide production"
cop_t2@units="pptv/h" cop_t2@units="pptv/h"
...@@ -139,7 +139,7 @@ printMinMax(rc_t1,0) ...@@ -139,7 +139,7 @@ printMinMax(rc_t1,0)
cop_t2@lon2d=lon2d 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=col_t2*1e9*3600
col_t2@long_name="carbon monoxide loss" col_t2@long_name="carbon monoxide loss"
col_t2@units="ppbv/h" col_t2@units="ppbv/h"
...@@ -167,7 +167,7 @@ zhat= a2->ZHAT(1:kmax+1) ...@@ -167,7 +167,7 @@ zhat= a2->ZHAT(1:kmax+1)
do i=0,imax-1 do i=0,imax-1
do j=0,jmax-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
end do end do
...@@ -250,23 +250,23 @@ alt@lon2d = lon2d ...@@ -250,23 +250,23 @@ alt@lon2d = lon2d
; PLOT ; PLOT
;=================================================; ;=================================================;
; interpolation des champs a 1250 m ; 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(rc_t1_plane,0)
printMinMax(alt,0) printMinMax(alt,0)
rc_t2_plane = wrf_user_intrp3d(rc_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,alt,"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,alt,"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,alt,"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,alt,"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,alt,"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,alt,"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,alt,"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,alt,"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,alt,"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,alt,"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,alt,"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,alt,"h",1250,0.,False) o3p_t2_plane= wrf_user_intrp3d(o3p_t2(0,:,:,:),alt(0,:,:,:),"h",1250,0.,False)
figname ="zsection_1250" figname ="zsection_1250"
......
MY_RUN/KTEST/009_ICARTT/006_ncl/plot_ICARTT_budget.ncl
+ 11
− 8
View file @ 19088a1e
...@@ -27,41 +27,44 @@ begin ...@@ -27,41 +27,44 @@ begin
;-------------------------------------------------; ;-------------------------------------------------;
; Liste de toutes les réactions impliquant O3 ; 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 ; Tableau 4D (reac,Z,Y,X) regroupant les termes de prod. et destr. de O3
; Niveau k=13 approx. 1250m ; 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/)) 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=o3_bud_t1*1e9*3600
o3_bud_t1@units="ppbv/h" o3_bud_t1@units="ppbv/h"
; Liste de toutes les réactions impliquant CO ; 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 ; Tableau 4D (reac,Z,Y,X) regroupant les termes de prod. et destr. de CO
; Niveau k=13 approx. 1250m ; 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=dim_avg_n(co_budget_t1,(/1,2/))
co_bud_t1=co_bud_t1*1e9*3600 co_bud_t1=co_bud_t1*1e9*3600
co_bud_t1@units="ppbv/h" co_bud_t1@units="ppbv/h"
; Liste de toutes les réactions impliquant O3 ; 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 ; Tableau 4D (reac,Z,Y,X) regroupant les termes de prod. et destr. de O3
; Niveau k=13 approx. 1250m ; 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=dim_avg_n(o3_budget_t2,(/1,2/))
o3_bud_t2=o3_bud_t2*1e9*3600 o3_bud_t2=o3_bud_t2*1e9*3600
o3_bud_t2@units="ppbv/h" o3_bud_t2@units="ppbv/h"
; Liste de toutes les réactions impliquant CO ; 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 ; Tableau 4D (reac,Z,Y,X) regroupant les termes de prod. et destr. de CO
; Niveau k=13 approx. 1250m ; 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=dim_avg_n(co_budget_t2,(/1,2/))
co_bud_t2=co_bud_t2*1e9*3600 co_bud_t2=co_bud_t2*1e9*3600
co_bud_t2@units="ppbv/h" co_bud_t2@units="ppbv/h"
......
MY_RUN/KTEST/011_KW78CHEM/005_ncl/plot_KW78_chem.ncl
+ 29
− 29
View file @ 19088a1e
...@@ -33,94 +33,94 @@ zhat = fichier1->ZHAT(1:kmax+1) ; ZHAT ...@@ -33,94 +33,94 @@ zhat = fichier1->ZHAT(1:kmax+1) ; ZHAT
yhat = fichier1->YHAT(1:jmax+jphext-1+1) ; YHAT yhat = fichier1->YHAT(1:jmax+jphext-1+1) ; YHAT
xhat = fichier1->XHAT(1:imax+jphext-1+1) ; XHAT xhat = fichier1->XHAT(1:imax+jphext-1+1) ; XHAT
wt= fichier1->WT(1:kmax+1,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; WT wt= fichier1->WT(0,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 vt= fichier1->VT(0,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 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@long_name="Potential Temperature"
tht1@units="K" 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@long_name="Large SCale Potential Temperature"
lsthm1@units="K" 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@long_name="Pressure"
pabst1@units="Pa" 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@long_name="Vapor mixing ratio"
rvt1@units="kg/kg" 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@long_name="Cloud mixing ratio"
rct1@units="kg/kg" 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@long_name="Rain mixing ratio"
rrt1@units="kg/kg" 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@long_name="Snow mixing ratio"
rst1@units="kg/kg" 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@long_name="Graupel mixing ratio"
rgt1@units="kg/kg" 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@long_name="Instantaneous Precipitation Rain Rate"
inprr@units="mm/h" 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@long_name="Accumulated Precipitation Rain Rate"
acprr@units="mm" 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@long_name="carbon monoxide"
cot1@units="ppp" 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@long_name="Ozone"
o3t1@units="ppp" 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@long_name="HCHO"
hchot1@units="ppp" 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@long_name="WC_HCHO"
wc_hchot1@units="ppp" 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@long_name="WR_HCHO"
wr_hchot1@units="ppp" 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@long_name="IC_HCHO"
ic_hchot1@units="ppp" 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@long_name="SO2"
so2t1@units="ppp" 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@long_name="WC_SO2"
wc_so2t1@units="ppp" 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@long_name="WR_SO2"
wr_so2t1@units="ppp" 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@long_name="IC_SO2"
ic_so2t1@units="ppp" 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@long_name="HNO3"
hno3t1@units="ppp" 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@long_name="WC_HNO3"
wc_hno3t1@units="ppp" 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@long_name="WR_HNO3"
wr_hno3t1@units="ppp" 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@long_name="IC_HNO3"
ic_hno3t1@units="ppp" 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@long_name="PHC"
phc1@units=" " 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@long_name="PHR"
phr1@units=" " phr1@units=" "
......
MY_RUN/KTEST/012_dust/007_ncl/plot_dust.ncl
+ 10
− 10
View file @ 19088a1e
...@@ -111,11 +111,11 @@ BETA = fichier1->BETA ...@@ -111,11 +111,11 @@ BETA = fichier1->BETA
;=================================================; ;=================================================;
; lecture des différents champs ; lecture des différents champs
;=================================================; ;=================================================;
wt= fichier1->WT(1:kmax+1,jphext:jmax+jphext-1,jphext:imax+jphext-1) ; WT wt= fichier1->WT(0,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 vt= fichier1->VT(0,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 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@long_name="Potential Temperature"
tht@units="K" tht@units="K"
tht@lat2d = lat2d tht@lat2d = lat2d
...@@ -123,37 +123,37 @@ tht@lon2d = lon2d ...@@ -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@long_name="DSTM03T"
dstm03t@units="ppp" dstm03t@units="ppp"
dstm03t@lat2d = lat2d dstm03t@lat2d = lat2d
dstm03t@lon2d = lon2d 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@long_name="DSTM33T"
dstm33t@units="ppp" dstm33t@units="ppp"
dstm33t@lat2d = lat2d dstm33t@lat2d = lat2d
dstm33t@lon2d = lon2d 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@long_name="DSTM02T"
dstm02t@units="ppp" dstm02t@units="ppp"
dstm02t@lat2d = lat2d dstm02t@lat2d = lat2d
dstm02t@lon2d = lon2d 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@long_name="DSTM32T"
dstm32t@units="ppp" dstm32t@units="ppp"
dstm32t@lat2d = lat2d dstm32t@lat2d = lat2d
dstm32t@lon2d = lon2d 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@long_name="DSTM01T"
dstm01t@units="ppp" dstm01t@units="ppp"
dstm01t@lat2d = lat2d dstm01t@lat2d = lat2d
dstm01t@lon2d = lon2d 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@long_name="DSTM31T"
dstm31t@units="ppp" dstm31t@units="ppp"
dstm31t@lat2d = lat2d dstm31t@lat2d = lat2d
......
MY_RUN/KTEST/014_LIMA/003_ncl/plot_LIMA.ncl
+ 6
− 6
View file @ 19088a1e
...@@ -16,32 +16,32 @@ begin ...@@ -16,32 +16,32 @@ begin
;=================================================; ;=================================================;
; lecture des différents champs ; 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= fichier1->$vname$(0,0,:,0,:)
RI___0001_P3@long_name="ice water content" RI___0001_P3@long_name="ice water content"
RI___0001_P3@units="g/kg" 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= fichier1->$vname2$(0,0,:,0,:)
NI___0001_P3@long_name="ice concentration" NI___0001_P3@long_name="ice concentration"
NI___0001_P3@units="/kg" 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= fichier1->$vname3$(0,0,:,0,:)
RS___0001_P3@long_name="snow water content" RS___0001_P3@long_name="snow water content"
RS___0001_P3@units="g/kg" 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= fichier1->$vname4$(0,0,:,0,:)
RG___0001_P3@long_name="graupel water content" RG___0001_P3@long_name="graupel water content"
RG___0001_P3@units="g/kg" 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= fichier1->$vname5$(0,0,:,0,:)
N_IFN_FREE_P3@long_name="concentration of free IFN" N_IFN_FREE_P3@long_name="concentration of free IFN"
N_IFN_FREE_P3@units="/kg" 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= fichier1->$vname6$(0,0,:,0,:)
N_IFN_NUCL_P3@long_name="concentration of nucleated IFN" N_IFN_NUCL_P3@long_name="concentration of nucleated IFN"
N_IFN_NUCL_P3@units="/kg" N_IFN_NUCL_P3@units="/kg"
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment

Documentation | Mentions légales | CGU | Accessibilité : non conforme