*
* comment:
* the thrid time step in fig8.bin is in fact January 13 at 9UT, and not January 12 at 9UT.
* the TDEF information in fig8.ctl is not accurate in this respected.
*

'reinit'
'set display color white'
'clear'
'set xlopts 1 4 0.15'
'set ylopts 1 4 0.15'
'set clopts 1 4 0.14'
'run ede.col'
'set annot 1 4 '
'set display color white'
'run ede.col'

'open  mowa.hepi.ctl'
'define mask1 = mask'
'close 1'
'set gxout vector'
col = 1

'open fig8.ctl'

'set mproj nps'
'set lat 30 60 '
'set lon 225 285'
'set mpvals 225 285 30 60'

****************************************************
**  first row: 10,11,13 JAN at 9 UT and z = 45 km ** 
****************************************************
'set lev 45  '
vel=120

* 10 JAN, 9UT, 45 km
'set t 1 '
'set annot 1 4'
'set vpage  0.5   6.2     4.30 8.00'
'set parea  0.55  5.65    3.84 7.04'
'set gxout shaded'
'set ccols  9   14     4   11   5    13   3  10   7   12  8   34  36  6'
'set clevs   -2.0 -1.6 -1.2 -0.8 -0.4 -0.2  0 0.2  0.4 0.8 1.2 1.6 2.0 '
'set grads off'
'd w'
xl = 5.74 ; xwid = 0.185 ; yb = 3.92  ; ywid= 0.217 ; sz=0.12
dummy = vcbar (xl,xwid,yb,ywid,sz)
dummy = drawgrid()
dummy = plotarrows(col,vel)
x0 = 5.15
y0 = 6.60
rc = arrow(x0,y0,vel)

* 11 JAN, 9UT, 45 km
'set t 2 '
'set annot 1 4 '
'set vpage  2.7   8.4     4.3  8.0'
'set parea  2.75  7.85    3.84 7.04'
'set ccols  9   14     4   11   5    13   3  10   7   12  8   34  36  6'
'set clevs   -2.0 -1.6 -1.2 -0.8 -0.4 -0.2  0 0.2  0.4 0.8 1.2 1.6 2.0 '
'set grads off'
'd w'
xl = 7.94 ; xwid = 0.18 ; yb = 3.92 ; ywid= 0.217 ; sz=0.12
dummy = vcbar (xl,xwid,yb,ywid,sz)
dummy = drawgrid()
dummy = plotarrows(col,vel)
x0 = 7.35
y0 = 6.60
rc = arrow(x0,y0,vel)

* 13 JAN, 9UT, 45 km; 3*48+19=163
'set t 3 '  
'set vpage  4.9   10.6    4.3  8.0'
'set parea  4.95  10.05   3.84 7.04'
'set ccols  9   14     4   11   5    13   3  10   7   12  8   34  36  6'
'set clevs   -2.0 -1.6 -1.2 -0.8 -0.4 -0.2  0 0.2  0.4 0.8 1.2 1.6 2.0 '
'set grads off'
'd w'
xl = 10.15 ; xwid = 0.18 ; yb = 3.92 ; ywid= 0.217 ; sz=0.12
dummy = vcbar (xl,xwid,yb,ywid,sz)
dummy = drawgrid()
dummy = plotarrows(col,vel)
x0 = 9.55
y0 = 6.60
rc = arrow(x0,y0,vel)


*****************************************************
**  second row: 10,11,13 JAN at 9 UT and z = 85 km ** 
*****************************************************
'set lev 85 '
vel= 70

* 10 JAN, 9UT, 85 km
'set t 1 '  
'set annot 1 4'
'set vpage  0.5   6.2     2.80  6.50'
'set parea  0.55  5.65    2.90  6.10'
'set ccols  9   14     4   11   5    13   3  10   7   12  8   34  36  6'
'set clevs   -2.5 -2.0 -1.5 -1.0 -0.5 -0.25 0 0.25 0.5 1.0 1.5 2.0 2.5 '
'set grads off'
'd w'
xl = 5.74 ; xwid = 0.185 ; yb = 2.98 ; ywid= 0.217 ; sz=0.12
dummy = vcbar (xl,xwid,yb,ywid,sz)
dummy = drawgrid()
dummy = plotarrows(col,vel)
x0 = 5.15
y0 = 5.65
rc = arrow(x0,y0,vel)


* 11 JAN, 9UT, 85 km
'set t 2 ' 
'set annot 1 4 '
'set vpage  2.7   8.4     2.8  6.5'
'set parea  2.75  7.85    2.9  6.1'
'set gxout shaded'
'set ccols  9   14     4   11   5    13   3  10   7   12  8   34  36  6'
'set clevs   -2.5 -2.0 -1.5 -1.0 -0.5 -0.25 0 0.25 0.5 1.0 1.5 2.0 2.5 '
'set grads off'
'd w'
xl = 7.94 ; xwid = 0.18 ; yb = 2.98  ; ywid= 0.217 ; sz=0.12
dummy = vcbar (xl,xwid,yb,ywid,sz)
dummy = drawgrid()
dummy = plotarrows(col,vel)
x0 = 7.35
y0 = 5.65
rc = arrow(x0,y0,vel)

* 13 JAN, 9UT, 85 km
'set t 3 ' 
'set vpage  4.9   10.6    2.8  6.5'
'set parea  4.95  10.05   2.9  6.1'
'set ccols  9   14     4   11   5    13   3  10   7   12  8   34  36  6'
'set clevs   -2.5 -2.0 -1.5 -1.0 -0.5 -0.25 0 0.25 0.5 1.0 1.5 2.0 2.5 '
'set grads off'
'd w'
xl = 10.15 ; xwid = 0.18 ; yb = 2.98   ; ywid= 0.217 ; sz=0.12
dummy = vcbar (xl,xwid,yb,ywid,sz)
dummy = drawgrid()
dummy = plotarrows(col,vel)
x0 = 9.55
y0 = 5.65
rc = arrow(x0,y0,vel)



****************************************************
**  third row: 10,11,13 JAN at 9 UT and z = 85 km ** 
****************************************************
'set lev 110'
vel=170

* 10 JAN, 9UT, 110 km
'set t 1'  
'set vpage  0.5   6.2     1.50  5.10'
'set parea  0.55  5.65    1.60  4.80'
'set ccols  9   14     4   11   5    13   3  10   7   12  8   34  36  6'
'set clevs   -5   -4   -3   -2   -1   -0.5  0 0.5  1   2   3   4   5   '
'set grads off'
'd w'
xl = 5.74 ; xwid = 0.185 ; yb = 1.68  ; ywid= 0.217 ; sz=0.12
dummy = vcbar (xl,xwid,yb,ywid,sz)
dummy = drawgrid()
dummy = plotarrows(col,vel)
x0 = 5.15
y0 = 4.40
rc = arrow(x0,y0,vel)

* 11 JAN, 19UT, 110 km
'set t 2 ' 
'set vpage  2.7   8.4     1.5  5.1'
'set parea  2.75  7.85    1.6  4.8'
'set ccols  9   14     4   11   5    13   3  10   7   12  8   34  36  6'
'set clevs   -5   -4   -3   -2   -1   -0.5  0 0.5  1   2   3   4   5   '
'set grads off'
'd w'
xl = 7.94 ; xwid = 0.18 ; yb = 1.68 ; ywid= 0.217 ; sz=0.12
dummy = vcbar (xl,xwid,yb,ywid,sz)
dummy = drawgrid()
dummy = plotarrows(col,vel)
x0 = 7.35
y0 = 4.40
rc = arrow(x0,y0,vel)

* 11 JAN, 23UT, 110 km
'set t 3 '
'set vpage  4.9   10.6    1.5  5.1'
'set parea  4.95  10.05   1.6  4.8'
'set ccols  9   14     4   11   5    13   3  10   7   12  8   34  36  6'
'set clevs   -5   -4   -3   -2   -1   -0.5  0 0.5  1   2   3   4   5   '
'set grads off'
'd w'
xl = 10.15 ; xwid = 0.18 ; yb = 1.68   ; ywid= 0.217 ; sz=0.12
dummy = vcbar (xl,xwid,yb,ywid,sz)
dummy = drawgrid()
dummy = plotarrows(col,vel)
x0 = 9.55
y0 = 4.40 
rc = arrow(x0,y0,vel)


*****************************************************
**  fourth row: 10,11,13 JAN at 9 UT and z = 85 km ** 
*****************************************************
'set lev 260'
vel = 160

* 10 JAN, 9UT, 260 km
'set t 1'  
'set vpage  0.5   6.2     0.2  3.9'
'set parea  0.55  5.65    0.3  3.5'
'set gxout shaded'
'set ccols  9   14     4   11   5    13   3  10   7   12  8   34  36  6'
'set clevs   -15  -12  -9   -6   -3   -1.5  0 1.5  3   6   9  12  15   '
'set grads off'
'd w'
xl = 5.74 ; xwid = 0.185 ; yb = 0.385 ; ywid= 0.217 ; sz=0.12
dummy = vcbar (xl,xwid,yb,ywid,sz)
dummy = drawgrid()
dummy = plotarrows(col,vel)
x0 = 5.15
y0 = 3.07
rc = arrow(x0,y0,vel)

* 11 JAN, 9UT, 260 km
'set t 2 '
'set vpage  2.7   8.4     0.2  3.9'
'set parea  2.75  7.85    0.3  3.5'
'set ccols  9   14     4   11   5    13   3  10   7   12  8   34  36  6'
'set clevs   -15  -12  -9   -6   -3   -1.5  0 1.5  3   6   9  12  15   '
'set grads off'
'd w'
xl = 7.94 ; xwid = 0.18 ; yb = 0.385; ywid= 0.217 ; sz=0.12
dummy = vcbar (xl,xwid,yb,ywid,sz)
dummy = drawgrid()
dummy = plotarrows(col,vel)
x0 = 7.35
y0 = 3.07
rc = arrow(x0,y0,vel)


* 13 JAN, 9UT, 260 km
'set t 3 '
'set vpage  4.9   10.6    0.2 3.9 '
'set parea  4.95  10.05   0.3 3.5 '
'set ccols  9   14     4   11   5    13   3  10   7   12  8   34  36  6'
'set clevs   -15  -12  -9   -6   -3   -1.5  0 1.5  3   6   9  12  15   '
'set grads off'
'd w'
xl = 10.15 ; xwid = 0.18 ; yb = 0.385  ; ywid= 0.217 ; sz=0.12
dummy = vcbar (xl,xwid,yb,ywid,sz)
dummy = drawgrid()
dummy = plotarrows(col,vel)
x0 = 9.55
y0 = 3.07
rc = arrow(x0,y0,vel)

'set parea off'
'set vpage off'

'set string 1 c 4 0'
'set strsiz 0.115'


"draw string 5.46  8.32  w & (U,V) (m/s), 09:00`a `nUT"

'set strsiz 0.100'
"draw string 2.10  8.05  (a) z=45`a `nkm,`a `n10JAN2017"
"draw string 5.45  8.05  (b) z=45`a `nkm,`a `n11JAN2017"
"draw string 8.80  8.05  (c) z=45`a `nkm,`a `n13JAN2017"

"draw string 2.10  6.06  (d) z=85`a `nkm,`a `n10JAN2017"
"draw string 5.45  6.06  (e) z=85`a `nkm,`a `n11JAN2017"
"draw string 8.80  6.06  (f) z=85`a `nkm,`a `n13JAN2017"

"draw string 2.10  4.09  (g) z=110`a `nkm,`a `n10JAN2017"
"draw string 5.45  4.09  (h) z=110`a `nkm,`a `n11JAN2017"
"draw string 8.80  4.09  (i) z=110`a `nkm,`a `n13JAN2017"

"draw string 2.10  2.12  (j) z=260`a `nkm,`a `n10JAN2017"
"draw string 5.45  2.12  (k) z=260`a `nkm,`a `n11JAN2017"
"draw string 8.80  2.12  (l) z=260`a `nkm,`a `n13JAN2017"

'gxprint  fig8.eps'

exit

********************************************************************
function vcbar (xl,xwid,yb,ywid,sz)  (E. Becker 14.10.2004)
*  Check shading information
'query shades'
shdinfo = result
*say shdinfo
if (subwrd(shdinfo,1)='None')
  say 'Cannot plot color bar: No shading information'
  return
endif
cnum = subwrd(shdinfo,5)
*  Plot colorbar
'set string 1 c 4'
*'set strsiz 0.120 0.120'
'set strsiz 'sz ' 'sz 
num   = 0
x1 = xl
y1 = yb
x2 = xl+xwid
while (num<cnum)
  rec = sublin(shdinfo,num+2)
  col = subwrd(rec,1)
  hi  = subwrd(rec,3)
  'set line 'col
  y2 = y1 + ywid
  'draw recf 'x1' 'y1' 'x2' 'y2
  if (num<cnum-1 )
*    if( num=0 | num=2 | num=4 | num=6 | num=8 | num=10 | num=12 | num=14 )
*      'draw string '%(xl+xwid+0.3)%' 'y2' 'hi
*    endif
    'draw string '%(xl+xwid+0.25)%' 'y2' 'hi
  endif
  num = num + 1
  y1 = y2
endwhile
'set line 1'
x1=xl; y1=yb
'draw line 'x1' 'y1' 'x1' 'y2
'draw line 'x1' 'y2' 'x2' 'y2
'draw line 'x2' 'y1' 'x2' 'y2
'draw line 'x1' 'y1' 'x2' 'y1
return
*******************************************************************


*******************************************************************
function plotarrows(col,vel)
'set annot 1 0'
'set gxout vector'
'set arrlab off '
'set arrowhead 0.11'
'set grid off'
len = 0.5
'set arrscl 'len' 'vel
'set ccolor 'col
'set cthick 1'
'd maskout(ul,mask1);maskout(vl,mask1)'
'set annot 1 4'
return
*******************************************************************
              
*******************************************************************
function arrow(x,y,vel)
'set line 1 1 5'
len = 0.5
'draw line 'x-len/2.' 'y' 'x+len/2.' 'y
'draw line 'x+len/2.-0.1' 'y+0.05' 'x+len/2.' 'y
'draw line 'x+len/2.-0.1' 'y-0.05' 'x+len/2.' 'y
'set string 1 c'
'set strsiz 0.13'
'draw string 'x' 'y-0.22' 'vel
return
*******************************************************************
function drawgrid()
'set annot 1 0 '
'set gxout contour'
'set cthick 1'
'set clab off'
'set ccolor 1'
'set cstyle 2'
'set clevs 30 40 50 60'
'set grads off'
'd lat'
'set ccolor 1'
'set cstyle 2'
'set clevs 225 240 255 270 284.9999'
'set grads off'
'd lon'
'set annot 1 4 '
'set gxout shaded'
return
*******************************************************************