怎么去掉右边纵坐标的title

601337653

请教一下各位，我用gsn_csm_pres_hgt画图，出来有两个纵坐标，左边是气压，右边是高度，图如下，我想把右边的height（km）去掉要怎么弄呢？试过几个命令，把高度坐标去掉了，可是title还在那边，束手无策啊！！

601337653

问题解决了，找到官网上相关的脚本了！

JHONG

楼主好，我也在画这种时间/高度图，花了很长时间还是没有把图画出来。
用的官网的案例 http://www.ncl.ucar.edu/Applications/Scripts/h_time_4.ncl
在其基础上进行修改的，在打印确认lev存在的情况下，仍报如下错误：
(0)        --------------------------------------
(0)        gsn_csm_pres_hgt: Fatal: The first dimension of the input data must
(0)        have a coordinate variable called 'lev.'
(0)        Cannot create plot.
fatal:Illegal right-hand side type for assignment

不知道您是否也遇到过一样的错误，求指导，谢谢了~~~

601337653

JHONG 发表于 2017-5-23 22:39
楼主好，我也在画这种时间/高度图，花了很长时间还是没有把图画出来。
用的官网的案例 http://www.ncl.uca ...

出现这种问题的原因还蛮多的，你可以把脚本贴一下

JHONG

601337653 发表于 2017-5-24 11:54
出现这种问题的原因还蛮多的，你可以把脚本贴一下

脚本本身比较简单。您看下呢~

; ===============================================
; h_time_4.ncl              
; ===============================================
;
; 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/wrf/WRFUserARW.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"

begin
  fils1 = systemfunc("ls /home/h/wrfchemoutput/tmp/wrfout_d02_*")
; fils1 = systemfunc("ls /home/h/wrfchemoutput/tmp/wrfout_d02_*")
  a  = addfiles(fils1,"r")   ;wrfout file
;-----------------------------------------------------------------------
  ListSetType(a,"cat")

  v    = wrf_user_getvar(a,"PM2_5_DRY",-1)    ; get tc for all times


printVarSummary(v)

dim = dimsizes(v)
data = v(:,:,:,50)

printVarSummary(data)

  slice   = data(south_north|0, bottom_top|:,Time|:)                  ; trick to copy meta data

  wks   = gsn_open_wks ("png", "h_time")           ; send graphics to PNG file
  res                      = True                  ; plot mods desired

  res@cnFillOn             = True                  ; turn on color
; res@cnFillPalette        = "BlWhRe"              ; set color map
  res@cnLinesOn            = False                 ; no contour lines

  res@cnLevelSelectionMode  = "ExplicitLevels"
  tres@cnLevels     = (/10,25, 50, 75, 100, 125, 150,200/)      
  tres@cnFillPalette         = "MPL_jet"
res@tmYRTickSpacingF     = 1                    ;
  
slice!0 = "lev"
slice!1 = "time"


;  print(slice)
  printVarSummary(slice)

  print("--------------------------------------")

  print(slice(lev|:, time|0))

  print("--------------------------------------")

;  terra = gsn_csm_contour (wks,slice,res)   

  plot = gsn_csm_pres_hgt (wks,slice,res)

end

601337653

JHONG 发表于 2017-5-24 13:06
脚本本身比较简单。您看下呢~

; ===============================================

你并没有给lev赋值呀，你只是把高度维用“lev”表示，但系统并不知道你的lev是哪些数

JHONG

601337653 发表于 2017-5-25 13:25
你并没有给lev赋值呀，你只是把高度维用“lev”表示，但系统并不知道你的lev是哪些数

slice!0 = "lev"
仅仅只是将原数组改了个名，数据都还在。

601337653

JHONG 发表于 2017-5-26 14:43
slice!0 = "lev"
仅仅只是将原数组改了个名，数据都还在。

你在把data复制到slice后，第一维还是高度维么

JHONG

601337653 发表于 2017-5-28 10:24
你在把data复制到slice后，第一维还是高度维么

谢谢，我再看下

老司机

请问是怎么去掉啊？我也遇到了这个问题。