Re-commited various bug-fixes etc.
Merging Ryan's bug "re-fixes"...
Merging Ryan's code
small changes to fast-timestep output to include variables for NACP intercomparison
merge mainline revision 2
Added an ANNUAL output option.  Currently this writes the states of slow state variables down to the cohort level, but does not yet have annual averages of important ecosystem or demographic rates

Also added a number of new variables output to the ANALYSIS (fast time step) files, which were added for the NACP (North American Carbon Plan)  model-data intercomparison.
added yearly output HDF5 file.  ED2IN flag is IYOUTPUT.  focuses currently on system state of slow, low-level variables (e.g. cohort dbh, patch soil carbon).  could easily add annual averages, demography (i.e. all the variables that were previously in the misc. annual ascii files)
Merging coupled ED with BRAMS-4.0.6.6, the revision right before I started changing 
    cumulus/microphysics/thermodynamics. Probably need some work at the initialisation...
Including the call to ED namelist when running ED, always a good thing...
Also, fixing some legacy on therm_lib.f90. Removed thermlibsupp.f90 because these 
functions will eventually be part of BRAMS, for microphysics.
Including the old cumulus scheme files back to bzr...
1. Adding back the old way to compute qsw (ed_params.f90) for mid-latitude PFTs.
2. Yeonjoo's bug fix in DBH loop (forestry.f90)
Bringing QSW and forestry changes.
1. Including temporary "bug inclusion" for mid-latitude PFTs when computing qsw;
2. Bug fix on forestry.f90
3. Switching back to old cumulus parameterization on BRAMS for now.
added time blocking on fast output, set up basic variables for blocking on history.

still need to:
 * check if it works for output frequencies other than 1 day
 * check if output is correct
 * implement history blocking
merge
A few bug fixes.
   1) Changed the order of some routines in the coupled driver so that
      memory was allocated at the correct time.
   2) Modified the grid search routine for type 1 ed restart files, such      that if no files were found on a given node, then it would search       the whole earth as a last resort.  This is necessary in large para      llel runs
   3) Modified the h5 io routines to prevent panicking on nodes that hav      e no polygons, ie nodes with only ocean.
   4) One other thing, but I forget right now
Fixed an soil/landcover search bug.
Merging Ryan's code.
Merging Ryan's code
1. BRAMS changes:
  a. Major re-structuring of BRAMS thermodynamics, including:
    - Saturation pressure is now computed using Murphy-Koop (2005) equations. Their
      equations put the triple point temperature at 273.16K rather than 273.00K;
    - Removed duplicated functions (for example, there used to be two functions to
      compute saturation mixing ratio, rs and rs?f. Now only rs?f remained). Changed
      all arguments to SI. This also affects the density functions!
    - Removed various "magic numbers" from the code, switching by rconstants.
    - The triple point is now defined by t3ple and es3ple, and t3ple has a different mean-
      ing from t00 (including internal energy calculations).
    - Ice-liquid potential temperature is now solved with the exponential form (Tripoli
      and Cotton 1981, but without expanding in Taylor series).
    - Brunt-Vaisala frequency subroutine (bruvais) was simplified since we now have the
      virtual temperature function. Also Nakanishi-Niino was adjusted to handle
      condensation properly;
    - Theta-E and Theta-IL are now solved through traditional numeric methods.
    - Virtual temperature considers condensation inside the run (only during makevfile it
      is assumed single phase).
  b. Converting therm_lib.f90 into a module, so it can handle optional arguments properly.
  c. Removed some TEB interfaces that were causing trouble. They are now passed as "normal"
     arguments rather than targets.

2. ED changes:
  a. Switching ED thermodynamics to be consistent with BRAMS. It is also a module now.
  b. Small fix on ed_state_vars, the old code was causing troubles on openmpi;

3. Coupled model is running in serial mode now.

4. Makefile from ED and EDBRAMS was changed. F_COMMAND and C_COMMAND were switched by
   F90_COMMAND, FPP_COMMAND and CXX_COMMAND. CXX is exactly the same as the old one, but
   now we distinguish F90_COMMAND for .f90 files and FPP_COMMAND for .F90 files, which 
   require the pre-processor. The -D variables are passed only in the FPP command.