- 积分
- 86
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2013-5-16
- 最后登录
- 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"
begin
fileout="avg1jtwc.txt"
;system("rm -f "+fileout)
system(fileout)
filein="cutdata1_8013_all_jtwc.txt"
nline=numAsciiCol(filein)
print (nline)
in=asciiread(filein,(/nline/),"string")
tc_time=stringtoint(str_get_field(in,1," "))
tc_lat=stringtofloat(str_get_field(in,3," "))
tc_lon=stringtofloat(str_get_field(in,4," "))
tc_wind=stringtofloat(str_get_field(in,5," "))
tc_name=str_get_field(in,8," ")
index=ind(tc_time.ne."")
n=dimsizes(index)
do i=0,n-1
t1=index(i)
if(tc_wind(t1).gt.33.4)then
cutlat=tc_lat(ind(tc_wind(t1)))
cutlon=tc_lon(ind(tc_wind(t1)))
if(cutlat.ge.0)then
system("echo '"+sprinti("%10i",tc_time(t1)) \
+sprintf("%9.1f",cutlat) \
+sprintf("%9.1f",cutlon) \
+"' >> "+fileout)
end if
end if
end do
end
|
|