Subroutine vom_add_hourly

1707 : 
1708 : !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1709 : !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1710 : !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1711 : 
1712 :       subroutine vom_add_hourly ()
1713 :       use vom_vegwat_mod
1714 :       implicit none
1715 : 
1716 :       REAL*8 :: asst__(3,3)
1717 :       REAL*8 :: assg__(3,3,3)
1718 :       INTEGER:: ii
1719 : 
1720 :     do ii = 1,3 !loop for LAI values
1721 :       asst__(:,ii) = (4.d0 * ca_h(th_) * gstomt + 8.d0 * gammastar        &
1722 :         &          * gstomt + jactt(:,ii) - 4.d0 * rlt_h(:,ii) - SQRT((-4.d0    &
1723 :         &          * ca_h(th_) * gstomt + 8.d0 * gammastar * gstomt       &
1724 :         &          + jactt(:,ii) - 4.d0 * rlt_h(:,ii)) ** 2.d0 + 16.d0          &
1725 :         &          * gammastar * gstomt * (8.d0 * ca_h(th_) * gstomt      &
1726 :         &          + jactt(:,ii) + 8.d0 * rlt_h(:,ii)))) / 8.d0  ! (3.22) ; (Out[319])
1727 :     end do
1728 :         asst_h(:,:) = asst_h(:,:) + asst__(:,:) * dt_
1729 : 
1730 :     do ii = 1,3 !loop for LAI values
1731 :         assg__(:,:,ii) = (4.d0 * ca_h(th_) * gstomg(:,:,ii) + 8.d0 * gammastar &
1732 :         &         * gstomg(:,:,ii) + jactg(:,:,ii) - 4.d0 * rlg_h(:,:,ii)       &
1733 :         &         - SQRT((-4.d0 * ca_h(th_) * gstomg(:,:,ii) + 8.d0       &
1734 :         &         * gammastar * gstomg(:,:,ii) + jactg(:,:,ii) - 4.d0        &
1735 :         &         * rlg_h(:,:,ii)) ** 2.d0 + 16.d0 * gammastar            &
1736 :         &         * gstomg(:,:,ii) * (8.d0 * ca_h(th_) * gstomg(:,:,ii)      &
1737 :         &         + jactg(:,:,ii) + 8.d0 * rlg_h(:,:,ii)))) / 8.d0  ! (3.22); (Out[319])
1738 :     end do
1739 :       assg_h(:,:,:) = assg_h(:,:,:) + assg__(:,:,:) * dt_
1740 :       ruptkt_h(:) = ruptkt_h(:) + ruptkt__(:) * dt_
1741 :       ruptkg_h(:) = ruptkg_h(:) + ruptkg__(:) * dt_
1742 :       if (optmode .eq. 0) then
1743 :         spgfcf_h    = spgfcf_h    + dt_ * spgfcf__
1744 :         infx_h      = infx_h      + dt_ * infx__
1745 :         io_h        = io_h        + dt_ * io__
1746 :         esoil_h     = esoil_h     + dt_ * esoil__
1747 :         etmt_h      = etmt_h      + dt_ * etmt__
1748 :         etmg_h      = etmg_h      + dt_ * etmg__(2,2,2)
1749 :         sumruptkt_h = sumruptkt_h + dt_ * SUM(ruptkt__(:))
1750 :       endif
1751 : 
1752 :       return
1753 :       end subroutine vom_add_hourly