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