含泪求助！！NCL绘图遇到问题，各位大神求你们进来看看吧QAQ!

hlyahs

    感谢大家
    lz在用cdo处理数据，使用 seltime 将daily资料的4个时次分别提取出来，想使用ncl绘图，但在绘图时ncl报错显示：
  fatal:NclMalloc Failed:[errno=12]
   Segmentation fault (核心已转储)
   
    楼主顿时懵逼了QAQ，然后泪流满面的查看资料，发现cdo处理后，time变成了double格式……不知道是不是这个原因造成的

Variable: f
Type: file
filename:       6wsgd0000
path:   /cygdrive/d/ncldate/best/20years/daily/6wsgd0000.nc
   file global attributes:
      CDI : Climate Data Interface version 1.6.9 (http://mpimet.mpg.de/cdi)
      Conventions : CF-1.6
      history : Mon Nov 13 11:29:06 2017: cdo selmon,6 wsgd0000.nc 6wsgd0000.nc
Mon Nov 13 01:26:57 2017: cdo selhour,06:00:00 9715wsgd.nc wsgd0000.nc
Mon Nov 13 01:21:14 2017: cdo selyear,1997/2015 wsgd_daily678.nc 9715wsgd.nc
2017-10-02 05:06:57 GMT by grib_to_netcdf-2.4.0: grib_to_netcdf /data/data02/scratch/_mars-atls09-a82bacafb5c306db76464bc7e824bb75-UOoXIW.grib -o /data/data04/scratch/_grib2netcdf-atls09-a562cefde8a29a7288fa0b8b7f9413f7-grWqND.nc -utime
      CDO : Climate Data Operators version 1.6.9 (http://mpimet.mpg.de/cdo)
   dimensions:
      longitude = 480
      latitude = 241
      time = 570  // unlimited
   variables:
      float longitude ( longitude )
         standard_name :        longitude
         long_name :    longitude
         units :        degrees_east
         axis : X

      float latitude ( latitude )
         standard_name :        latitude
         long_name :    latitude
         units :        degrees_north
         axis : Y

     double time ( time )
         standard_name :        time
         long_name :    time
         units :        hours since 1900-01-01 00:00:00
         calendar :     standard

      short z ( time, latitude, longitude )
         standard_name :        geopotential
         long_name :    Geopotential
         units :        m**2 s**-2
         add_offset :   44209.41894082705
         scale_factor : 0.2363370959096943
         _FillValue :   -32767
         missing_value :        -32767

        请问有大神知道解决方法么，明天要给老板汇报，现在弄不出来快急哭了，学渣分分钟哭出来的节奏嘤嘤嘤QAQ!
       万分感谢各位的帮助！！！！！

hlyahs

楼主看了一下之前用cdo转换的资料，time格式也是double的，但是时次没有那么多，
但因为这次数据是逐日的，所以时次比较多，刚才重新修改了一下脚本，还是报错了：fatal:NclMalloc Failed:[errno=12]
fatal:dim_avg_n: Could not coerce input data to double, can't continue
fatal:["Execute.c":8575]:Execute: Error occurred at or near line 2352 in file $NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl


脚本是最基础最简单的脚本，就是画一个填色图而已，连overlay都没用上的，不知道到底哪出了问题，有大神愿意站出来么{:eb303:}{:eb303:}{:eb303:}{:eb303:}

    shpfn2 = "D:/ncldate/better/tibetshape/tibetan.shp"

    f = addfile("D:/ncldate/best/20years/daily/wsgd0600.nc", "r")
   

    wsgd5= short2flt(f->z(0:1747,:,:))
   
    wsgdave5 = dim_avg_n_Wrap(wsgd5,0)
    wsgdave5 = wsgdave5/9.8
  
    wks = gsn_open_wks("png","D:/ncldate/best/daily0600")


    res = True
    res@gsnDraw = False ; don't draw
    res@gsnFrame = False ; don't advance frame
    res@gsnMaximize    = True
   
    res@pmTickMarkDisplayMode = "Always"

    res@cnFillOn                = True
    res@cnLinesOn               = False
    res@cnLineLabelsOn = False
    res@cnFillPalette = "BlueYellowRed" ;  BlueWhiteOrangeRed  BlueYellowRed

   
    res@mpFillOn              =  False
    res@mpMaxLatF = 40
    res@mpMinLatF = 25
    res@mpMaxLonF = 105
    res@mpMinLonF = 70

    res@gsnAddCyclic             = False
    res@gsnLeftString = ""
    res@gsnRightString= ""
    res@cnLevelSelectionMode = "ManualLevels"
    res@cnMinLevelValF = -10
    tres@cnMaxLevelValF = 10

    res@mpFillOn              =  False          ; turn off map fill
    res@mpOutlineDrawOrder    = "PostDraw"
    res@mpPerimOn             = True
   
    res@mpDataSetName              = "Earth..4"
    res@mpDataBaseVersion          = "MediumRes" ; or "Ncarg4_1"
    res@mpAreaMaskingOn            = True
    res@mpMaskAreaSpecifiers       = (/"China"/)
    res@mpOutlineSpecifiers        = (/"China","China:Provinces"/)
    res@mpGeophysicalLineThicknessF = 1.5

    plot=  gsn_csm_contour_map(wks,wsgdave5,res)
  

    plres = True
    plres@gsLineColor = "black"
    plres@gsLineThicknessF = 3

    shp_plot     = gsn_add_shapefile_polylines(wks,tplot,shpfn2,plres)

    draw(tplot)
    frame(wks)

hlyahs

啊……各位大神不用理会lz了，楼主脑子有毛病
突然意识到其实只要再用cdo时间平均后，出来的数据就小了呀
泪崩，让这帖子沉了吧
还是感谢各位点进来看的大神们，鞠躬！

sdhjdfjsb

什么意思，我也遇到类似问题

hlyahs

sdhjdfjsb 发表于 2018-4-26 16:21
什么意思，我也遇到类似问题

我这个是数据太大照成了，ncl要求计算的数据在2G内

nwafuhcz

怎么解决的呀

yeaihua

同问，怎么解决的啊？数据就是很大