- 积分
- 794
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-10-11
- 最后登录
- 1970-1-1

|

楼主 |
发表于 2016-12-21 18:57:25
|
显示全部楼层
begin
nt=101
nlev=81
nlat=161
nlon=161
q=new((/nt,nlev,nlat,nlon/),"float")
u=q
v=q
w=q
th=q
p=q
th2D=new((/nlon,nlat/),"float")
x=new(nlon,"float")
y=new(nlat,"float")
do gg=0,nt-1
q(gg,:,:,:)=fbinrecread("./123.dat",0+6*gg,(/nlev,nlat,nlon/),"float")
u(gg,:,:,:)=fbinrecread("./123.dat",1+6*gg,(/nlev,nlat,nlon/),"float")
v(gg,:,:,:)=fbinrecread("./123.dat",2+6*gg,(/nlev,nlat,nlon/),"float")
w(gg,:,:,:)=fbinrecread("./123.dat",3+6*gg,(/nlev,nlat,nlon/),"float")
th(gg,:,:,:)=fbinrecread("./123.dat",4+6*gg,(/nlev,nlat,nlon/),"float")
p(gg,:,:,:)=fbinrecread("./123.dat",5+6*gg,(/nlev,nlat,nlon/),"float")
end do
do j=0,nlat-1
do i=0,nlon-1
th2D(i,j)=th(nt-1,11,j,i)
x(i)=0.0625*i
y(j)=0.0625*j
end do
end do
;th2D=th2D(:,:)
;th2D@lon2d=x(:,:)
;th2D@lat2d=y(:,:)
th2D!0="lon"
th2D!1="lat"
th2D&lon=x
th2D&lat=y
printVarSummary(th2D)
wks=gsn_open_wks("x11","plot")
res=True
res@tmXBMode="Explicit"
res@tmXBValues=ispan(3,8,1)
res@tmXBLabels=(/3,4,5,6,7,8/)
res@tmLabelAutoStride=True
res@tmYLMode="Explicit"
res@tmYLValues=ispan(3,8,1)
res@tmYLLabels=(/3,4,5,6,7,8/)
res@tmLabelAutoStride=True
plot=gsn_csm_contour(wks,th2D,res)
end |
|