c -----------------------------------------------------
c 
c TAO Project Office,  NOAA/PMEL/OCRD
c
c Programmed by Nancy Soreide, 1991. Modified by Nancy, 
c   Dai McClurg, and Willa Zhu to initiate a subsequent
c   plot using stack_met.ppc.
c 
c sensor.ppc
c
C COMMAND FILE FOR WIND,TEMP FROM REAL-TIME DATA
c  pressures removed summer 1994
c
c parameters
c
c  p1 latitude, eg, 2n (up1 is 2N)
c  p2 longitude, eg, 110w (up2 is 110W)
c  p3 is path, eg, /home/heron/taohome
c  p4 is Process id, eg, PID
c  p5 ATLAS, CM, or none
c  p6 met file name or cm file name
c  p7 data subdirectory
c  p8 starttime, e.g. 910101
c  p9 endtime, e.g. 940831
c
c -----------------------------------------------------*
c
c set prflag to 1 to also plot presures, else 0.
c
c The exported version of this script MUST have prflag = 0.
c
set prflag 0
c
tics,0.05,0.15,0.05,0.15
c
set *starttime "'p8'0000"
set *endtime "'p9'0000"
c
window,on
echo off
set up1 $edit(p1,UPCASE)
set up2 $edit(p2,UPCASE)
@'p3'/ppc/color.ppc SUN
c
pen,1,'red'
BOX,OFF
WINDOW,ON
CROSS,2
size 8 10.7
c xgeom 550 750
window,on
dfltfnt,CR
c
set xx 6.0
set yy 1.1
set ytemps 3.5
set yyp 0.8
set ysp 0
c
set ypos 8.7
set xpos 1.2
origin 'xpos' 'ypos'
axlen,'xx','yy'
c 
CONSET,,,,,,0.,0
set xt 'xx' / 2.
set yt 'yy' + .5
set xl 'xx' - .1
c
c  Read climatology for this location into variables for u,v,air,sst
c
e openr "'p3'/data/buoy/c'p1''p2'.cdf"
e nextr
e slab=[i=1,j=1,k=1,t='starttime':'endtime']
e cv1=eps423[slab]
e cu1=eps422[slab]
e cr1=eps910[slab]
e ca1=eps21[slab]
e cs1=eps25[slab]
e close epsin
c *********************** Open wind netCDF file, get u, v, and rh; close it.
e openr "'p3'/'p7'/w'p1''p2'_dy.cdf"
e nextr
e yt=[y='up1',t='starttime':'endtime']
e status u
c 
if eps$fld_tmax .lt. starttime then
else
  e u1=u[yt]
  e v1=v[yt]
  e r1=eps910[yt]
  e u1
  e status u1
  line,'ppl$line_count',,0
  pen,'ppl$line_count','red'
endif
e close
yaxis,-10,10,5
e cu1
line,'ppl$line_count',,0
pen,'ppl$line_count','magenta'
c
taxis,on             
time,w'starttime',w'endtime',w'eps$fld_tmin'
set yl 'yy' - .2
labs:nouser,1,'xl','yl',+1,"@p'red'U"
pltnme,sensor'up1''up2'_'p4'.plt
c  eg, sensor2202.plt001 where 2202 is PID
pltype,0
%opnplt
C *********************** Plot U
c labs:nouser,2,'xt','yt',0,"'up1' 'up2'"
set ytop 'yt' + .3
hlabs,2,.18
labs:nouser,2,'xt','yt',0,"@trTAO Daily Averaged Data"
set ytop 'yt' + .6
c labs:nouser,3,'xx','ytop',1,"'date'"
YAXIS,   -10.00,    10.00,    2.50
axlint,1,2
YFOR,(i5)
YLAB,W-(m/s)-E
DATPT,2
TXTYPE,YR,MON3
txlabp,0,1
title
plot:nowait
C
C *********************** Now plot V
c
set ypos 'ypos' - 'yy'
origin 'xpos' 'ypos'
labs
labs:nouser,1,'xl','yl',+1,"@p'blue'V"
TAXIS,ON
TXTYPE,YR,MON3
txlabp,0
if eps$fld_tmax .lt. starttime then
else
  e v1
  line,'ppl$line_count',,0
  pen,'ppl$line_count','blue'
endif
e cv1
c line,'ppl$line_count',20,1
line,'ppl$line_count',,0
pen,'ppl$line_count','cyan'
YLAB,S-(m/s)-N
DATPT,2
title
PLOT:nowait
C
C *********************** Now plot RH
c
set ypos 'ypos' - 'yy'
origin 'xpos' 'ypos'
c
yaxis,55,115,15
c
labs
labs:nouser,1,'xl','yl',+1,"@p'black'RH"
TAXIS,ON
TXTYPE,YR,MON3
txlabp,0
if eps$fld_tmax .lt. starttime then
else
  e r1
  line,'ppl$line_count',,0
  pen,'ppl$line_count','black'
endif
e cr1
c line,'ppl$line_count',20,1
line,'ppl$line_count',,0
pen,'ppl$line_count','green'
YLAB,"@AS%"
DATPT,2
title
PLOT:nowait
c
cccccccccccccccccccccccccc Temperatures cccccccccccccccccccccccccccccccc
c
nlines
set ypos 'ypos' - 'ytemps'
origin 'xpos' 'ypos'
axlen,'xx','ytemps'
c *****************  Open netCDF file
e openr "'p3'/'p7'/t'p1''p2'_dy.cdf"
e nextr
C *****************  Now plot temperatures from each sensor
c
limits,0,YLE,ON
labs
YAXIS,   7.00,    31.00,    2.00
YLAB,TEMPERATURE
AXLINT,1,1
DATPT,2
c *********************** Read air temperature, load into pplus buffers
e yt2=[y=1,t='starttime':'endtime']
e status eps21
taxis,on             
time,w'starttime',w'endtime',w'eps$fld_tmin'
if eps$fld_tmax .lt. starttime then
else
  e at1=eps21[yt2]
  e at1
  e status at1
  pen,'ppl$line_count' 'black'
  line,'ppl$line_count',,0
endif
c yaxis,20,35,5
c ********************** Read water temperatures, load into pplus buffers
set count 1
set pcount 2
set numt 11
c
e status eps20
c set ppl$xlast2 1.e34
c set ppl$ylast2 1.e34
if eps$fld_tmax .lt. starttime then
  enter
  1,1.e35
  end
else
   while count .le. 'numt' then
        e t1=eps20[yt2,k='count']
	e status t1
	if eps$fld_tmax .eq. eps$fld_tmin then
	else
	  e t1
        endif
    set cnt 'count' + 1
    line,'cnt',,0
    pen,'cnt','pcount'
    inc pcount
    set pcount $integer(pcount)
    if pcount .eq. 9 then
         set pcount 1
    endif
    inc count
    set count $integer(count)
    endw
endif
c if ppl$ylast(2) .gt. 1.e30 then
c     set ppl$ylast(2) 1.e34
c     set ppl$xlast(2) 1.e34
c endif
c if ppl$ylast(1) .gt. 1.e30 then
c     set ppl$ylast(1) 1.e34
c     set ppl$xlast(1) 1.e34
c endif
c
e xloc = 'ppl$xlen' - 0.4
e load(xloc)
set xloc 'eps$scaler'
c
e yloc1 = 'ppl$ylen' - 0.25
e load(yloc1)
set yloc1 'eps$scaler'
c
e yloc2 = 'ppl$ylen' - 0.5
e load(yloc2)
set yloc2 'eps$scaler'
c
yaxis,5,34,5
labs:nouser,1,'xloc','yloc1',-1,"@p2SST"
labs:nouser,2,'xloc','yloc2',-1,"@p1AIR"
c
e cs1
line,'ppl$line_count',,0
pen,'ppl$line_count','magenta'
c
c AT climatology
c
c e ca1
c line,'ppl$line_count',,0
c pen,'ppl$line_count','cyan'
c
title
YLAB,TEMPERATURE (C)
if prflag .eq. 1 then
  txlabp,0
  plot:nowait
else
  txlabp,-1
  set realpos "   ('eps$latitude'  'eps$longitude')"
  plot:nowait "'up1' 'up2' 'realpos'"
endif
c
e close epsin
c
yaxis
labs
set ypos 'ypos' - 'yyp'
origin 'xpos' 'ypos'
axlen,'xx','yyp'
c
if prflag .eq. 1 then
  e openr "/home/tao2/taodata/alltao/t'p1''p2'.cdf"
  e nextr
C *****************PLOT P300
  pen,1,'magenta'
  e p1=eps1[yt2,k=1]
  e status p1
  if eps$fld_tmax .lt. starttime then
    enter
    1,1.e35
    end
  else
    e p1
  endif
  YAXIS 320,270,10
  YLAB,P (db)
  AXLINT,2,2
  title
  PLOT:nowait
C *****************PLOT P500
  txlabp,-1
  set ypos 'ypos' - 'yyp'
  origin 'xpos' 'ypos'
  limits,520,YGE,on
  pen,1,'aqua'
  e p1=eps1[yt2,k=2]
  e status p1
  if eps$fld_tmax .lt. starttime then
    enter
    1,1.e35
    end
  else
    e p1
    set dplot "YES"
  endif
  e close epsin
  TAXIS,ON
  yaxis,520,470,10
  YLAB,P (db)
  txlabp,-1
  title
  set realpos "   ('eps$latitude'  'eps$longitude')"
  plot:nowait "'up1' 'up2' 'realpos'"
endif
c
limits,yge,off
limits,yle,off
c
yaxis
nlines
c
@'p3'/ppc/util/proj.ppc
c
@'p3'/ppc/util/date_stamp.ppc lr
c
%clsplt
c
c ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c  high resolution data from the met files
c ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c
c parameters
c
c  p1 latitude, eg, 2n (up1 is 2N)
c  p2 longitude, eg, 110w (up2 is 110W)
c  p3 is path, eg, /home/heron/taohome
c  p4 is Process id, eg, PID
c  p5 ATLAS, CM, or none
c  p6 met or CM file name
c
if p6 .ne. "none" then
  if p5 .eq. "ATLAS" then
     @'p3'/ppc/atlas_oper/stack_met.ppc 'p3' 'p4' 'p6' 'p1' 'p2'
  endif
  if p5 .eq. "CM" then
     @'p3'/ppc/cm_oper/stack_cm.ppc 'p3' 'p4' 'p6' 'p1' 'p2'
  endif
endif
c