- 积分
- 237
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2016-4-12
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
以下是我的脚本
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
begin
list_geo="d:/cloudsat/1507_geoprof.lst"
listname=asciiread(list_geo,-1,"string")
nl=dimsizes(listname);
do i=0,nl-1
str=stringtochar(listname(i));
if(str(0:1).ne."20")
yr=str(0:1)
mn=str(2:3)
dy=str(4:5)
str2="d:/cloudsat/"+yr+mn+"/"+yr+mn+dy+"/"
delete(str);
else
filename="d:/cloudsat/"+str(0:54);
a=addfile(filename,"r")
var=a->Radar_Reflectivity;
l=dimsizes(var);
radar=new((/l(1),l(0)/),"float")
radar=transpose(var)/100.0;
delete(var);
var=a->CPR_Cloud_mask;
cldmask=new((/l(1),l(0)/),"integer")
cldmask=transpose(var);
delete(var);
fileradar=str2+"radar."+str(14:18)+".grd"
filemask=str2+"cldmask."+str(14:18)+".grd"
fbindirwrite(fileradar,radar);
fbindirwrite(filemask,cldmask);
delete(l)
delete(radar);
delete(cldmask);
delete(str);
delete(a);
end if
end do
end
line 35 是红色的那行,
小白想知道为什么出问题了(脚本是老师给的,所以我也不知道哪里错了)
|
|