- 积分
- 6892
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2018-4-8
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
代码如下
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
begin
dir= "D:/tibet_plateau/code/"
filepath = dir+"190702.txt"
nrow0=numAsciiRow(filepath)
Lightning0=asciiread(filepath,-1,"string")
lon=stringtocharacter(str_get_field(Lightning0,1," "))
lat=stringtocharacter(str_get_field(Lightning0,2," "))
lon1=tofloat(chartostring(lon(:,0:4)))
lat1=tofloat(chartostring(lat(:,0:4)))
a=dimsizes(lat1)-1
xi = ispan(73,104,1)
yi = ispan(26,40,1)
n = 31
m = 14
md=new((/n+1,m+1/),float)
do i = 0,n
do j = 0,m
md(i,j)=0
end do
end do
do i = 0,n
do j = 0,m
do k = 0,a
if ((lat1(k)-26).ge.j.and.(lat1(k)-26).lt.(j+1).and.(lon1(k)-73).ge.i.and.(lon1(k)-73).lt.(i+1)) then
md(i,j)=md(i,j)+1
end if
end do
end do
end do
data=new((/480,3/),float)
num1=reshape(md,(/480/))
do i = 0,n
do j = 0,m
do k = 0,479
data(k,0)=yi(i)
data(k,1)=xi(j)
data(k,2)=num1(k)
end do
end do
end do
printVarSummary(data)
exit
|
|