已解决：NCL绘制wrfout 垂直剖面横坐标转换和设置

lsucky

小白一枚，最近在学习ncl，遇到的两个问题在大家的帮助下已解决，总结如下：
（1）将横坐标由格点数转换为经度表示，类似于下图：

添加经纬度信息，补全的脚本如下：
  ; Interpolate to pressure
    av_plane = wrf_user_intrp3d(avo,p,"v",plane,0.,False)
    p_plane = wrf_user_intrp3d(p, p,"v",plane, 0.,False)
   lon_plane = wrf_user_intrp2d(xlon,plane,90.,False)
  
    av_plane!0="level"
    av_plane!1="lon"
    av_plane&level = p_plane(:, 0)
    av_plane&lon = lon_plane
    ; MAKE PLOTS                                       
     plot  = gsn_csm_pres_hgt(wks,av_plane,res )

（2）坐标间隔设置问题,指定起点终点和间隔：
      res@tmXBMode = "Explicit"
      res@tmXBValues = fspan(xxx,xxx,4)
      res@tmXBLabels = (/"xxx", "xxx",  "xxx", "xxx"/)

最后再次感谢家园的乐于助人的小伙伴！比心~

lsucky

pandasp1213 发表于 2017-1-7 09:36
把经度坐标信息补上就行了
av_plane!1="lon"
av_plane&lon = xxx

按照您说的尝试了下，显示错误如下：
fatal:Coordinate variables must be the same dimension as their dimension
fatal:["Execute.c":8575]:Execute: Error occurred at or near line 74 in file p-av.ncl

修改代码如下：

pandasp1213

把经度坐标信息补上就行了
av_plane!1="lon"
av_plane&lon = xxx

pandasp1213

lsucky 发表于 2017-1-8 10:56
按照您说的尝试了下，显示错误如下：
fatal:Coordinate variables must be the same dimension as the ...

从你图上看经度方向应该有900个点左右，113.63,113.66两个肯定是不够的呀，所以提示must be the same dimension，得给变量900个经度

lsucky

修改了维度已经初步得到解决了，非常感谢！

lsucky

pandasp1213 发表于 2017-1-8 14:36
从你图上看经度方向应该有900个点左右，113.63,113.66两个肯定是不够的呀，所以提示must be the same dim ...

试过修改维度，问题已经基本上解决了，非常感谢！
现在的图的坐标已经是经度了，
顺便再求教一下，横坐标的精度怎么调节？
像上图中只出现一个113，我想让它的间隔为0.1显示。
试过            
      ;res@tmXBValues    = fspan(xxx, xxx, xx)   但并没有成功

愛、星晴

lsucky 发表于 2017-1-8 18:20
试过修改维度，问题已经基本上解决了，非常感谢！
现在的图的坐标已经是经度了，
顺便再求教一下，横 ...

添加以下语句：
res@tmXBMode                = "Explicit"
res@tmXBValues              = fspan(101,105,5)                    ; Create tick marks
res@tmXBLabels              = (/"101E","102E","103E","104E","105E"/)  ; Create labels

lsucky

愛、星晴 发表于 2017-1-8 19:26
添加以下语句：
res@tmXBMode                = "Explicit"
res@tmXBValues              = fspan(101, ...

谢谢！但是是这样的，资料的精度很小，经度范围在112-113.8之间，
function fspan （        start  [1] : numeric,        
                finish [1] : numeric,        
                npts   [1] : integral type   )
第三项必须是整数，所以应该是不能用这个

愛、星晴

lsucky 发表于 2017-1-8 20:24
谢谢！但是是这样的，资料的精度很小，经度范围在112-113.8之间，
function fspan （        start  [1] : nume ...

怎么不可以用，查查这个函数的意思你就知道怎么回事，当然，还有一种方法就是你把想要的刻度都列出来。

lsucky

愛、星晴 发表于 2017-1-8 21:05
怎么不可以用，查查这个函数的意思你就知道怎么回事，当然，还有一种方法就是你把想要的刻度都列出来。

fatal:fspan: number of elements parameter must be of an integral type, cannot continue
刚刚试了一下，确实是不可以呢