- 积分
- 21
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2017-8-22
- 最后登录
- 1970-1-1
|

楼主 |
发表于 2019-8-30 13:27:27
|
显示全部楼层
原代码如下:
begin
file1 = systemfunc("ls ../out_data/temp_diff_150y.nc")
in = addfile(file1,"r")
temp =in ->TEMP_DIFF ;(150,50)
zt =in ->depth
zt2 =zt(::-1)
temp1 = transpose(temp)
wks = gsn_open_wks("png","diff_21_170y")
res = True
res@xyYStyle = "Irregular"
res@xyYIrregularPoints = (/200,400,600,800,1000,2000,3000,4000,5000/)
res@gsnMaximize = True
res@cnFillOn = True ;turn on color fill
res@cnFillPalette = "BlueYellowRed" ;color map
res@cnLineLabelsOn = False ;contour line labels
res@cnLinesOn = True ;contour lines
res@cnInfoLabelOn = False ;contour info label
res@lbOrientation = "vertical" ;vertical label bar
res@sfXArray = temp&year ;xlabel
;res@sfYArray = temp&depth ;ylabel
;res@tiXAxisString = "Model Year";xaxis name
;res@tiYAxisString = "Depth(meters)" ;yaxis name
;res@trYReverse = True
plot = gsn_csm_contour(wks,temp1,res) ; create PLOT
end
|
|