- 积分
- 42
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2019-1-21
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2020-4-26 13:43:28
|
显示全部楼层
begin
y=addfile("/Users/ms/Documents/albedo8/tavg.02025.01.01.nc","r")
O_temp=y->O_temp
time1=y->time
latitude1=y->latitude
longitude1=y->longitude
temp1=O_temp(15,0,:,:) ;2100年,试验组的0层海洋温度,0至-50m
temp2=O_temp(15,4,:,:) ;2100年,试验组的4层海洋温度,-380至-550m
temp3=O_temp(15,7,:,:) ;2100年,试验组的7层海洋温度,-980至-1240m
yy=addfile("/Users/ms/Documents/albedo8/rcp85_emi/tavg.01805.01.01.nc", "r")
O_temp1=yy->O_temp
time2=yy->time
latitude2=yy->latitude
longitude2=yy->longitude
temp11=O_temp1(59,0,:,:) ; 2100年,对照组的0层海洋温度
temp12=O_temp1(59,4,:,:)
temp13=O_temp1(59,7,:,:)
temp21=O_temp(15,0,:,:)-O_temp1(59,0,:,:)
temp22=O_temp(15,4,:,:)-O_temp1(15,4,:,:)
temp23=O_temp(15,7,:,:)-O_temp1(15,7,:,:)
wks=gsn_open_wks("png", "OARCP_oceantemp")
gsn_define_colormap(wks, "NCV_blue_red")
res=True
res@gsnDraw = False
res@gsnFrame = False
res@gsnLeftString=" "
res@gsnRightString=" "
res@cnFillOn = True
res@cnLinesOn = False ;等值线不显示
res@cnLineLabelsOn = False ;
res@lbLabelBarOn=False
res@cnLevelSelectionMode="AutomaticLevels"
res@cnSmoothingOn = True
res@cnSmoothingTensionF = 0
res@cnSmoothingDistanceF=0.0000000000000000000000000000000000000000000000000000000000000000000000000001
res@lbTitleString="~R~[~S~o~N~C~R~]"
res@lbTitleFontHeightF=0.02
res@pmLabelBarOrthogonalPosF=0.15
res@cnFillDrawOrder = "PreDraw"
pres=True
pres@gsnMaximize = True ; Maximize plots in panel.
pres@gsnPanelLabelBar = True ; Turn on panel labelbar.
pres@pmLabelBarWidthF = 0.8 ; Change labelbar width.
pres@lbTitleString="~R~[~S~o~N~C~R~]"
pres@lbTitleFontHeightF=0.01
pres@pmLabelBarOrthogonalPosF=0.00000000001
pres@gsnPanelYWhiteSpacePercent=2
plot=new(3,graphic)
plot(0)=gsn_csm_contour_map(wks, temp21, res)
plot(1)=gsn_csm_contour_map(wks, temp22, res)
plot(2)=gsn_csm_contour_map(wks, temp23, res)
gsn_panel(wks,plot, (/1,3/), pres)
end |
|