- 积分
- 13
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2020-4-29
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2020-5-9 20:35:59
|
显示全部楼层
FNAME = OLOC+'UofTentoSixty.ps'
!P.FONT = 0
ENTRY_DEVICE = !D.NAME
SET_PLOT, 'PS'
DEVICE, FILENAME=FNAME, /COLOR, /HELVETICA, /BOLD, BITS_PER_PIXEL=8, /LANDSCAPE
lon0=0
lon1=360
lat0=5
lat1=65
;map_set,0,180,/mercator,limit=[lat0,lon0,lat1,lon1],position = [0.1, 0.35, 0.9, 0.6];,/noerase,/noborder;,$
index_valid = WHERE(u1 GE 20, count_valid)
min_value = MIN(u1[index_valid])
levels=[19.99,20,25,30]
nlevels=n_elements(levels)
ncolors=nlevels
bottom=1
c_labels=[0,1,1,1]
;c_labels=[0,replicate(1,nlevels)]
c_colors=indgen(ncolors)+bottom
loadct,22,ncolors=ncolors,bottom=bottom
map_set,0,180,/mercator,limit=[lat0,lon0,lat1,lon1],position = [0.1, 0.35, 0.9, 0.6];,/noerase,/noborder;,$
contour,u1,longitude,latitude,levels=levels,c_colors=c_colors,/fill,/overplot
contour,u1,longitude,latitude,levels=levels,c_labels=c_labels,c_thick=1,/overplot
;levels=[-30,-25,-20]
;nlevels=n_elements(levels)
;ncolors=nlevels
;c_labels=[1,1,1]
;c_colors=indgen(ncolors)
;loadct,21,ncolors=ncolors
;c_labels=[1,1,1]
;c_colors=250
;contour,u1,longitude,latitude,levels=levels,c_colors=c_colors,/overplot,/fill
;contour,u1,longitude,latitude,levels=levels,c_labels=c_labels,c_thick=1,/overplot
map_continents,/continents,/countries
;print,max(u1),min(u1)
map_grid,lons =[0,30,60,90,120,150,180,210,240,270,300,330,360 ],lats = [10, 25, 40, 55 ],londel = 13, latdel = 3,/box;,charsize = 1.2
device, /close_file
set_plot, entry_device
end |
|