- 积分
- 42203
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-3-5
- 最后登录
- 1970-1-1
|
发表于 2017-8-8 04:57:06
|
显示全部楼层
本帖最后由 royalosyin 于 2017-8-9 05:13 编辑
用NCL就可以的,极地投影。
下面是主要的code,可以参照。
;---Start the graphics
wks_name = "Arctic_SIC_" + token + "_" + tostring(thisyear) + sprinti("%0.2i", thismonth)
wks_type = "png"
wks_type@wkWidth = 960
wks_type@wkHeight = 960
wks = gsn_open_wks(wks_type, wks_name)
gsn_define_colormap(wks,colorBounds)
res = True
res@gsnPolar = "NH"
res@gsnMaximize = True ; maximize plot in frame
res@gsnDraw = False ; don't draw plot yet
res@gsnFrame = False ; don't advance frame yet
res@cnFillOn = True
res@cnFillMode = "RasterFill" ; Raster Mode
res@cnRasterSmoothingOn = True
;res@mpGridMaskMode = "MaskNotOcean"
res@cnLineLabelsOn = False
res@cnLinesOn = False
;---Make sure both plots have same contour levels
res@cnLevelSelectionMode = "ExplicitLevels"
res@cnLevels = (/-50.0, -40.0, -30.0, -20.0, -10.0, -5.0, 0.0, 5.0, 10.0, 20.0, 30.0, 40.0, 50.0/)
res@lbLabelBarOn = False
res@lbOrientation = "vertical"
res@lbLabelJust = "CenterLeft" ; left justify labels
res@gsnAddCyclic = True
res@mpPerimOn = False
res@gsnTickMarksOn = True ; dow show box of the draw zone
res@mpGridAndLimbOn = True
res@gsnStringFontHeightF = 0.010
res@gsnPolarLabelFontHeightF = 0.007
res@gsnBoxMargin = 0.0;
res@gsnRightString = ""
res@gsnLeftString = ""
res@gsnCenterString = strmons(thismonth-1) + " " + tostring(thisyear) + " Sea Ice Concentration Anomaly"
res@gsnCenterStringFontHeightF = 0.010
res@mpDataBaseVersion = "MediumRes" ; Better map resolution
res@mpDataSetName = "Earth..4"
res@mpLandFillColor = "tan"
res@mpOceanFillColor = "White"
res@mpInlandWaterFillColor = "White"
res@mpMinLatF = 55.5 ; Specify min lat ; Specify min lat!
;---Create plot
map_data = gsn_csm_contour_map_polar(wks,vals_One,res)
draw(map_data);
frame(wks)
|
-
|