爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 5879|回复: 0

[作图] 用GSview打开ps文件,在一个工作站怎么让它翻页?

[复制链接]

新浪微博达人勋

发表于 2015-10-27 11:37:04 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
原数据来自不同的nc文件,日期不同,但数据的格式一致,想知道其中的变量名是不是都需要修改?load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"

begin
  dir = "E:/gekco/"
  fin1  = addfile(dir+"v_gekco2_19980601.nc","r")
;  print(fin)
;  uf = addfile(dir+"/cdf/Ustorm.cdf","r")  ; Open two netCDF files.
;  vf = addfile(dir+"/cdf/Vstorm.cdf","r")

  u1   = fin1->u
  v1   = fin1->v
  lat1 = fin1->lat
  lon1 = fin1->lon

  wks = gsn_open_wks("ps",dir+"1998Jun1-5") ; Open a workstation.

  u1@_FillValue = 0
  v1@_FillValue = 0

  res1 = True
  nlon1 = dimsizes(lon1)
  nlat1 = dimsizes(lat1)
  res1@pmTickMarkDisplayMode = "Always"; use NCL default lat/lon labels
  res1@gsnAddCyclic = False
  ;res@mpLimitMode = "LatLon"
  res1@mpMinLonF   = 120
  res1@mpMaxLonF     = 270
  res1@mpMinLatF     = 0
  res1@mpMaxLatF     = 70
  res1@mpCenterLonF  = 195
  res1@vfXCStartV  = lon1(0)      ; Define lat/lon corners
  res1@vfXCEndV    = lon1(nlon1-1) ; for vector plot.
  res1@vfYCStartV  = lat1(0)
  res1@vfYCEndV    = lat1(nlat1-1)

  res1@gsnFrame   = False            ; so we can draw time stamp
  res1@vcRefAnnoOrthogonalPosF = -1.0             ; move ref vector up
  res1@vcMinDistanceF = 0.017
  res1@vcRefMagnitudeF = 20          ; define vector ref mag
  res1@vcRefLengthF    = 0.01           ; define length of vec ref
  res1@vcGlyphStyle  = "CurlyVector"    ; turn on curly vectors

  res1@tiMainString  = "Ocean Surface Currents Velocity"
  res1@gsnLeftString    = "19980601"
  ;res@gsnRightString   = "units: meter/sec"
  ;res@vcGlyphStyle  ="LineArrow"  

  map1 = gsn_csm_vector_map_ce(wks,u1(0,:,:),v1(0,:,:),res1)  ; Draw a vector plot of
                                                              ; u and v over a map.
frame(wks)

fin2   = addfile(dir+"v_gekco2_19980602.nc","r")

  u2   = fin2->u
  v2   = fin2->v
  lat2 = fin2->lat
  lon2 = fin2->lon

  u2@_FillValue = 0
  v2@_FillValue = 0

  res2 = True
  nlon2 = dimsizes(lon2)
  nlat2 = dimsizes(lat2)
  res2@pmTickMarkDisplayMode = "Always"; use NCL default lat/lon labels
  res2@gsnAddCyclic = False
  ;res@mpLimitMode = "LatLon"
  res2@mpMinLonF   = 120
  res2@mpMaxLonF     = 270
  res2@mpMinLatF     = 0
  res2@mpMaxLatF     = 70
  res2@mpCenterLonF  = 195
  res2@vfXCStartV  = lon2(0)      ; Define lat/lon corners
  res2@vfXCEndV    = lon2(nlon2-1) ; for vector plot.
  res2@vfYCStartV  = lat2(0)
  res2@vfYCEndV    = lat2(nlat2-1)

  res2@gsnFrame   = False            ; so we can draw time stamp
  res2@vcRefAnnoOrthogonalPosF = -1.0             ; move ref vector up
  res2@vcMinDistanceF = 0.017
  res2@vcRefMagnitudeF = 20          ; define vector ref mag
  res2@vcRefLengthF    = 0.01           ; define length of vec ref
  res2@vcGlyphStyle  = "CurlyVector"    ; turn on curly vectors

  res2@tiMainString  = "Ocean Surface Currents Velocity"
  res2@gsnLeftString    = "19980602"
  ;res@gsnRightString   = "units: meter/sec"
  ;res@vcGlyphStyle  ="LineArrow"  

  map2 = gsn_csm_vector_map_ce(wks,u2(0,:,:),v2(0,:,:),res2)  ; Draw a vector plot of
                                                              ; u and v over a map.
frame(wks)

end
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表