SURFEX  V8_0
Surfex V8_0 release
 All Classes Files Functions Variables
modd_diag_evap_isban.F90
Go to the documentation of this file.
1 !SFX_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
2 !SFX_LIC This is part of the SURFEX software governed by the CeCILL-C licence
3 !SFX_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
4 !SFX_LIC for details. version 1.
5 !########################
7 !########################
8 !
9 !!**** *MODD_DIAG_ISBA - declaration of packed surface parameters for ISBA scheme
10 !!
11 !! PURPOSE
12 !! -------
13 !
14 !!
15 !!** IMPLICIT ARGUMENTS
16 !! ------------------
17 !! None
18 !!
19 !! REFERENCE
20 !! ---------
21 !!
22 !! AUTHOR
23 !! ------
24 !! P. Le Moigne *Meteo France*
25 !!
26 !! MODIFICATIONS
27 !! -------------
28 !! Original 07/11/03
29 !! P. Samuelsson 04/2012 MEB
30 !
31 !* 0. DECLARATIONS
32 ! ------------
33 !
34 !
35 USE yomhook ,ONLY : lhook, dr_hook
36 USE parkind1 ,ONLY : jprb
37 !
38 IMPLICIT NONE
39 
41 !------------------------------------------------------------------------------
42 !
43  LOGICAL :: LSURF_EVAP_BUDGET ! flag for all terms of evaporation
44  LOGICAL :: LSURF_BUDGETC ! flag for surface cumulated energy budget
45  LOGICAL :: LRESET_BUDGETC ! flag for surface cumulated energy budget
46  LOGICAL :: LWATER_BUDGET ! flag for isba water budget including input
47  ! fluxes (rain and snow) and reservoir tendencies
48 !
49 !* variables for each patch
50 !
51  REAL, POINTER, DIMENSION(:,:) :: XLEG ! latent heat of evaporation over the ground (W/m2)
52  REAL, POINTER, DIMENSION(:,:) :: XLEGI ! surface soil ice sublimation (W/m2)
53  REAL, POINTER, DIMENSION(:,:) :: XLEV ! latent heat of evaporation over vegetation (W/m2)
54  REAL, POINTER, DIMENSION(:,:) :: XLES ! latent heat of sublimation over the snow (W/m2)
55  REAL, POINTER, DIMENSION(:,:) :: XLESL ! latent heat of evaporation over the snow (W/m2)
56  REAL, POINTER, DIMENSION(:,:) :: XLER ! evaporation from canopy water interception (W/m2)
57  REAL, POINTER, DIMENSION(:,:) :: XLETR ! evapotranspiration of the vegetation (W/m2)
58  REAL, POINTER, DIMENSION(:,:) :: XEVAP ! evapotranspiration (kg/m2/s)
59  REAL, POINTER, DIMENSION(:,:) :: XSUBL ! sublimation (kg/m2/s)
60  REAL, POINTER, DIMENSION(:,:) :: XSNDRIFT ! blowing snow sublimation (ES or Crocus) (kg/m2/s)
61  REAL, POINTER, DIMENSION(:,:) :: XDRAIN ! soil drainage flux (kg/m2/s)
62  REAL, POINTER, DIMENSION(:,:) :: XQSB ! lateral subsurface flux (dif option) (kg/m2/s)
63  REAL, POINTER, DIMENSION(:,:) :: XRUNOFF ! sub-grid and supersaturation runoff (kg/m2/s)
64  REAL, POINTER, DIMENSION(:,:) :: XHORT ! sub-grid Horton runoff from the SGH scheme (kg/m2/s)
65  REAL, POINTER, DIMENSION(:,:) :: XRRVEG ! precipitation intercepted by the vegetation (kg/m2/s)
66  REAL, POINTER, DIMENSION(:,:) :: XMELT ! snow melt (kg/m2/s)
67  REAL, POINTER, DIMENSION(:,:) :: XIFLOOD ! Floodplains infiltration (kg/m2/s)
68  REAL, POINTER, DIMENSION(:,:) :: XPFLOOD ! Precipitation intercepted by the floodplains (kg/m2/s)
69  REAL, POINTER, DIMENSION(:,:) :: XLE_FLOOD ! Floodplains evapotration (W/m2)
70  REAL, POINTER, DIMENSION(:,:) :: XLEI_FLOOD ! Floodplains evapotration (W/m2)
71  REAL, POINTER, DIMENSION(:,:) :: XDRIP ! dripping from the vegetation reservoir (kg/m2/s)
72  REAL, POINTER, DIMENSION(:,:) :: XIRRIG_FLUX ! irrigation rate (as soil input) (kg/m2/s)
73 !
74  REAL, POINTER, DIMENSION(:,:) :: XGPP ! Gross Primary Production (kgCO2/m2/s)
75  REAL, POINTER, DIMENSION(:,:) :: XRESP_AUTO ! Autotrophic respiration (kgCO2/m2/s)
76  REAL, POINTER, DIMENSION(:,:) :: XRESP_ECO ! Ecosystem respiration (kgCO2/m2/s)
77 !
78  REAL, POINTER, DIMENSION(:,:) :: XLEVCV ! MEB: total evapotranspiration from vegetation canopy overstory [W/m2]
79  REAL, POINTER, DIMENSION(:,:) :: XLESC ! MEB: total snow sublimation from vegetation canopy overstory [W/m2]
80  REAL, POINTER, DIMENSION(:,:) :: XLETRGV ! MEB: transpiration from understory vegetation [W/m2]
81  REAL, POINTER, DIMENSION(:,:) :: XLETRCV ! MEB: transpiration from overstory canopy vegetation [W/m2]
82  REAL, POINTER, DIMENSION(:,:) :: XLERGV ! MEB: interception evaporation from understory vegetation [W/m2]
83  REAL, POINTER, DIMENSION(:,:) :: XLELITTER ! MEB: interception evaporation from understory vegetation [W/m2]
84  REAL, POINTER, DIMENSION(:,:) :: XLELITTERI ! MEB: interception evaporation from understory vegetation [W/m2]
85  REAL, POINTER, DIMENSION(:,:) :: XDRIPLIT !
86  REAL, POINTER, DIMENSION(:,:) :: XRRLIT !
87  REAL, POINTER, DIMENSION(:,:) :: XLERCV ! MEB: interception evaporation from overstory canopy vegetation [W/m2]
88  REAL, POINTER, DIMENSION(:,:) :: XLE_V_C ! MEB: latent heat flux from vegetation canopy overstory [W/m2]
89  REAL, POINTER, DIMENSION(:,:) :: XLE_G_C ! MEB: latent heat flux from understory [W/m2]
90  REAL, POINTER, DIMENSION(:,:) :: XLE_C_A ! MEB: latent heat flux from canopy air space to the atmosphere [W/m2]
91  ! NOTE total latent heat flux to the atmosphere also possibly
92  ! includes a contribution from snow covering the canopy
93  REAL, POINTER, DIMENSION(:,:) :: XLE_N_C ! MEB: latent heat flux from the snow on the ground [W/m2]
94  ! NOTE total latent heat flux from the snowpack
95  ! possibly includes a contribution from snow covering the canopy
96  REAL, POINTER, DIMENSION(:,:) :: XSWNET_V ! MEB: net vegetation canopy shortwave radiation [W/m2]
97  REAL, POINTER, DIMENSION(:,:) :: XSWNET_G ! MEB: net ground shortwave radiation [W/m2]
98  REAL, POINTER, DIMENSION(:,:) :: XSWNET_N ! MEB: net snow shortwave radiation [W/m2]
99  REAL, POINTER, DIMENSION(:,:) :: XSWNET_NS ! MEB: net snow shortwave radiation for *surface* layer
100  ! (i.e. net snow shortwave radiation less absorbed radiation) [W/m2]
101  REAL, POINTER, DIMENSION(:,:) :: XLWNET_V ! MEB: net vegetation canopy longwave radiation [W/m2]
102  REAL, POINTER, DIMENSION(:,:) :: XLWNET_G ! MEB: net ground longwave radiation [W/m2]
103  REAL, POINTER, DIMENSION(:,:) :: XLWNET_N ! MEB: net snow longwave radiation [W/m2]
104  REAL, POINTER, DIMENSION(:,:) :: XH_V_C ! MEB: sensible heat flux from vegetation canopy overstory [W/m2]
105  REAL, POINTER, DIMENSION(:,:) :: XH_G_C ! MEB: sensible heat flux from understory [W/m2]
106  REAL, POINTER, DIMENSION(:,:) :: XH_C_A ! MEB: sensible heat flux from canopy air space to the atmosphere [W/m2]
107  ! NOTE total sensible heat flux to the atmosphere also possibly
108  ! includes a contribution from snow covering the canopy
109  REAL, POINTER, DIMENSION(:,:) :: XH_N_C ! MEB: sensible heat flux from the snow on the ground [W/m2]
110  ! NOTE total sensible heat flux from the snowpack
111  ! possibly includes a contribution from snow covering the canopy
112  REAL, POINTER, DIMENSION(:,:) :: XSR_GN ! MEB: snow unloading rate from the overstory reservoir [kg/m2/s]
113  REAL, POINTER, DIMENSION(:,:) :: XMELTCV ! MEB: snow melt rate from the overstory snow reservoir [kg/m2/s]
114  REAL, POINTER, DIMENSION(:,:) :: XFRZCV ! MEB: snow refreeze rate from the overstory snow reservoir [kg/m2/s]
115  REAL, POINTER, DIMENSION(:,:) :: XSWDOWN_GN ! MEB: total shortwave radiation transmitted through the canopy
116  ! reaching the snowpack/ground understory [W/m2]
117  REAL, POINTER, DIMENSION(:,:) :: XLWDOWN_GN ! MEB: total shortwave radiation transmitted through and emitted by the canopy
118  ! reaching the snowpack/ground understory (explicit part) [W/m2]
119 !
120  REAL, POINTER, DIMENSION(:,:) :: XDWG ! liquid soil moisture time tendencies (kg/m2/s)
121  REAL, POINTER, DIMENSION(:,:) :: XDWGI ! solid soil moisture time tendencies (kg/m2/s)
122  REAL, POINTER, DIMENSION(:,:) :: XDWR ! canopy water time tendencies (kg/m2/s)
123  REAL, POINTER, DIMENSION(:,:) :: XDSWE ! snow water equivalent time tendencies (kg/m2/s)
124  REAL, POINTER, DIMENSION(:,:) :: XWATBUD ! ISBA water budget (kg/m2/s)
125 !
126 !* average variables
127 !
128  REAL, POINTER, DIMENSION(:) :: XAVG_LEG ! latent heat of evaporation over the ground (W/m2)
129  REAL, POINTER, DIMENSION(:) :: XAVG_LEGI ! surface soil ice sublimation (W/m2)
130  REAL, POINTER, DIMENSION(:) :: XAVG_LEV ! latent heat of evaporation over vegetation (W/m2)
131  REAL, POINTER, DIMENSION(:) :: XAVG_LES ! latent heat of sublimation over the snow (W/m2)
132  REAL, POINTER, DIMENSION(:) :: XAVG_LESL ! latent heat of evaporation over the snow (W/m2)
133  REAL, POINTER, DIMENSION(:) :: XAVG_LER ! evaporation from canopy water interception (W/m2)
134  REAL, POINTER, DIMENSION(:) :: XAVG_LETR ! evapotranspiration of the vegetation (W/m2)
135  REAL, POINTER, DIMENSION(:) :: XAVG_EVAP ! evapotranspiration (kg/m2/s)
136  REAL, POINTER, DIMENSION(:) :: XAVG_SUBL ! sublimation (kg/m2/s)
137  REAL, POINTER, DIMENSION(:) :: XAVG_SNDRIFT ! blowing snow sublimation (ES or Crocus) (kg/m2/s)
138  REAL, POINTER, DIMENSION(:) :: XAVG_DRAIN ! soil drainage flux (kg/m2/s)
139  REAL, POINTER, DIMENSION(:) :: XAVG_QSB ! lateral subsurface flux (dif option) (kg/m2/s)
140  REAL, POINTER, DIMENSION(:) :: XAVG_RUNOFF ! sub-grid and supersaturation runoff (kg/m2/s)
141  REAL, POINTER, DIMENSION(:) :: XAVG_HORT ! sub-grid Horton runoff from the SGH scheme (kg/m2/s)
142  REAL, POINTER, DIMENSION(:) :: XAVG_DRIP ! dripping from the vegetation reservoir (kg/m2/s)
143  REAL, POINTER, DIMENSION(:) :: XAVG_MELT ! snow melt (kg/m2/s)
144  REAL, POINTER, DIMENSION(:) :: XAVG_IFLOOD ! Floodplains infiltration (kg/m2/s)
145  REAL, POINTER, DIMENSION(:) :: XAVG_PFLOOD ! Precipitation intercepted by the floodplains (kg/m2/s)
146  REAL, POINTER, DIMENSION(:) :: XAVG_LE_FLOOD ! Floodplains evapotration (W/m2)
147  REAL, POINTER, DIMENSION(:) :: XAVG_LEI_FLOOD! Floodplains evapotration (W/m2)
148  REAL, POINTER, DIMENSION(:) :: XAVG_RRVEG ! precipitation intercepted by the vegetation (kg/m2/s)
149  REAL, POINTER, DIMENSION(:) :: XAVG_IRRIG_FLUX! irrigation rate (as soil input) (kg/m2/s)
150 !
151  REAL, POINTER, DIMENSION(:) :: XAVG_GPP ! Gross Primary Production (kgCO2/m2/s)
152  REAL, POINTER, DIMENSION(:) :: XAVG_RESP_AUTO! Autotrophic respiration (kgCO2/m2/s)
153  REAL, POINTER, DIMENSION(:) :: XAVG_RESP_ECO ! Ecosystem respiration (kgCO2/m2/s)
154 !
155  REAL, POINTER, DIMENSION(:) :: XAVG_LEVCV ! MEB: total evapotranspiration from vegetation canopy overstory [W/m2]
156  REAL, POINTER, DIMENSION(:) :: XAVG_LESC ! MEB: total snow sublimation from vegetation canopy overstory [W/m2]
157  REAL, POINTER, DIMENSION(:) :: XAVG_LETRGV ! MEB: transpiration from understory vegetation [W/m2]
158  REAL, POINTER, DIMENSION(:) :: XAVG_LETRCV ! MEB: transpiration from overstory canopy vegetation [W/m2]
159  REAL, POINTER, DIMENSION(:) :: XAVG_LERGV ! MEB: interception evaporation from understory vegetation [W/m2]
160  REAL, POINTER, DIMENSION(:) :: XAVG_LELITTER ! MEB: interception evaporation from understory vegetation [W/m2]
161  REAL, POINTER, DIMENSION(:) :: XAVG_LELITTERI ! MEB: interception evaporation from understory vegetation [W/m2]
162  REAL, POINTER, DIMENSION(:) :: XAVG_DRIPLIT !
163  REAL, POINTER, DIMENSION(:) :: XAVG_RRLIT !
164  REAL, POINTER, DIMENSION(:) :: XAVG_LERCV ! MEB: interception evaporation from overstory canopy vegetation [W/m2]
165  REAL, POINTER, DIMENSION(:) :: XAVG_LE_V_C ! MEB: latent heat flux from vegetation canopy overstory [W/m2]
166  REAL, POINTER, DIMENSION(:) :: XAVG_LE_G_C ! MEB: latent heat flux from understory [W/m2]
167  REAL, POINTER, DIMENSION(:) :: XAVG_LE_C_A ! MEB: latent heat flux from canopy air space to the atmosphere [W/m2]
168  ! NOTE total latent heat flux to the atmosphere also possibly
169  ! includes a contribution from snow covering the canopy
170  REAL, POINTER, DIMENSION(:) :: XAVG_LE_N_C ! MEB: latent heat flux from the snow on the ground [W/m2]
171  ! NOTE total latent heat flux from the snowpack
172  ! possibly includes a contribution from snow covering the canopy
173  REAL, POINTER, DIMENSION(:) :: XAVG_SWNET_V ! MEB: net vegetation canopy shortwave radiation [W/m2]
174  REAL, POINTER, DIMENSION(:) :: XAVG_SWNET_G ! MEB: net ground shortwave radiation [W/m2]
175  REAL, POINTER, DIMENSION(:) :: XAVG_SWNET_N ! MEB: net snow shortwave radiation [W/m2]
176  REAL, POINTER, DIMENSION(:) :: XAVG_SWNET_NS ! MEB: net snow shortwave radiation for *surface* layer
177  ! (i.e. net snow shortwave radiation less absorbed radiation) [W/m2]
178  REAL, POINTER, DIMENSION(:) :: XAVG_LWNET_V ! MEB: net vegetation canopy longwave radiation [W/m2]
179  REAL, POINTER, DIMENSION(:) :: XAVG_LWNET_G ! MEB: net ground longwave radiation [W/m2]
180  REAL, POINTER, DIMENSION(:) :: XAVG_LWNET_N ! MEB: net snow longwave radiation [W/m2]
181  REAL, POINTER, DIMENSION(:) :: XAVG_H_V_C ! MEB: sensible heat flux from vegetation canopy overstory [W/m2]
182  REAL, POINTER, DIMENSION(:) :: XAVG_H_G_C ! MEB: sensible heat flux from understory [W/m2]
183  REAL, POINTER, DIMENSION(:) :: XAVG_H_C_A ! MEB: sensible heat flux from canopy air space to the atmosphere [W/m2]
184  ! NOTE total sensible heat flux to the atmosphere also possibly
185  ! includes a contribution from snow covering the canopy
186  REAL, POINTER, DIMENSION(:) :: XAVG_H_N_C ! MEB: sensible heat flux from the snow on the ground [W/m2]
187  ! NOTE total sensible heat flux from the snowpack
188  ! possibly includes a contribution from snow covering the canopy
189  REAL, POINTER, DIMENSION(:) :: XAVG_SR_GN ! MEB: snow unloading rate from the overstory reservoir [kg/m2/s]
190  REAL, POINTER, DIMENSION(:) :: XAVG_MELTCV ! MEB: snow melt rate from the overstory snow reservoir [kg/m2/s]
191  REAL, POINTER, DIMENSION(:) :: XAVG_FRZCV ! MEB: snow refreeze rate from the overstory snow reservoir [kg/m2/s]
192  REAL, POINTER, DIMENSION(:) :: XAVG_SWDOWN_GN ! MEB: total shortwave radiation transmitted through the canopy
193  ! reaching the snowpack/ground understory [W/m2]
194  REAL, POINTER, DIMENSION(:) :: XAVG_LWDOWN_GN ! MEB: total shortwave radiation transmitted through and emitted by the canopy
195  ! reaching the snowpack/ground understory (explicit part) [W/m2]
196 !
197  REAL, POINTER, DIMENSION(:) :: XRAINFALL ! input rainfall rate for LWATER_BUDGET (kg/m2/s)
198  REAL, POINTER, DIMENSION(:) :: XSNOWFALL ! input snowfall rate for LWATER_BUDGET (kg/m2/s)
199  REAL, POINTER, DIMENSION(:) :: XAVG_DWG ! liquid soil moisture time tendencies (kg/m2/s)
200  REAL, POINTER, DIMENSION(:) :: XAVG_DWGI ! solid soil moisture time tendencies (kg/m2/s)
201  REAL, POINTER, DIMENSION(:) :: XAVG_DWR ! canopy water time tendencies (kg/m2/s)
202  REAL, POINTER, DIMENSION(:) :: XAVG_DSWE ! snow water equivalent time tendencies (kg/m2/s)
203  REAL, POINTER, DIMENSION(:) :: XAVG_WATBUD ! ISBA water budget (kg/m2/s)
204 !
205 !* budget summation variables for each patch
206 !
207  REAL, POINTER, DIMENSION(:,:) :: XRNC ! net radiation at surface (J/m2)
208  REAL, POINTER, DIMENSION(:,:) :: XHC ! sensible heat flux (J/m2)
209  REAL, POINTER, DIMENSION(:,:) :: XLEC ! total latent heat flux (J/m2)
210  REAL, POINTER, DIMENSION(:,:) :: XLEIC ! sublimation latent heat flux (J/m2)
211  REAL, POINTER, DIMENSION(:,:) :: XGFLUXC ! net soil-vegetation flux (J/m2)
212  REAL, POINTER, DIMENSION(:,:) :: XLEGC ! latent heat of evaporation over the ground (J/m2)
213  REAL, POINTER, DIMENSION(:,:) :: XLEGIC ! surface soil ice sublimation (J/m2)
214  REAL, POINTER, DIMENSION(:,:) :: XLEVC ! latent heat of evaporation over vegetation (J/m2)
215  REAL, POINTER, DIMENSION(:,:) :: XLESAC ! latent heat of sublimation over the snow (J/m2)
216  REAL, POINTER, DIMENSION(:,:) :: XLESLC ! latent heat of evaporation over the snow (J/m2)
217  REAL, POINTER, DIMENSION(:,:) :: XLERC ! evaporation from canopy water interception (J/m2)
218  REAL, POINTER, DIMENSION(:,:) :: XLETRC ! evapotranspiration of the vegetation (J/m2)
219  REAL, POINTER, DIMENSION(:,:) :: XEVAPC ! evapotranspiration (kg/m2)
220  REAL, POINTER, DIMENSION(:,:) :: XSUBLC ! sublimation (kg/m2)
221  REAL, POINTER, DIMENSION(:,:) :: XSNDRIFTC ! blowing snow sublimation (ES or Crocus) (kg/m2)
222  REAL, POINTER, DIMENSION(:,:) :: XDRAINC ! soil drainage flux (kg/m2)
223  REAL, POINTER, DIMENSION(:,:) :: XQSBC ! lateral subsurface flux (dif option) (kg/m2)
224  REAL, POINTER, DIMENSION(:,:) :: XRUNOFFC ! sub-grid and supersaturation runoff (kg/m2)
225  REAL, POINTER, DIMENSION(:,:) :: XHORTC ! sub-grid Horton runoff from the SGH scheme (kg/m2)
226  REAL, POINTER, DIMENSION(:,:) :: XDRIPC ! dripping from the vegetation reservoir (kg/m2/s)
227  REAL, POINTER, DIMENSION(:,:) :: XMELTC ! snow melt (kg/m2)
228  REAL, POINTER, DIMENSION(:,:) :: XIFLOODC ! Floodplains infiltration (kg/m2)
229  REAL, POINTER, DIMENSION(:,:) :: XPFLOODC ! Precipitation intercepted by the floodplains (kg/m2)
230  REAL, POINTER, DIMENSION(:,:) :: XLE_FLOODC ! Floodplains evapotration (J/m2)
231  REAL, POINTER, DIMENSION(:,:) :: XLEI_FLOODC ! Floodplains evapotration (J/m2)
232  REAL, POINTER, DIMENSION(:,:) :: XICEFLUXC ! ice calving flux (kg/m2)
233  REAL, POINTER, DIMENSION(:,:) :: XRRVEGC ! precipitation intercepted by the vegetation (kg/m2)
234  REAL, POINTER, DIMENSION(:,:) :: XIRRIG_FLUXC ! irrigation rate (as soil input) (kg/m2)
235 !
236  REAL, POINTER, DIMENSION(:,:) :: XGPPC ! Gross Primary Production (kgCO2/m2)
237  REAL, POINTER, DIMENSION(:,:) :: XRESPC_AUTO ! Autotrophic respiration (kgCO2/m2)
238  REAL, POINTER, DIMENSION(:,:) :: XRESPC_ECO ! Ecosystem respiration (kgCO2/m2)
239 !
240  REAL, POINTER, DIMENSION(:,:) :: XLEVCVC ! MEB: total evapotranspiration from vegetation canopy overstory [J/m2]
241  REAL, POINTER, DIMENSION(:,:) :: XLESCC ! MEB: total snow sublimation from vegetation canopy overstory [J/m2]
242  REAL, POINTER, DIMENSION(:,:) :: XLETRGVC ! MEB: transpiration from understory vegetation [J/m2]
243  REAL, POINTER, DIMENSION(:,:) :: XLETRCVC ! MEB: transpiration from overstory canopy vegetation [J/m2]
244  REAL, POINTER, DIMENSION(:,:) :: XLERGVC ! MEB: interception evaporation from understory vegetation [J/m2]
245  REAL, POINTER, DIMENSION(:,:) :: XLERCVC ! MEB: interception evaporation from overstory canopy vegetation [J/m2]
246  REAL, POINTER, DIMENSION(:,:) :: XLE_V_CC ! MEB: latent heat flux from vegetation canopy overstory [J/m2]
247  REAL, POINTER, DIMENSION(:,:) :: XLE_G_CC ! MEB: latent heat flux from understory [J/m2]
248  REAL, POINTER, DIMENSION(:,:) :: XLE_C_AC ! MEB: latent heat flux from canopy air space to the atmosphere [J/m2]
249  ! NOTE total latent heat flux to the atmosphere also possibly
250  ! includes a contribution from snow covering the canopy
251  REAL, POINTER, DIMENSION(:,:) :: XLE_N_CC ! MEB: latent heat flux from the snow on the ground [J/m2]
252  ! NOTE total latent heat flux from the snowpack
253  ! possibly includes a contribution from snow covering the canopy
254  REAL, POINTER, DIMENSION(:,:) :: XSWNET_VC ! MEB: net vegetation canopy shortwave radiation [W/m2]
255  REAL, POINTER, DIMENSION(:,:) :: XSWNET_GC ! MEB: net ground shortwave radiation [W/m2]
256  REAL, POINTER, DIMENSION(:,:) :: XSWNET_NC ! MEB: net snow shortwave radiation [W/m2]
257  REAL, POINTER, DIMENSION(:,:) :: XSWNET_NSC ! MEB: net snow shortwave radiation for *surface* layer
258  ! (i.e. net snow shortwave radiation less absorbed radiation) [W/m2]
259  REAL, POINTER, DIMENSION(:,:) :: XLWNET_VC ! MEB: net vegetation canopy longwave radiation [W/m2]
260  REAL, POINTER, DIMENSION(:,:) :: XLWNET_GC ! MEB: net ground longwave radiation [W/m2]
261  REAL, POINTER, DIMENSION(:,:) :: XLWNET_NC ! MEB: net snow longwave radiation [W/m2]
262  REAL, POINTER, DIMENSION(:,:) :: XH_V_CC ! MEB: sensible heat flux from vegetation canopy overstory [W/m2]
263  REAL, POINTER, DIMENSION(:,:) :: XH_G_CC ! MEB: sensible heat flux from understory [W/m2]
264  REAL, POINTER, DIMENSION(:,:) :: XH_C_AC ! MEB: sensible heat flux from canopy air space to the atmosphere [W/m2]
265  ! NOTE total sensible heat flux to the atmosphere also possibly
266  ! includes a contribution from snow covering the canopy
267  REAL, POINTER, DIMENSION(:,:) :: XH_N_CC ! MEB: sensible heat flux from the snow on the ground [W/m2]
268  ! NOTE total sensible heat flux from the snowpack
269  ! possibly includes a contribution from snow covering the canopy
270  REAL, POINTER, DIMENSION(:,:) :: XSR_GNC ! MEB: snow unloading rate from the overstory reservoir [kg/m2/s]
271  REAL, POINTER, DIMENSION(:,:) :: XMELTCVC ! MEB: snow melt rate from the overstory snow reservoir [kg/m2/s]
272  REAL, POINTER, DIMENSION(:,:) :: XFRZCVC ! MEB: snow refreeze rate from the overstory snow reservoir [kg/m2/s]
273  REAL, POINTER, DIMENSION(:,:) :: XSWDOWN_GNC ! MEB: total shortwave radiation transmitted through the canopy
274  ! reaching the snowpack/ground understory [W/m2]
275  REAL, POINTER, DIMENSION(:,:) :: XLWDOWN_GNC ! MEB: total shortwave radiation transmitted through and emitted by the canopy
276  ! reaching the snowpack/ground understory (explicit part) [W/m2]
277 !
278  REAL, POINTER, DIMENSION(:,:) :: XDWGC ! liquid soil moisture time tendencies (kg/m2)
279  REAL, POINTER, DIMENSION(:,:) :: XDWGIC ! solid soil moisture time tendencies (kg/m2)
280  REAL, POINTER, DIMENSION(:,:) :: XDWRC ! canopy water time tendencies (kg/m2)
281  REAL, POINTER, DIMENSION(:,:) :: XDSWEC ! snow water equivalent time tendencies (kg/m2)
282  REAL, POINTER, DIMENSION(:,:) :: XWATBUDC ! ISBA water budget (kg/m2)
283 !
284 !* average budget summation variables
285 !
286  REAL, POINTER, DIMENSION(:) :: XAVG_RNC ! net radiation at surface (J/m2)
287  REAL, POINTER, DIMENSION(:) :: XAVG_HC ! sensible heat flux (J/m2)
288  REAL, POINTER, DIMENSION(:) :: XAVG_LEC ! total latent heat flux (J/m2)
289  REAL, POINTER, DIMENSION(:) :: XAVG_LEIC ! sublimation latent heat flux (J/m2)
290  REAL, POINTER, DIMENSION(:) :: XAVG_GFLUXC ! net soil-vegetation flux (J/m2)
291  REAL, POINTER, DIMENSION(:) :: XAVG_LEGC ! latent heat of evaporation over the ground (J/m2)
292  REAL, POINTER, DIMENSION(:) :: XAVG_LEGIC ! surface soil ice sublimation (J/m2)
293  REAL, POINTER, DIMENSION(:) :: XAVG_LEVC ! latent heat of evaporation over vegetation (J/m2)
294  REAL, POINTER, DIMENSION(:) :: XAVG_LESAC ! latent heat of sublimation over the snow (J/m2)
295  REAL, POINTER, DIMENSION(:) :: XAVG_LESLC ! latent heat of evaporation over the snow (J/m2)
296  REAL, POINTER, DIMENSION(:) :: XAVG_LERC ! evaporation from canopy water interception (J/m2)
297  REAL, POINTER, DIMENSION(:) :: XAVG_LETRC ! evapotranspiration of the vegetation (J/m2)
298  REAL, POINTER, DIMENSION(:) :: XAVG_EVAPC ! evapotranspiration (kg/m2)
299  REAL, POINTER, DIMENSION(:) :: XAVG_SUBLC ! sublimation (kg/m2)
300  REAL, POINTER, DIMENSION(:) :: XAVG_SNDRIFTC ! blowing snow sublimation (ES or Crocus) (kg/m2)
301  REAL, POINTER, DIMENSION(:) :: XAVG_DRAINC ! soil drainage flux (kg/m2)
302  REAL, POINTER, DIMENSION(:) :: XAVG_QSBC ! lateral subsurface flux (dif option) (kg/m2)
303  REAL, POINTER, DIMENSION(:) :: XAVG_RUNOFFC ! sub-grid and supersaturation runoff (kg/m2)
304  REAL, POINTER, DIMENSION(:) :: XAVG_HORTC ! sub-grid Horton runoff from the SGH scheme (kg/m2)
305  REAL, POINTER, DIMENSION(:) :: XAVG_DRIPC ! dripping from the vegetation reservoir (kg/m2/s)
306  REAL, POINTER, DIMENSION(:) :: XAVG_MELTC ! snow melt (kg/m2)
307  REAL, POINTER, DIMENSION(:) :: XAVG_IFLOODC ! Floodplains infiltration (kg/m2)
308  REAL, POINTER, DIMENSION(:) :: XAVG_PFLOODC ! Precipitation intercepted by the floodplains (kg/m2)
309  REAL, POINTER, DIMENSION(:) :: XAVG_LE_FLOODC ! Floodplains evapotration (J/m2)
310  REAL, POINTER, DIMENSION(:) :: XAVG_LEI_FLOODC! Floodplains evapotration (J/m2)
311  REAL, POINTER, DIMENSION(:) :: XAVG_ICEFLUXC ! ice calving flux (kg/m2)
312  REAL, POINTER, DIMENSION(:) :: XAVG_RRVEGC ! precipitation intercepted by the vegetation (kg/m2)
313  REAL, POINTER, DIMENSION(:) :: XAVG_IRRIG_FLUXC! irrigation rate (as soil input) (kg/m2)
314 !
315  REAL, POINTER, DIMENSION(:) :: XAVG_GPPC ! Gross Primary Production (kgCO2/m2)
316  REAL, POINTER, DIMENSION(:) :: XAVG_RESPC_AUTO! Autotrophic respiration (kgCO2/m2)
317  REAL, POINTER, DIMENSION(:) :: XAVG_RESPC_ECO ! Ecosystem respiration (kgCO2/m2)
318 !
319  REAL, POINTER, DIMENSION(:) :: XAVG_LEVCVC ! MEB: total evapotranspiration from vegetation canopy overstory [J/m2]
320  REAL, POINTER, DIMENSION(:) :: XAVG_LESCC ! MEB: total snow sublimation from vegetation canopy overstory [J/m2]
321  REAL, POINTER, DIMENSION(:) :: XAVG_LETRGVC ! MEB: transpiration from understory vegetation [J/m2]
322  REAL, POINTER, DIMENSION(:) :: XAVG_LETRCVC ! MEB: transpiration from overstory canopy vegetation [J/m2]
323  REAL, POINTER, DIMENSION(:) :: XAVG_LERGVC ! MEB: interception evaporation from understory vegetation [J/m2]
324  REAL, POINTER, DIMENSION(:) :: XAVG_LERCVC ! MEB: interception evaporation from overstory canopy vegetation [J/m2]
325  REAL, POINTER, DIMENSION(:) :: XAVG_LE_V_CC ! MEB: latent heat flux from vegetation canopy overstory [J/m2]
326  REAL, POINTER, DIMENSION(:) :: XAVG_LE_G_CC ! MEB: latent heat flux from understory [J/m2]
327  REAL, POINTER, DIMENSION(:) :: XAVG_LE_C_AC ! MEB: latent heat flux from canopy air space to the atmosphere [J/m2]
328  ! NOTE total latent heat flux to the atmosphere also possibly
329  ! includes a contribution from snow covering the canopy
330  REAL, POINTER, DIMENSION(:) :: XAVG_LE_N_CC ! MEB: latent heat flux from the snow on the ground [J/m2]
331  ! NOTE total latent heat flux from the snowpack
332  ! possibly includes a contribution from snow covering the canopy
333  REAL, POINTER, DIMENSION(:) :: XAVG_SWNET_VC ! MEB: net vegetation canopy shortwave radiation [W/m2]
334  REAL, POINTER, DIMENSION(:) :: XAVG_SWNET_GC ! MEB: net ground shortwave radiation [W/m2]
335  REAL, POINTER, DIMENSION(:) :: XAVG_SWNET_NC ! MEB: net snow shortwave radiation [W/m2]
336  REAL, POINTER, DIMENSION(:) :: XAVG_SWNET_NSC ! MEB: net snow shortwave radiation for *surface* layer
337  ! (i.e. net snow shortwave radiation less absorbed radiation) [W/m2]
338  REAL, POINTER, DIMENSION(:) :: XAVG_LWNET_VC ! MEB: net vegetation canopy longwave radiation [W/m2]
339  REAL, POINTER, DIMENSION(:) :: XAVG_LWNET_GC ! MEB: net ground longwave radiation [W/m2]
340  REAL, POINTER, DIMENSION(:) :: XAVG_LWNET_NC ! MEB: net snow longwave radiation [W/m2]
341  REAL, POINTER, DIMENSION(:) :: XAVG_H_V_CC ! MEB: sensible heat flux from vegetation canopy overstory [W/m2]
342  REAL, POINTER, DIMENSION(:) :: XAVG_H_G_CC ! MEB: sensible heat flux from understory [W/m2]
343  REAL, POINTER, DIMENSION(:) :: XAVG_H_C_AC ! MEB: sensible heat flux from canopy air space to the atmosphere [W/m2]
344  ! NOTE total sensible heat flux to the atmosphere also possibly
345  ! includes a contribution from snow covering the canopy
346  REAL, POINTER, DIMENSION(:) :: XAVG_H_N_CC ! MEB: sensible heat flux from the snow on the ground [W/m2]
347  ! NOTE total sensible heat flux from the snowpack
348  ! possibly includes a contribution from snow covering the canopy
349  REAL, POINTER, DIMENSION(:) :: XAVG_SR_GNC ! MEB: snow unloading rate from the overstory reservoir [kg/m2/s]
350  REAL, POINTER, DIMENSION(:) :: XAVG_MELTCVC ! MEB: snow melt rate from the overstory snow reservoir [kg/m2/s]
351  REAL, POINTER, DIMENSION(:) :: XAVG_FRZCVC ! MEB: snow refreeze rate from the overstory snow reservoir [kg/m2/s]
352  REAL, POINTER, DIMENSION(:) :: XAVG_SWDOWN_GNC ! MEB: total shortwave radiation transmitted through the canopy
353  ! reaching the snowpack/ground understory [W/m2]
354  REAL, POINTER, DIMENSION(:) :: XAVG_LWDOWN_GNC ! MEB: total shortwave radiation transmitted through and emitted by the canopy
355  ! reaching the snowpack/ground understory (explicit part) [W/m2]
356 !
357  REAL, POINTER, DIMENSION(:) :: XRAINFALLC ! input rainfall rate for LWATER_BUDGET (kg/m2)
358  REAL, POINTER, DIMENSION(:) :: XSNOWFALLC ! input snowfall rate for LWATER_BUDGET (kg/m2)
359  REAL, POINTER, DIMENSION(:) :: XAVG_DWGC ! liquid soil moisture time tendencies (kg/m2)
360  REAL, POINTER, DIMENSION(:) :: XAVG_DWGIC ! solid soil moisture time tendencies (kg/m2)
361  REAL, POINTER, DIMENSION(:) :: XAVG_DWRC ! canopy water time tendencies (kg/m2)
362  REAL, POINTER, DIMENSION(:) :: XAVG_DSWEC ! snow water equivalent time tendencies (kg/m2)
363  REAL, POINTER, DIMENSION(:) :: XAVG_WATBUDC ! ISBA water budget (kg/m2)
364 !
365 !------------------------------------------------------------------------------
366 !
367 
368 END TYPE diag_evap_isba_t
369 
370 
371 
372  CONTAINS
373 !
374 
375 
376 
377 
378 SUBROUTINE diag_evap_isba_init(YDIAG_EVAP_ISBA)
379 TYPE(diag_evap_isba_t), INTENT(INOUT) :: ydiag_evap_isba
380 REAL(KIND=JPRB) :: zhook_handle
381 IF (lhook) CALL dr_hook("MODD_DIAG_EVAP_ISBA_N:DIAG_EVAP_ISBA_INIT",0,zhook_handle)
382  nullify(ydiag_evap_isba%XLEG)
383  nullify(ydiag_evap_isba%XLEGI)
384  nullify(ydiag_evap_isba%XLEV)
385  nullify(ydiag_evap_isba%XLES)
386  nullify(ydiag_evap_isba%XLESL)
387  nullify(ydiag_evap_isba%XLER)
388  nullify(ydiag_evap_isba%XLETR)
389  nullify(ydiag_evap_isba%XEVAP)
390  nullify(ydiag_evap_isba%XSUBL)
391  nullify(ydiag_evap_isba%XSNDRIFT)
392  nullify(ydiag_evap_isba%XDRAIN)
393  nullify(ydiag_evap_isba%XQSB)
394  nullify(ydiag_evap_isba%XRUNOFF)
395  nullify(ydiag_evap_isba%XHORT)
396  nullify(ydiag_evap_isba%XRRVEG)
397  nullify(ydiag_evap_isba%XMELT)
398  nullify(ydiag_evap_isba%XIFLOOD)
399  nullify(ydiag_evap_isba%XPFLOOD)
400  nullify(ydiag_evap_isba%XLE_FLOOD)
401  nullify(ydiag_evap_isba%XLEI_FLOOD)
402 !
403  nullify(ydiag_evap_isba%XLEVCV)
404  nullify(ydiag_evap_isba%XLESC)
405  nullify(ydiag_evap_isba%XLETRGV)
406  nullify(ydiag_evap_isba%XLETRCV)
407  nullify(ydiag_evap_isba%XLERGV)
408  nullify(ydiag_evap_isba%XLELITTER)
409  nullify(ydiag_evap_isba%XLELITTERI)
410  nullify(ydiag_evap_isba%XDRIPLIT)
411  nullify(ydiag_evap_isba%XRRLIT)
412  nullify(ydiag_evap_isba%XLERCV)
413  nullify(ydiag_evap_isba%XLE_V_C)
414  nullify(ydiag_evap_isba%XLE_G_C)
415  nullify(ydiag_evap_isba%XLE_C_A)
416  nullify(ydiag_evap_isba%XLE_N_C)
417 !
418  nullify(ydiag_evap_isba%XSWNET_V)
419  nullify(ydiag_evap_isba%XSWNET_G)
420  nullify(ydiag_evap_isba%XSWNET_N)
421  nullify(ydiag_evap_isba%XSWNET_NS)
422  nullify(ydiag_evap_isba%XLWNET_V)
423  nullify(ydiag_evap_isba%XLWNET_G)
424  nullify(ydiag_evap_isba%XLWNET_N)
425  nullify(ydiag_evap_isba%XSWDOWN_GN)
426  nullify(ydiag_evap_isba%XLWDOWN_GN)
427  nullify(ydiag_evap_isba%XH_V_C)
428  nullify(ydiag_evap_isba%XH_G_C)
429  nullify(ydiag_evap_isba%XH_C_A)
430  nullify(ydiag_evap_isba%XH_N_C)
431  nullify(ydiag_evap_isba%XSR_GN)
432  nullify(ydiag_evap_isba%XMELTCV)
433  nullify(ydiag_evap_isba%XFRZCV)
434 !
435  nullify(ydiag_evap_isba%XDRIP)
436  nullify(ydiag_evap_isba%XIRRIG_FLUX)
437  nullify(ydiag_evap_isba%XGPP)
438  nullify(ydiag_evap_isba%XRESP_AUTO)
439  nullify(ydiag_evap_isba%XRESP_ECO)
440  nullify(ydiag_evap_isba%XDWG)
441  nullify(ydiag_evap_isba%XDWGI)
442  nullify(ydiag_evap_isba%XDWR)
443  nullify(ydiag_evap_isba%XDSWE)
444  nullify(ydiag_evap_isba%XWATBUD)
445  nullify(ydiag_evap_isba%XAVG_LEG)
446  nullify(ydiag_evap_isba%XAVG_LEGI)
447  nullify(ydiag_evap_isba%XAVG_LEV)
448  nullify(ydiag_evap_isba%XAVG_LES)
449  nullify(ydiag_evap_isba%XAVG_LESL)
450  nullify(ydiag_evap_isba%XAVG_LER)
451  nullify(ydiag_evap_isba%XAVG_LETR)
452  nullify(ydiag_evap_isba%XAVG_EVAP)
453  nullify(ydiag_evap_isba%XAVG_SUBL)
454  nullify(ydiag_evap_isba%XAVG_SNDRIFT)
455  nullify(ydiag_evap_isba%XAVG_DRAIN)
456  nullify(ydiag_evap_isba%XAVG_QSB)
457  nullify(ydiag_evap_isba%XAVG_RUNOFF)
458  nullify(ydiag_evap_isba%XAVG_HORT)
459  nullify(ydiag_evap_isba%XAVG_DRIP)
460  nullify(ydiag_evap_isba%XAVG_MELT)
461  nullify(ydiag_evap_isba%XAVG_IFLOOD)
462  nullify(ydiag_evap_isba%XAVG_PFLOOD)
463  nullify(ydiag_evap_isba%XAVG_LE_FLOOD)
464  nullify(ydiag_evap_isba%XAVG_LEI_FLOOD)
465  nullify(ydiag_evap_isba%XAVG_RRVEG)
466 !
467  nullify(ydiag_evap_isba%XAVG_LEVCV)
468  nullify(ydiag_evap_isba%XAVG_LESC)
469  nullify(ydiag_evap_isba%XAVG_LETRGV)
470  nullify(ydiag_evap_isba%XAVG_LETRCV)
471  nullify(ydiag_evap_isba%XAVG_LERGV)
472  nullify(ydiag_evap_isba%XAVG_LELITTER)
473  nullify(ydiag_evap_isba%XAVG_LELITTERI)
474  nullify(ydiag_evap_isba%XAVG_DRIPLIT)
475  nullify(ydiag_evap_isba%XAVG_RRLIT)
476  nullify(ydiag_evap_isba%XAVG_LERCV)
477  nullify(ydiag_evap_isba%XAVG_LE_V_C)
478  nullify(ydiag_evap_isba%XAVG_LE_G_C)
479  nullify(ydiag_evap_isba%XAVG_LE_C_A)
480  nullify(ydiag_evap_isba%XAVG_LE_N_C)
481 !
482  nullify(ydiag_evap_isba%XAVG_SWNET_V)
483  nullify(ydiag_evap_isba%XAVG_SWNET_G)
484  nullify(ydiag_evap_isba%XAVG_SWNET_N)
485  nullify(ydiag_evap_isba%XAVG_SWNET_NS)
486  nullify(ydiag_evap_isba%XAVG_LWNET_V)
487  nullify(ydiag_evap_isba%XAVG_LWNET_G)
488  nullify(ydiag_evap_isba%XAVG_LWNET_N)
489  nullify(ydiag_evap_isba%XAVG_SWDOWN_GN)
490  nullify(ydiag_evap_isba%XAVG_LWDOWN_GN)
491  nullify(ydiag_evap_isba%XAVG_H_V_C)
492  nullify(ydiag_evap_isba%XAVG_H_G_C)
493  nullify(ydiag_evap_isba%XAVG_H_C_A)
494  nullify(ydiag_evap_isba%XAVG_H_N_C)
495  nullify(ydiag_evap_isba%XAVG_SR_GN)
496  nullify(ydiag_evap_isba%XAVG_MELTCV)
497  nullify(ydiag_evap_isba%XAVG_FRZCV)
498 !
499  nullify(ydiag_evap_isba%XAVG_IRRIG_FLUX)
500  nullify(ydiag_evap_isba%XAVG_GPP)
501  nullify(ydiag_evap_isba%XAVG_RESP_AUTO)
502  nullify(ydiag_evap_isba%XAVG_RESP_ECO)
503  nullify(ydiag_evap_isba%XRAINFALL)
504  nullify(ydiag_evap_isba%XSNOWFALL)
505  nullify(ydiag_evap_isba%XAVG_DWG)
506  nullify(ydiag_evap_isba%XAVG_DWGI)
507  nullify(ydiag_evap_isba%XAVG_DWR)
508  nullify(ydiag_evap_isba%XAVG_DSWE)
509  nullify(ydiag_evap_isba%XAVG_WATBUD)
510  nullify(ydiag_evap_isba%XRNC)
511  nullify(ydiag_evap_isba%XHC)
512  nullify(ydiag_evap_isba%XLEC)
513  nullify(ydiag_evap_isba%XLEIC)
514  nullify(ydiag_evap_isba%XGFLUXC)
515  nullify(ydiag_evap_isba%XLEGC)
516  nullify(ydiag_evap_isba%XLEGIC)
517  nullify(ydiag_evap_isba%XLEVC)
518  nullify(ydiag_evap_isba%XLESAC)
519  nullify(ydiag_evap_isba%XLESLC)
520  nullify(ydiag_evap_isba%XLERC)
521  nullify(ydiag_evap_isba%XLETRC)
522  nullify(ydiag_evap_isba%XEVAPC)
523  nullify(ydiag_evap_isba%XSUBLC)
524  nullify(ydiag_evap_isba%XSNDRIFTC)
525  nullify(ydiag_evap_isba%XDRAINC)
526  nullify(ydiag_evap_isba%XQSBC)
527  nullify(ydiag_evap_isba%XRUNOFFC)
528  nullify(ydiag_evap_isba%XHORTC)
529  nullify(ydiag_evap_isba%XDRIPC)
530  nullify(ydiag_evap_isba%XMELTC)
531  nullify(ydiag_evap_isba%XIFLOODC)
532  nullify(ydiag_evap_isba%XPFLOODC)
533  nullify(ydiag_evap_isba%XLE_FLOODC)
534  nullify(ydiag_evap_isba%XLEI_FLOODC)
535  nullify(ydiag_evap_isba%XICEFLUXC)
536  nullify(ydiag_evap_isba%XRRVEGC)
537 !
538  nullify(ydiag_evap_isba%XLEVCVC)
539  nullify(ydiag_evap_isba%XLESCC)
540  nullify(ydiag_evap_isba%XLETRGVC)
541  nullify(ydiag_evap_isba%XLETRCVC)
542  nullify(ydiag_evap_isba%XLERGVC)
543  nullify(ydiag_evap_isba%XLERCVC)
544  nullify(ydiag_evap_isba%XLE_V_CC)
545  nullify(ydiag_evap_isba%XLE_G_CC)
546  nullify(ydiag_evap_isba%XLE_C_AC)
547  nullify(ydiag_evap_isba%XLE_N_CC)
548 !
549  nullify(ydiag_evap_isba%XSWNET_VC)
550  nullify(ydiag_evap_isba%XSWNET_GC)
551  nullify(ydiag_evap_isba%XSWNET_NC)
552  nullify(ydiag_evap_isba%XSWNET_NSC)
553  nullify(ydiag_evap_isba%XLWNET_VC)
554  nullify(ydiag_evap_isba%XLWNET_GC)
555  nullify(ydiag_evap_isba%XLWNET_NC)
556  nullify(ydiag_evap_isba%XSWDOWN_GNC)
557  nullify(ydiag_evap_isba%XLWDOWN_GNC)
558  nullify(ydiag_evap_isba%XH_V_CC)
559  nullify(ydiag_evap_isba%XH_G_CC)
560  nullify(ydiag_evap_isba%XH_C_AC)
561  nullify(ydiag_evap_isba%XH_N_CC)
562  nullify(ydiag_evap_isba%XSR_GNC)
563  nullify(ydiag_evap_isba%XMELTCVC)
564  nullify(ydiag_evap_isba%XFRZCVC)
565 !
566  nullify(ydiag_evap_isba%XIRRIG_FLUXC)
567  nullify(ydiag_evap_isba%XGPPC)
568  nullify(ydiag_evap_isba%XRESPC_AUTO)
569  nullify(ydiag_evap_isba%XRESPC_ECO)
570  nullify(ydiag_evap_isba%XDWGC)
571  nullify(ydiag_evap_isba%XDWGIC)
572  nullify(ydiag_evap_isba%XDWRC)
573  nullify(ydiag_evap_isba%XDSWEC)
574  nullify(ydiag_evap_isba%XWATBUDC)
575  nullify(ydiag_evap_isba%XAVG_RNC)
576  nullify(ydiag_evap_isba%XAVG_HC)
577  nullify(ydiag_evap_isba%XAVG_LEC)
578  nullify(ydiag_evap_isba%XAVG_LEIC)
579  nullify(ydiag_evap_isba%XAVG_GFLUXC)
580  nullify(ydiag_evap_isba%XAVG_LEGC)
581  nullify(ydiag_evap_isba%XAVG_LEGIC)
582  nullify(ydiag_evap_isba%XAVG_LEVC)
583  nullify(ydiag_evap_isba%XAVG_LESAC)
584  nullify(ydiag_evap_isba%XAVG_LESLC)
585  nullify(ydiag_evap_isba%XAVG_LERC)
586  nullify(ydiag_evap_isba%XAVG_LETRC)
587  nullify(ydiag_evap_isba%XAVG_EVAPC)
588  nullify(ydiag_evap_isba%XAVG_SUBLC)
589  nullify(ydiag_evap_isba%XAVG_SNDRIFTC)
590  nullify(ydiag_evap_isba%XAVG_DRAINC)
591  nullify(ydiag_evap_isba%XAVG_QSBC)
592  nullify(ydiag_evap_isba%XAVG_RUNOFFC)
593  nullify(ydiag_evap_isba%XAVG_HORTC)
594  nullify(ydiag_evap_isba%XAVG_DRIPC)
595  nullify(ydiag_evap_isba%XAVG_MELTC)
596  nullify(ydiag_evap_isba%XAVG_IFLOODC)
597  nullify(ydiag_evap_isba%XAVG_PFLOODC)
598  nullify(ydiag_evap_isba%XAVG_LE_FLOODC)
599  nullify(ydiag_evap_isba%XAVG_LEI_FLOODC)
600  nullify(ydiag_evap_isba%XAVG_ICEFLUXC)
601  nullify(ydiag_evap_isba%XAVG_RRVEGC)
602 !
603  nullify(ydiag_evap_isba%XAVG_LEVCVC)
604  nullify(ydiag_evap_isba%XAVG_LESCC)
605  nullify(ydiag_evap_isba%XAVG_LETRGVC)
606  nullify(ydiag_evap_isba%XAVG_LETRCVC)
607  nullify(ydiag_evap_isba%XAVG_LERGVC)
608  nullify(ydiag_evap_isba%XAVG_LERCVC)
609  nullify(ydiag_evap_isba%XAVG_LE_V_CC)
610  nullify(ydiag_evap_isba%XAVG_LE_G_CC)
611  nullify(ydiag_evap_isba%XAVG_LE_C_AC)
612  nullify(ydiag_evap_isba%XAVG_LE_N_CC)
613 !
614  nullify(ydiag_evap_isba%XAVG_SWNET_VC)
615  nullify(ydiag_evap_isba%XAVG_SWNET_GC)
616  nullify(ydiag_evap_isba%XAVG_SWNET_NC)
617  nullify(ydiag_evap_isba%XAVG_SWNET_NSC)
618  nullify(ydiag_evap_isba%XAVG_LWNET_VC)
619  nullify(ydiag_evap_isba%XAVG_LWNET_GC)
620  nullify(ydiag_evap_isba%XAVG_LWNET_NC)
621  nullify(ydiag_evap_isba%XAVG_SWDOWN_GNC)
622  nullify(ydiag_evap_isba%XAVG_LWDOWN_GNC)
623  nullify(ydiag_evap_isba%XAVG_H_V_CC)
624  nullify(ydiag_evap_isba%XAVG_H_G_CC)
625  nullify(ydiag_evap_isba%XAVG_H_C_AC)
626  nullify(ydiag_evap_isba%XAVG_H_N_CC)
627  nullify(ydiag_evap_isba%XAVG_SR_GNC)
628  nullify(ydiag_evap_isba%XAVG_MELTCVC)
629  nullify(ydiag_evap_isba%XAVG_FRZCVC)
630 !
631  nullify(ydiag_evap_isba%XAVG_IRRIG_FLUXC)
632  nullify(ydiag_evap_isba%XAVG_GPPC)
633  nullify(ydiag_evap_isba%XAVG_RESPC_AUTO)
634  nullify(ydiag_evap_isba%XAVG_RESPC_ECO)
635  nullify(ydiag_evap_isba%XRAINFALLC)
636  nullify(ydiag_evap_isba%XSNOWFALLC)
637  nullify(ydiag_evap_isba%XAVG_DWGC)
638  nullify(ydiag_evap_isba%XAVG_DWGIC)
639  nullify(ydiag_evap_isba%XAVG_DWRC)
640  nullify(ydiag_evap_isba%XAVG_DSWEC)
641  nullify(ydiag_evap_isba%XAVG_WATBUDC)
642 ydiag_evap_isba%LSURF_EVAP_BUDGET=.false.
643 ydiag_evap_isba%LSURF_BUDGETC=.false.
644 ydiag_evap_isba%LRESET_BUDGETC=.false.
645 ydiag_evap_isba%LWATER_BUDGET=.false.
646 IF (lhook) CALL dr_hook("MODD_DIAG_EVAP_ISBA_N:DIAG_EVAP_ISBA_INIT",1,zhook_handle)
647 END SUBROUTINE diag_evap_isba_init
648 
649 
650 END MODULE modd_diag_evap_isba_n
subroutine diag_evap_isba_init(YDIAG_EVAP_ISBA)