| 
 
	积分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
 | 
 |