⽔汽通量的经向输送分析
摘要：使⽤2012年NCEP/NCAR再分析经向风资料和⽐湿资料，计算总的经向输送(TOTAL)、定常经向环流输送(MMC)、定常波输送(SE)以及瞬变部分输送(TE)。利⽤NCL绘制各项的垂直剖⾯图，分析⽔汽通量的经向输送。分析发现MMC对于⽔汽的输送作⽤最为明显，主要体现为Hadley环流将⽔汽向低纬地区输送以及Ferrel环流将⽔汽向⾼纬地区输送，⽽TE与SE都主要将⽔汽向⾼纬地区输送，且TE的作⽤⽐SE稍⼤。 关键词：⽔汽通量，经向输送，时空分解
1、引⾔：⽔汽的输送，直接决定着⼤⽓的各个物理过程，决定着全球的⽓候，以及天⽓状况，影响着全球经济⽣产活动。
2、使⽤资料：本⽂使⽤了2012年NCEP/NCAR再分析资料中的经向风场和⽐湿场资料，垂直⽅向取8层(1000、900、850、
700、600、500、400、300hPa)。⽔平分辨率均为2.5度。
3、⽔汽经向输送分析：计算总的经向输送(TOTAL)、定常经向环流输送(MMC)定常波输送(SE)，以及瞬变经向环流与瞬变涡动输送之和(TE)。分别见下图
从图中可以发现⽔汽主要集中在850hPa以下的低层⼤⽓。⽔汽的定常经向环流输送(MMC)量级最⼤，接近总的⽔汽输送(TOTAL)。它的主要作⽤范围在低纬地区。北半球低纬地区(5-30N)是个强的负值区；南半球⾚道附近地区(30S-5N)存在⼀个强的正值区，这些区域正是Hadley环流所在区，因此，MMC反映了⽔汽在经向的输送，主要是Hadley环流将⽔汽向低纬输送，此外还有中纬度的Ferrel环流将⽔汽向⾼纬输送。定常波输送(SE)的作⽤则相对较⼩。图中南半球⾚道附近及低纬地区是负值中⼼，⽽北半球中纬附近地区为正值中⼼。因此，在南北半球，SE都起到了向极地的⽔汽输送作⽤。瞬变部分输送(TE)主要在中⾼纬地区，作⽤⽐SE稍⼤。输送主要来⾃于向极地移动的暖空⽓和向低纬移动的冷空⽓。图中北半球地区为正值区，南半球地区为负值区，都表⽰向极地输送⽔汽。
4、总结：本⽂得到的结果与实际结果相符，TOTAL=MMC+SE+TE，并且反应出了各项对⽔汽输送的重要性。
附：NCL数据处理脚本：
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/l"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/l"
load "$NCARG_ROOT/lib/ncarg/nclscripts/l"
begin
a=addfile("./","r")
b=addfile("./","r")
v1=short2flt(a->vwnd)
q=short2flt(b->shum)*1000
lat=b->lat
lev_p=b->level
v=v1(:,0:7,:,:)
Avgtime_v=dim_avg_n(v,0)
timetb_v=v
do z=0,7,1
do t=0,365,1
timetb_v(t,z,y,x)=v(t,z,y,x) - Avgtime_v(z,y,x) end do
end do
end do
end do
Avgzon_v=dim_avg_n(v,3)
zontb_v=v
do t =0,365,1
do z =0,7,1
do y=0,72,1
do x=0,143,1
zontb_v(t,z,y,x)=v(t,z,y,x) - Avgzon_v(t,z,y) end do
end do
end do
end do
Avgtime_q=dim_avg_n(q,0)
timetb_q=q
do z=0,7,1
do y=0,72,1
do x=0,143,1
do t=0,365,1
timetb_q(t,z,y,x)=q(t,z,y,x) - Avgtime_q(z,y,x) end do
end do
end do
end do
Avgzon_q=dim_avg_n(q,3)
zontb_q=q
do t=0,365,1
do z=0,7,1
zontb_q(t,z,y,x)=q(t,z,y,x) - Avgzon_q(t,z,y)
end do
end do
end do
end do
AvgzonAvgtime_v=dim_avg_n(Avgtime_v,2)
网络辅导
AvgzonAvgtimetb_v=Avgtime_v
do z=0,7,1
do y=0,72,1
do x=0,143,1
AvgzonAvgtimetb_v(z,y,x)=Avgtime_v(z,y,x) - AvgzonAvgtime_v(z,y) end do
end do
end do
AvgzonAvgtime_q=dim_avg_n(Avgtime_q,2)
AvgzonAvgtimetb_q=Avgtime_q
do z=0,7,1
do y=0,72,1
do x=0,143,1
AvgzonAvgtimetb_q(z,y,x)=Avgtime_q(z,y,x) - AvgzonAvgtime_q(z,y) end do
end do
end do
探秘古墓
;Total
vq=v
do t=0,365,1
do z=0,7,1
do y=0,72,1
do x=0,143,1
vq(t,z,y,x)=v(t,z,y,x) * q(t,z,y,x)
end do
end do
end do
vq_Avgtime=dim_avg_n(vq,0)
Total=dim_avg_n(vq_Avgtime,2)
;MMC
MMC=AvgzonAvgtime_v
吉尔达鞋业do z=0,7,1
do y=0,72,1
MMC(z,y)=AvgzonAvgtime_v(z,y) * AvgzonAvgtime_q(z,y)
end do
end do
;SE
Totaltb=AvgzonAvgtimetb_v
do z=0,7,1
do y=0,72,1
do x=0,143,1
Totaltb(z,y,x)=AvgzonAvgtimetb_v(z,y,x) * AvgzonAvgtimetb_q(z,y,x) end do
大连破网售假药案end do
end do
成像设计SE=dim_avg_n(Totaltb,2)
;TE
Totaltb1=timetb_v
do t=0,365,1
do z=0,7,1
do y=0,72,1
do x=0,143,1
Totaltb1(t,z,y,x)=timetb_v(t,z,y,x) * timetb_q(t,z,y,x)
end do
end do
end do
end do
功能高分子材料TE1=dim_avg_n(Totaltb1,0)
TE=dim_avg_n(TE1,2)
Total!0="lev"
Total!1="lat"
Total&lev=lev_p
Total&lat=lat
MMC!0 = "lev"
MMC!1 = "lat"
MMC&lev = lev_p
MMC&lat = lat
SE!0 = "lev"
SE!1 = "lat"
SE&lev = lev_p
SE&lat = lat
TE!0 = "lev"
TE!1 = "lat"
TE&lev = lev_p
TE&lat = lat
;output
system("bin/rm -f ./")
d = addfile("./","c")
d->Total = Total
system("bin/rm -f ./")
e = addfile("./","c")
e->MMC = MMC
system("bin/rm -f ./SE.nc")
f = addfile("./SE.nc","c")
f->SE = SE
system("bin/rm -f ./TE.nc")
h = addfile("./TE.nc","c")
h->TE = TE
end
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议