;
pro plot_osc1
;
phycon
;
;save, filename = sav_filename, aTkP, aUP, aTkA, azP, asec, adUpdt, gamma_A, TkAo, ho
sav_filename = 'osc1.sav'
;
restore, filename = sav_filename, /verbose
;
az_plot    = [ 0.,     1500.] ; height relative to base
aTkA_plot  = [TkAo + 1000.*gamma_A, TkAo - 500.*gamma_A] ; 1000 m below to 500 m above
;
;..set up the colors
;
tvlct, [   0, 255,   0,   0, 255], $
  [   0,   0,   0, 255,   0], $
  [   0,   0, 255,   0, 255]
red     = 1
blue    = 2
magenta = 4
green   = 3
black   = 0
;
;..setup the postscript file
;
xsize_inch            =  8.5
ysize_inch            =  11.
;
width_inch            =  1.25
width_normal          = width_inch / xsize_inch
height_inch           =  4.0
height_normal         = height_inch / ysize_inch
height_spacer_inch    =  0.25
height_spacer_normal  = height_spacer_inch / ysize_inch
width_spacer_inch     =  0.45
width_spacer_normal   = width_spacer_inch / xsize_inch
yoffset_inch          =  0.5
yoffset_normal        = yoffset_inch / ysize_inch
xoffset_inch          =  1.0
xoffset_normal        = xoffset_inch / xsize_inch
;
set_plot, 'ps'
!p.font    =    0
device, /times, /inches, /color, font_size = 12, filename = 'osc1.ps', $
 yoffset = 0, xoffset = 0., xsize = xsize_inch, ysize = ysize_inch
;
ix                    = 1
iy                    = 1
x_coord_min_normal    = xoffset_normal + (ix-1)*(width_spacer_normal + width_normal)
x_coord_max_normal    = x_coord_min_normal + width_normal
y_coord_min_normal    = 1.0 - yoffset_normal - iy*height_normal - (iy-1)*height_spacer_normal
y_coord_max_normal    = y_coord_min_normal + height_normal
!p.position           = [x_coord_min_normal,y_coord_min_normal,x_coord_max_normal,y_coord_max_normal]
;
xyouts, x_coord_min_normal, y_coord_max_normal + 0.01, sav_filename, /normal
;
xrange = [7., 17.]
xticks = 2
xminor = 5
yrange = [0., 1500.]
yticks = 3
yminor = 5
;
plot, /nodata, /noerase, [0], ytitle='Height Relative to Base, m', xstyle=1, xrange=xrange, $
 xticks=xticks, xminor=xminor, xtitle='T!DP!N and T!DA!N, !Uo!NC', ystyle=1, yrange=yrange, yticks=yticks, $
 yminor=yminor, xticklen = -0.02, yticklen = -0.04
for i = 0, n_elements(aTkP) - 1, 10 do begin
 plots, aTkP[i] - !Tko, azP[i], psym = 4, symsize = 0.10, color = red
endfor
oplot, aTkA_plot - !Tko, az_plot, linestyle = 1
;
ix                    = 2
iy                    = 1
x_coord_min_normal    = xoffset_normal + (ix-1)*(width_spacer_normal + width_normal)
x_coord_max_normal    = x_coord_min_normal + width_normal
y_coord_min_normal    = 1.0 - yoffset_normal - iy*height_normal - (iy-1)*height_spacer_normal
y_coord_max_normal    = y_coord_min_normal + height_normal
!p.position           = [x_coord_min_normal,y_coord_min_normal,x_coord_max_normal,y_coord_max_normal]
;
xrange = [-3., 3.]
xticks = 2
xminor = 3
yrange = [0., 1500.]
yticks = 3
yminor = 5
;
plot, /nodata, /noerase, [0], xstyle=1, xrange=xrange, $
  xticks=xticks, xminor=xminor, xtitle='U!DP!N, m s!U-1!N', ystyle=1, yrange=yrange, yticks=yticks, $
  yminor=yminor, xticklen = -0.02, yticklen = -0.04, ytickname = replicate(' ', yticks + 1)
for i = 0, n_elements(aTkP) - 1, 10 do begin
 plots, aUP[i], azP[i], psym = 4, symsize = 0.10, color = red
endfor
;
ix                    = 3
iy                    = 1
x_coord_min_normal    = xoffset_normal + (ix-1)*(width_spacer_normal + width_normal)
x_coord_max_normal    = x_coord_min_normal + width_normal
y_coord_min_normal    = 1.0 - yoffset_normal - iy*height_normal - (iy-1)*height_spacer_normal
y_coord_max_normal    = y_coord_min_normal + height_normal
!p.position           = [x_coord_min_normal,y_coord_min_normal,x_coord_max_normal,y_coord_max_normal]
;
xrange = [-0.03, 0.03]
xticks = 2
xminor = 3
yrange = [0., 1500.]
yticks = 3
yminor = 5
;
plot, /nodata, /noerase, [0], xstyle=1, xrange=xrange, $
  xticks=xticks, xminor=xminor, xtitle='dU!DP!N/dt, m s!U-2!N', ystyle=1, yrange=yrange, yticks=yticks, $
  yminor=yminor, xticklen = -0.02, yticklen = -0.04, ytickname = replicate(' ', yticks + 1)
for i = 0, n_elements(aTkP) - 1, 10 do begin
 plots, adUpdt[i], azP[i], psym = 4, symsize = 0.10, color = red
endfor
;
width_inch            =  7.0
width_normal          = width_inch / xsize_inch
height_inch           =  1.0
height_normal         = height_inch / ysize_inch
height_spacer_inch    =  0.25
height_spacer_normal  = height_spacer_inch / ysize_inch
width_spacer_inch     =  0.5
width_spacer_normal   = width_spacer_inch / xsize_inch
yoffset_inch          =  5.5
yoffset_normal        = yoffset_inch / ysize_inch
xoffset_inch          =  1.0
xoffset_normal        = xoffset_inch / xsize_inch
;
ix                    = 1
iy                    = 1
x_coord_min_normal    = xoffset_normal + (ix-1)*(width_spacer_normal + width_normal)
x_coord_max_normal    = x_coord_min_normal + width_normal
y_coord_min_normal    = 1.0 - yoffset_normal - iy*height_normal - (iy-1)*height_spacer_normal
y_coord_max_normal    = y_coord_min_normal + height_normal
!p.position           = [x_coord_min_normal,y_coord_min_normal,x_coord_max_normal,y_coord_max_normal]
;
xrange = [0., 3600.]
xticks = 4
xminor = 5
yrange = [7., 17.]
yticks = 2
yminor = 5
;
plot, /nodata, /noerase, [0], xstyle=1, xrange=xrange, $
  xticks=xticks, xminor=xminor, ytitle='T!DP!N, !Uo!NC', ystyle=1, yrange=yrange, yticks=yticks, $
  yminor=yminor, xticklen = -0.05, yticklen = -0.01, xtickname = replicate(' ', xticks + 1)
for i = 0, n_elements(aTkP) - 1, 5 do begin
 plots, asec[i], aTkp[i] - !Tko, psym = 4, symsize = 0.20, color = red
endfor
;
ix                    = 1
iy                    = 2
x_coord_min_normal    = xoffset_normal + (ix-1)*(width_spacer_normal + width_normal)
x_coord_max_normal    = x_coord_min_normal + width_normal
y_coord_min_normal    = 1.0 - yoffset_normal - iy*height_normal - (iy-1)*height_spacer_normal
y_coord_max_normal    = y_coord_min_normal + height_normal
!p.position           = [x_coord_min_normal,y_coord_min_normal,x_coord_max_normal,y_coord_max_normal]
;
xrange = [0., 3600.]
xticks = 4
xminor = 5
yrange = [0., 1500.]
yticks = 3
yminor = 5
;
plot, /nodata, /noerase, [0], xstyle=1, xrange=xrange, $
  xticks=xticks, xminor=xminor, ytitle='Z!DP!N, m', ystyle=1, yrange=yrange, yticks=yticks, $
  yminor=yminor, xticklen = -0.05, yticklen = -0.01, xtickname = replicate(' ', xticks + 1)
for i = 0, n_elements(azP) - 1, 5 do begin
 plots, asec[i], azp[i], psym = 4, symsize = 0.20, color = red
endfor
;
ix                    = 1
iy                    = 3
x_coord_min_normal    = xoffset_normal + (ix-1)*(width_spacer_normal + width_normal)
x_coord_max_normal    = x_coord_min_normal + width_normal
y_coord_min_normal    = 1.0 - yoffset_normal - iy*height_normal - (iy-1)*height_spacer_normal
y_coord_max_normal    = y_coord_min_normal + height_normal
!p.position           = [x_coord_min_normal,y_coord_min_normal,x_coord_max_normal,y_coord_max_normal]
;
xrange = [0., 3600.]
xticks = 4
xminor = 5
yrange = [-3., 3.]
yticks = 2
yminor = 3
;
plot, /nodata, /noerase, [0], xstyle=1, xrange=xrange, $
  xticks=xticks, xminor=xminor, ytitle='U!DP!N, m s!U-1!N', ystyle=1, yrange=yrange, yticks=yticks, $
  yminor=yminor, xticklen = -0.05, yticklen = -0.01, xtickname = replicate(' ', xticks + 1)
for i = 0, n_elements(aUP) - 1, 5 do begin
 plots, asec[i], aUp[i], psym = 4, symsize = 0.20, color = red
endfor
;
ix                    = 1
iy                    = 4
x_coord_min_normal    = xoffset_normal + (ix-1)*(width_spacer_normal + width_normal)
x_coord_max_normal    = x_coord_min_normal + width_normal
y_coord_min_normal    = 1.0 - yoffset_normal - iy*height_normal - (iy-1)*height_spacer_normal
y_coord_max_normal    = y_coord_min_normal + height_normal
!p.position           = [x_coord_min_normal,y_coord_min_normal,x_coord_max_normal,y_coord_max_normal]
;
xrange = [0., 3600.]
xticks = 4
xminor = 5
yrange = [-0.03, 0.03]
yticks = 2
yminor = 3
;
plot, /nodata, /noerase, [0], xstyle=1, xrange=xrange, $
  xticks=xticks, xminor=xminor, ytitle='dU!DP!N/dt, m s!U-2!N', ystyle=1, yrange=yrange, yticks=yticks, $
  yminor=yminor, xticklen = -0.05, yticklen = -0.01, xtitle = 'Time, s'
for i = 0, n_elements(aUP) - 1, 5 do begin
 plots, asec[i], adUpdt[i], psym = 4, symsize = 0.20, color = red
endfor
;
device, /close
stop
;
end