ncl trmm卫星格式转换

风之牧语

太感谢了，没有一点问题！{:eb502:}{:eb502:}{:eb502:}

风之牧语

转成nc格式后进行发图发现时次多谢的话就会出现错误
fatal:NclMalloc Failed:[errno=12]
fatal:dim_avg_n: Could not coerce input data to double, can't continue
fatal:["Execute.c":8128]:Execute: Error occurred at or near line 2410 in file $NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl

fatal:["Execute.c":8128]:Execute: Error occurred at or near line 11 in file test.ncl
在nc读取是这样的
      files=systemfunc("ls -1 3B42RT.201305*.7.nc")
      f=addfiles(files,"r")
      ListSetType(f,"cat")
      printVarSummary(f)
      p=f[:]->PRC(:,:,:)
      t_p=dim_avg_n_Wrap(p,0)
      printVarSummary(t_p)
      t_p1=t_p*2
      copy_VarMeta(t_p,t_p1)
      printVarSummary(t_p1)

风之牧语

是不是需要在转nc格式的时候要全部tofloat或者todouble啊？

longlivehj

风之牧语 发表于 2014-7-18 09:26
转成nc格式后进行发图发现时次多谢的话就会出现错误
fatal:NclMalloc Failed:[errno=12]
fatal:dim_avg_n ...

NclMalloc 内存分配的问题，就是说运算的数据量太大了。
dim_avg_n之类的数学函数内部运算时，一般都是先转成double类型的。因此，写程序的时候，要适当考虑一下。

风之牧语

那像这种数据量比较大的怎么办呢，我上次在做模式图的时候也出现内存不够的情况

longlivehj

风之牧语 发表于 2014-7-18 09:52
那像这种数据量比较大的怎么办呢，我上次在做模式图的时候也出现内存不够的情况

p=f[:]->PRC(:,:,:)，对数据的空间范围进行限制。

风之牧语

没有明白，你是说少读取数据的维数吗？还是用其他的语句对数据进行限制？像刚才那个问题，我使用todouble或者tofloat同样就会出现这样的错误

longlivehj

风之牧语 发表于 2014-7-18 10:07
没有明白，你是说少读取数据的维数吗？还是用其他的语句对数据进行限制？像刚才那个问题，我使用todouble或 ...

点回复，否则不能及时看到。
如果是对局部地区进行研究，就只去除相应空间范围的数据就好。
不用显式的去转换数据类型。

风之牧语

longlivehj 发表于 2014-7-18 10:48
点回复，否则不能及时看到。
如果是对局部地区进行研究，就只去除相应空间范围的数据就好。
不用显式的 ...

恩。好的！我前面在画模式图的时候模拟中国区域之后就需要画整个区域的，也出现了内存不够，你说的这个我到可以在这个卫星数据里面用！

爽歪歪！

神贴加神人