Subroutine vom_add_yearly

2068 : 
2069 : !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
2070 : !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
2071 : !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
2072 : 
2073 :       subroutine vom_add_yearly ()
2074 :       use vom_vegwat_mod
2075 :       implicit none
2076 : 
2077 :       if (fyear(nday) .eq. nyear) then
2078 :         rain_y   = rain_y  + rain_d(nday)  ! in [mm]
2079 :         par_y    = par_y   + par_d(nday)
2080 :         srad_y   = srad_y  + srad_d(nday)  ! srad originally in MJ/day
2081 :         vd_y     = vd_y    + vd_d / 24.d0
2082 :         etm_y    = etm_y   + (etmt_d + etmg_d) * 1000.d0  ! in[mm]
2083 :         esoil_y  = esoil_y + esoil_d           * 1000.d0  ! in [mm]
2084 : !       * for grasses
2085 :         etmg_y   = etmg_y  + etmg_d * 1000.d0 ! in [mm]
2086 :         assg_y   = assg_y  + assg_d(2,2,2)
2087 :         rlg_y    = rlg_y   + rlg_d
2088 :         rrg_y    = rrg_y   + rrg_d      * 3600.d0 * 24.d0
2089 :         cpccg_y  = cpccg_y + cpccg_d(2) * 3600.d0 * 24.d0
2090 :         tcg_y    = tcg_y   + tcg_d(2,2)   * 3600.d0 * 24.d0
2091 : !       * for trees
2092 :         etmt_y   = etmt_y  + etmt_d * 1000.d0  ! in [mm]
2093 :         asst_y   = asst_y  + asst_d(2,2)
2094 :         rlt_y    = rlt_y   + rlt_d
2095 :         rrt_y    = rrt_y   + rrt_d     * 3600.d0 * 24.d0
2096 :         cpcct_y  = cpcct_y + q_cpcct_d * 3600.d0 * 24.d0
2097 :         tct_y    = tct_y   + q_tct_d(2)   * 3600.d0 * 24.d0
2098 :       else
2099 :         nyear    = fyear(nday)
2100 :         rain_y   = rain_d(nday)
2101 :         par_y    = par_d(nday)
2102 :         srad_y   = srad_d(nday)              ! srad originally in MJ/day
2103 :         vd_y     = vd_d / 24.d0
2104 :         etm_y    = (etmg_d + etmt_d) * 1000.d0
2105 :         esoil_y  = esoil_d * 1000.d0
2106 : !       * for grasses
2107 :         etmg_y   = etmg_d * 1000.d0
2108 :         assg_y   = assg_d(2,2,2)
2109 :         rlg_y    = rlg_d
2110 :         rrg_y    = rrg_d      * 3600.d0 * 24.d0
2111 :         cpccg_y  = cpccg_d(2) * 3600.d0 * 24.d0
2112 :         tcg_y    = tcg_d(2,2)   * 3600.d0 * 24.d0
2113 : !       * for trees
2114 :         etmt_y = etmt_d * 1000.d0
2115 :         asst_y   = asst_d(2,2)
2116 :         rlt_y    = rlt_d
2117 :         rrt_y    = rrt_d     * 3600.d0 * 24.d0
2118 :         cpcct_y  = q_cpcct_d * 3600.d0 * 24.d0
2119 :         tct_y    = q_tct_d(2)   * 3600.d0 * 24.d0
2120 :       endif
2121 : 
2122 :       return
2123 :       end subroutine vom_add_yearly