pro perf,ps=ps
;
;
if (keyword_set(ps)) then begin
  set_plot,'PS'
endif else begin
  set_plot,'X'
endelse
;
proc         = [1.  , 4. , 8. , 16., 24., 32.] ; use with secs/step
nproc = n_elements(proc)
;
; Bluevista:
;
time_bv = [81.8, 3.89, 2.38, 1.66, 1.51, 1.48] ; sec/step
ytitle="Ave sec/step"
yrange=[1.,7.]
plot,proc,time_bv,xtitle="ntask (1,4,8,16,24,32)",ytitle=ytitle,$
  yrange=yrange,position=[.1,.2,.9,.9],/norm,$
  title="TIEGCM 1-day simulation (hres=2.5, vres=0.25, step=90s)"
xyouts,proc[nproc-1]+1.,time_bv[nproc-1],"Bluevista"
;
; Bluesky:
;
time_bs = [16.7, 6.46, 4.25, 3.18, 3.25, 2.30] ; sec/step
oplot,proc,time_bs,linestyle=2
xyouts,proc[nproc-1]+1.,time_bs[nproc-1],"Bluesky"
;
; Lightning:
;
time_ln = [16.4, 5.94, 4.59, 3.67, 4.01, 4.32] ; sec/step
oplot,proc,time_ln,linestyle=3
xyouts,proc[nproc-1]+1.,time_ln[nproc-1],"Lightning"
;
if (keyword_set(ps)) then begin
  device,/close
endif
;
end