| 
 
	积分7697贡献 精华在线时间 小时注册时间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
 
 
 | 
 |