; ------------------------------------------------ ;
; ----- Cloud_Layer_Conditional_multi_base_RL ----- ;
; ------------------------------------------------ ;

;;;if (flag_Cloud_Layer_Conditional_multi_base_RL eq 1) then begin

; now plot the cond probs for multi layer

set_plot, 'z'
device, set_resolution=[1000,1000]
;window, 1, xsize=900, ysize=800
;loadct, 4;5
!p.background=!d.n_colors-1
!p.multi=[0,16,6]
pos_kl=[0.05,0.20,0.35,0.50,0.65] & pos_kl=pos_kl+0.05
pos_ku=[0.15,0.30,0.45,0.60,0.75] & pos_ku=pos_ku+0.05
pos_jl=[0.05,0.30,0.55,0.80] ;& pos_jl=pos_jl+0.05
pos_jr=[0.20,0.45,0.70,0.95] ;& pos_jr=pos_jr+0.05

;pos_jl=[0.05,0.35,0.65] & pos_jl=pos_jl+0.05
;pos_jr=[0.25,0.55,0.85] & pos_jr=pos_jr+0.05

plot_array=fltarr(n_elements(dz_bin_mid), n_elements(cbtz_bin_mid))
;plot_array=fltarr(5,8)
;plot_array=fltarr(7,8)

xstr_ind=0
for j=0,3 do begin

ystr_ind=0

for k=0,4 do begin

plot_array[*,*]=cbz_dz_cond_freq[j,k,*,*]

cbz_dz_freq_string=strtrim(string((fix(cbz_dz_freq[j,k]*100.)),format='(i2)'))
cbz_dz_multi_freq_string=strtrim(string((fix(cbz_dz_multi_freq[j,k]*100.)),format='(i2)'))

;print, "cbz_dz_freq_string: ", cbz_dz_freq_string
;print, "cbz_dz_multi_freq_string: ", cbz_dz_multi_freq_string

!p.charsize=1.5
title_string=' base:'+yaxis_string[ystr_ind]+'-'+yaxis_string[ystr_ind+2]+$
      ' thick:'+xaxis_string[xstr_ind]+'-'+xaxis_string[xstr_ind+2]+ ' '+   $
      cbz_dz_freq_string + ' '+cbz_dz_multi_freq_string
      
xFormat = 2

xl=pos_jl[j] & yl=pos_kl[k] & xr=pos_jr[j] & yr=pos_ku[k]   ; upper right

plot_2d_image_999_maxminin2_axes, plot_array, $
     (dz_bin_mid-dz_dbin)/1000., (cbtz_bin_mid-cbtz_dbin)/1000.,    $
     n_elements(cbz_dz_multi_freq[*,0]), n_elements(cbz_dz_freq[0,*]), $
     title_string, 'Layer Thickness [km]', 'Base Height [km]', [xl,yl,xr,yr], $
    [xr+0.035,yl,+xr+0.04,yr], ' ',max(plot_array)+0.1,0.,xaxis_string,yaxis_string, $
    xFormat

lat_string=strtrim(string(location_lat,format='(f5.1)'),2)
lon_string=strtrim(string(location_lon,format='(f6.1)'),2)
;title_string='Cloudsat Multilayer Conditional Layer Occurrence for '+lat_string+' Latitude '+lon_string+' Longitude'
title_string='Cloudsat Multilayer Conditional Layer Occurrence' 
xyouts, 0.5,0.98, title_string,color=0,alignment=0.5, /normal,charsize=1.5

;xyouts, 0.5,0.94, 'All Cloud Layers',color=0,alignment=0.5, /normal,charsize=1.5

astring='Period Covered: '+datelist[0]+' to '+datelist[n_elements(datelist)-1]
xyouts, 0.5,0.96, astring,color=0,alignment=0.5, /normal,charsize=1.5

;astring='Averaging Region: '+strtrim(string(2.*dlat,format='(f5.1)'),2)+' Degrees Latitude by '$
;+strtrim(string(2.*dlon,format='(f5.1)'),2)+' Degrees Longitude'
astring='Averaging Region: Latitude (' + latStart_string + ', ' + latEnd_string + '), ' + $
                          'Longitude (' + lonStart_string + ', ' + lonEnd_string + ')' 
xyouts, 0.5,0.94, astring,color=0,alignment=0.5, /normal,charsize=1.5

astring='Total Cloud Coverage: '+strtrim(string(total_coverage,format='(f5.2)'),2);+$
;'  Estimated Column Attenuation Frequency: '+strtrim(string(atten_freq,format='(f5.2)'),2)
xyouts, 0.01,0.90, astring,color=0,alignment=0.0, /normal,charsize=1.0


caldat, systime(/julian), cur_month, cur_day, cur_year, cur_hour, cur_minute

a_string='Plot produced by Jay Mace on '+strtrim(string(fix(cur_month), format='(i2)'),2)+'/'+strtrim(string(cur_day),2)+$
'/'+strtrim(string(cur_year),2)

xyouts, 0.01,0.001, a_string,color=0,alignment=0.0, /normal,charsize=0.75

tvlct, red, green, blue, /get
;write_gif, output_path+'Cloud_Layer_Conditional_multi_base_RL'+'_'+datelist[0]+'-'+$
;datelist[n_elements(datelist)-1]+$
;'_'+lat_string+'_'+lon_string+'_'+strtrim(string(2.*dlat,format='(f5.1)'),2)+'by'$
;+strtrim(string(2.*dlon,format='(f5.1)'),2)+'.gif',tvrd();, red, green, blue
write_gif, output_path + 'Cloud_Layer_Conditional_multi_base_RL' + $
          '_' + datelist[0] + '-' + datelist[n_elements(datelist)-1] + $
          '_lat_' + latStart_string + '_' + latEnd_string + $
          '_lon_' + lonStart_string + '_' + lonEnd_string + $
          '_' + string_day_night + '.gif', $
          tvrd();, red, green, blue
          
ystr_ind=ystr_ind+2

endfor  ; for k=0,n_elements(dz_bin_mid)-1 do begin

xstr_ind=xstr_ind+2

endfor ;for j=0,n_elements(cbtz_bin_mid)-1 do begin

;;;endif      ; if (flag_Cloud_Layer_Conditional_multi_base_RL eq 1)