NCL 时间-经度 剖面图 时间转换

朝夕木兰

ncl 画时间-经度图，大家可能会和我一样碰到时间坐标转换的问题。将时间转换，并且显示在图上，主要需要两部分命令。
一是加载  load "$NCARG_ROOT/lib/ncarg/nclscripts/contrib/time_axis_labels.ncl"
二是坐标轴转换
restick = True
restick@ttmFormat = "%Y"
restick@ttmAxis   = "YL"
restick@ttmMajorStride = 30
time_axis_labels( var&time,  res, restick )
代码贴出来供大家参考。

1. ; ===========================================
   
2. ; hov.ncl
   
3. ; ===========================================
   
4. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
   
5. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
   
6. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
   
7. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
   
8. load "$NCARG_ROOT/lib/ncarg/nclscripts/contrib/time_axis_labels.ncl"
   
9. ; ===========================================
   
10. begin
    
11. f     = addfile("sst_ff.nc", "r")
    
12. var   = f->sst(:,10,:)         
    
13. ;time  =  cd_calendar(f->time,-1)
    
14. ;print (time)
    
15. ;==============================================
    
16. ; create color plot
    
17. ;=============================================
    
18. wks  = gsn_open_wks ("x11", "sst" )            ; open ps file
    
19. gsn_define_colormap(wks,"BlWhRe")              ; choose colormap
    
20. res                      = True                ; plot mods desired
    
21. res@cnFillOn             = True               ; turn on color fill
    
22. res@trYReverse           = True                ; reverses y-axis
    
23. ;------------------time convert for axis-----------------
    
24. restick = True
    
25. restick@ttmFormat = "%Y"
    
26. restick@ttmAxis   = "YL"
    
27. restick@ttmMajorStride = 30
    
28. time_axis_labels( var&time,  res, restick )
    
29. ;---------------- modulate details -----------------------
    
30. res@gsnStringFont = "helvetica"
    
31. res@tmXBLabelFont = "helvetica"
    
32. res@tmYLLabelFont = "helvetica"
    
33. 
    
34. res@gsnSpreadColors      = True               ; use full range of colors
    
35. 
    
36. res@lbOrientation        = "vertical"
    
37. res@lbLabelFont          = "helvetica"
    
38. res@lbLabelFontHeightF   = 0.02
    
39. 
    
40. res@tiYAxisString = "Year"
    
41. res@tiYAxisFont   = "helvetica"
    
42. ;-----------------------plot-----------------------------
    
43. plot = gsn_csm_hov(wks, var, res)
    
44. end
    

复制代码

朝夕木兰

tsl08

把时间序列倒过来是不是更好一些？

朝夕木兰

tsl08 发表于 2015-10-27 09:20
把时间序列倒过来是不是更好一些？

我自己的分析需要，就用res@trYReverse           = True                ; reverses y-axis   把时间倒过来了。

jiaohe1990

赞  赞  正需要  

清柳雨菡

谢谢楼主呀~~~~~~

novelbean

不错，学习学习

subtropical

学习学习谢谢楼主

subtropical

楼主，如果纬度不是确定的一个，而是一段范围怎么设置呢？求教

朝夕木兰

subtropical 发表于 2016-1-17 14:40
楼主，如果纬度不是确定的一个，而是一段范围怎么设置呢？求教

抱歉，现在才看到。
一般时间－经度剖面图是针对某一个纬度的。比如想要南纬45度的，就把十二行改成：
var1  = f->sst(:,{-45},:)
你说的一个纬度范围，应当是要对这个范围内的数据作处理再呈现其剖面图（比如之后再平均）
var2  = f->sst(:,{-45:-60},:)选定了南纬45～60度范围