求教，ncl求几十年温度场平均的问题

橙子鲜

请问下，我用大模式输出资料画45年的2m温度平均，资料为nc,每十年放在一个nc之中，分辨率1.875*1.875，只是想求平均温度，就有如下两个问题

1.fatal：dim_avg_n：could not coerce input data to double,can't continue
2.warning: non-monotonic coordinate array generated--check validity of the aggregated  (leftmost dimesion)

脚本：
begin
diri="e:/LMDZ/data/historical/tas/"
fils=systemfunc ("ls " + diri + "tas_day_MPI-ESM-MR_historical_r1i1p1*")
data=addfiles(fils,"r")
tas=data[:]->tas;;;;;;;;;;;            tas(3653,96,192)             1.875*1.875

lat=data[:]->lat
lon=data[:]->lon
time=data[:]->time


t=tas-273.15
copy_VarCoords(tas,t)
； printVarSummary(t)
tavg= dim_avg_n_Wrap (t,0)


麻烦各位帮忙看一下这些问题，问题意思我大概知道，但是不知道如何改，还请提出宝贵意见，谢谢

longlivehj

1. 数据类型转换出错，有可能是数组太大，内存不够
2. waring问题不大

andrewsoong

你怎么不试试CDO

橙子鲜

longlivehj 发表于 2015-6-28 08:36
1. 数据类型转换出错，有可能是数组太大，内存不够
2. waring问题不大

恩恩，对是数组的问题，只能换个方法了，多谢

橙子鲜

andrewsoong 发表于 2015-6-28 15:19
你怎么不试试CDO

这个从来没用过呢，我学习试试

橙子鲜

longlivehj 发表于 2015-6-28 08:36
1. 数据类型转换出错，有可能是数组太大，内存不够
2. waring问题不大

可不可以顺便问个问题，我今天做双线性插值，插值结果出现好多9.9621e+36,只有零零星星的数值，什么问题呢，脚本在下面
f1=new((/80,77/),float)
filelat="e:/LMDZ/data/1.txt"
f1 = asciiread(filelat,(/80,77/),"float")
  ;print(f1)
f1!0="lat"
  f1&lat=lat
f1!1="lon"
  f1&lon=lon;;;;;;;;;;;;;;;;;;;;;;;;将f2的经纬度赋值给f1
; printVarSummary(f1)

    ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;做双线性插值
  LON = fspan(86.15,134.15,49);;;;0.6
  LAT = fspan(53.81723,5.81723,81);;;0.6


  fo = linint2_Wrap (lon,lat,f3, True, LON,LAT,0)
   print(fo)

zhangkf_7

你好，我也遇到这个问题，还有其他方法可以对这个海温场求平均吗？

橙子鲜

zhangkf_7 发表于 2017-9-23 12:57
你好，我也遇到这个问题，还有其他方法可以对这个海温场求平均吗？

数组太大，无法一次新读入进去，可以用循环，将小维度的数组一条一条循环计算