pro Mean_Nd_Capricorn22

  fid=ncdf_open('Mean_Nd_Capricorn22.cdf')
  xid=ncdf_varid(fid,'julian_day') & ncdf_varget,fid,xid,julian_day
  xid=ncdf_varid(fid,'nd_mean') & ncdf_varget,fid,xid,nd_mean
  xid=ncdf_varid(fid,'nd_sdv') & ncdf_varget,fid,xid,nd_sdv
  xid=ncdf_varid(fid,'lat_vec') & ncdf_varget,fid,xid,lat_vec
  xid=ncdf_varid(fid,'min_lat') & ncdf_varget,fid,xid,min_lat
  xid=ncdf_varid(fid,'max_lat') & ncdf_varget,fid,xid,max_lat
  xid=ncdf_varid(fid,'nd_num') & ncdf_varget,fid,xid,nd_num      
  ncdf_close,fid

do_plot='yes'
if do_plot eq 'yes' then begin
  
;  Plot size
pxdim=800 & pydim=800
;  Position the plots
xl=0.12  &  xr=0.95
yb=0.06  &   yt=0.915
sx=0.08
sy=0.08
numplots_x=1
numplots_y=3
position_plots,xl,xr,yb,yt,sx,sy,numplots_x,numplots_y,pos
;  Set up the plot space
pnum=0
p0=plot([0,1],[0,1],position=pos[pnum,*],axis_style=4,/nodata,/buffer,$
  dimensions=[pxdim,pydim])

dummy=label_date(date_format=['%M %D'])
fs1=12
xinterval=5
dx=0.01 & dy=0.01

;*****  Plot 0 **********  
pnum=0
p0=plot(julian_day,nd_num,/current,position=pos[pnum,*],$
  xtickunits='days',xtickinterval=xinterval,$
  xrange=[julian_day[0]-2D,julian_day[-1]+2D],$
  title='Capricorn ll Number of Retrievals',$
  /hist,xtickformat='label_date',ytitle='Number of Retrievals')
t1=text(pos[pnum,0]-9.0*dx,pos[pnum,3]+2*dy,'c)',font_size=14)

;*****  Plot 1 ********** 
pnum=1
numtimes=n_elements(lat_vec)
yerror=make_array(2,numtimes,/float)
yerror[0,*]=lat_vec-min_lat
yerror[1,*]=max_lat-lat_vec
;yerror=max_lat-lat_vec
p0=errorplot(julian_day,lat_vec,yerror,/current,position=pos[pnum,*],$
;p0=plot(julian_day,lat_vec,/current,position=pos[pnum,*],$
  xtickunits='days',xtickinterval=xinterval,$
  xrange=[julian_day[0]-2D,julian_day[-1]+2D],$
  linestyle=6,symbol='o',sym_size=0.5,/sym_filled,$
  title='Capricorn ll Daily Mean Ship Latitude',$
  xtickformat='label_date',ytitle='Mean Ship Latitude')
;p0=plot(julian_day,min_lat,/overplot,color='blue',linestyle=6,symbol='.')
;p0=plot(julian_day,max_lat,/overplot,color='red',linestyle=6,symbol='.')    
t1=text(pos[pnum,0]-9.0*dx,pos[pnum,3]+2*dy,'b)',font_size=14)

;*****  Plot 3 **********  
pnum=2
;numtimes=n_elements(lat_vec)
;yerror=make_array(2,numtimes,/float)
;yerror[0,*]=lat_vec-min_lat
;yerror[1,*]=max_lat-lat_vec
;yerror=nd_mean-nd_sdv
yerror=nd_sdv
p0=errorplot(julian_day,nd_mean,yerror,/current,position=pos[pnum,*],$
;p0=plot(julian_day,nd_mean,/current,position=pos[pnum,*],$
  xtickunits='days',xtickinterval=xinterval,$
  xrange=[julian_day[0]-2D,julian_day[-1]+2D],$
  linestyle=6,symbol='o',/sym_filled,$  
  title='Capricorn ll Mean Cloud Droplet Number (per cubic cm)',$
  xtickformat='label_date',ytitle='Cloud Nd (per cubic cm)')
t1=text(pos[pnum,0]-9.0*dx,pos[pnum,3]+2*dy,'a)',font_size=14)
      
p0.save,'mean_nd_capricorn.png',height=pydim   
p0.save,'mean_nd_capricorn.eps',height=pydim
stop
endif  ;end of do_plot

end