begin

  filename = "Fire_C_1997-2006_monthly_0.9.nc"
  
  in_file1 = addfile("surfdata_0.9x1.25_c070501.nc","r")
  in_file2 = addfile("lnd_0.9.nc","r")
  in_file3 = addfile("Fire_C_1997-2006_monthly_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_PFT
  lsmpftdatamask = in_file1->PFTDATA_MASK
  lsmpft = in_file1->PCT_PFT
  
  lndlandmask = in_file2->landmask
  lndarea = in_file2->area
  bigfire = in_file3->FIRE_C
  bigdate = in_file3->date
    
  latlondims = dimsizes(lsmlongxy)
  latdim = latlondims(0)
  londim = latlondims(1)
  
  latlonbigdims = dimsizes(bigfire)
  yeardim = latlonbigdims(0)
  monthdim = latlonbigdims(1)
  latbigdim = latlonbigdims(2)
  lonbigdim = latlonbigdims(3)

  dimNames = (/ "lon",  "lat", "month",  "year"   /)
  dimSizes = (/ londim, latdim, monthdim, yeardim /)
  dimUnlim = (/ False,  False,  False,    True  /)
  filedimdef( out_file, dimNames, dimSizes, dimUnlim )

  filevardef(out_file, "FIRE_C",typeof(lsmpft), (/ "year", "month", "lat", "lon" /))
  filevardef(out_file, "lat",typeof(lsmpft), (/ "lat" /))
  filevardef(out_file, "lon",typeof(lsmpft), (/ "lon" /))
  filevardef(out_file, "date",typeof(bigdate), (/ "year", "month" /))
  filevardef(out_file, "area",typeof(lsmpft), (/ "lat", "lon" /))

  newfire = new((/ yeardim, monthdim, latdim, londim /),float)
    
  do latbigindex = 0,latbigdim-1
      do lonbigindex = 0,lonbigdim-1
          newfire(:,:,latbigindex*2,lonbigindex*2) = bigfire(:,:,latbigindex,lonbigindex)
          newfire(:,:,latbigindex*2,lonbigindex*2+1) = bigfire(:,:,latbigindex,lonbigindex)
          newfire(:,:,latbigindex*2+1,lonbigindex*2) = bigfire(:,:,latbigindex,lonbigindex)
          newfire(:,:,latbigindex*2+1,lonbigindex*2+1) = bigfire(:,:,latbigindex,lonbigindex)
      end do
  end do
  
  out_file->FIRE_C@long_name = "Fire Emission C"
  out_file->FIRE_C@units = "gC/m2/month"
  out_file->FIRE_C@_FillValue = 1.0e36
  out_file->FIRE_C = (/ newfire /)
  out_file->lat = in_file2->lat
  out_file->lon = in_file2->lon
  out_file->date = in_file3->date
  out_file->area = in_file2->area
          
end