pro analyze_re_liq
 
;**************************************
;  Commons
;**************************************

;  Color common
common colors,r_orig,g_orig,b_orig,r_curr,g_curr,b_curr
 
;**************************************
;  How does this function behave
;**************************************

num=10000 
ze=(findgen(num)*0.01)-90.0
lwc=1.e-4
re_liq=make_array(/float,num,value=-9999)
 ;10.*alog10(ze_input[i]) gt 20.

for j=0,9999 do begin 
;  Calculate liquid effective radius
if 10.*alog10(ze[j]) lt 0.0 and ze[j] gt -8888. and lwc ge 1.e-4  then begin
   re_liq[j]=19.5*exp(0.0384*(10.*alog10(ze[j])))
   re_liq[j]=re_liq[j]+(0.20*re_liq[j])
endif else begin
   if ze[j] gt -8888. and lwc ge 1.e-4 then re_liq[j]=15.   ; was 8.4 6/29/05
endelse
endfor   
 
 
 ;********************************
;  Set up the plot device
;********************************

;  Set up the device
set_plot,'z'

;  xsize of the plot window
xwin=1000  

;  ysize of the plot window
ywin=700  

;  Size the window
device,set_resolution=[xwin,ywin]

;********************************
;  Set up the colors
;  0-253 rainbow
;  254   gray
;  255   white
;********************************

;  Load the color table
loadct,39   

;  Make gray for missing data in color 253
r_curr[254]=200
g_curr[254]=200
b_curr[254]=200

;  Make the background white
!p.background=!d.n_colors-1  

erase

;********************************
;  Position vectors
;********************************

xl=0.06  &  xr=0.99
yb=0.06  &  yt=0.97
sx=0.01        ;x buffer space
sy=0.06        ;y buffer space
numplots_x=1  ;number of plots in x direction
numplots_y=2  ;number of plots in y direction
position_plots,xl,xr,yb,yt,sx,sy,numplots_x,numplots_y,pos

plot,ze,re_liq,color=0,position=pos[0,*],xstyle=1,ystyle=1,psym=-4,$
  xtitle='ze',ytitle='re_liq',xrange=[-1,2]
  
plot,10.*alog10(ze),re_liq,color=0,position=pos[1,*],xstyle=1,ystyle=1,psym=-4,$
  xtitle='10*alog10(ze)',ytitle='re_liq',/noerase;,xrange=[-5,5]
  

 write_gif,'/Users/u0079358/Documents/fifth_cloudret/flux_closureL/analyze_re_liq.gif',tvrd()
 
 
    
plot,ze,re_liq,color=0
stop


end