Subroutine vom_write_hourly
1774 :
1775 : !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1776 : !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1777 : !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1778 :
1779 : subroutine vom_write_hourly ()
1780 : use vom_vegwat_mod
1781 : implicit none
1782 :
1783 : CHARACTER(60) :: hourlyformat
1784 : CHARACTER(3) :: str
1785 :
1786 : if (fyear(nday) .ge. i_firstyear .and. fyear(nday) .le. i_lastyear) then
1787 : ! * internal write to convert from number to string
1788 : write(str,'(i3)') wlayer_
1789 : ! * includes a column for each sublayer
1790 : hourlyformat = '(I6,I6,I4,I7,I5,'//str//'E14.6)'
1791 :
1792 : write(kfile_resultshourly,'(I6,I7,I7,I7,I7,22E15.5)') &
1793 : & fyear(nday), fmonth(nday), fday(nday), nday, nhour, &
1794 : & rain_h(th_), tair_h(th_), par_h(th_), vd_h(th_), esoil_h, &
1795 : & o_cai + pcg_d(2), jmax25t_d(2), jmax25g_d(2), mqt_, &
1796 : & rlt_h(2,2) + rlg_h(2,2,2), lambdat_d, lambdag_d, rrt_d + rrg_d, &
1797 : & asst_h(2,2), assg_h(2,2,2), etmt_h, etmg_h, su__(1), zw_, wsnew, &
1798 : & spgfcf_h, infx_h
1799 :
1800 : write(kfile_delzhourly,hourlyformat) fyear(nday), &
1801 : & fmonth(nday), fday(nday), nday, nhour, s_delz(1:wlayer_)
1802 :
1803 : write(kfile_ruptkthourly,hourlyformat) fyear(nday), &
1804 : & fmonth(nday), fday(nday), nday, nhour, ruptkt_h(1:wlayer_)
1805 :
1806 : write(kfile_suhourly,hourlyformat) fyear(nday), &
1807 : & fmonth(nday), fday(nday), nday, nhour, su__(1:wlayer_)
1808 : endif
1809 :
1810 : return
1811 : end subroutine vom_write_hourly