pro plot_merged_sounding_main

common  merged_sounding_data, num_times2, hrfrac2, jday2, temp, $
	height2, num_heights2, press, rh, theta, mxrat, svp, smxrat

common merged_sounding_data2, sflag

common mmcr_merged_data_dbz, num_times33, hrfrac33, jday33, dbz33, $
	height33, num_heights33


common mmcr_merged_data_laser, laser_base33

common mmcr_merged_data, num_times, hrfrac, jday, dbz, $
	height, num_heights, vel, width

common mmcr_merged_data2, snr

common mmcr_merged_data3, laser_base


; construct the filename
gname_prefix='Dsgpcombined.C1.'
gname_suffix='.000000.merged_sounding.gif'

fname=''
openr, 1, '/data/mace4/arm_data/drivers/merged_sounding_files.dat'
readf, 1, fname
close, 1
sepname=str_sep(fname,'.')
i=n_elements(sepname)-4
cur_day=sepname[i]
six_day=strmid(cur_day,2,6)

date_string_fm_adtg, cur_day, date_string

loadct, 5 ;  This is the stern color scheme
gamma_ct, 1.0

!p.charsize=2.00
!p.font=-1

; Get the data

get_merged_sounding, fname, 0.0, 23.99, 0., 15000.


if num_times gt 1 then begin

WINDOW, 0, /pixmap, XSIZE=800, YSIZE=800, TITLE='Merged Sounding '
!p.multi = [0,3,3]


; Temperature

height[*]=height[*]/1000.
height2[*]=height2[*]/1000.

xl=0.075
yl=0.40
xr=0.75
yr=0.65


plot_2d_image_999_maxminin_test, float(temp), hrfrac, $
	height2, num_times, num_heights, $
	' ',' ','Heght (km)',[xl,yl,xr,yr], [xr+0.080,yl,xr+0.095,yr], $
	'Temperature (k)',330.,190., 200.

a_string='Temperature'
xyouts, (800*(xr+0.095)+800)/2., 800*((yl+yr)/2.), $
alignment=0.5,a_string, /device, color=0, charsize=1.25
a_string=''
xyouts, (800*(xr+0.095)+800)/2., 800*((yl+yr)/2.)-20, $
alignment=0.5,a_string, /device, color=0, charsize=1.25

; RH

xl=0.075
yl=0.70
xr=0.75
yr=0.95

plot_2d_image_999_maxminin_test, float(rh), hrfrac, $
	height2, num_times, num_heights, $
	' ',' ','Heght (km)',[xl,yl,xr,yr], [xr+0.080,yl,xr+0.095,yr], $
	'RH (%) ',100.,0., 200.


a_string='Relative'
xyouts, (800*(xr+0.095)+800)/2., 800*((yl+yr)/2.), $
alignment=0.5,a_string, /device, color=0, charsize=1.25
a_string='Humidity'
xyouts, (800*(xr+0.095)+800)/2., 800*((yl+yr)/2.)-20, $
alignment=0.5,a_string, /device, color=0, charsize=1.25

; Theta

xl=0.075
yl=0.10
xr=0.75
yr=0.35

plot_2d_image_999_maxminin_test, float(theta), hrfrac, $
	height2, num_times, num_heights, $
	' ','Time (UTC hours) ','Heght (km)',[xl,yl,xr,yr], [xr+0.080,yl,xr+0.095,yr], $
	'Theta (k)',390.,250., 200.


a_string='Potential'
xyouts, (800*(xr+0.095)+800)/2., 800*((yl+yr)/2.), $
alignment=0.5,a_string, /device, color=0, charsize=1.25
a_string='Temperature'
xyouts, (800*(xr+0.095)+800)/2., 800*((yl+yr)/2.)-20, $
alignment=0.5,a_string, /device, color=0, charsize=1.25


a_string='SGP Merged Sounding, '+date_string

xyouts, 400, 5, alignment=0.5,a_string, /device, color=0, charsize=2

path='/data/mace4/arm_data/temp_image/'
gname=path+gname_prefix+cur_day+gname_suffix
;if timeinc eq 0 then gname=path+gname_prefix+cur_day+'.Moments.00-12UTC'+fname_suffix
;if timeinc eq 1 then gname=path+gname_prefix+cur_day+'.Moments.12-00UTC'+fname_suffix


write_gif, gname, tvrd()

window, 0, /pixmap,  xsize=500, ysize=600, title='Merged Sounding Status Flag'
!p.multi=[0,0,0]

plot, hrfrac, sflag, background=!d.n_colors-1, color=0, charsize=1.5, $
xtitle="Time (UTC)", xrange=[0,24], xstyle=1, xticks=12, psym=6, $
yrange=[0,18], ystyle=1, yticks=18,  pos=[.075,.4,.925,.88]

xyouts,250, 575, alignment=0.5, date_string, /device, color=0, charsize=2

a_string='Merged Sounding Status Flag'
xyouts,250, 550, alignment=0.5, a_string, /device, color=0, charsize=2

a_string='1: CART sonde, mwr vapor, aeri profile'
xyouts, 25, 160, alignment=0.0, a_string, /device, color=0, charsize=1.25

a_string='2: CART sonde, mwr vapor'
xyouts, 25, 145, alignment=0.0, a_string, /device, color=0, charsize=1.25

a_string='3: CART sonde, aeri profile'
xyouts, 25, 130, alignment=0.0, a_string, /device, color=0, charsize=1.25

a_string='4: CART sonde'
xyouts, 25, 115, alignment=0.0, a_string, /device, color=0, charsize=1.25

a_string='5: OUN Sonde, mwr vapor, aeri profile'
xyouts, 25, 100, alignment=0.0, a_string, /device, color=0, charsize=1.25

a_string='6: OUN Sonde, mwr vapor'
xyouts, 25,85, alignment=0.0, a_string, /device, color=0, charsize=1.25

a_string='7: OUN Sonde, aeri profile'
xyouts, 25, 70, alignment=0.0, a_string, /device, color=0, charsize=1.25

a_string='8: MOLTS, mwr vapor, aeri profile'
xyouts, 25, 55, alignment=0.0, a_string, /device, color=0, charsize=1.25

a_string='9: MOLTS, mwr vapor'
xyouts, 275, 160, alignment=0.0, a_string, /device, color=0, charsize=1.25

a_string='10: MOLTS, AERI Profile'
xyouts, 275, 145, alignment=0.0, a_string, /device, color=0, charsize=1.25

a_string='11: ECMWF, mwr vapor, aeri profile'
xyouts, 275, 130, alignment=0.0, a_string, /device, color=0, charsize=1.25

a_string='12: ECMWF, mwr vapor'
xyouts, 275, 115, alignment=0.0, a_string, /device, color=0, charsize=1.25

a_string='13: ECMWF, aeri profile'
xyouts, 275, 100, alignment=0.0, a_string, /device, color=0, charsize=1.25

a_string='14: RUC, mwr vapor, aeri profile'
xyouts, 275, 85, alignment=0.0, a_string, /device, color=0, charsize=1.25

a_string='15: RUC, mwr vapor'
xyouts, 275, 70, alignment=0.0, a_string, /device, color=0, charsize=1.25

a_string='17: No Sounding'
xyouts, 275, 55, alignment=0.0, a_string, /device, color=0, charsize=1.25

path='/data/mace4/arm_data/temp_image/'
gname_prefix='Dsgpstatusflag.C1.'
gname=path+gname_prefix+cur_day+gname_suffix

write_gif, gname, tvrd()

endif	; check for num_times

;endfor	; time increment

end_program: print, 'Done'
end