SURFEX  V8_0
Surfex V8_0 release
 All Classes Files Functions Variables
modd_diag_idealn.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_IDEAL - declaration of diagnostics for IDEAL 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 04/2009
29 !! P. Le Moigne 03/2015: add diagnostics IDEAL case
30 !
31 !* 0. DECLARATIONS
32 ! ------------
33 !
34 !
35 !
36 USE yomhook ,ONLY : lhook, dr_hook
37 USE parkind1 ,ONLY : jprb
38 !
39 IMPLICIT NONE
40 
42 !------------------------------------------------------------------------------
43 !
44  REAL :: XDIAG_TSTEP ! time step for diagnostics writing
45 !
46  INTEGER :: N2M ! flag for 2 meters (and 10 meters) quantities
47  LOGICAL :: L2M_MIN_ZS ! flag for 2 meters quantities evaluated on
48 ! ! the minimum orographyy of the grid
49  LOGICAL :: LSURF_BUDGET ! flag for surface energy budget
50  LOGICAL :: LRAD_BUDGET ! flag for radiative energy budget
51  LOGICAL :: LCOEF ! flag for transfer coefficients
52  LOGICAL :: LSURF_VARS ! flag for surface variables
53  LOGICAL :: LSURF_BUDGETC ! flag for surface cumulated energy budget
54  LOGICAL :: LRESET_BUDGETC ! flag for surface cumulated energy budget
55 
56 !* averaged variables
57 !
58  REAL, POINTER, DIMENSION(:) :: XRI ! Bulk-Richardson number (-)
59  REAL, POINTER, DIMENSION(:) :: XCD ! drag coefficient for wind (W/s2)
60  REAL, POINTER, DIMENSION(:) :: XCH ! drag coefficient for heat (W/s)
61  REAL, POINTER, DIMENSION(:) :: XCE ! drag coefficient for vapor (W/s/K)
62  REAL, POINTER, DIMENSION(:) :: XZ0 ! roughness length for momentum (m)
63  REAL, POINTER, DIMENSION(:) :: XZ0H ! roughness length for heat (m)
64  REAL, POINTER, DIMENSION(:) :: XRN ! net radiation at surface (W/m2)
65  REAL, POINTER, DIMENSION(:) :: XH ! sensible heat flux (W/m2)
66  REAL, POINTER, DIMENSION(:) :: XLE ! latent heat flux (W/m2)
67  REAL, POINTER, DIMENSION(:) :: XLEI ! sublimation latent heat flux (W/m2)
68  REAL, POINTER, DIMENSION(:) :: XGFLUX ! net soil-vegetation flux (W/m2)
69  REAL, POINTER, DIMENSION(:) :: XEVAP ! total evaporation (kg/m2/s)
70  REAL, POINTER, DIMENSION(:) :: XSUBL ! sublimation (kg/m2/s)
71  REAL, POINTER, DIMENSION(:) :: XT2M ! air temperature at 2 meters (K)
72  REAL, POINTER, DIMENSION(:) :: XT2M_MIN ! Minimum air temperature at 2 meters (K)
73  REAL, POINTER, DIMENSION(:) :: XT2M_MAX ! Maximum air temperature at 2 meters (K)
74  REAL, POINTER, DIMENSION(:) :: XQ2M ! air humidity at 2 meters (kg/kg)
75  REAL, POINTER, DIMENSION(:) :: XHU2M ! relative humidity at 2 meters (-)
76  REAL, POINTER, DIMENSION(:) :: XHU2M_MIN! Minimum relative humidity at 2 meters (-)
77  REAL, POINTER, DIMENSION(:) :: XHU2M_MAX! Maximum relative humidity at 2 meters (-)
78  REAL, POINTER, DIMENSION(:) :: XQS ! humidity at surface (kg/kg)
79  REAL, POINTER, DIMENSION(:) :: XZON10M ! zonal wind at 10 meters (m/s)
80  REAL, POINTER, DIMENSION(:) :: XMER10M ! meridian wind at 10 meters (m/s)
81  REAL, POINTER, DIMENSION(:) :: XWIND10M ! wind at 10 meters (m/s)
82  REAL, POINTER, DIMENSION(:) :: XWIND10M_MAX! Maximum wind at 10 meters (m/s)
83  REAL, POINTER, DIMENSION(:) :: XLWD ! downward long wave radiation (W/m2)
84  REAL, POINTER, DIMENSION(:) :: XLWU ! upward long wave radiation (W/m2)
85  REAL, POINTER, DIMENSION(:) :: XSWD ! downward short wave radiation (W/m2)
86  REAL, POINTER, DIMENSION(:) :: XSWU ! upward short wave radiation (W/m2)
87  REAL, POINTER, DIMENSION(:,:) :: XSWBD ! downward short wave radiation by spectral band (W/m2)
88  REAL, POINTER, DIMENSION(:,:) :: XSWBU ! upward short wave radiation by spectral band (W/m2)
89  REAL, POINTER, DIMENSION(:) :: XFMU ! horizontal momentum flux zonal (m2/s2)
90  REAL, POINTER, DIMENSION(:) :: XFMV ! horizontal momentum flux meridian (m2/s2)
91  REAL, POINTER, DIMENSION(:) :: XDIAG_TS ! water surface temperature (K)
92  REAL, POINTER, DIMENSION(:) :: XALBT ! Total Albedo
93  REAL, POINTER, DIMENSION(:) :: XSWE ! snow water equivalent (kg/m2)
94 !
95 !* cumulated averaged variables
96 !
97  REAL, POINTER, DIMENSION(:) :: XRNC ! net radiation at surface (J/m2)
98  REAL, POINTER, DIMENSION(:) :: XHC ! sensible heat flux (J/m2)
99  REAL, POINTER, DIMENSION(:) :: XLEC ! total latent heat flux (J/m2)
100  REAL, POINTER, DIMENSION(:) :: XLEIC ! sublimation latent heat flux (J/m2)
101  REAL, POINTER, DIMENSION(:) :: XGFLUXC ! net soil-vegetation flux (J/m2)
102  REAL, POINTER, DIMENSION(:) :: XEVAPC ! total evaporation (kg/m2)
103  REAL, POINTER, DIMENSION(:) :: XSUBLC ! sublimation (kg/m2)
104  REAL, POINTER, DIMENSION(:) :: XLWDC ! downward long wave radiation (J/m2)
105  REAL, POINTER, DIMENSION(:) :: XLWUC ! upward long wave radiation (J/m2)
106  REAL, POINTER, DIMENSION(:) :: XSWDC ! downward short wave radiation (J/m2)
107  REAL, POINTER, DIMENSION(:) :: XSWUC ! upward short wave radiation (J/m2)
108  REAL, POINTER, DIMENSION(:) :: XFMUC ! horizontal momentum flux zonal (kg/ms)
109  REAL, POINTER, DIMENSION(:) :: XFMVC ! horizontal momentum flux meridian (kg/ms)
110 !
111 !------------------------------------------------------------------------------
112 !
113 
114 END TYPE diag_ideal_t
115 
116 
117 
118  CONTAINS
119 
120 !
121 
122 
123 
124 
125 SUBROUTINE diag_ideal_init(YDIAG_IDEAL)
126 TYPE(diag_ideal_t), INTENT(INOUT) :: ydiag_ideal
127 REAL(KIND=JPRB) :: zhook_handle
128 IF (lhook) CALL dr_hook("MODD_DIAG_IDEAL_N:DIAG_IDEAL_INIT",0,zhook_handle)
129  nullify(ydiag_ideal%XRI)
130  nullify(ydiag_ideal%XCD)
131  nullify(ydiag_ideal%XCH)
132  nullify(ydiag_ideal%XCE)
133  nullify(ydiag_ideal%XZ0)
134  nullify(ydiag_ideal%XZ0H)
135  nullify(ydiag_ideal%XRN)
136  nullify(ydiag_ideal%XH)
137  nullify(ydiag_ideal%XLE)
138  nullify(ydiag_ideal%XLEI)
139  nullify(ydiag_ideal%XGFLUX)
140  nullify(ydiag_ideal%XEVAP)
141  nullify(ydiag_ideal%XSUBL)
142  nullify(ydiag_ideal%XT2M)
143  nullify(ydiag_ideal%XT2M_MIN)
144  nullify(ydiag_ideal%XT2M_MAX)
145  nullify(ydiag_ideal%XQ2M)
146  nullify(ydiag_ideal%XHU2M)
147  nullify(ydiag_ideal%XHU2M_MIN)
148  nullify(ydiag_ideal%XHU2M_MAX)
149  nullify(ydiag_ideal%XQS)
150  nullify(ydiag_ideal%XZON10M)
151  nullify(ydiag_ideal%XMER10M)
152  nullify(ydiag_ideal%XWIND10M)
153  nullify(ydiag_ideal%XWIND10M_MAX)
154  nullify(ydiag_ideal%XLWD)
155  nullify(ydiag_ideal%XLWU)
156  nullify(ydiag_ideal%XSWD)
157  nullify(ydiag_ideal%XSWU)
158  nullify(ydiag_ideal%XSWBD)
159  nullify(ydiag_ideal%XSWBU)
160  nullify(ydiag_ideal%XFMU)
161  nullify(ydiag_ideal%XFMV)
162  nullify(ydiag_ideal%XDIAG_TS)
163  nullify(ydiag_ideal%XALBT)
164  nullify(ydiag_ideal%XSWE)
165  nullify(ydiag_ideal%XRNC)
166  nullify(ydiag_ideal%XHC)
167  nullify(ydiag_ideal%XLEC)
168  nullify(ydiag_ideal%XLEIC)
169  nullify(ydiag_ideal%XGFLUXC)
170  nullify(ydiag_ideal%XEVAPC)
171  nullify(ydiag_ideal%XSUBLC)
172  nullify(ydiag_ideal%XLWDC)
173  nullify(ydiag_ideal%XLWUC)
174  nullify(ydiag_ideal%XSWDC)
175  nullify(ydiag_ideal%XSWUC)
176  nullify(ydiag_ideal%XFMUC)
177  nullify(ydiag_ideal%XFMVC)
178 ydiag_ideal%XDIAG_TSTEP=0.
179 ydiag_ideal%N2M=0
180 ydiag_ideal%L2M_MIN_ZS=.false.
181 ydiag_ideal%LSURF_BUDGET=.false.
182 ydiag_ideal%LRAD_BUDGET=.false.
183 ydiag_ideal%LCOEF=.false.
184 ydiag_ideal%LSURF_VARS=.false.
185 ydiag_ideal%LSURF_BUDGETC=.false.
186 ydiag_ideal%LRESET_BUDGETC=.false.
187 IF (lhook) CALL dr_hook("MODD_DIAG_IDEAL_N:DIAG_IDEAL_INIT",1,zhook_handle)
188 END SUBROUTINE diag_ideal_init
189 
190 
191 END MODULE modd_diag_ideal_n
subroutine diag_ideal_init(YDIAG_IDEAL)