Ncl结果出现这个的原因是内存不足导致的吗？

王大琳

图片里的line19是我构建工作台的命令  19    wks = gsn_open_wks("png","S4_ET_Sum_annual")

风之牧语

这种情况内存不够几率是比较高的。你这个是应该数据量太大在调用NCL自带函数时候内存不够返回的错误

御风者

信息太少了

男紫汗

line 19.。。

xigua

内存不足好像会显示memory字样

王大琳

御风者 发表于 2017-12-3 22:49
信息太少了

1 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
  2 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
  3 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
  4 begin
  5    f_S2= addfile("CABLE-POP_S2_evapotrans.nc","r")
  6    ET_S2=f_S2->evapotrans(:,:,:)
  7    ET_sum_S2=dim_sum_n_Wrap(ET_S2,0)
  8    ET_sum_annual_S2=ET_sum_S2/156
  9    copy_VarMeta(ET_sum_S2,ET_sum_annual_S2)
10    f_S3= addfile("CABLE-POP_S3_evapotrans.nc","r")
11    ET_S3= f_S3->evapotrans(:,:,:)
12    ET_sum_S3=dim_sum_n_Wrap(ET_S3,0)
13    ET_sum_annual_S3=ET_sum_S3/156
14    copy_VarMeta(ET_sum_S3,ET_sum_annual_S3)
15    ET_sum_annual_S4=ET_sum_annual_S3-ET_sum_annual_S2
16    copy_VarMeta(ET_sum_annual_S3,ET_sum_annual_S4)
17    printVarSummary(ET_sum_annual_S4)
18    delete([/ET_sum_annual_S3,ET_sum_annual_S2/])
19   print(ET_sum_annual_S4)
20     wks = gsn_open_wks("png","S4_ET_Sum_annual")
21 res = True
22    res@cnLevelSelectionMode = "ManualLevels"
23    res@cnFillOn = True ; Turn on contour fill
24    res@lbOrientation = "Vertical" ; Move labelbar
25    res@cnFillPalette = "BlueYellowRed"
26    res@cnLinesOn = False
27   plot = gsn_csm_contour_map(wks,ET_sum_annual_S4,res)
28  end
上面是我的代码，我写这段代码是想读取两个.nc文件的ET（time，x，y），由于这个time是月份，我想求的是蒸发量ET的年平均值，然后我用dim_sum_n_Wrap 得到这么多月（1872个月，156年）总的蒸发量   ，然后再除以156得到年均值，这就得到两个.nc文件的ET（） 的年平均值，接着我想让他两相减，得到的差值就可以知道由CO2单独引起的蒸散量的变化是多少。但是在画图的时候出现上图的错误。