begin

  filename = "LAI_2000-2005_MONS_0.5.nc"
  
  in_file1 = addfile("surfdata_360x720_mcrop2000.nc","r")
  in_file2 = addfile("lnd_0.5.nc","r")
  in_file3 = addfile("LAI_2000-2005_MONS_1.9.nc","r")
  out_file = addfile(filename,"c")

  fileAttr = True
  fileAttr@title = filename
  fileAttr@creation_date = systemfunc("date")
  fileAttr@source = "Peter Lawrence, NCAR"
  fileattdef(out_file,fileAttr)
  
  lsmedgen = in_file1->EDGEN
  lsmedgee = in_file1->EDGEE
  lsmedges = in_file1->EDGES
  lsmedgew = in_file1->EDGEW
  lsmlatn = in_file1->LATN
  lsmlone = in_file1->LONE
  lsmlats = in_file1->LATS
  lsmlonw = in_file1->LONW
  lsmarea = in_file1->AREA
  lsmnumlon = in_file1->NUMLON
  lsmlongxy = in_file1->LONGXY
  lsmlatixy = in_file1->LATIXY
  lsmlandfracpft = in_file1->LANDFRAC
  lsmpftdatamask = in_file1->LANDMASK
  lsmpft = in_file1->PCT_PFT
  
  lndlandmask = in_file2->landmask
  monthlai = in_file3->LAI
  monthlaitime = in_file3->time 
    
  latlondims = dimsizes(lsmlongxy)
  latdim = latlondims(0)
  londim = latlondims(1)
  
  latlonbigdims = dimsizes(monthlai)
  timedim = latlonbigdims(0)
  latbigdim = latlonbigdims(1)
  lonbigdim = latlonbigdims(2)

  dimNames = (/ "lon",  "lat", "time"      /)
  dimSizes = (/ londim, latdim, timedim /)
  dimUnlim = (/ False,  False, False  /)
  filedimdef( out_file, dimNames, dimSizes, dimUnlim )

  filevardef(out_file, "LAI",typeof(lsmpft), (/ "time", "lat", "lon" /))
  filevardef(out_file, "time", typeof(monthlaitime), (/ "time" /))
  filevardef(out_file, "lat",typeof(lsmpft), (/ "lat" /))
  filevardef(out_file, "lon",typeof(lsmpft), (/ "lon" /))

  cleanmonthlai = where(monthlai .ne. monthlai@_FillValue,monthlai,0.0)
  cleanmonthlai@missing_value = cleanmonthlai@_FillValue
  
  newmonthlai = new((/ timedim, latdim, londim /),float)
  newmonthlai = (/ 0.0 /)
  
  
  do latindex = 0,latdim-1
      do lonindex = 0,londim-1
      
          latbigindex = floattointeger(latindex * 0.5 / 1.875)
          lonbigindex = floattointeger(lonindex * 0.5 / 2.5)

          if (.not. ismissing(cleanmonthlai(0,latbigindex,lonbigindex))) then
              newmonthlai(:,latindex,lonindex) = cleanmonthlai(:,latbigindex,lonbigindex)
	  end if
      end do
  end do
  
  out_file->LAI@long_name = "Leaf Area Index"
  out_file->LAI@units = "none"
  out_file->LAI@_FillValue = 1.0e36
  do timeindex = 0,timedim-1
  out_file->LAI(timeindex,:,:) = (/ where(lndlandmask.eq.1,newmonthlai(timeindex,:,:),out_file->LAI@_FillValue) /)
  end do
  out_file->time = in_file3->time
  out_file->lat = in_file2->lat
  out_file->lon = in_file2->lon
          
end