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