ncl读观测数据，有缺省值，影响数据读取，怎么办

黄小仙儿

ncl读观测数据，有缺省值，影响数据读取，怎么办？
怎么设置缺省值啊？缺省值是PPC，直接读float是不是有影响


ascii_filename = "/public/home/huanglei/data/ob.txt"
  seismic = asciiread(ascii_filename,(/725,6/),"float")   
  
     time_o= seismic(:,1)  ; Column 1 of file contains X values.
      p_o = seismic(:,2)  ; Column 2 of file contains Y values.
     t2_o= seismic(:,3)  ; Column 3 of file contains Z values.

这是数据的一部分，缺测值显示为PPC
57131      2015050120  0             161        255        9  
57131      2015050121  0             151        164        9  
57131      2015050122  0             155        142        16
57131      2015050123  0             148        PPC        0  
57131      2015050200  0             142        297        3  


网站上没找到怎么办，求大神帮忙指点

黄小仙儿

刚刚设置了：ob@_FillValue  = PPC

提示有错：
warning:asciiread: End of file reached and only (4343) elements were read from the file, filling remaining elements with the default missing value for the requested type
fatal:Variable (PPC) is undefined
fatal:["Execute.c":8567]:Execute: Error occurred at or near line 116 in file cmaq_xy_2.ncl

longlivehj

当作字符串去处理
str_get_field stringtofloat

liuyanfeng4220

fdsaaaaaaaaaaaaaaaa

黄小仙儿

pangzi_xu 发表于 2015-8-21 17:04
竟然还用cmaq，同行啊，不过数据缺测的问题我一直没有好的办法。看看别人怎么说

我只是想检验一下，画个满意的图好难啊