Newer
Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
!MNH_LIC Copyright 1995-2021 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
!-----------------------------------------------------------------
! ########################
MODULE MODI_PHYS_PARAM_n
! ########################
!
!
INTERFACE
!
SUBROUTINE PHYS_PARAM_n( KTCOUNT, TPFILE, &
PRAD, PSHADOWS, PKAFR, PGROUND, PMAFL, PDRAG,PEOL, PTURB, &
PTRACER, PTIME_BU, PWETDEPAER, OMASKkids, OCLOUD_ONLY )
!
USE MODD_IO, ONLY: TFILEDATA
use modd_precision, only: MNHTIME
!
INTEGER, INTENT(IN) :: KTCOUNT ! temporal iteration count
TYPE(TFILEDATA), INTENT(IN) :: TPFILE ! Synchronous output file
! advection schemes
REAL(kind=MNHTIME), DIMENSION(2), INTENT(INOUT) :: PRAD,PSHADOWS,PKAFR,PGROUND,PTURB,PMAFL,PDRAG,PTRACER,PEOL ! to store CPU
! time for computing time
REAL(kind=MNHTIME), DIMENSION(2), INTENT(INOUT) :: PTIME_BU ! time used in budget&LES budgets statistics
REAL, DIMENSION(:,:,:,:), INTENT(INOUT) :: PWETDEPAER
LOGICAL, DIMENSION(:,:), INTENT(IN) :: OMASKkids ! kids domains mask
LOGICAL, INTENT(OUT) :: OCLOUD_ONLY ! conditionnal radiation computations for
! the only cloudy columns
!
END SUBROUTINE PHYS_PARAM_n
!
END INTERFACE
!
END MODULE MODI_PHYS_PARAM_n
!
! ########################################################################################
SUBROUTINE PHYS_PARAM_n( KTCOUNT, TPFILE, &
PRAD, PSHADOWS, PKAFR, PGROUND, PMAFL, PEOL, PDRAG, PTURB, &
PTRACER, PTIME_BU, PWETDEPAER, OMASKkids, OCLOUD_ONLY )
! ########################################################################################
!
!!**** *PHYS_PARAM_n * -monitor of the parameterizations used by model _n
!!
!! PURPOSE
!! -------
! The purpose of this routine is to update the sources by adding the
! parameterized terms. This is realized by sequentially calling the
! specialized routines.
!
!!** METHOD
!! ------
!! The first parametrization is the radiation scheme:
!! ----------------
!! * CRAD = 'FIXE'
!! In this case, a temporal interpolation is performed for the downward
!! surface fluxes XFLALWD and XFLASWD.
!! * CRAD = 'ECMWF'
!! Several tests are performed before calling the radiation computations
!! interface with the ECMWF radiation scheme code. A control is made to
!! ensure that:
!! - the full radiation code is called at the first model timestep
!! - there is a priority for calling the full radiation instead of the
!! cloud-only approximation if both must be called at the current
!! timestep
!! - the cloud-only option (approximation) is coherent with the
!! occurence of one cloudy vertical column at least
!! If all the above conditions are fulfilled (GRAD is .TRUE.) then the
!! position of the sun is computed in routine SUNPOS_n and the interfacing
!! routine RADIATIONS is called to update the radiative tendency XDTHRAD
!! and the downward surface fluxes XFLALWD and XFLASWD. Finally, the
!! radiative tendency is integrated as a source term in the THETA prognostic
!! equation.
!!
!! The second parameterization is the soil scheme:
!! -----------
!!
!! externalized surface
!!
!! The third parameterization is the turbulence scheme:
!! -----------------
!! * CTURB='NONE'
!! no turbulent mixing is taken into account
!! * CTURB='TKEL'
!! The turbulent fluxes are computed according to a one and half order
!! closure of the hydrodynamical equations. This scheme is based on a
!! prognostic for the turbulent kinetic energy and a mixing length
!! computation ( the mesh size or a physically based length). Other
!! turbulent moments are diagnosed according to a stationarization of the
!! second order turbulent moments. This turbulent scheme forecasts
!! either a purely vertical turbulent mixing or 3-dimensional mixing
!! according to its internal degrees of freedom.
!!
!!
!! The LAST parameterization is the chemistry scheme:
!! -----------------
!! The chemistry part of MesoNH has two namelists, NAM_SOLVER for the
!! parameters concerning the stiff solver, and NAM_MNHCn concerning the
!! configuration and options of the chemistry module itself.
!! The switch LUSECHEM in NAM_CONF acitvates or deactivates the chemistry.
!! The only variables of MesoNH that are modified by chemistry are the
!! scalar variables. If calculation of chemical surface fluxes is
!! requested, those fluxes are calculated before
!! entering the turbulence scheme, since those fluxes are taken into
!! account by TURB as surface boundary conditions.
!! CAUTION: chemistry has allways to be called AFTER ALL OTHER TERMS
!! that affect the scalar variables (dynamical terms, forcing,
!! parameterizations (like TURB, CONVECTION), since it uses the variables
!! XRSVS as input in case of the time-split option.
!!
!! EXTERNAL
!! --------
!! Subroutine SUNPOS_n : computes the position of the sun
!! Subroutine RADIATIONS : computes the radiative tendency and fluxes
!! Subroutine TSZ0 : computes the surface from temporally
!! interpolated Ts and given z0
!! Subroutine ISBA : computes the surface fluxes from a soil scheme
!! Subroutine TURB : computes the turbulence source terms
!! Subroutine CONVECTION : computes the convection source term
!! Subroutine CH_SURFACE_FLUX_n: computes the surface flux for chemical
!! species
!! Subroutine CH_MONITOR_n : computes the chemistry source terms
!! that are applied to the scalar variables
!!
!! IMPLICIT ARGUMENTS
!! ------------------
!! USE MODD_DYN
!! USE MODD_CONF
!! USE MODD_CONF_n
!! USE MODD_CURVCOR_n
!! USE MODD_DYN_n
!! USE MODD_FIELD_n
!! USE MODD_GR_FIELD_n
!! USE MODD_LSFIELD_n
!! USE MODD_GRID_n
!! USE MODD_LBC_n
!! USE MODD_PARAM_RAD_n
!! USE MODD_RADIATIONS_n
!! USE MODD_REF_n
!! USE MODD_LUNIT_n
!! USE MODD_TIME_n
!! USE MODD_CH_MNHC_n
!!
!! REFERENCE
!! ---------
!! None
!!
!! AUTHOR
!! ------
!! J. Stein * Meteo-France *
!!
!! MODIFICATIONS
!! -------------
!! Original 05/01/95
!! Modifications Feb 14, 1995 (J.Cuxart) add the I/O arguments,
!! the director cosinus and change the names of the surface fluxes
!! Modifications March 21, 1995 (J.M.Carriere) take into account liquid
!! water
!! June 30,1995 (J.Stein) initialize at 0 the surf. fluxes
!! Modifications Sept. 1, 1995 (S.Belair) ISBA scheme
!! Modifications Sept.25, 1995 (J.Stein) switch on the radiation scheme
!! Modifications Sept. 11, 1995 (J.-P. Pinty) radiation scheme
!! Nov. 15, 1995 (J.Stein) cleaning + change the temporal
!! algorithm for the soil scheme-turbulence
!! Jan. 23, 1996 (J.Stein) add a new option for the surface
!! fluxes where Ts and z0 are given
!! March 18, 1996 (J.Stein) add the cloud fraction
!! March 28, 1996 (J.Stein) the soil scheme gives energy
!! fluxes + cleaning
!! June 17, 1996 (Lafore) statistics of computing time
!! August 4, 1996 (K. Suhre) add chemistry
!! Oct. 12, 1996 (J.Stein) use XSRCM in the turbulence
!! scheme
!! Nov. 18, 1996 (J.-P. Pinty) add domain translation
!! change arg. in radiations
!! Fev. 4, 1997 (J.Viviand) change isba's calling for ice
!! Jun. 22, 1997 (J.Stein) change the equation system and use
!! the absolute pressure
!! Jul. 09, 1997 (V.Masson) add directional z0
!! Jan. 24, 1998 (P.Bechtold) add convective transport for tracers
!! Jan. 24, 1998 (J.-P. Pinty) split SW and LW part for radiation
!! Mai. 10, 1999 (P.Bechtold) shallow convection
!! Oct. 20, 1999 (P.Jabouille) domain translation for turbulence
!! Jan. 04, 2000 (V.Masson) removes TSZ0 case
!! Jan. 04, 2000 (V.Masson) modifies albedo computation
! Jul 02, 2000 (F.Solmon/V.Masson) adaptation for patch approach
!! Nov. 15, 2000 (V.Masson) LES routines
!! Nov. 15, 2000 (V.Masson) effect of slopes on surface fluxes
!! Feb. 02, 2001 (P.Tulet) add friction velocities and aerodynamical
!! resistance (patch approach)
!! Jan. 04, 2000 (V.Masson) modify surf_rad_modif computation
!! Mar. 04, 2002 (F.Solmon) new interface for radiation call
!! Nov. 06, 2002 (V.Masson) LES budgets & budget time counters
!! Jan. 2004 (V.Masson) surface externalization
!! Jan. 13, 2004 (J.Escobar) bug correction : compute "GRAD" in parallel
!! Jan. 20, 2005 (P. Tulet) add dust sedimentation
!! Jan. 20, 2005 (P. Tulet) climatologic SSA
!! Jan. 20, 2005 (P. Tulet) add aerosol / dust scavenging
!! Jul. 2005 (N. Asencio) use the two-way result-fields
!! before ground_param call
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
Loading
Loading full blame...