SURFEX  V8_0
Surfex V8_0 release
 All Classes Files Functions Variables
modd_pack_diag_isba.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_PACK_DIAG_ISBA - declaration of packed diagnostics for ISBA scheme
10 !!
11 !! PURPOSE
12 !! -------
13 !
14 !!
15 !!** IMPLICIT ARGUMENTS
16 !! ------------------
17 !! None
18 !!
19 !! REFERENCE
20 !! ---------
21 !!
22 !! AUTHOR
23 !! ------
24 !! V. Masson *Meteo France*
25 !!
26 !! MODIFICATIONS
27 !! -------------
28 !! Original 01/2004
29 !! Modified 10/2004 by P. Le Moigne: add Halstead coefficient
30 !! Modified 11/2009 by S. Senesi: add precipitation intercepted by the vegetation (XP_RRVEG)
31 !! Modified 04-09 by A.L. Gibelin : Add carbon diagnostics
32 !! Modified 10-14 by P. Samuelsson: MEB
33 !
34 !* 0. DECLARATIONS
35 ! ------------
36 !
37 USE yomhook ,ONLY : lhook, dr_hook
38 USE parkind1 ,ONLY : jprb
39 !
40 IMPLICIT NONE
41 !
42 !------------------------------------------------------------------------------
43 !
45 
46  INTEGER :: NSIZE_SIMPLE
47  INTEGER :: NSIZE_GROUND
48  INTEGER :: NSIZE_SNOW
49  INTEGER :: NSIZE_KSW
50  INTEGER :: NSIZE_ABC
51  INTEGER :: NSIZE_0
52  INTEGER :: NSIZE_00
53  REAL, POINTER, DIMENSION(:,:) :: XBLOCK_SIMPLE
54  REAL, POINTER, DIMENSION(:,:,:) :: XBLOCK_GROUND
55  REAL, POINTER, DIMENSION(:,:,:) :: XBLOCK_SNOW
56  REAL, POINTER, DIMENSION(:,:,:) :: XBLOCK_KSW
57  REAL, POINTER, DIMENSION(:,:,:) :: XBLOCK_ABC
58  REAL, POINTER, DIMENSION(:,:) :: XBLOCK_0
59  REAL, POINTER, DIMENSION(:,:,:) :: XBLOCK_00
60 !
61  REAL, POINTER, DIMENSION(:) :: XP_RNSNOW ! net radiative flux from snow (ISBA-ES:3-L) (W/m2)
62  REAL, POINTER, DIMENSION(:) :: XP_HSNOW ! sensible heat flux from snow (ISBA-ES:3-L) (W/m2)
63  REAL, POINTER, DIMENSION(:) :: XP_HPSNOW ! heat release from rainfall (ISBA-ES:3-L) (W/m2)
64  REAL, POINTER, DIMENSION(:) :: XP_GFLUXSNOW ! net surface energy flux into snowpack
65 ! ! (ISBA-ES:3-L) (W/m2)
66  REAL, POINTER, DIMENSION(:) :: XP_USTARSNOW ! friction velocity over snow
67 ! ! (ISBA-ES:3-L) (m/s)
68  REAL, POINTER, DIMENSION(:) :: XP_GRNDFLUX ! soil/snow interface heat flux (ISBA-ES:3-L) (W/m2)
69  REAL, POINTER, DIMENSION(:) :: XP_SRSFC ! snowfall over snowpack (ISBA-ES:3-L) (kg/m2/s)
70  REAL, POINTER, DIMENSION(:) :: XP_RRSFC ! rainfall over snowpack (ISBA-ES:3-L) (kg/m2/s)
71  REAL, POINTER, DIMENSION(:) :: XP_LESL ! snowpack evaporation (ISBA-ES:3-L) (W/m2)
72  REAL, POINTER, DIMENSION(:) :: XP_CDSNOW ! snow drag coefficient (ISBA-ES:3-L) (-)
73  REAL, POINTER, DIMENSION(:) :: XP_CHSNOW ! heat turbulent transfer coefficient
74 ! ! (ISBA-ES:3-L) (-)
75  REAL, POINTER, DIMENSION(:,:)::XP_SNOWTEMP ! snow temperature profile (ISBA-ES:3-L) (K)
76  REAL, POINTER, DIMENSION(:,:)::XP_SNOWLIQ ! snow liquid water profile (ISBA-ES:3-L) (m)
77  REAL, POINTER, DIMENSION(:,:)::XP_SNOWDZ ! snow layer thicknesses (m)
78  REAL, POINTER, DIMENSION(:) :: XP_SNOWHMASS ! heat content change due to mass
79 ! ! changes in snowpack: for budget
80 ! ! calculations only. (ISBA-ES:3-L) (J/m2)
81  REAL, POINTER, DIMENSION(:) :: XP_RN_ISBA ! net radiative flux from snow-free surface
82 ! ! (ISBA-ES:3-L) (W/m2)
83  REAL, POINTER, DIMENSION(:) :: XP_H_ISBA ! sensible heat flux from snow-free surface
84 ! ! (ISBA-ES:3-L) (W/m2)
85  REAL, POINTER, DIMENSION(:) :: XP_LEG_ISBA ! baresoil evaporation from snow-free surface
86 ! ! (ISBA-ES:3-L) (W/m2)
87  REAL, POINTER, DIMENSION(:) :: XP_LEGI_ISBA ! baresoil sublimation from snow-free surface
88 ! ! (ISBA-ES:3-L) (W/m2)
89  REAL, POINTER, DIMENSION(:) :: XP_LEV_ISBA ! total evapotranspiration from vegetation over
90 ! ! snow-free surface (ISBA-ES:3-L) (W/m2)
91  REAL, POINTER, DIMENSION(:) :: XP_LETR_ISBA ! transpiration from snow-free surface
92 ! ! (ISBA-ES:3-L) (W/m2)
93  REAL, POINTER, DIMENSION(:) :: XP_USTAR_ISBA! friction velocity from snow-free
94 ! ! surface (ISBA-ES:3-L) (m/s)
95  REAL, POINTER, DIMENSION(:) :: XP_LER_ISBA ! evaporation from canopy water interception
96 ! ! store over snow-free surface (ISBA-ES:3-L) (W/m2)
97  REAL, POINTER, DIMENSION(:) :: XP_LE_ISBA ! total latent heat flux from snow-free surface
98  REAL, POINTER, DIMENSION(:) :: XP_LEI_ISBA ! sublimation latent heat flux from snow-free surface
99 ! ! (ISBA-ES:3-L) (W/m2)
100  REAL, POINTER, DIMENSION(:) :: XP_GFLUX_ISBA! net energy flux into the snow-free surface
101 ! ! (ISBA-ES:3-L) (W/m2)
102  REAL, POINTER, DIMENSION(:) :: XP_MELTADV ! advective energy from snow melt water
103 ! ! (ISBA-ES:3-L) (W/m2)
104  REAL, POINTER, DIMENSION(:) :: XP_CH ! thermal diffusion coefficient (W/s)
105  REAL, POINTER, DIMENSION(:) :: XP_CE ! transfer coefficient for vapor (W/s/K)
106  REAL, POINTER, DIMENSION(:) :: XP_CD ! drag coefficient (-)
107  REAL, POINTER, DIMENSION(:) :: XP_CDN ! neutral drag coefficient (-)
108  REAL, POINTER, DIMENSION(:) :: XP_RI ! Bulk-Richardson number (-)
109  REAL, POINTER, DIMENSION(:) :: XP_HU ! area averaged surface humidity coefficient (-)
110  REAL, POINTER, DIMENSION(:) :: XP_HUG ! baresoil surface humidity coefficient (-)
111  REAL, POINTER, DIMENSION(:) :: XP_HV ! Halstead coefficient (-)
112 !
113  REAL, POINTER, DIMENSION(:) :: XP_ALBT ! Total Albedo (-)
114 !
115  REAL, POINTER, DIMENSION(:) :: XP_RN ! net radiation at surface (W/m2)
116  REAL, POINTER, DIMENSION(:) :: XP_H ! sensible heat flux (W/m2)
117  REAL, POINTER, DIMENSION(:) :: XP_LEG ! baresoil evaporation (W/m2)
118  REAL, POINTER, DIMENSION(:) :: XP_LEGI ! baresoil sublimation (W/m2)
119  REAL, POINTER, DIMENSION(:) :: XP_LEV ! total evapotranspiration from vegetation (W/m2)
120  REAL, POINTER, DIMENSION(:) :: XP_LES ! snow sublimation (W/m2)
121  REAL, POINTER, DIMENSION(:) :: XP_LER ! evaporation from canopy water interception
122 ! ! store (W/m2)
123  REAL, POINTER, DIMENSION(:) :: XP_LETR ! transpiration (W/m2)
124  REAL, POINTER, DIMENSION(:) :: XP_EVAP ! evapotranspiration (kg/m2/s)
125  REAL, POINTER, DIMENSION(:) :: XP_SUBL ! sublimation (kg/m2/s)
126  REAL, POINTER, DIMENSION(:) :: XP_SNDRIFT ! blowing snow sublimation (kg/m2/s)
127  REAL, POINTER, DIMENSION(:) :: XP_LEI ! sublimation latent heat flux (W/m2)
128  REAL, POINTER, DIMENSION(:) :: XP_GFLUX ! net soil-vegetation flux (W/m2)
129  REAL, POINTER, DIMENSION(:) :: XP_RESTORE ! surface energy budget restore term (W/m2)
130  REAL, POINTER, DIMENSION(:) :: XP_DRAIN ! soil drainage flux (kg/m2/s)
131  REAL, POINTER, DIMENSION(:) :: XP_QSB ! lateral subsurface flux (kg/m2/s)
132  REAL, POINTER, DIMENSION(:) :: XP_RUNOFF ! sub-grid and supersaturation runoff (kg/m2/s)
133  REAL, POINTER, DIMENSION(:) :: XP_MELT ! melting rate of the snow (kg/m2/s)
134  REAL, POINTER, DIMENSION(:) :: XP_SNOWFREE_ALB ! snow-free global albedo (-)
135  REAL, POINTER, DIMENSION(:) :: XP_SNOWFREE_ALB_VEG ! snow-free global albedo of vegetation
136  REAL, POINTER, DIMENSION(:) :: XP_SNOWFREE_ALB_SOIL! snow-free soil albedo
137  REAL, POINTER, DIMENSION(:) :: XP_Z0_WITH_SNOW ! total roughness length (including snow) (m)
138  REAL, POINTER, DIMENSION(:) :: XP_Z0H_WITH_SNOW! roughness length for heat (including snow) (m)
139  REAL, POINTER, DIMENSION(:) :: XP_Z0EFF ! effective roughness length (with relief added)(m)
140 !
141  REAL, POINTER, DIMENSION(:,:)::XP_IACAN ! PAR in the canopy at different gauss level (micmolphot/m2/s)
142 !
143  REAL, POINTER, DIMENSION(:) :: XP_CG ! heat capacity of the ground
144  REAL, POINTER, DIMENSION(:) :: XP_C1 ! coefficients for the moisure
145  REAL, POINTER, DIMENSION(:) :: XP_C2 ! equation.
146  REAL, POINTER, DIMENSION(:) :: XP_WGEQ ! equilibrium volumetric water
147 ! ! content
148  REAL, POINTER, DIMENSION(:) :: XP_CT ! area-averaged heat capacity
149  REAL, POINTER, DIMENSION(:) :: XP_RS ! stomatal resistance (s/m)
150 !
151 !------------------------------------------------------------------------------
152 !
153  REAL, POINTER, DIMENSION(:) :: XP_TS ! Surface temperature (K)
154  REAL, POINTER, DIMENSION(:) :: XP_TSRAD ! Radiative surface temperature (K)
155  REAL, POINTER, DIMENSION(:) :: XP_T2M ! Air temperature at 2 meters (K)
156  REAL, POINTER, DIMENSION(:) :: XP_Q2M ! Air spec. humidity at 2 meters (kg/kg)
157  REAL, POINTER, DIMENSION(:) :: XP_HU2M ! Air rela. humidity at 2 meters (-)
158  REAL, POINTER, DIMENSION(:) :: XP_ZON10M ! zonal Wind at 10 meters (m/s)
159  REAL, POINTER, DIMENSION(:) :: XP_MER10M ! meridian Wind at 10 meters (m/s)
160 !
161 !------------------------------------------------------------------------------
162 !
163  REAL, POINTER, DIMENSION(:) :: XP_QS ! humidity at surface (Kg/kg)
164  REAL, POINTER, DIMENSION(:,:) :: XP_SWI ! soil wetness index profile (-)
165  REAL, POINTER, DIMENSION(:,:) :: XP_TSWI ! total soil wetness index profile (-)
166 !
167 !------------------------------------------------------------------------------
168 !
169  REAL, POINTER, DIMENSION(:) :: XP_TWSNOW ! total snow reservoir (kg/m2)
170  REAL, POINTER, DIMENSION(:) :: XP_TDSNOW ! total snow height (m)
171 !
172 !------------------------------------------------------------------------------
173 !
174  REAL, POINTER, DIMENSION(:) :: XP_SWD ! downward short wave radiation (W/m2)
175  REAL, POINTER, DIMENSION(:) :: XP_SWU ! upward short wave radiation (W/m2)
176  REAL, POINTER, DIMENSION(:,:) :: XP_SWBD ! downward short wave radiation by spectral band (W/m2)
177  REAL, POINTER, DIMENSION(:,:) :: XP_SWBU ! upward short wave radiation by spectral band (W/m2)
178  REAL, POINTER, DIMENSION(:) :: XP_LWD ! downward long wave radiation (W/m2)
179  REAL, POINTER, DIMENSION(:) :: XP_LWU ! upward long wave radiation (W/m2)
180  REAL, POINTER, DIMENSION(:) :: XP_FMU ! horizontal momentum flux zonal (m2/s2)
181  REAL, POINTER, DIMENSION(:) :: XP_FMV ! horizontal momentum flux meridian (m2/s2)
182 !
183 !------------------------------------------------------------------------------
184 !
185  REAL, POINTER, DIMENSION(:) :: XP_HORT ! sub-grid Horton runoff from the SGH scheme (kg/m2/s)
186  REAL, POINTER, DIMENSION(:) :: XP_DRIP ! dripping from the vegetation reservoir (kg/m2/s)
187  REAL, POINTER, DIMENSION(:) :: XP_IFLOOD ! flood infiltration (kg/m2/s)
188  REAL, POINTER, DIMENSION(:) :: XP_PFLOOD ! precipitation intercepted by the floodplains (kg/m2/s)
189  REAL, POINTER, DIMENSION(:) :: XP_LE_FLOOD ! flood evaporation (W/m2)
190  REAL, POINTER, DIMENSION(:) :: XP_LEI_FLOOD ! frozen flood evaporation (W/m2)
191  REAL, POINTER, DIMENSION(:) :: XP_ICEFLUX
192  REAL, POINTER, DIMENSION(:) :: XP_RRVEG ! precipitation intercepted by the vegetation (kg/m2/s)
193  REAL, POINTER, DIMENSION(:) :: XP_IRRIG_FLUX! irrigation rate (kg/m2/s)
194 !
195 !------------------------------------------------------------------------------
196 !
197  REAL, POINTER, DIMENSION(:) :: XP_GPP ! Gross primary production (kgCO2/m2/s)
198  REAL, POINTER, DIMENSION(:) :: XP_RESP_AUTO ! Autotrophic respiration (kgCO2/m2/s)
199  REAL, POINTER, DIMENSION(:) :: XP_RESP_ECO ! Ecosystem respiration (kgCO2/m2/s)
200  REAL, POINTER, DIMENSION(:) :: XP_FAPAR ! Fapar of vegetation
201  REAL, POINTER, DIMENSION(:) :: XP_FAPIR ! Fapir of vegetation
202  REAL, POINTER, DIMENSION(:) :: XP_FAPAR_BS ! Fapar of bare soil
203  REAL, POINTER, DIMENSION(:) :: XP_FAPIR_BS ! Fapir of bare soil
204 !
205 !------------------------------------------------------------------------------
206 !
207 !* diagnostic variables for multi-energy balance (MEB)
208 ! ---------------------------------------------------
209 !
210  REAL, POINTER, DIMENSION(:) :: XP_SWUP ! MEB: net *total* (surface) upwelling shortwave radiation to atmosphere [W/m2]
211  REAL, POINTER, DIMENSION(:) :: XP_SWNET_V ! MEB: net vegetation canopy shortwave radiation [W/m2]
212  REAL, POINTER, DIMENSION(:) :: XP_SWNET_G ! MEB: net ground shortwave radiation [W/m2]
213  REAL, POINTER, DIMENSION(:) :: XP_SWNET_N ! MEB: net snow shortwave radiation [W/m2]
214  REAL, POINTER, DIMENSION(:) :: XP_SWNET_NS ! MEB: net snow shortwave radiation for *surface* layer
215  ! (i.e. net snow shortwave radiation less absorbed radiation) [W/m2]
216  REAL, POINTER, DIMENSION(:) :: XP_LWUP ! MEB: net *total* (surface) upwelling longwave radiation to atmosphere [W/m2]
217  REAL, POINTER, DIMENSION(:) :: XP_LWNET_V ! MEB: net vegetation canopy longwave radiation [W/m2]
218  REAL, POINTER, DIMENSION(:) :: XP_LWNET_G ! MEB: net ground longwave radiation [W/m2]
219  REAL, POINTER, DIMENSION(:) :: XP_LWNET_N ! MEB: net snow longwave radiation [W/m2]
220  REAL, POINTER, DIMENSION(:) :: XP_LEVCV ! MEB: total evapotranspiration from vegetation canopy overstory [W/m2]
221  REAL, POINTER, DIMENSION(:) :: XP_LESC ! MEB: total snow sublimation from vegetation canopy overstory [W/m2]
222  REAL, POINTER, DIMENSION(:) :: XP_H_V_C ! MEB: sensible heat flux from vegetation canopy overstory [W/m2]
223  REAL, POINTER, DIMENSION(:) :: XP_H_G_C ! MEB: sensible heat flux from understory [W/m2]
224  REAL, POINTER, DIMENSION(:) :: XP_LETRGV ! MEB: transpiration from understory vegetation [W/m2]
225  REAL, POINTER, DIMENSION(:) :: XP_LETRCV ! MEB: transpiration from overstory canopy vegetation [W/m2]
226  REAL, POINTER, DIMENSION(:) :: XP_LERGV ! MEB: interception evaporation from understory vegetation [W/m2
227  REAL, POINTER, DIMENSION(:) :: XP_LELITTER ! MEB: water evaporation from litter[W/m2]
228  REAL, POINTER, DIMENSION(:) :: XP_LELITTERI ! MEB: ice sublimation from litter [W/m2]
229  REAL, POINTER, DIMENSION(:) :: XP_DRIPLIT !
230  REAL, POINTER, DIMENSION(:) :: XP_RRLIT !
231  REAL, POINTER, DIMENSION(:) :: XP_LERCV ! MEB: interception evaporation from overstory canopy vegetation [W/m2]
232  REAL, POINTER, DIMENSION(:) :: XP_H_C_A ! MEB: sensible heat flux from canopy air space to the atmosphere [W/m2]
233  REAL, POINTER, DIMENSION(:) :: XP_H_N_C ! MEB: sensible heat flux from the snow on the ground [W/m2]
234  ! NOTE total sensible heat flux to the atmosphere also possibly
235  ! includes a contribution from snow covering the canopy
236  REAL, POINTER, DIMENSION(:) :: XP_LE_V_C ! MEB: latent heat flux from vegetation canopy overstory [W/m2]
237  REAL, POINTER, DIMENSION(:) :: XP_LE_G_C ! MEB: latent heat flux from understory [W/m2]
238  REAL, POINTER, DIMENSION(:) :: XP_LE_C_A ! MEB: latent heat flux from canopy air space to the atmosphere [W/m2]
239  ! NOTE total latent heat flux to the atmosphere also possibly
240  ! includes a contribution from snow covering the canopy
241  REAL, POINTER, DIMENSION(:) :: XP_LE_N_C ! MEB: latent heat flux from the snow on the ground [W/m2]
242  ! NOTE total latent heat flux from the snowpack
243  ! possibly includes a contribution from snow covering the canopy
244  REAL, POINTER, DIMENSION(:) :: XP_EVAP_N_C ! MEB: Total evap from snow on the ground to canopy air space [kg/m2/s]
245  REAL, POINTER, DIMENSION(:) :: XP_EVAP_G_C ! MEB: Total evap from ground to canopy air space [kg/m2/s]
246  REAL, POINTER, DIMENSION(:) :: XP_SR_GN ! MEB: total snow reacing the ground snow [kg/m2/s]
247  REAL, POINTER, DIMENSION(:) :: XP_MELTCV ! MEB: snow melt rate from the overstory snow reservoir [kg/m2/s]
248  REAL, POINTER, DIMENSION(:) :: XP_FRZCV ! MEB: snow refreeze rate from the overstory snow reservoir [kg/m2/s]
249  REAL, POINTER, DIMENSION(:) :: XP_SWDOWN_GN ! MEB: total shortwave radiation transmitted through the canopy
250  ! reaching the snowpack/ground understory [W/m2]
251  REAL, POINTER, DIMENSION(:) :: XP_LWDOWN_GN ! MEB: total shortwave radiation transmitted through and emitted by the canopy
252 ! ! reaching the snowpack/ground understory (explicit part) [W/m2]
253 !
254 !------------------------------------------------------------------------------
255 !
256  REAL, POINTER, DIMENSION(:) :: XP_DWG ! liquid soil moisture time tendencies (kg/m2/s)
257  REAL, POINTER, DIMENSION(:) :: XP_DWGI ! solid soil moisture time tendencies (kg/m2/s)
258  REAL, POINTER, DIMENSION(:) :: XP_DWR ! canopy water time tendencies (kg/m2/s)
259  REAL, POINTER, DIMENSION(:) :: XP_DSWE ! snow water equivalent time tendencies (kg/m2/s)
260  REAL, POINTER, DIMENSION(:) :: XP_WATBUD ! ISBA water budget (kg/m2/s)
261 
262 END TYPE pack_diag_isba_t
263 !
264 !-------------------------------------------------------------------------------
265 !
266 
267 
268  CONTAINS
269 
270 !
271 !
272 
273 
274 
275 
276 SUBROUTINE pack_diag_isba_init(YPACK_DIAG_ISBA)
277 TYPE(pack_diag_isba_t), INTENT(INOUT) :: ypack_diag_isba
278 REAL(KIND=JPRB) :: zhook_handle
279 IF (lhook) CALL dr_hook("MODD_PACK_DIAG_ISBA_N:PACK_DIAG_ISBA_INIT",0,zhook_handle)
280  nullify(ypack_diag_isba%XBLOCK_SIMPLE)
281  nullify(ypack_diag_isba%XBLOCK_GROUND)
282  nullify(ypack_diag_isba%XBLOCK_SNOW)
283  nullify(ypack_diag_isba%XBLOCK_KSW)
284  nullify(ypack_diag_isba%XBLOCK_ABC)
285  nullify(ypack_diag_isba%XBLOCK_0)
286  nullify(ypack_diag_isba%XBLOCK_00)
287  nullify(ypack_diag_isba%XP_RNSNOW)
288  nullify(ypack_diag_isba%XP_HSNOW)
289  nullify(ypack_diag_isba%XP_HPSNOW)
290  nullify(ypack_diag_isba%XP_GFLUXSNOW)
291  nullify(ypack_diag_isba%XP_USTARSNOW)
292  nullify(ypack_diag_isba%XP_GRNDFLUX)
293  nullify(ypack_diag_isba%XP_SRSFC)
294  nullify(ypack_diag_isba%XP_RRSFC)
295  nullify(ypack_diag_isba%XP_LESL)
296  nullify(ypack_diag_isba%XP_CDSNOW)
297  nullify(ypack_diag_isba%XP_CHSNOW)
298  nullify(ypack_diag_isba%XP_SNOWTEMP)
299  nullify(ypack_diag_isba%XP_SNOWLIQ)
300  nullify(ypack_diag_isba%XP_SNOWDZ)
301  nullify(ypack_diag_isba%XP_SNOWHMASS)
302  nullify(ypack_diag_isba%XP_RN_ISBA)
303  nullify(ypack_diag_isba%XP_H_ISBA)
304  nullify(ypack_diag_isba%XP_LEG_ISBA)
305  nullify(ypack_diag_isba%XP_LEGI_ISBA)
306  nullify(ypack_diag_isba%XP_LEV_ISBA)
307  nullify(ypack_diag_isba%XP_LETR_ISBA)
308  nullify(ypack_diag_isba%XP_USTAR_ISBA)
309  nullify(ypack_diag_isba%XP_LER_ISBA)
310  nullify(ypack_diag_isba%XP_LE_ISBA)
311  nullify(ypack_diag_isba%XP_LEI_ISBA)
312  nullify(ypack_diag_isba%XP_GFLUX_ISBA)
313  nullify(ypack_diag_isba%XP_MELTADV)
314  nullify(ypack_diag_isba%XP_CH)
315  nullify(ypack_diag_isba%XP_CE)
316  nullify(ypack_diag_isba%XP_CD)
317  nullify(ypack_diag_isba%XP_CDN)
318  nullify(ypack_diag_isba%XP_RI)
319  nullify(ypack_diag_isba%XP_HU)
320  nullify(ypack_diag_isba%XP_HUG)
321  nullify(ypack_diag_isba%XP_HV)
322  nullify(ypack_diag_isba%XP_ALBT)
323  nullify(ypack_diag_isba%XP_RN)
324  nullify(ypack_diag_isba%XP_H)
325  nullify(ypack_diag_isba%XP_LEG)
326  nullify(ypack_diag_isba%XP_LEGI)
327  nullify(ypack_diag_isba%XP_LEV)
328  nullify(ypack_diag_isba%XP_LES)
329  nullify(ypack_diag_isba%XP_LER)
330  nullify(ypack_diag_isba%XP_LETR)
331  nullify(ypack_diag_isba%XP_EVAP)
332  nullify(ypack_diag_isba%XP_SUBL)
333  nullify(ypack_diag_isba%XP_SNDRIFT)
334  nullify(ypack_diag_isba%XP_LEI)
335  nullify(ypack_diag_isba%XP_GFLUX)
336  nullify(ypack_diag_isba%XP_RESTORE)
337  nullify(ypack_diag_isba%XP_DRAIN)
338  nullify(ypack_diag_isba%XP_QSB)
339  nullify(ypack_diag_isba%XP_RUNOFF)
340  nullify(ypack_diag_isba%XP_MELT)
341  nullify(ypack_diag_isba%XP_SNOWFREE_ALB)
342  nullify(ypack_diag_isba%XP_SNOWFREE_ALB_VEG)
343  nullify(ypack_diag_isba%XP_SNOWFREE_ALB_SOIL)
344  nullify(ypack_diag_isba%XP_Z0_WITH_SNOW)
345  nullify(ypack_diag_isba%XP_Z0H_WITH_SNOW)
346  nullify(ypack_diag_isba%XP_Z0EFF)
347  nullify(ypack_diag_isba%XP_IACAN)
348  nullify(ypack_diag_isba%XP_CG)
349  nullify(ypack_diag_isba%XP_C1)
350  nullify(ypack_diag_isba%XP_C2)
351  nullify(ypack_diag_isba%XP_WGEQ)
352  nullify(ypack_diag_isba%XP_CT)
353  nullify(ypack_diag_isba%XP_RS)
354  nullify(ypack_diag_isba%XP_TS)
355  nullify(ypack_diag_isba%XP_TSRAD)
356  nullify(ypack_diag_isba%XP_T2M)
357  nullify(ypack_diag_isba%XP_Q2M)
358  nullify(ypack_diag_isba%XP_HU2M)
359  nullify(ypack_diag_isba%XP_ZON10M)
360  nullify(ypack_diag_isba%XP_MER10M)
361  nullify(ypack_diag_isba%XP_QS)
362  nullify(ypack_diag_isba%XP_SWI)
363  nullify(ypack_diag_isba%XP_TSWI)
364  nullify(ypack_diag_isba%XP_TWSNOW)
365  nullify(ypack_diag_isba%XP_TDSNOW)
366  nullify(ypack_diag_isba%XP_SWD)
367  nullify(ypack_diag_isba%XP_SWU)
368  nullify(ypack_diag_isba%XP_SWBD)
369  nullify(ypack_diag_isba%XP_SWBU)
370  nullify(ypack_diag_isba%XP_LWD)
371  nullify(ypack_diag_isba%XP_LWU)
372  nullify(ypack_diag_isba%XP_FMU)
373  nullify(ypack_diag_isba%XP_FMV)
374  nullify(ypack_diag_isba%XP_HORT)
375  nullify(ypack_diag_isba%XP_DRIP)
376  nullify(ypack_diag_isba%XP_IFLOOD)
377  nullify(ypack_diag_isba%XP_PFLOOD)
378  nullify(ypack_diag_isba%XP_LE_FLOOD)
379  nullify(ypack_diag_isba%XP_LEI_FLOOD)
380  nullify(ypack_diag_isba%XP_ICEFLUX)
381  nullify(ypack_diag_isba%XP_RRVEG)
382  nullify(ypack_diag_isba%XP_IRRIG_FLUX)
383  nullify(ypack_diag_isba%XP_GPP)
384  nullify(ypack_diag_isba%XP_RESP_AUTO)
385  nullify(ypack_diag_isba%XP_RESP_ECO)
386  nullify(ypack_diag_isba%XP_FAPAR)
387  nullify(ypack_diag_isba%XP_FAPIR)
388  nullify(ypack_diag_isba%XP_FAPAR_BS)
389  nullify(ypack_diag_isba%XP_FAPIR_BS)
390  nullify(ypack_diag_isba%XP_SWUP)
391  nullify(ypack_diag_isba%XP_SWNET_V)
392  nullify(ypack_diag_isba%XP_SWNET_G)
393  nullify(ypack_diag_isba%XP_SWNET_N)
394  nullify(ypack_diag_isba%XP_SWNET_NS)
395  nullify(ypack_diag_isba%XP_LWUP)
396  nullify(ypack_diag_isba%XP_LWNET_V)
397  nullify(ypack_diag_isba%XP_LWNET_G)
398  nullify(ypack_diag_isba%XP_LWNET_N)
399  nullify(ypack_diag_isba%XP_LEVCV)
400  nullify(ypack_diag_isba%XP_LESC)
401  nullify(ypack_diag_isba%XP_H_V_C)
402  nullify(ypack_diag_isba%XP_H_G_C)
403  nullify(ypack_diag_isba%XP_LETRGV)
404  nullify(ypack_diag_isba%XP_LETRCV)
405  nullify(ypack_diag_isba%XP_LERGV)
406  nullify(ypack_diag_isba%XP_LELITTER)
407  nullify(ypack_diag_isba%XP_LELITTERI)
408  nullify(ypack_diag_isba%XP_DRIPLIT)
409  nullify(ypack_diag_isba%XP_RRLIT)
410  nullify(ypack_diag_isba%XP_LERCV)
411  nullify(ypack_diag_isba%XP_H_C_A)
412  nullify(ypack_diag_isba%XP_H_N_C)
413  nullify(ypack_diag_isba%XP_LE_V_C)
414  nullify(ypack_diag_isba%XP_LE_G_C)
415  nullify(ypack_diag_isba%XP_LE_C_A)
416  nullify(ypack_diag_isba%XP_LE_N_C)
417  nullify(ypack_diag_isba%XP_EVAP_N_C)
418  nullify(ypack_diag_isba%XP_EVAP_G_C)
419  nullify(ypack_diag_isba%XP_SR_GN)
420  nullify(ypack_diag_isba%XP_MELTCV)
421  nullify(ypack_diag_isba%XP_FRZCV)
422  nullify(ypack_diag_isba%XP_SWDOWN_GN)
423  nullify(ypack_diag_isba%XP_LWDOWN_GN)
424  nullify(ypack_diag_isba%XP_DWG)
425  nullify(ypack_diag_isba%XP_DWGI)
426  nullify(ypack_diag_isba%XP_DWR)
427  nullify(ypack_diag_isba%XP_DSWE)
428  nullify(ypack_diag_isba%XP_WATBUD)
429 IF (lhook) CALL dr_hook("MODD_PACK_DIAG_ISBA_N:PACK_DIAG_ISBA_INIT",1,zhook_handle)
430 END SUBROUTINE pack_diag_isba_init
431 
432 
433 END MODULE modd_pack_diag_isba
subroutine pack_diag_isba_init(YPACK_DIAG_ISBA)