PRO PLOT_SONDE_RH_WINDS_5KM

; plots the name relative humidity with winds overlain
; at 1 mb pressure resolution
printflag=1
CASE (printflag) OF
0:	begin
	set_plot,'win'
	device,decompose=0
	window,/free,xsize=900,ysize=300,title='pacs03 rh'
	white=!d.n_colors-1
	!p.background=white & !p.color=0
	!p.font=0
	device,set_font='Helvetica*bold*proof*16'
	end
1:  begin
    set_plot,'ps'
    !p.font=0
    device,/times,font_size=12,font_index=7
    device,/inches,ysize=8.0,yoffset=2.0
    device,filename='/Users/pz/rico/rhuv_5km.ps'
    device,/color,bits=8
    device,/landscape
    end
ENDCASE
loadct,39
;!p.multi=[0,1,2,0,0]
!p.multi=[0,0,0,0,0]
RESTORE,'/Users/pz/rico/'+'rico_sonde_1p_leg2.idlsave'
	;pgrid,zgrid,tnew,rhnew,filenames,unew,vnew

z=zgrid/1000.
; tack in a vector of null values for 17/19
dates=[dates[0:6],'05011720',dates[7:nsonde-1]]
empty=fltarr(1,1)+!values.f_nan
empty2=fltarr(1,1010)+!values.f_nan
tnew=[tnew,empty2] & rhnew=[rhnew,empty2]
unew=[unew,empty2] & vnew=[vnew,empty2]
FOR i=0,1009 DO begin
	tnew[0:nsonde,i]=[tnew[0:6,i],empty,tnew[7:nsonde-1,i]]
	rhnew[0:nsonde,i]=[rhnew[0:6,i],empty,rhnew[7:nsonde-1,i]]
	unew[0:nsonde,i]=[unew[0:6,i],empty,unew[7:nsonde-1,i]]
	vnew[0:nsonde,i]=[vnew[0:6,i],empty,vnew[7:nsonde-1,i]]
ENDFOR
nsonde=nsonde+1
if printflag eq 0 then !p.position=[0.1,0.1,0.9,0.9]
;if printflag eq 1 then !p.position=[0.1,0.5,0.9,0.9]

mylevels=findgen(21)*5.
mycolors=[findgen(21)*12.]
mylevels=findgen(41)*2.5
mycolors=[findgen(41)*6.]
day=float(strmid(dates[1:nsonde-1],4,2))
hr=float(strmid(dates[1:nsonde-1],6,2))
hr=hr+1
b=where(hr eq 24)
hr[b]=0 & day[b]=day[b]+1
b=where(hr eq 1,count) & if count gt 0 then hr[b]=0
b=where(hr eq 13,count)  & if count gt 0 then hr[b]=12
b=where(hr eq 9,count)  & if count gt 0 then hr[b]=8
daynames=string(day,format='(i2)')  ;+':'+string(hr,format='(i2)')
dontuse=where(hr ne 0.) ;leg1
;dontuse=where(hr ne 12.) ;leg2
daynames[dontuse]=' '
;daynames=daynames[[indgen(nsonde/2)*2]]
CONTOUR,rhnew[1:nsonde-1,10:456],indgen(nsonde-1),pgrid[10:456], $
	xrange=[0,nsonde-1],/xstyle, $
	yrange=[1000.0,555.0],/ystyle,ytit='mb', $
	xtit='January 2005 (day label at 0 UTC)', $
	tit='RICO Leg 2 relative humidity & u,v winds', levels=mylevels, $
	c_colors=mycolors, xticks=n_elements(daynames)-1,xtickname=daynames, $
	max_value=110.0,min_value=0.0,/cell_fill, $
	position=[0.05,0.1,0.9,0.9],/nodata,yticks=3

data=rhnew[1:nsonde-1,10:456]
minval=0.0 & maxval=100.0
good=where(data ge minval and data le maxval)
data[good]=data[good]/(maxval-minval)*254.
bad=where(data eq -999.0,count)
if count gt 0 then data[bad]=255

FIT_IN_WINDOW,1,data,top=255,/noscale
axis,/yaxis,yrange=[1000.0,555.0],/ystyle,ytit='km',yticks=5, $
	ytickname=['0.1','1','2','3','4','5']
axis,xaxis=1,xticklen=0.0,xtickname=[' ',' ',' ',' ',' ',' ']

label=string(mylevels,format='(i3)')
for k=1,40,4 do label[k:k+2]=' '
if printflag eq 0 then $
MAKE_KEY, 0.93,0.1,0.03,0.8,/normal,orientation=1,colors=mycolors, $
	labels=label,charsize=0.8
if printflag eq 1 then $
MAKE_KEY, 0.95,0.55,0.03,0.35,/normal,orientation=1,colors=mycolors, $
	labels=label,charsize=0.8

unew2=unew[1:nsonde-1,10:455] & vnew2=vnew[1:nsonde-1,10:455]
b=where(unew2 eq -999.0,count) & if count gt 0 then unew2[b]=!values.f_nan
b=where(vnew2 eq -999.0,count) & if count gt 0 then vnew2[b]=!values.f_nan
unew2=congrid(unew2,nsonde-1,44)
vnew2=congrid(vnew2,nsonde-1,44)
p2=congrid(pgrid[10:455],44)
loadct,42,file='/Users/pz/myidl_lib/colors1.tbl'
OVECTOR,unew2,vnew2,indgen(nsonde-1),p2,length=4,thick=2, $
	hthick=2,color=3,hsize=2,/center,standard=10, $
	/legend,leg_loc=[nsonde-8,695]
XYOUTS,nsonde-8,545.0,/data,'10 m s!U-1!N',charsize=1.0
XYOUTS,nsonde+8,750.0,/data,'% RH'
if printflag eq 0 then begin
fpng='C:\paquita\pacs03\pacs03_sondes\sonde_rh_uv.png'
write_png,fpng,tvrd(true=1)
endif
if printflag eq 1 then begin
	timestamp,'/Users/pz/rico/plot_sonde_rh_winds_5km.pro'
	device,/close_file
endif
print,'dun'

END