爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 5727|回复: 1

急求ncl程序问题原因

[复制链接]

新浪微博达人勋

发表于 2015-1-6 20:25:09 | 显示全部楼层 |阅读模式

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

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

x
本人编了一段ncl程序,用于将netcdf转为bin,程序运行显示segmentation fault,求大神指教!!!
;  Open pressure data
    if (isfilepresent(Source_monthly_File_PATH_pressure+"air.mon.mean.nc")) then
      f01 = addfile(Source_monthly_File_PATH_pressure+"air.mon.mean.nc","r")
    else
      print("File "+Source_monthly_File_PATH_pressure+"air.mon.mean.nc NOT exist!")
    end if

;
;   Extract the variable data from NetCDF file
;
    air  = f01->air

;
;  Calculate monthly data
;
    data_monthly_air  = air  * air@scale_factor  + air@add_offset
    delete(air)
;
    f_time = f01->time
    n_time = dimsizes(f_time)-1
    utc_date = ut_calendar(f_time,0)
    delete(f_time)
    do i = n_time-4, n_time
      data_year  = utc_date(i,0)
      data_month = utc_date(i,1)
      c_file_year = flt2string(data_year)
      if ( data_month.lt.10 ) then
         c_file_month = "0" + flt2string(data_month)
      else
         c_file_month = flt2string(data_month)
      end if
      print(c_file_year+c_file_month)
;
;  Output monthly data
;
       file_name = monthly_File_PATH_pressure+"mon_"+c_file_year+c_file_month+".bin"
       if (isfilepresent(file_name)) then
           system("rm -f "+file_name)
       end if
       do k = 0, 23
          fbindirwrite(file_name, data_monthly_air(i,k,::-1,:))
       end do

    end do
end
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-1-6 21:14:08 | 显示全部楼层
最好能用print之类的找找看是运行到哪一行出现的问题。
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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