- 积分
- 25737
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2014-3-31
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
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
|
|