ncl画图官网实例zonal——4.ncl实现

wjjjjyyyyy

已于 2023-04-13 16:20:26 修改

阅读量421

点赞数 1

文章标签：笔记

于 2023-04-13 16:10:27 首次发布

本文链接：https://blog.csdn.net/m0_70443361/article/details/130132616

版权

文章是关于如何处理NCL(zonal_4.ncl)代码运行时因缺少atmos.nc文件导致的问题。解决方案是到指定官网下载文件并确保与代码在同一目录，或者正确指定文件路径。此外，文章还提到了如何修改代码以在本地显示图形结果，以及代码涉及的变量处理和等压面计算。最后，展示了如何使用NCL进行等区平均图的绘制。

摘要由CSDN通过智能技术生成

这是是官方代码官网zonal_4.ncl代码

出现的问题：

atmos.nc打不开

解决方法：由于atmos.nc文件缺失导致的，

去官网这个地方Data files for some of the NCL examples (ucar.edu)找到所需要的文件下载即可

下载的所需文件一定要与ncl在同一个文件夹下，否则要写清楚路径

运行成功但是没有出来图像，生成了png类型的文件，可以去文件夹里找到

要是想直接看到画出的结果可以将官方代码中的

wks = gsn_open_wks("png","zonal")替换为了

wks = gsn_open_wks("x11","zonal")

;*************************************************
; zonal_4.ncl
;
; Concepts illustrated:
;   - Drawing a zonal means plot
;   - Customzing the axes of a zonal means plot
;************************************************
;
; These files are loaded by default in NCL V6.2.0 and newer
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
;************************************************
begin
;************************************************
; variable and file handling
;************************************************
  in  = addfile("atmos.nc","r")
  T   = in->T                                     ; select variable to ave
;************************************************
;  calculate zonal average of 3d variable
;************************************************
  zave = dim_avg_Wrap(T)                             ; calculate zonal ave
;************************************************
;  calculate eta
;************************************************
  a=in->hyam                                    ; select hyam 
  b=in->hybm                                    ; select hybm
  p=in->P0                                      ; select P0
  eta = (a+b)*p                                 ; calc eta
  eta = eta/100                                 ; scale eta by 100
;************************************************
; plotting
;************************************************
  wks  = gsn_open_wks("x11","zonal")              ; send graphics to PNG file
  
  res               = True
  res@tiXAxisString = "latitude"                  ; x-axis label
  res@tiYAxisString = "eta x 1000"                ; y-axis label
  res@sfXArray      = zave&lat                    ; put lat on x axis
  res@sfYArray      = eta                         ; put eta on y axis
  res@trYReverse    = True                        ; reverses y axis
  res@trXReverse    = True                        ; reverse X axis
  
  plot=gsn_contour(wks,zave(0,:,:),res)           ; plot zonal ave
;***********************************************
end

运行成功