c
c -----------------------------------------------------
c
c TAO Project Office,  NOAA/PMEL/OCRD
c
c  Programmed by Dai McClurg, Spring 1994.
c
c  Plot ocean current vectors
c
c vec.ppc
c
c  p1 latitude, eg, 2n (up1 is 2N)
c  p2 longitude, eg, 110w (up2 is 110W)
c  p3 path, eg, /home/heron/taohome
c  p4 frequency identifier: hr (hourly); dy (daily); hf (5-day); lf (monthly)
c  p5 starttime (e.g. 920415)
c  p6 endtime   (e.g. 930415)
c  p7 min depth
c  p8 max depth
c  p9 vector length (inches)
c  p10 vector magnitude
c  p11 data subdirectory
c
c -----------------------------------------------------
c
window,on
c
set up1 $edit(p1,UPCASE)
set up2 $edit(p2,UPCASE)
c
c Following command sets tx1 and tx2
c
set *starttime "'p5'0000"
set *endtime   "'p6'0000"
c
cross,0
c
e xl2 = 'xx' - 0.1
e load(xl2)
set xl2 'eps$scaler'
c
set lsz 0.12
c
c tics,0.08,0.13,0.08,0.13,-1,-1
c
origin 'xpos' 'ypos'
taxis,on
axlint,1,2
yfor,(i4)
c
c cccccccccccccccccccccccccc Current Vectors ccccccccccccccccccccccccccccccc
c
e openr "'p3'/'p11'/cur'p1''p2'_'p4'.cdf"
e nextr
c
@'p3'/ppc/util/buoy_tim_set.ppc:quiet 'p5' 'p6' eps320
c
if error .ne. "none" then
  taxis,off
  e close epsin
  set xtemp 0
  set ytemp 'yy' / 2
  %label:nouser,'xtemp','ytemp',-1,0,0.13,"@p'red''error'"
  set ytemp 'ytemp' - 0.3
  %label:nouser,'xtemp','ytemp',-1,0,0.13,"@p'red''eps$fld_name' 'up1' 'up2'"
c
  datpt,0
c
  txlabp,-1
  axset,1,1,1,1
  line,1,,0
  cross,0
c
  e ypos = 'ypos' + 'yup'
  e load(ypos)
  set *ypos 'eps$scaler'
c
  return
endif
c
e xlp = 'ppl$xlen' / 2
e load(xlp)
set xlp 'eps$scaler'
c
e ylp = 'ppl$ylen' + 0.05
e load(ylp)
set ylp 'eps$scaler'
c
%label:nouser,'xlp','ylp',0,0,0.13,"Ocean Current (cm s^-^1) 'up1' 'up2'"
c
vecset,'p9','p10'
c
e ylp = 'ppl$ylen' + 0.1
e load(ylp)
set ylp 'eps$scaler'
c
e xlp = 'ppl$xlen' - 1.15
e load(xlp)
set xlp 'eps$scaler'
c
veckey:nouser,'xlp','ylp',,(f5.0)
c
e r = abs('p8' - 'p7')
e z1 = 'p8' - 0.1 * r
e load(z1)
set z1 'eps$scaler'
c
%range,'p7','z1',8
c yaxis,'ppl$range_high','ppl$range_low','ppl$range_inc'
yaxis,'p7','z1','ppl$range_inc'
datpt,2
c
e sl=[t='starttime':'endtime',z='p7':'p8']
e u1=eps320[sl]
e v1=eps321[sl]
e status u1
c
e fm = abs('eps$fld_max')
e load(fm)
set fm 'eps$scaler'
c
if fm .eq. 1e+35 then
  e xtemp = 'ppl$xlen' / 2
  e load(xtemp)
  set xtemp 'eps$scaler'
  e ytemp = 'ppl$ylen' / 2
  e load(ytemp)
  set ytemp 'eps$scaler'
  %label:nouser,'xtemp','ytemp',0,'lsz',-
    "@p'red'ALL DATA MISSING IN THIS TIME WINDOW"
endif
c
time,w'tx1',w'tx2',w'eps$fld_tmin'
c
e u1
e v1
e vector(u1,v1)
c
set dec 1
c
if p4 .eq. "hf" then
  e epsin.vector_interval
  set dec 'eps$attr'
endif
c
if p4 .eq. "lf" then
  e epsin.vector_interval
  set dec 'eps$attr'
endif
c
e close epsin
c
ylab "Depth (m)"
c
if year .eq. "null" then
  title,
else
  title,"'year'"
endif
c
pen,1,'blue'
c
vector:nowait,'dec'
c
e xtemp = 'ppl$xlen' + 0.2
e load(xtemp)
set xtemp 'eps$scaler'
c
e ytemp = 'ppl$ylen' / 2.0
e load(ytemp)
set ytemp 'eps$scaler'
c
@'p3'/ppc/util/vector_key.ppc 'p3' 'xtemp' 'ytemp' 1 merup
c
datpt,0
c
txlabp,0
set *year "null"
axset,1,1,1,1
line,1,,0
cross,0
c
e ypos = 'ypos' + 'yup'
e load(ypos)
set *ypos 'eps$scaler'
c