SURFEX  V8_0
Surfex V8_0 release
 All Classes Files Functions Variables
coupling_flake_sbln.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 ! ###############################################################################
6 SUBROUTINE coupling_flake_sbl_n (FM, DST, SLT, &
7  hprogram, hcoupling, ptstep, kyear, kmonth, kday, ptime, &
8  ki, ksv, ksw, ptsun, pzenith, pzenith2, pazim, pzref, puref, pu, pv, &
9  pqa, pta, prhoa, psv, pco2, hsv,prain, psnow, plw, pdir_sw, psca_sw, &
10  psw_bands, pps, ppa, psftq, psfth, psfts, psfco2, psfu, psfv, &
11  ptrad, pdir_alb, psca_alb, pemis, ptsurf, pz0, pz0h, pqsurf, &
12  ppew_a_coef, ppew_b_coef, ppet_a_coef, ppeq_a_coef, ppet_b_coef, &
13  ppeq_b_coef, htest )
14 ! ###############################################################################
15 !
16 !!**** *COUPLING_FLAKE_SBL_n * - Adds a SBL into SEAFLUX
17 !!
18 !! PURPOSE
19 !! -------
20 !
21 !!** METHOD
22 !! ------
23 !!
24 !! REFERENCE
25 !! ---------
26 !!
27 !!
28 !! AUTHOR
29 !! ------
30 !! V. Masson
31 !!
32 !! MODIFICATIONS
33 !! -------------
34 !! Original 09/2007
35 !! V. Masson 05/2009 Implicitation of momentum fluxes
36 !! S. Riette 06/2009 Initialisation of XT, PQ, XU and XTKE on canopy levels
37 !! S. Riette 10/2009 Iterative computation of XZ0
38 !! S. Riette 01/2010 Use of interpol_sbl to compute 10m wind diagnostic
39 !! B. Decharme 04/2013 new coupling variables
40 !----------------------------------------------------------------
41 !
42 !
43 !
44 USE modd_surfex_n, ONLY : flake_model_t
45 !
46 USE modd_dst_n, ONLY : dst_t
47 USE modd_slt_n, ONLY : slt_t
48 !
49 USE modd_surf_par, ONLY : xundef
50 USE modd_csts, ONLY : xcpd
51 !
52 USE mode_coupling_canopy
53 !
54 USE modi_init_water_sbl
55 !
56 USE modi_canopy_evol
57 USE modi_canopy_grid_update
58 !
59 USE modi_coupling_flake_n
60 !
61 USE yomhook ,ONLY : lhook, dr_hook
62 USE parkind1 ,ONLY : jprb
63 !
64 IMPLICIT NONE
65 !
66 !* 0.1 declarations of arguments
67 !
68 !
69 TYPE(flake_model_t), INTENT(INOUT) :: fm
70 TYPE(dst_t), INTENT(INOUT) :: dst
71 TYPE(slt_t), INTENT(INOUT) :: slt
72 !
73  CHARACTER(LEN=6), INTENT(IN) :: hprogram ! program calling surf. schemes
74  CHARACTER(LEN=1), INTENT(IN) :: hcoupling ! type of coupling
75  ! 'E' : explicit
76  ! 'I' : implicit
77 INTEGER, INTENT(IN) :: kyear ! current year (UTC)
78 INTEGER, INTENT(IN) :: kmonth ! current month (UTC)
79 INTEGER, INTENT(IN) :: kday ! current day (UTC)
80 REAL, INTENT(IN) :: ptime ! current time since midnight (UTC, s)
81 INTEGER, INTENT(IN) :: ki ! number of points
82 INTEGER, INTENT(IN) :: ksv ! number of scalars
83 INTEGER, INTENT(IN) :: ksw ! number of short-wave spectral bands
84 REAL, DIMENSION(KI), INTENT(IN) :: ptsun ! solar time (s from midnight)
85 REAL, INTENT(IN) :: ptstep ! atmospheric time-step (s)
86 REAL, DIMENSION(KI), INTENT(IN) :: pzref ! height of T,q forcing (m)
87 REAL, DIMENSION(KI), INTENT(IN) :: puref ! height of wind forcing (m)
88 !
89 REAL, DIMENSION(KI), INTENT(IN) :: pta ! air temperature forcing (K)
90 REAL, DIMENSION(KI), INTENT(IN) :: pqa ! air humidity forcing (kg/m3)
91 REAL, DIMENSION(KI), INTENT(IN) :: prhoa ! air density (kg/m3)
92 REAL, DIMENSION(KI,KSV),INTENT(IN) :: psv ! scalar variables
93 ! ! chemistry: first char. in HSV: '#' (molecule/m3)
94 ! !
95  CHARACTER(LEN=6), DIMENSION(KSV),INTENT(IN):: hsv ! name of all scalar variables
96 REAL, DIMENSION(KI), INTENT(IN) :: pu ! zonal wind (m/s)
97 REAL, DIMENSION(KI), INTENT(IN) :: pv ! meridian wind (m/s)
98 REAL, DIMENSION(KI,KSW),INTENT(IN) :: pdir_sw ! direct solar radiation (on horizontal surf.)
99 ! ! (W/m2)
100 REAL, DIMENSION(KI,KSW),INTENT(IN) :: psca_sw ! diffuse solar radiation (on horizontal surf.)
101 ! ! (W/m2)
102 REAL, DIMENSION(KSW),INTENT(IN) :: psw_bands ! mean wavelength of each shortwave band (m)
103 REAL, DIMENSION(KI), INTENT(IN) :: pzenith ! zenithal angle at t (radian from the vertical)
104 REAL, DIMENSION(KI), INTENT(IN) :: pzenith2 ! zenithal angle at t+1(radian from the vertical)
105 REAL, DIMENSION(KI), INTENT(IN) :: pazim ! azimuthal angle (radian from North, clockwise)
106 REAL, DIMENSION(KI), INTENT(IN) :: plw ! longwave radiation (on horizontal surf.)
107 ! ! (W/m2)
108 REAL, DIMENSION(KI), INTENT(IN) :: pps ! pressure at atmospheric model surface (Pa)
109 REAL, DIMENSION(KI), INTENT(IN) :: ppa ! pressure at forcing level (Pa)
110 REAL, DIMENSION(KI), INTENT(IN) :: pco2 ! CO2 concentration in the air (kg/m3)
111 REAL, DIMENSION(KI), INTENT(IN) :: psnow ! snow precipitation (kg/m2/s)
112 REAL, DIMENSION(KI), INTENT(IN) :: prain ! liquid precipitation (kg/m2/s)
113 !
114 !
115 REAL, DIMENSION(KI), INTENT(OUT) :: psfth ! flux of heat (W/m2)
116 REAL, DIMENSION(KI), INTENT(OUT) :: psftq ! flux of water vapor (kg/m2/s)
117 REAL, DIMENSION(KI), INTENT(OUT) :: psfu ! zonal momentum flux (Pa)
118 REAL, DIMENSION(KI), INTENT(OUT) :: psfv ! meridian momentum flux (Pa)
119 REAL, DIMENSION(KI), INTENT(OUT) :: psfco2 ! flux of CO2 (m/s*kg_CO2/kg_air)
120 REAL, DIMENSION(KI,KSV),INTENT(OUT):: psfts ! flux of scalar var. (kg/m2/s)
121 !
122 REAL, DIMENSION(KI), INTENT(OUT) :: ptrad ! radiative temperature (K)
123 REAL, DIMENSION(KI,KSW),INTENT(OUT):: pdir_alb! direct albedo for each spectral band (-)
124 REAL, DIMENSION(KI,KSW),INTENT(OUT):: psca_alb! diffuse albedo for each spectral band (-)
125 REAL, DIMENSION(KI), INTENT(OUT) :: pemis ! emissivity (-)
126 !
127 REAL, DIMENSION(KI), INTENT(OUT) :: ptsurf ! surface effective temperature (K)
128 REAL, DIMENSION(KI), INTENT(OUT) :: pz0 ! roughness length for momentum (m)
129 REAL, DIMENSION(KI), INTENT(OUT) :: pz0h ! roughness length for heat (m)
130 REAL, DIMENSION(KI), INTENT(OUT) :: pqsurf ! specific humidity at surface (kg/kg)
131 !
132 REAL, DIMENSION(KI), INTENT(IN) :: ppew_a_coef! implicit coefficients
133 REAL, DIMENSION(KI), INTENT(IN) :: ppew_b_coef! needed if HCOUPLING='I'
134 REAL, DIMENSION(KI), INTENT(IN) :: ppet_a_coef
135 REAL, DIMENSION(KI), INTENT(IN) :: ppeq_a_coef
136 REAL, DIMENSION(KI), INTENT(IN) :: ppet_b_coef
137 REAL, DIMENSION(KI), INTENT(IN) :: ppeq_b_coef
138  CHARACTER(LEN=2), INTENT(IN) :: htest ! must be equal to 'OK'
139 !
140 !* 0.2 declarations of local variables
141 !
142 !* forcing variables
143 !
144 REAL, DIMENSION(KI) :: zwind ! lowest atmospheric level wind speed (m/s)
145 REAL, DIMENSION(KI) :: zexna ! Exner function at lowest SBL scheme level (-)
146 REAL, DIMENSION(KI) :: zta ! temperature (K)
147 REAL, DIMENSION(KI) :: zpa ! pressure (Pa)
148 REAL, DIMENSION(KI) :: zzref ! temperature forcing level (m)
149 REAL, DIMENSION(KI) :: zuref ! wind forcing level (m)
150 REAL, DIMENSION(KI) :: zu ! zonal wind (m/s)
151 REAL, DIMENSION(KI) :: zv ! meridian wind (m/s)
152 REAL, DIMENSION(KI) :: zqa ! specific humidity (kg/m3)
153 REAL, DIMENSION(KI) :: zpeq_a_coef ! specific humidity implicit
154 REAL, DIMENSION(KI) :: zpeq_b_coef ! coefficients (hum. in kg/kg)
155 !
156 !
157 ! SBL turbulence scheme
158 !
159 REAL, DIMENSION(KI) :: zsflux_u ! Surface flux u'w' (m2/s2)
160 REAL, DIMENSION(KI) :: zsflux_t ! Surface flux w'T' (mK/s)
161 REAL, DIMENSION(KI) :: zsflux_q ! Surface flux w'q' (kgm2/s)
162 REAL, DIMENSION(KI,FM%FSB%NLVL) :: zforc_u ! tendency due to drag force for wind
163 REAL, DIMENSION(KI,FM%FSB%NLVL) :: zdforc_udu! formal derivative of
164 ! ! tendency due to drag force for wind
165 REAL, DIMENSION(KI,FM%FSB%NLVL) :: zforc_e ! tendency due to drag force for TKE
166 REAL, DIMENSION(KI,FM%FSB%NLVL) :: zdforc_ede! formal derivative of
167 ! ! tendency due to drag force for TKE
168 REAL, DIMENSION(KI,FM%FSB%NLVL) :: zforc_t ! tendency due to drag force for Temp
169 REAL, DIMENSION(KI,FM%FSB%NLVL) :: zdforc_tdt! formal derivative of
170 ! ! tendency due to drag force for Temp
171 REAL, DIMENSION(KI,FM%FSB%NLVL) :: zforc_q ! tendency due to drag force for Temp
172 REAL, DIMENSION(KI,FM%FSB%NLVL) :: zdforc_qdq! formal derivative of
173 ! ! tendency due to drag force for hum.
174 REAL, DIMENSION(KI,FM%FSB%NLVL) :: zlmo ! MO length
175 REAL, DIMENSION(KI,FM%FSB%NLVL) :: zlm ! mixing length
176 REAL, DIMENSION(KI,FM%FSB%NLVL) :: zleps ! dissipative length
177 REAL, DIMENSION(KI) :: zh ! canopy height (m)
178 REAL, DIMENSION(KI) :: zustar ! friction velocity (m/s)
179 !
180 REAL, DIMENSION(KI) :: zpet_a_coef ! temperature implicit
181 REAL, DIMENSION(KI) :: zpet_b_coef ! coefficients (K)
182 REAL, DIMENSION(KI) :: zpew_a_coef ! wind implicit
183 REAL, DIMENSION(KI) :: zpew_b_coef ! coefficients (m/s)
184 
185 REAL, DIMENSION(KI) :: zalfau ! V+(1) = - alfa rho u'w'(1) + beta
186 REAL, DIMENSION(KI) :: zbetau ! V+(1) = - alfa rho u'w'(1) + beta
187 REAL, DIMENSION(KI) :: zalfath ! Th+(1) = - alfa rho w'th'(1) + beta
188 REAL, DIMENSION(KI) :: zbetath ! Th+(1) = - alfa rho w'th'(1) + beta
189 REAL, DIMENSION(KI) :: zalfaq ! Q+(1) = - alfa rho w'q'(1) + beta
190 REAL, DIMENSION(KI) :: zbetaq ! Q+(1) = - alfa rho w'q'(1) + beta
191 !
192  CHARACTER(LEN=1) :: gcoupling
193 REAL(KIND=JPRB) :: zhook_handle
194 !-------------------------------------------------------------------------------------
195 !
196 !
197 !* 1. Preliminary computations of the SBL scheme
198 ! ------------------------------------------
199 !
200 IF (lhook) CALL dr_hook('COUPLING_FLAKE_SBL_N',0,zhook_handle)
201 IF (fm%F%LSBL) THEN
202 !
203 !* 1.1 Updates SBL vertical grid as a function of forcing height
204 ! ---------------------------------------------------------
205 !
206 !* determines where is the forcing level and modifies the upper levels of the canopy grid
207 !
208  zh = 0.
209  CALL canopy_grid_update(ki,fm%FSB%NLVL,zh,puref,&
210  fm%FSB%XZ,fm%FSB%XZF,fm%FSB%XDZ,fm%FSB%XDZF)
211 !
212 !
213 !
214 !* 1.2 Initialisation at first time step
215 ! ---------------------------------
216 !
217  IF(any(fm%FSB%XT(:,:) == xundef)) THEN
218  CALL init_water_sbl(fm%FSB%NLVL, ppa, pps, pta, pqa, prhoa, pu, pv, prain, psnow, &
219  psfth, psftq, pzref, puref, fm%F%XTS, fm%F%XZ0, fm%FSB%XZ, &
220  fm%FSB%XT, fm%FSB%XQ, fm%FSB%XU, fm%FSB%XTKE, fm%FSB%XP)
221  ENDIF
222 !
223 !
224 !* 1.3 Allocations
225 ! -----------
226 !
227  CALL init_forc( zforc_u, zdforc_udu, zforc_e, zdforc_ede, &
228  zforc_t, zdforc_tdt, zforc_q, zdforc_qdq )
229 !
230  zsflux_u = 0.
231  zsflux_t = 0.
232  zsflux_q = 0.
233 !
234  zlmo = spread(fm%FSB%XLMO,2,fm%FSB%NLVL)
235 !
236 !* 1.3 Computes coefficients for implicitation
237 ! ---------------------------------------
238 !
239  zwind = sqrt(pu**2+pv**2)
240  CALL canopy_evol(ki,fm%FSB%NLVL,ptstep,1,fm%FSB%XZ,zwind,pta,pqa,ppa,prhoa, &
241  zsflux_u,zsflux_t,zsflux_q, &
242  zforc_u,zdforc_udu,zforc_e,zdforc_ede, &
243  zforc_t,zdforc_tdt,zforc_q,zdforc_qdq, &
244  fm%FSB%XZ,fm%FSB%XZF,fm%FSB%XDZ,fm%FSB%XDZF, &
245  fm%FSB%XU,fm%FSB%XTKE,fm%FSB%XT, &
246  fm%FSB%XQ,zlmo,zlm,zleps,fm%FSB%XP,zustar, &
247  zalfau,zbetau,zalfath,zbetath,zalfaq,zbetaq )
248 
249 !
250 !* 1.5 Goes from atmospheric forcing to canopy forcing height
251 ! ------------------------------------------------------
252 !
253  gcoupling = 'E'
254 !
255  CALL init_coupling_canopy( fm%FSB%XP(:,1), ppa, fm%FSB%XT(:,1), &
256  fm%FSB%XQ(:,1), pu, pv, fm%FSB%XZ(:,1), &
257  fm%FSB%XU(:,1), &
258  prhoa, zalfau, zbetau, zalfath, &
259  zbetath, zalfaq, zbetaq, &
260  zpa, zta, zqa, zu, zv, &
261  zuref, zzref, zexna, &
262  zpew_a_coef, zpew_b_coef, &
263  zpet_a_coef, zpet_b_coef, &
264  zpeq_a_coef, zpeq_b_coef )
265 !
266 !-------------------------------------------------------------------------------------
267 ELSE
268 !-------------------------------------------------------------------------------------
269 !
270 !* 2. If no SBL scheme is used, forcing is not modified
271 ! -------------------------------------------------
272 !
273  gcoupling = hcoupling
274 !
275  CALL init_coupling( hcoupling, &
276  pps, ppa, pta, pqa, pu, pv, &
277  puref, pzref, &
278  ppew_a_coef, ppew_b_coef, &
279  ppet_a_coef, ppet_b_coef, &
280  ppeq_a_coef, ppeq_b_coef, &
281  zpa, zta, zqa, zu, zv, &
282  zuref, zzref, &
283  zpew_a_coef, zpew_b_coef, &
284  zpet_a_coef, zpet_b_coef, &
285  zpeq_a_coef, zpeq_b_coef )
286 !
287 END IF
288 !
289 !-------------------------------------------------------------------------------------
290 !
291 !* 2. Call of SEAFLUX
292 ! ------------
293 !
294  CALL coupling_flake_n(fm, dst, slt, &
295  hprogram, gcoupling, &
296  ptstep, kyear, kmonth, kday, ptime, &
297  ki, ksv, ksw, &
298  ptsun, pzenith, pzenith2, pazim, &
299  zzref, zuref, zu, zv, zqa, zta, prhoa, psv, pco2, hsv, &
300  prain, psnow, plw, pdir_sw, psca_sw, psw_bands, pps, zpa, &
301  psftq, psfth, psfts, psfco2, psfu, psfv, &
302  ptrad, pdir_alb, psca_alb, pemis, ptsurf, pz0, pz0h, pqsurf, &
303  zpew_a_coef, zpew_b_coef, &
304  zpet_a_coef, zpeq_a_coef, zpet_b_coef, zpeq_b_coef, &
305  htest )
306 !
307 !-------------------------------------------------------------------------------------
308 !
309 !* 3. End if no SBL is used
310 ! ---------------------
311 !
312 IF (.NOT. fm%F%LSBL .AND. lhook) CALL dr_hook('COUPLING_FLAKE_SBL_N',1,zhook_handle)
313 IF (.NOT. fm%F%LSBL) RETURN
314 !
315 !-------------------------------------------------------------------------------------
316 !
317 !* 4. Computes the impact of canopy and surfaces on air
318 ! -------------------------------------------------
319 !
320  CALL init_forc( zforc_u, zdforc_udu, zforc_e, zdforc_ede, &
321  zforc_t, zdforc_tdt, zforc_q, zdforc_qdq )
322 !
323 zsflux_u = - sqrt(psfu(:)**2+psfv(:)**2) / prhoa(:)
324 zsflux_t(:) = psfth(:) / xcpd * zexna(:) / prhoa(:)
325 zsflux_q(:) = psftq(:)
326 !
327 !-------------------------------------------------------------------------------------
328 !
329 !* 6. Evolution of canopy air due to these impacts
330 ! --------------------------------------------
331 !
332 zwind = sqrt(pu**2+pv**2)
333  CALL canopy_evol(ki,fm%FSB%NLVL,ptstep,2,fm%FSB%XZ,zwind,pta,pqa,ppa,prhoa, &
334  zsflux_u,zsflux_t,zsflux_q, &
335  zforc_u,zdforc_udu,zforc_e,zdforc_ede, &
336  zforc_t,zdforc_tdt,zforc_q,zdforc_qdq, &
337  fm%FSB%XZ,fm%FSB%XZF,fm%FSB%XDZ,fm%FSB%XDZF, &
338  fm%FSB%XU,fm%FSB%XTKE,fm%FSB%XT, &
339  fm%FSB%XQ,zlmo,zlm,zleps,fm%FSB%XP,zustar, &
340  zalfau,zbetau,zalfath,zbetath,zalfaq,zbetaq )
341 !
342 fm%FSB%XLMO(:) = zlmo(:,fm%FSB%NLVL)
343 !
344 !-------------------------------------------------------------------------------------
345 !
346 !* 7. 2m and 10m diagnostics if canopy is used
347 ! ----------------------------------------
348 !
349 !
350 IF (fm%F%LSBL .AND. fm%DGF%N2M>=1) CALL init_2m_10m( fm%FSB%XP(:,2), fm%FSB%XT(:,2), &
351  fm%FSB%XQ(:,2), fm%FSB%XU, fm%FSB%XZ, &
352  pu, pv, zwind, prhoa, &
353  fm%DGF%XT2M, fm%DGF%XQ2M, fm%DGF%XHU2M, &
354  fm%DGF%XZON10M, fm%DGF%XMER10M,&
355  fm%DGF%XWIND10M, fm%DGF%XWIND10M_MAX, fm%DGF%XT2M_MIN, &
356  fm%DGF%XT2M_MAX, fm%DGF%XHU2M_MIN, fm%DGF%XHU2M_MAX )
357 !
358 IF (lhook) CALL dr_hook('COUPLING_FLAKE_SBL_N',1,zhook_handle)
359 !
360 !-------------------------------------------------------------------------------------
361 !
362 END SUBROUTINE coupling_flake_sbl_n
modd_slt_n
Definition: modd_sltn.F90:5
modd_slt_n::slt_t
Definition: modd_sltn.F90:18
mode_coupling_canopy::init_coupling_canopy
subroutine init_coupling_canopy(PP, PPA, PT, PQ, PU, PV, PZ, PXU, PRHOA, PALFAU, PBETAU, PALFATH, PBETATH, PALFAQ, PBETAQ, PPPA, PTTA, PQQA, PUU, PVV, PUUREF, PZZREF, PEXNA, PPEW_AA_COEF, PPEW_BB_COEF, PPET_AA_COEF, PPET_BB_COEF, PPEQ_AA_COEF, PPEQ_BB_COEF)
Definition: mode_coupling_canopy.F90:51
modd_surfex_n
Definition: modd_surfexn.F90:5
modd_surfex_n::flake_model_t
Definition: modd_surfexn.F90:95
coupling_flake_sbl_n
subroutine coupling_flake_sbl_n(FM, DST, SLT, HPROGRAM, HCOUPLING, PTSTEP, KYEAR, KMONTH, KDAY, PTIME, KI, KSV, KSW, PTSUN, PZENITH, PZENITH2, PAZIM, PZREF, PUREF, PU, PV, PQA, PTA, PRHOA, PSV, PCO2, HSV, PRAIN, PSNOW, PLW, PDIR_SW, PSCA_SW, PSW_BANDS, PPS, PPA, PSFTQ, PSFTH, PSFTS, PSFCO2, PSFU, PSFV, PTRAD, PDIR_ALB, PSCA_ALB, PEMIS, PTSURF, PZ0, PZ0H, PQSURF, PPEW_A_COEF, PPEW_B_COEF, PPET_A_COEF, PPEQ_A_COEF, PPET_B_COEF, PPEQ_B_COEF, HTEST)
Definition: coupling_flake_sbln.F90:6
canopy_grid_update
subroutine canopy_grid_update(KI, KLVL, PH, PZFORC, PZ, PZF, PDZ, PDZF)
Definition: canopy_grid_update.F90:6
mode_coupling_canopy::init_coupling
subroutine init_coupling(HCOUPLING, PPS, PPA, PTA, PQA, PU, PV, PUREF, PZREF, PPEW_A_COEF, PPEW_B_COEF, PPET_A_COEF, PPET_B_COEF, PPEQ_A_COEF, PPEQ_B_COEF, PPPA, PTTA, PQQA, PUU, PVV, PUUREF, PZZREF, PPEW_AA_COEF, PPEW_BB_COEF, PPET_AA_COEF, PPET_BB_COEF, PPEQ_AA_COEF, PPEQ_BB_COEF)
Definition: mode_coupling_canopy.F90:132
modd_dst_n::dst_t
Definition: modd_dstn.F90:19
coupling_flake_n
subroutine coupling_flake_n(FM, DST, SLT, HPROGRAM, HCOUPLING, PTSTEP, KYEAR, KMONTH, KDAY, PTIME, KI, KSV, KSW, PTSUN, PZENITH, PZENITH2, PAZIM, PZREF, PUREF, PU, PV, PQA, PTA, PRHOA, PSV, PCO2, HSV, PRAIN, PSNOW, PLW, PDIR_SW, PSCA_SW, PSW_BANDS, PPS, PPA, PSFTQ, PSFTH, PSFTS, PSFCO2, PSFU, PSFV, PTSRAD, PDIR_ALB, PSCA_ALB, PEMIS, PTSURF, PZ0, PZ0H, PQSURF, PPEW_A_COEF, PPEW_B_COEF, PPET_A_COEF, PPEQ_A_COEF, PPET_B_COEF, PPEQ_B_COEF, HTEST)
Definition: coupling_flaken.F90:6
mode_coupling_canopy::init_2m_10m
subroutine init_2m_10m(PP, PT, PQ, PXU, PXZ, PU, PV, PWIND, PRHOA, PT2M, PQ2M, PHU2M, PZON10M, PMER10M, PWIND10M, PWIND10M_MAX, PT2M_MIN, PT2M_MAX, PHU2M_MIN, PHU2M_MAX)
Definition: mode_coupling_canopy.F90:210
mode_coupling_canopy
Definition: mode_coupling_canopy.F90:6
modd_surf_par
Definition: modd_surf_par.F90:6
mode_coupling_canopy::init_forc
subroutine init_forc(PFORC_U, PDFORC_UDU, PFORC_E, PDFORC_EDE, PFORC_T, PDFORC_TDT, PFORC_Q, PDFORC_QDQ)
Definition: mode_coupling_canopy.F90:15
init_water_sbl
subroutine init_water_sbl(KLVL, PPA, PPS, PTA, PQA, PRHOA, PU, PV, PRAIN, PSNOW, PSFTH, PSFTQ, PZREF, PUREF, PTS, PZ0, PZ, PT, PQ, PWIND, PTKE, PP)
Definition: init_water_sbl.F90:6
canopy_evol
subroutine canopy_evol(KI, KLVL, PTSTEP, KIMPL, PZZ, PWIND, PTA, PQA, PPA, PRHOA, PSFLUX_U, PSFLUX_T, PSFLUX_Q, PFORC_U, PDFORC_UDU, PFORC_E, PDFORC_EDE, PFORC_T, PDFORC_TDT, PFORC_Q, PDFORC_QDQ, PZ, PZF, PDZ, PDZF, PU, PTKE, PT, PQ, PLMO, PLM, PLEPS, PP, PUSTAR, PALFAU, PBETAU, PALFATH, PBETATH, PALFAQ, PBETAQ, ONEUTRAL)
Definition: canopy_evol.F90:6
modd_dst_n
Definition: modd_dstn.F90:5
modd_csts
Definition: modd_csts.F90:6
Generated on Fri May 13 2016 09:31:01 for SURFEX by   doxygen 1.8.5