;$Id: pgravt.pro,v 1.1 2024/12/16 10:55:21 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=12,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
end
;
!p.charsize=1.7
!x.margin=[10.8,.5]
!y.margin=[3.2,1.9]
;
default,lread,1
if lread then begin
  power,'hel_GWs','hel_GWh',k=k,spec1=grav2,spec2=grah2,i=n,tt=t,/noplot,/lks
  lread=0
endif
;
cwd,run
!x.title='!8t!6-!8t!6!d*!n'
!y.title="!8d/dt!6(!8A!dh'!n-!8k!6!u2!nA!d!8h!6!n)"
!p.title='!6'+run
nt=n_elements(t)
k2=spread(k^2,[1],[nt])
plot_oi,t-1,deriv(t,total(grav2-k2*grah2,1))
;oplot,t,total(grah2[*,*],1),col=122
;
print,"$mv idl.ps ~/Overleaf/Axel/Azi-pol/fig/pgravt_"+run+".eps"
END