[求助][已解决]xy plot图中，x轴自定义数组值怎么替换进去

xuebiz

画某点的时间序列图，数据文件中有一项定义好的时间“yyyymmddhh”类型是char

画图时我定义了一个x轴数组：  t = ispan(0,64,1)
用什么resource属性，可以让x轴的标注为time数组，而不是1,2，...呢

数据文件中的定义：
dimensions:
    time = UNLIMITED ; // (65 currently)
    lon = 720 ;
    lat = 361 ;
    char = 10 ;
variables:
    float lon(lon) ;
        lon:long_name = "Longitude" ;
        lon:units = "degrees_east" ;
     float lat(lat) ;
        lat:long_name = "Latitude" ;
        lat:units = "degrees_north" ;
     char time(time, char) ;
        time:long_name = "Time" ;
        time:units = "Gregorian calendar ,UTC" ;
        time:format = "yyyymmddhh" ;

喵星教授Mamonde

用这个res@tmXBMode， res@tmXBValues，res@tmXBLabels试试呢？

xuebiz

谢谢，查到的类似用法。。。

  res2@tmXBMode        = "Explicit"   
  res2@tmXBValues      = (/114,115,116,117,118,119,120/)  
  res2@tmXBLabels   = (/"114~S~o~N~E","115~S~o~N~E","116~S~o~N~E","117~S~o~N~E",\  
  "118~S~o~N~E","119~S~o~N~E","120~S~o~N~E"/)  

xuebiz

  ts = f->time
  t = ispan(0,64,1)

  res@tmXBMode = "Explicit"
  res@tmXBValue = t
  res@tmXBLables = ts

报错如下：
warning:tmXBValue is not a valid resource in gfs-wind-timeq_xy at this time
warning:tmXBLables is not a valid resource in gfs-wind-timeq_xy at this time

难道要让我把ts（0，*）自己组合成一个新数组么？

看其他变量的定义中time也正常使用：
short temp(time, lat, lon)

-------------------------------------------------------------------------

我ncdump -v time 的结果是这样子的：
data:
time =
  "2015031000",
  "2015031003",
  "2015031006",


但是我在NCL里printVarSummarg的结果是这样子的：
Variable: ts
Type: character
Total Size: 650 bytes
            650 values
Number of Dimensions: 2
Dimensions and sizes:        [time | 65] x [char | 10]
Coordinates:
Number Of Attributes: 3
  long_name :        Time
  units :        Gregorian calendar ,UTC
  format :        yyyymmddhh

xuebiz

我做了一个strT：

  strT = new(65,string)
  do n=0,64,1      
      strT(n)=str_concat(charactertostring(ts(n,:)))
  end do
  print(strT)

可是还是出这个错：
warning:tmXBValue is not a valid resource in gfs-wind-timeq_xy at this time
warning:tmXBLables is not a valid resource in gfs-wind-timeq_xy at this time
报错不怕，可X轴还是1,2，3。。。

print（strT）的结果：

Variable: strT
Type: string
Total Size: 520 bytes
            65 values
Number of Dimensions: 1
Dimensions and sizes:        [65]
Coordinates:
Number Of Attributes: 1
  _FillValue :        missing
(0)        2015031000
(1)        2015031003
(2)        2015031006

longlivehj

贴出全部脚本

xuebiz

longlivehj 发表于 2016-4-25 17:01
贴出全部脚本

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"

begin

  f = addfile("../data/wind.2015031000.nc", "r")

  wks = gsn_open_wks("png","wind-timeq")
  ts= f->time
  t = ispan(0,64,1)
  u = short2flt(f->uwnd(:,{30},{130}))
  v = short2flt(f->vwnd(:,{30},{130}))
;  printVarSummary(ts)

  wind = sqrt(u*u+v*v)
  
  strT = new(65,string)
  do n=0,64,1      
      strT(n)=str_concat(charactertostring(ts(n,:)))
  end do
  print(strT)

  res = True  
  res@tiMainString = "wind filetime=2015031000"
  res@tiXAxisString = "DateTime"
  res@tiYAxisString = "Wind"
  
  res@tmXBMode = "Explicit"
  res@tmXBValue = t
  res@tmXBLables = strT
  
  plot = gsn_xy(wks, t, wind, res)

end

longlivehj

tmXBValues
tmXBLabels
注意拼写！

xuebiz

longlivehj 发表于 2016-4-25 17:29
tmXBValues
tmXBLabels
注意拼写！

被自己蠢死~~
OK了
不过下面一团浆糊了，我要改标签间隔和斜体去了
谢谢~~大版主

xuebiz

longlivehj 发表于 2016-4-25 17:01
贴出全部脚本

  strT = new(65,string)
  do n=0,64,1      
      strT(n)=str_concat(charactertostring(ts(n,:)))
  end do
各位达人，版主大人，这里妥不妥，是不是有更好的办法~~

----------------------------------------------------------------------------
char time(time, char)
short temp(time, lat, lon)

我用ncdump出来time的结果是yyyymmddhh，IDL直接用time也可以，NCL就不行。但是temp等数组里用起来也没有问题。。。