C NCLFORTSTART
       subroutine bin_stats(sf,nlat,nlon,missing,tb_bins,ntb_bins,
     +     ntb_binsc,lat_bins,nlat_bins,nlat_binsc,gbin,lsm,
     +	    lat2d,outbin,outbin2)
       integer j,tb,la,jl,jk
       integer nlat_bins,ntb_bins,nlat_binsc,ntb_binsc
       integer tb_bins(ntb_bins)
       integer lat_bins(nlat_bins)
       integer nlat,nlon,sf
       real    missing
       real    gbin(nlat,nlon),lsm(nlat,nlon),lat2d(nlat,nlon)
       real    outbin(sf,ntb_binsc),
     + outbin2(sf,nlat_binsc,ntb_binsc)
C NCLEND