;$Id: pspec_save.pro,v 1.3 2020/12/28 16:03:37 brandenb Exp $
if !d.name eq 'PS' then begin
  ;device,xsize=18,ysize=8,yoffset=3
  device,xsize=18,ysize=11,yoffset=3
  !p.charthick=1.6 & !p.thick=1.6 & !x.thick=1.6 & !y.thick=1.6
end
;
;  mv idl.ps ~/tex/roper/GW/fig/pspecm_scl2.ps
;  mv idl.ps ~/tex/prop/computing/pdc/18/fig/pspecm_scl3n.ps
;
@parameters
!p.charsize=2.0
!x.margin=[8.2,0.5]
!y.margin=[3.1,0.4]
!p.multi=[0,1,2]
siz=1.
si2=.9
;
default,i3,1
default,regime,2
default,narrow,1
print,"$sed.csh data/param2.nml"
pc_read_param,obj=param
;
default,iread,0
if iread eq 0 then begin
  power,'_GWh','_GWs',k=k,spec1=grah1,spec2=grav1,i=n,tt=t,/noplot,/lks
  power,'_mag','_Str',k=k,spec1=spec1,spec2=stre1,i=n,tt=t,/noplot,/lks
  iread=1
endif
default,fact1,1.
default,fact2,1.
;
!x.title='!8k!6'
!y.title='!8E!6!dGW!n(!8k,t!6)'
if narrow then yr=[1e-27,1e-16]*fact1 else yr=[1e-37,1e-6]*fact1
plot_oo,k[1:*],k(1:*),yr=yr,/nodata
for i=i1,i2,i3 do begin
  oplot,k[1:*],grav1(1:*,i)
endfor
xx=[120.,1400.] & oplot,xx,1e-18*(xx/100.)^2
xx=[120.,1400.] & oplot,xx,1e-23*(xx/100.)^2
xx=[120.,1400.] & oplot,xx,1e-28*(xx/100.)^2
xx=[120.,1400.] & oplot,xx,1e-38*(xx/100.)^2
xx=[120.,7400.] & oplot,xx,1e-13/(xx/100.)^2
xx=[120.,7400.] & oplot,xx,1e-23/(xx/100.)^2
;
;  panel 4
;
!y.title='!8S!dh!6!n(!8k,t!6)'
if narrow then yr=[1d-34,1e-20]*fact2 else yr=[1d-44,1e-10]*fact2
plot_oo,k[1:*],k(1:*),yr=yr,/nodata
for i=i1,i2,i3 do begin
  oplot,k[1:*],1d0*grah1(1:*,i)
endfor
xx=[120.,1400.] & oplot,xx,1e-18*(xx/100.)^2
xx=[120.,1400.] & oplot,xx,1e-28*(xx/100.)^2
xx=[120.,1400.] & oplot,xx,1e-38*(xx/100.)^2
xx=[120.,7400.] & oplot,xx,1e-13/(xx/100.)^2
xx=[120.,7400.] & oplot,xx,1e-23/(xx/100.)^2
;
grav1_range=grav1(*,i1:i2)
grah1_range=grah1(*,i1:i2)
spec1_range=spec1(*,i1:i2)
stre1_range=stre1(*,i1:i2)
save,file='spec_range.sav',k,grav1_range,grah1_range, spec1_range, stre1_range
;cwd,run
;print,'$mv idl.ps ~/tex/rei/GW/fig/pspec_save_'+run+'.ps'
END