爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 6004|回复: 0

小白在做ncl提取数据时出了问题,求大神帮帮忙

[复制链接]
发表于 2018-4-8 23:14:43 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册

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 是红色的那行,
小白想知道为什么出问题了(脚本是老师给的,所以我也不知道哪里错了)
图片1.png
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表