c
c Arguments
c
c p1	pltype
c p2	4 digit date (e.g. 9506)
c p3	taopath
c p4	data subdirectory (e.g. data/lat_lon)
c p5	Process ID number
c p6	frequency identifier (e.g. lf or hf)
c p7	variable short name (e.g. array, windarray, etc)
c
c -----------------------------------------------------
c
c TAO Project Office,  NOAA/PMEL/OCRD
c
c Programmed by Dai McClurg July 1995
c
pltype,'p1'
c
dfltfnt,cr
window,on
c
pltnme,lat_lon_'p2'_'p7'_'p6'_'p5'.plt
c
axlint,1,1
c
set yaxmin -20
set yaxmax  20
c
set xaxmin 125
set xaxmax 275
c
yaxis,'yaxmin','yaxmax',10
xaxis,'xaxmin','xaxmax',20
c
set xx 7.1
e yy = 'xx' * ('yaxmax' - 'yaxmin') / ('xaxmax' - 'xaxmin')
e load(yy)
set yy 'eps$scaler'
c
axlen 'xx' 'yy'
c
e sx = 'xx' + 1.5
e load(sx)
set sx 'eps$scaler'
c
e sy = 'yy' + 2
e load(sy)
set sy 'eps$scaler'
c
size,'sx','sy'
c
e xg = 'sx' * 100
e load(xg)
set xg 'eps$scaler'
c
e yg = 'sy' * 100
e load(yg)
set yg 'eps$scaler'
c
xgeom,'xg','yg'
c
yfor (i3,''lat'')
xfor (i3,''lone'')
ylab,
Xlab,
axnmtc,3,2
tics,0.07,0.17,0.07,0.17,-1,-1
labset,,.15,.15
axlsze,.14,.14
c
c Note that x1 and y1 are used in lat_lon_sam
c
e x1 = 'xx' / 2
e load(x1)
set *x1 'eps$scaler'
c
e y1 = 'yy' + 0.3
e load(y1)
set *y1 'eps$scaler'
c
%opnplt
c
origin,0.85,1
c
@'p3'/ppc/util/month.ppc:quiet "'p2'15"
c
set ypos 'y1' + 0.35
c
%label:nouser 'x1','ypos',0,0,0.2,"@trTAO Array"
c
e openr "'p3'/'p4'/'p2'_'p7'_xy_lf.cdf"
e nextr
c
@'p3'/ppc/util/lat_lon_sam.ppc:quiet 'p3'
c
e close epsin
c
@'p3'/ppc/util/eqline.ppc:quiet 'xaxmin' 'xaxmax'
c
loadcmap,'p3'/ppc/gray.1
c
@'p3'/ppc/util/mask.ppc:quiet 'p3'
c
c refresh axis lines
c
enter
  1.0e35,1.0e35
end
c
plot
c
@'p3'/ppc/util/proj.ppc:quiet
c
%clsplt:wait