2076ca15b468c7865a5caec8eb7566b69a171e78
[MNH-git_open_source-lfs.git] / src / MNH / default_desfmn.f90
1 !MNH_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
2 !MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
3 !MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
4 !MNH_LIC for details. version 1.
5 ! _Source: /home/cvsroot/MNH-VX-Y-Z/src/MNH/default_desfmn.f90,v _
6 !-----------------------------------------------------------------
7 !     ###########################
8       MODULE MODI_DEFAULT_DESFM_n
9 !     ###########################
10 !
11 INTERFACE
12 !
13 SUBROUTINE DEFAULT_DESFM_n(KMI)
14 INTEGER,         INTENT(IN)  :: KMI       ! Model index
15 END SUBROUTINE DEFAULT_DESFM_n
16 !
17 END INTERFACE
18 !
19 END MODULE MODI_DEFAULT_DESFM_n
20 !
21 !
22 !
23 !     ###############################
24       SUBROUTINE DEFAULT_DESFM_n(KMI)
25 !     ###############################
26 !
27 !!****  *DEFAULT_DESFM_n * - set default values for descriptive variables of
28 !!                         model KMI
29 !!
30 !!    PURPOSE
31 !!    -------
32 !       The purpose of this routine is to set default values for the variables
33 !     in descriptor files by filling the corresponding variables which
34 !     are stored in modules.
35 !
36 !
37 !!**  METHOD
38 !!    ------
39 !!      Each variable in modules, which can be initialized by  reading its
40 !!    value in the descriptor file is set to a default value.
41 !!     When this routine is used during INIT, the modules of the first model
42 !!   are used to temporarily store  the variables associated with a nested
43 !!   model.
44 !!     When this routine is used during  SPAWNING, the modules of a second
45 !!   model must be initialized.
46 !!     Default values for variables common to all models are set only
47 !!   at the first call of DEFAULT_DESFM_n (i.e. when KMI=1)
48 !!
49 !!
50 !!    EXTERNAL
51 !!    --------
52 !!     NONE
53 !!
54 !!    IMPLICIT ARGUMENTS
55 !!    ------------------
56 !!      Module MODD_PARAMETERS : JPHEXT,JPVEXT
57 !!
58 !!      Module MODD_CONF       : CCONF,L2D,L1D,LFLAT,NMODEL,NVERB
59 !!
60 !!      Module MODD_DYN        : XSEGLEN,XASSELIN,LCORIO,LNUMDIFF
61 !!                               XALKTOP,XALZBOT
62 !!
63 !!      Module MODD_FMOUT      : XOUT
64 !!
65 !!      Module MODD_NESTING    : NDAD(m),NDTRATIO(m),XWAY(m)
66 !!
67 !!      Module MODD_CONF_n    : LUSERV,LUSERC,LUSERR,LUSERI,LUSERS
68 !!                              LUSERG,LUSERH,CSEG,CEXP
69 !!
70 !!      Module MODD_LUNIT_n   : CINIFILE,CCPLFILE
71 !!
72 !!
73 !!      Module MODD_DYN_n     : XTSTEP,CPRESOPT,NITR,XRELAX,LHO_RELAX
74 !!         LVE_RELAX,XRIMKMAX,NRIMX,NRIMY
75 !!
76 !!      Module MODD_ADV_n : CUVW_ADV_SCHEME,CMET_ADV_SCHEME,CSV_ADV_SCHEME,NLITER
77 !!
78 !!      Module MODD_PARAM_n : CTURB,CRAD,CDCONV,CSCONV
79 !!
80 !!      Module MODD_LBC_n : CLBCX, CLBCY,NLBLX,NLBLY,XCPHASE,XCPHASE_PBL,XPOND
81 !!
82 !!      Module MODD_TURB_n : XIMPL,CTURBLEN,CTURBDIM,LTURB_FLX,LTURB_DIAG,LSUBG_COND
83 !!                           LTGT_FLX
84 !!
85 !!
86 !!      Module MODD_PARAM_RAD_n:
87 !!          XDTRAD,XDTRAD_CLONLY,LCLEAR_SKY,NRAD_COLNBR, NRAD_DIAG
88 !!
89 !!      Module MODD_BUDGET : CBUTYPE,NBUMOD,XBULEN,NBUKL, NBUKH,LBU_KCP,XBUWRI
90 !!         NBUIL, NBUIH,NBUJL, NBUJH,LBU_ICP,LBU_JCP,NBUMASK
91 !!         LBU_RU,LBU_RV,LBU_RW,LBU_RTH,LBU_RTKE,LBU_RRV,LBU_RRC,LBU_RRR
92 !!         LBU_RRI,LBU_RRS,LBU_RRG,LBU_RRH,LBU_RSVx
93 !!         NADVXU, NADVYU, NADVZU, NCURVU, NCORU, NDIFU, NRELU, NHTURBU,
94 !!         NVTURBU, NPRESU
95 !!         NADVXV, NADVYV, NADVZV, NCURVV, NCORV, NDIFV, NRELV, NHTURBV,
96 !!         NVTURBV, NPRESV
97 !!         NADVXW, NADVYW, NADVZW, NCURVW, NCORW, NGRAVW, NDIFW, NRELW,
98 !!         NHTURBW, NVTURBW, NPRESW
99 !!         NADVXTH, NADVYTH, NADVZTH, NPREFTH, NDIFTH, NRELTH, NHTURBTH,
100 !!         NVTURBTH, NREVATH, NCONDTH
101 !!         NADVXTKE, NADVYTKE, NADVZTKE, NDIFTKE, NDPTKE, NTPTKE, NDISSTKE,
102 !!         NTRTKE
103 !!         NADVXRV, NADVYRV, NADVZRV, NDIFRV, NRELRV, NHTURBRV, NVTURBRV,
104 !!         NREVARV, NCONDRV
105 !!         NADVXRC, NADVYRC, NADVZRC, NDIFRC, NHTURBRC, NVTURBRC, NACCRRC,
106 !!         NAUTORC, NCONDRC
107 !!         NADVXRR, NADVYRR, NADVZRR, NDIFRR, NACCRRR, NAUTORR, NREVARR,
108 !!         NSEDIRR
109 !!         NADVXRI, NADVYRI, NADVZRI, NDIFRI
110 !!         NADVXRS, NADVYRS, NADVZRS, NDIFRS
111 !!         NADVXRG, NADVYRG, NADVZRG, NDIFRG
112 !!         NADVXRH, NADVYRH, NADVZRH, NDIFRH
113 !!         NADVXSVx, NADVYSVx,  NADVZSVx, NDIFSVx, NHTURBSVx, NVTURBSVx
114 !!
115 !!      Module MODD_BLANK :
116 !!
117 !!          XDUMMYi, NDUMMYi, LDUMMYi, CDUMMYi
118 !!
119 !!      Module MODD_FRC :
120 !!
121 !!          LGEOST_UV_FRC,LGEOST_TH_FRC,LTEND_THRV_FRC
122 !!          LVERT_MOTION_FRC,LRELAX_THRV_FRC,LRELAX_UV_FRC,XRELAX_TIME_FRC
123 !!          XRELAX_HEIGHT_FRC,CRELAX_HEIGHT_TYPE,LTRANS,XUTRANS,XVTRANS,
124 !!          LPGROUND_FRC
125 !!
126 !!      Module MODD_PARAM_ICE :
127 !!
128 !!          LWARM,CPRISTINE_ICE
129 !!
130 !!      Module MODD_PARAM_KAFR_n :
131 !!
132 !!          XDTCONV,LREFRESH_ALL,LDOWN,NICE,LCHTRANS
133 !!
134 !!      Module MODD_PARAM_MFSHALL_n :
135 !!
136 !!         CMF_UPDRAFT,LMIXUV,CMF_CLOUD,XIMPL_MF,LMF_FLX              
137 !!
138 !!
139 !!
140 !!
141 !!    REFERENCE
142 !!    ---------
143 !!      Book2 of the documentation (routine DEFAULT_DESFM_n)
144 !!
145 !!
146 !!    AUTHOR
147 !!    ------
148 !!      V. Ducrocq       * Meteo France *
149 !!
150 !!    MODIFICATIONS
151 !!    -------------
152 !!      Original      02/06/94
153 !!      Modifications 17/10/94  (Stein)  For LCORIO
154 !!      Modifications 06/12/94  (Stein)  remove LBOUSS+add LABSLAYER, LNUMDIFF
155 !!                                       ,LSTEADYLS
156 !!      Modifications 06/12/94  (Stein)  remove LABSLAYER, add LHO_RELAX,
157 !!                                       LVE_RELAX, NRIMX, NRIMY, XRIMKMAX
158 !!      Modifications 09/01/95  (Lafore) add LSTEADY_DMASS
159 !!      Modifications 09/01/95  (Stein)  add the turbulence scheme namelist
160 !!      Modifications 09/01/95  (Stein)  add the 1D switch
161 !!      Modifications 10/03/95  (Mallet) add the coupling files
162 !!                    29/06/95  ( Stein, Nicolau, Hereil) add the budgets
163 !!      Modifications 25/09/95  ( Stein )add the LES tools
164 !!      Modifications 25/10/95  ( Stein )add the radiations
165 !!      Modifications 23/10/95  (Vila, lafore) new scalar advection scheme
166 !!      Modifications 24/02/96  (Stein)  change the default value for CCPLFILE
167 !!      Modifications 12/02/96  (Lafore) transformation to DEFAULT_DESFM_n for
168 !!                                       spawning
169 !!      Modifications 25/04/96  (Suhre)  add the blank module
170 !!      Modifications 29/07/96  (Pinty&Suhre) add module MODD_FRC
171 !!      Modifications 11/04/96  (Pinty)  add the rain-ice scheme and modify
172 !!                                       the splitted arrays in MODD_PARAM_RAD_n
173 !!      Modifications 11/01/97  (Pinty)  add the deep convection scheme
174 !!      Modifications 24/11/96  (Masson)  add LREFRESH_ALL in deep convection
175 !!      Modifications 12/02/96  (Lafore) transformation to DEFAULT_DESFM_n for spawning
176 !!      Modifications 22/07/96  (Lafore) gridnesting implementation
177 !!      Modifications 29/07/96  (Lafore) add the module MODD_FMOUT
178 !!      Modifications 23/06/97  (Stein)  add the equation system name
179 !!      Modifications 10/07/97  (Masson) add MODD_PARAM_GROUNDn : CROUGH
180 !!      Modifications 28/07/97  (Masson) remove LREFRESH_ALL and LSTEADY_DMASS
181 !!      Modifications 08/10/97  (Stein)  switch (_n=1) to initialize the
182 !!                                       parameters common to all models
183 !!      Modifications 24/01/98 (Bechtold) add LREFRESH_ALL, LCHTRANS,
184 !!                                         LTEND_THRV_FR and LSST_FRC
185 !!      Modifications 18/07/99  (Stein)  add LRAD_DIAG
186 !!      Modification  15/03/99 (Masson)  use of XUNDEF
187 !!      Modification  11/12/00 (Tomasini) Add CSEA_FLUX to MODD_PARAMn
188 !!      Modification  22/01/01 (Gazen) delete NSV and add LHORELAX_SVC2R2
189 !!                                     LHORELAX_SVCHEM,LHORELAX_SVLG
190 !!      Modification 15/03/02 (Solmon) radiation scheme: remove NSPOT and add
191 !!                                   default for aerosol and cloud rad. prop. control
192 !!      Modification 22/05/02 (Jabouille) put chimical default here
193 !!      Modification 01/2004  (Masson) removes surface (externalization)
194 !!                      09/04 (M. Tomasini) New namelist to modify the
195 !!                                             Cloud mixing length
196 !!                   07/05 (P.Tulet) New namelists for dust and aerosol
197 !!      Modification 01/2007  (Malardel, Pergaud) Add MODD_PARAM_MFSHALL_n
198 !!      Modification 10/2009  (Aumond) Add user multimasks for LES
199 !!      Modification 10/2009  (Aumond) Add MEAN_FIELD              
200 !!      Modification 12/04/07 (Leriche) add LUSECHAQ for aqueous chemistry
201 !!      Modification 30/05/07 (Leriche) add LCH_PH and XCH_PHINIT for pH
202 !!      Modification 25/04/08 (Leriche) add XRTMIN_AQ LWC threshold for aq. chemistry
203 !!                   16/07/10           add LHORELAX_SVIC
204 !!                   16/09/10           add LUSECHIC
205 !!                   13/01/11           add LCH_RET_ICE
206 !!                   01/07/11 (F.Couvreux) Add CONDSAMP
207 !!                   01/07/11 (B.Aouizerats) Add CAOP    
208 !!                   07/2013  (C.Lac) add WENO, LCHECK              
209 !!                   07/2013  (Bosseur & Filippi) adds Forefire
210 !!                   08/2015  (Redelsperger & Pianezze) add XPOND coefficient for LBC
211 !-------------------------------------------------------------------------------
212 !
213 !*       0.    DECLARATIONS
214 !              ------------
215 USE MODD_PARAMETERS
216 USE MODD_CONF             !        For INIT only DEFAULT_DESFM1
217 USE MODD_CONFZ
218 USE MODD_DYN
219 USE MODD_NESTING
220 USE MODD_FMOUT
221 USE MODD_SERIES
222 USE MODD_CONF_n           ! modules used to set the default values is only
223 USE MODD_LUNIT_n          ! the one corresponding to model 1. These memory
224 USE MODD_DIM_n            ! addresses will then be filled by the values read in
225 USE MODD_DYN_n            ! the DESFM corresponding to model n which may have
226 USE MODD_ADV_n            ! missing values. This is why we affect default values.
227 USE MODD_PARAM_n          !      For SPAWNING DEFAULT_DESFM2 is also used
228 USE MODD_LBC_n
229 USE MODD_OUT_n
230 USE MODD_TURB_n
231 USE MODD_BUDGET
232 USE MODD_LES
233 USE MODD_PARAM_RAD_n
234 USE MODD_BLANK
235 USE MODD_FRC
236 USE MODD_PARAM_ICE
237 USE MODD_PARAM_C2R2
238 USE MODD_TURB_CLOUD
239 USE MODD_PARAM_KAFR_n
240 USE MODD_PARAM_MFSHALL_n
241 USE MODD_CH_MNHC_n
242 USE MODD_SERIES_n
243 USE MODD_NUDGING_n
244 USE MODD_CH_AEROSOL
245 USE MODD_DUST
246 USE MODD_SALT
247 USE MODD_PASPOL
248 USE MODD_CONDSAMP
249 USE MODD_MEAN_FIELD
250 USE MODD_DRAGTREE
251 USE MODD_LATZ_EDFLX
252 USE MODD_2D_FRC
253 #ifdef MNH_FOREFIRE
254 USE MODD_FOREFIRE
255 #endif
256 !
257 IMPLICIT NONE
258 !
259 !*       0.1   declarations of arguments
260 !
261 INTEGER,         INTENT(IN)  :: KMI       ! Model index
262 !
263 !*       0.2   declaration of local variables
264 !
265 INTEGER             :: JM      ! loop index
266 !
267 !-------------------------------------------------------------------------------
268 !
269 !*      1.    SET DEFAULT VALUES FOR MODD_LUNIT_n :
270 !             ----------------------------------
271 !
272 CINIFILE='INIFILE'
273 CINIFILEPGD=''
274 CCPLFILE(:)='    '
275 !
276 !-------------------------------------------------------------------------------
277 !
278 !*      2.    SET DEFAULT VALUES FOR MODD_CONF AND MODD_CONF_n :
279 !             ------------------------------------------------
280 !
281 IF (KMI == 1) THEN
282   CCONF      ='START'
283   LTHINSHELL = .FALSE.
284   L2D        = .FALSE.
285   L1D        = .FALSE.
286   LFLAT      = .FALSE.
287   NMODEL     = 1
288   CEQNSYS    = 'DUR'
289   NVERB      = 5
290   CEXP       = 'EXP01'
291   CSEG       = 'SEG01'
292   LFORCING   = .FALSE.
293   L2D_ADV_FRC= .FALSE.
294   L2D_REL_FRC= .FALSE. 
295   XRELAX_HEIGHT_BOT = 0.
296   XRELAX_HEIGHT_TOP = 30000.
297   XRELAX_TIME = 864000.
298   LPACK      = .TRUE.
299   NHALO      = 3
300 #ifdef MNH_SX5
301   CSPLIT     ='YSPLITTING'   ! NEC vectoriel architecture , low number of PROC
302 #else
303   CSPLIT     ='BSPLITTING'   ! Scalaire architecture , high number of PROC
304 #endif
305   NZ_PROC    = 0               !JUAN Z_SPLITTING :: number of proc in Z splitting
306   NZ_SPLITTING = 10            !JUAN Z_SPLITTING :: for debug NZ=1=flat_inv;  NZ=10=flat_invz; NZ=1+2 the two
307   LLG        = .FALSE.
308   LINIT_LG   = .FALSE.
309   CINIT_LG   = 'FMOUT'
310   LNOMIXLG   = .FALSE.
311   LCHECK     = .FALSE.
312 END IF
313 !
314 CCLOUD    = 'NONE'
315 LUSERV    = .TRUE.
316 LUSERC    = .FALSE.
317 LUSERR    = .FALSE.
318 LUSERI    = .FALSE.
319 LUSERS    = .FALSE.
320 LUSERG    = .FALSE.
321 LUSERH    = .FALSE.
322 !NSV      = 0
323 !NSV_USER = 0
324 LUSECI    = .FALSE.
325 !
326 !-------------------------------------------------------------------------------
327 !
328 !*      3.    SET DEFAULT VALUES FOR MODD_DYN AND MODD_DYN_n :
329 !             -----------------------------------------------
330 !
331 IF (KMI == 1) THEN
332   XSEGLEN   = 43200.
333   XASSELIN  = 0.2
334   XASSELIN_SV = 0.02
335   LCORIO    = .TRUE.
336   LNUMDIFU  = .FALSE.
337   LNUMDIFTH = .FALSE.
338   LNUMDIFSV = .FALSE.
339   XALZBOT   = 4000.
340   XALKTOP   = 0.01
341   XALKGRD   = 0.01
342   XALZBAS   = 0.01
343 END IF
344 !
345 XTSTEP    = 60.
346 CPRESOPT  = 'CRESI'
347 NITR      = 4
348 LITRADJ   = .TRUE.
349 LRES      = .FALSE.
350 XRES      = 1.E-07
351 XRELAX    = 1.
352 LVE_RELAX = .FALSE.
353 LVE_RELAX_GRD = .FALSE.
354 XRIMKMAX  = 0.01 / XTSTEP
355 XT4DIFU   = 1800.          
356 XT4DIFTH  = 1800.           
357 XT4DIFSV  = 1800.           
358 !
359 IF (KMI == 1) THEN   ! for model 1 we have a Large scale information
360   NRIMX = JPRIMMAX   ! for U,V,W,TH,Rv used for the hor. relaxation
361   NRIMY = JPRIMMAX
362 ELSE
363   NRIMX = 0          ! for inner models we use only surfacic fields to
364   NRIMY = 0          ! give the lbc and no hor. relaxation is used
365 END IF
366 !
367 LHORELAX_UVWTH = .FALSE.
368 LHORELAX_RV = .FALSE.
369 LHORELAX_RC = .FALSE. ! for all these fields, no large scale is usally available
370 LHORELAX_RR = .FALSE. ! for model 1 and for inner models, we only use surfacic
371 LHORELAX_RS = .FALSE. ! fiels ( no hor. relax. )
372 LHORELAX_RI = .FALSE.
373 LHORELAX_RG = .FALSE.
374 LHORELAX_RH = .FALSE.
375 LHORELAX_TKE = .FALSE.
376 LHORELAX_SV(:) = .FALSE.
377 LHORELAX_SVC2R2 = .FALSE.
378 LHORELAX_SVC1R3 = .FALSE.
379 LHORELAX_SVELEC = .FALSE.
380 LHORELAX_SVLG   = .FALSE.
381 LHORELAX_SVCHEM = .FALSE.
382 LHORELAX_SVCHIC = .FALSE.
383 LHORELAX_SVDST  = .FALSE.
384 LHORELAX_SVSLT  = .FALSE.
385 LHORELAX_SVPP   = .FALSE.
386 LHORELAX_SVCS   = .FALSE.
387 LHORELAX_SVAER  = .FALSE.
388 #ifdef MNH_FOREFIRE
389 LHORELAX_SVFF   = .FALSE.
390 #endif
391 !
392 !-------------------------------------------------------------------------------
393 !
394 !*      4.    SET DEFAULT VALUES FOR MODD_NESTING :
395 !             -----------------------------------
396 !
397 IF (KMI == 1) THEN
398   NDAD(1)=1
399   DO JM=2,JPMODELMAX
400     NDAD(JM)  = JM - 1
401   END DO
402   NDTRATIO(:) = 1
403   XWAY(:)     = 2.      ! two-way interactive gridnesting
404   XWAY(1)     = 0.      ! except for model 1
405 END IF
406 !
407 !-------------------------------------------------------------------------------
408 !
409 !*      5.    SET DEFAULT VALUES FOR MODD_ADV_n :
410 !             ----------------------------------
411 !
412 CUVW_ADV_SCHEME =  'WENO_K'
413 CMET_ADV_SCHEME =  'PPM_01'
414 CSV_ADV_SCHEME  =  'PPM_01'
415 CTEMP_SCHEME    =  'RK53'        
416 NWENO_ORDER = 3
417 NSPLIT      = 1
418 LSPLIT_CFL  = .TRUE.
419 XSPLIT_CFL  = 0.8
420 LCFL_WRIT   = .FALSE.
421 !
422 !-------------------------------------------------------------------------------
423 !
424 !*      6.    SET DEFAULT VALUES FOR MODD_PARAM_n :
425 !             -----------------------------------
426 !
427 CTURB   = 'NONE'
428 CRAD    = 'NONE'
429 CDCONV  = 'NONE'
430 CSCONV  = 'NONE'
431 CELEC   = 'NONE'
432 CACTCCN = 'NONE'
433 !
434 !-------------------------------------------------------------------------------
435 !
436 !*      7.    SET DEFAULT VALUES FOR MODD_LBC_n :
437 !             ---------------------------------
438 !
439 CLBCX(1) ='CYCL'
440 CLBCX(2) ='CYCL'
441 CLBCY(1) ='CYCL'
442 CLBCY(2) ='CYCL'
443 NLBLX(:) = 1
444 NLBLY(:) = 1
445 XCPHASE = 20.
446 XCPHASE_PBL = 0.
447 XCARPKMAX = XUNDEF
448 XPOND = 1.0
449 !
450 !-------------------------------------------------------------------------------
451 !
452 !*      8.    SET DEFAULT VALUES FOR MODD_NUDGING_n :
453 !             ---------------------------------
454 !
455 LNUDGING = .FALSE.
456 XTNUDGING = 21600.
457 !
458 !-------------------------------------------------------------------------------
459 !
460 !*      9.    SET DEFAULT VALUES FOR MODD_FMOUT and MODD_OUT_n :
461 !             ------------------------------------------------
462 !
463 IF (KMI == 1) XFMOUT (:,:) = XUNDEF
464 !
465 !
466 !-------------------------------------------------------------------------------
467 !
468 !*      10.   SET DEFAULT VALUES FOR MODD_TURB_n :
469 !             ----------------------------------
470 !
471 XIMPL     = 1.
472 XKEMIN    = 0.01
473 CTURBLEN  = 'BL89'
474 CTURBDIM  = '1DIM'
475 LTURB_FLX =.FALSE.
476 LTURB_DIAG=.FALSE.
477 LSUBG_COND=.FALSE.
478 CSUBG_AUCV='NONE' 
479 LSIGMAS   =.TRUE.
480 LSIG_CONV =.FALSE.
481 LRMC01    =.FALSE.
482 CTOM      ='NONE'
483 VSIGQSAT  = 0.
484 !-------------------------------------------------------------------------------
485 !
486 !*      10b.   SET DEFAULT VALUES FOR MODD_DRAGTREE :
487 !             ----------------------------------
488 !
489 LDRAGTREE = .FALSE.
490 !
491 !-------------------------------------------------------------------------------
492 !
493 !*      11.   SET DEFAULT VALUES FOR MODD_BUDGET :
494 !             ------------------------------------
495 !
496 !       11.1 General  budget variables
497 !
498 IF (KMI == 1) THEN
499   CBUTYPE = 'NONE'
500   NBUMOD = 1
501   XBULEN = XSEGLEN
502   XBUWRI = XSEGLEN
503   NBUKL  = 1
504   NBUKH  = 0        
505   LBU_KCP = .TRUE.
506 !
507 !       11.2 Variables for the cartesian box
508 !
509   NBUIL = 1
510   NBUIH = 0      
511   NBUJL = 1
512   NBUJH = 0          
513   LBU_ICP = .TRUE.
514   LBU_JCP = .TRUE.
515 !
516 !       11.3 Variables for the mask
517 !
518   NBUMASK = 1
519 !
520 !       11.4 Variables for budget and processes choice
521   LBU_RU = .FALSE.
522   NASSEU  = 0
523   NNESTU  = 0
524   NADVU   = 0
525   NFRCU   = 0
526   NNUDU   = 0
527   NCURVU  = 0
528   NCORU   = 0
529   NDIFU   = 0
530   NRELU   = 0
531   NMAFLU  = 0
532   NDRAGU  = 0
533   NVTURBU = 0
534   NHTURBU = 0
535   NPRESU  = 0
536 !
537 !                    Budget of RV
538   LBU_RV = .FALSE.
539   NASSEV  = 0
540   NNESTV  = 0
541   NADVV   = 0
542   NFRCV   = 0
543   NNUDV   = 0
544   NCURVV  = 0
545   NCORV   = 0
546   NDIFV   = 0
547   NRELV   = 0
548   NMAFLV  = 0
549   NDRAGV  = 0
550   NVTURBV = 0
551   NHTURBV = 0
552   NPRESV  = 0
553 !
554 !                    Budget of RW
555   LBU_RW = .FALSE.
556   NASSEW  = 0
557   NNESTW  = 0
558   NADVW   = 0
559   NFRCW   = 0
560   NNUDW   = 0
561   NCURVW  = 0
562   NCORW   = 0
563   NGRAVW  = 0
564   NDIFW   = 0
565   NRELW   = 0
566   NVTURBW = 0
567   NHTURBW = 0
568   NPRESW  = 0
569 !
570 !                    Budget of RTH
571   LBU_RTH = .FALSE.
572   NASSETH  = 0
573   NNESTTH  = 0
574   NADVTH   = 0
575   NFRCTH   = 0
576   NNUDTH   = 0
577   NPREFTH  = 0
578   NDIFTH   = 0
579   NRELTH   = 0
580   NRADTH   = 0
581   NDCONVTH = 0
582   NMAFLTH  = 0
583   NVTURBTH = 0
584   NHTURBTH = 0
585   NDISSHTH = 0
586   NNEGATH  = 0
587   NREVATH  = 0
588   NCONDTH  = 0
589   NHENUTH  = 0
590   NHONTH   = 0
591   NSFRTH   = 0
592   NDEPSTH  = 0
593   NDEPGTH  = 0
594   NREVATH  = 0
595   NRIMTH   = 0
596   NACCTH   = 0
597   NCFRZTH  = 0
598   NWETGTH  = 0
599   NDRYGTH  = 0
600   NGMLTTH  = 0
601   NIMLTTH  = 0
602   NBERFITH = 0
603   NCDEPITH = 0
604   NWETHTH = 0
605   NHMLTTH = 0
606 !
607 !                    Budget of RTKE
608   LBU_RTKE = .FALSE.
609   NASSETKE = 0
610   NADVTKE  = 0
611   NFRCTKE  = 0
612   NDIFTKE  = 0
613   NRELTKE  = 0
614   NDRAGTKE = 0
615   NDPTKE   = 0
616   NTPTKE   = 0
617   NDISSTKE = 0
618   NTRTKE   = 0
619 !
620 !                    Budget of RRV
621   LBU_RRV = .FALSE.
622   NASSERV  = 0
623   NNESTRV  = 0
624   NADVRV   = 0
625   NFRCRV   = 0
626   NNUDRV   = 0
627   NDIFRV   = 0
628   NRELRV   = 0
629   NDCONVRV = 0
630   NMAFLRV  = 0
631   NVTURBRV = 0
632   NHTURBRV = 0
633   NNEGARV  = 0
634   NREVARV  = 0
635   NCONDRV  = 0
636   NHENURV  = 0
637   NDEPSRV  = 0
638   NDEPGRV  = 0
639   NREVARV  = 0
640   NCDEPIRV = 0
641 !
642 !                    Budget of RRC
643   LBU_RRC = .FALSE.
644   NASSERC  = 0
645   NNESTRC  = 0
646   NADVRC   = 0
647   NFRCRC   = 0
648   NDIFRC   = 0
649   NRELRC   = 0
650   NDCONVRC = 0
651   NVTURBRC = 0
652   NHTURBRC = 0
653   NNEGARC  = 0
654   NACCRRC  = 0
655   NAUTORC  = 0
656   NCONDRC  = 0
657   NAUTORC  = 0
658   NACCRRC  = 0
659   NHONRC   = 0
660   NRIMRC   = 0
661   NWETGRC  = 0
662   NDRYGRC  = 0
663   NIMLTRC  = 0
664   NBERFIRC = 0
665   NCDEPIRC = 0
666   NWETHRC = 0
667 !
668 !                    Budget of RRR
669   LBU_RRR = .FALSE.
670   NASSERR  = 0
671   NNESTRR  = 0
672   NADVRR   = 0
673   NFRCRR   = 0
674   NDIFRR   = 0
675   NRELRR   = 0
676   NNEGARR  = 0
677   NACCRRR  = 0
678   NAUTORR  = 0
679   NREVARR  = 0
680   NSEDIRR  = 0
681   NSFRRR   = 0
682   NACCRR   = 0
683   NCFRZRR  = 0
684   NWETGRR  = 0
685   NDRYGRR  = 0
686   NGMLTRR  = 0
687   NWETHRR  = 0
688   NHMLTRR  = 0
689 !
690 !                    Budget of RRI
691   LBU_RRI = .FALSE.
692   NASSERI  = 0
693   NNESTRI  = 0
694   NADVRI   = 0
695   NFRCRI   = 0
696   NDIFRI   = 0
697   NRELRI   = 0
698   NDCONVRI = 0
699   NVTURBRI = 0
700   NHTURBRI = 0
701   NNEGARI  = 0
702   NSEDIRI  = 0
703   NHENURI  = 0
704   NHONRI   = 0
705   NAGGSRI  = 0
706   NAUTSRI  = 0
707   NCFRZRI  = 0
708   NWETGRI  = 0
709   NDRYGRI  = 0
710   NIMLTRI  = 0
711   NBERFIRI = 0
712   NCDEPIRI = 0
713   NWETHRI = 0
714 !
715 !                    Budget of RRS
716   LBU_RRS = .FALSE.
717   NASSERS  = 0
718   NNESTRS  = 0
719   NADVRS   = 0
720   NFRCRS   = 0
721   NDIFRS   = 0
722   NRELRS   = 0
723   NNEGARS  = 0
724   NSEDIRS  = 0
725   NDEPSRS  = 0
726   NAGGSRS  = 0
727   NAUTSRS  = 0
728   NRIMRS   = 0
729   NACCRS   = 0
730   NCMELRS  = 0
731   NWETGRS  = 0
732   NDRYGRS  = 0
733   NWETHRS  = 0
734 !
735 !                    Budget of RRG
736   LBU_RRG = .FALSE.
737   NASSERG  = 0
738   NNESTRG  = 0
739   NADVRG   = 0
740   NFRCRG   = 0
741   NDIFRG   = 0
742   NRELRG   = 0
743   NNEGARG  = 0
744   NSEDIRG  = 0
745   NSFRRG   = 0
746   NDEPGRG  = 0
747   NRIMRG   = 0
748   NACCRG   = 0
749   NCMELRG  = 0
750   NCFRZRG  = 0
751   NWETGRG  = 0
752   NDRYGRG  = 0
753   NGMLTRG  = 0
754   NWETHRG  = 0
755   NCOHGRG  = 0
756 !
757 !                    Budget of RRH
758   LBU_RRH = .FALSE.
759   NASSERH  = 0
760   NNESTRH  = 0
761   NADVRH   = 0
762   NFRCRH   = 0
763   NDIFRH   = 0
764   NRELRH   = 0
765   NNEGARH  = 0
766   NWETGRH  = 0
767   NWETHRH  = 0
768   NCOHGRH  = 0
769   NHMLTRH  = 0
770 !
771 !                    Budget of RSVx
772   LBU_RSV = .FALSE.
773   NASSESV  = 0
774   NNESTSV  = 0
775   NADVSV   = 0
776   NFRCSV   = 0
777   NDIFSV   = 0
778   NRELSV   = 0
779   NDCONVSV = 0
780   NMAFLSV  = 0
781   NVTURBSV = 0
782   NHTURBSV = 0
783   NCHEMSV  = 0
784   NNEGASV  = 0
785 !
786 !
787 END IF
788 !
789 !-------------------------------------------------------------------------------
790 !
791 !*      12.    SET DEFAULT VALUES FOR MODD_LES :
792 !             ---------------------------------
793 !
794 IF (KMI == 1) THEN
795   LLES_MEAN               = .FALSE.
796   LLES_RESOLVED           = .FALSE.
797   LLES_SUBGRID            = .FALSE.
798   LLES_UPDRAFT            = .FALSE.
799   LLES_DOWNDRAFT          = .FALSE.
800   LLES_SPECTRA            = .FALSE.
801 !
802   NLES_LEVELS             = NUNDEF
803   XLES_ALTITUDES          = XUNDEF
804   NSPECTRA_LEVELS         = NUNDEF
805   XSPECTRA_ALTITUDES      = XUNDEF
806   NLES_TEMP_SERIE_I       = NUNDEF
807   NLES_TEMP_SERIE_J       = NUNDEF
808   NLES_TEMP_SERIE_Z       = NUNDEF
809   CLES_NORM_TYPE          = 'NONE'
810   CBL_HEIGHT_DEF          = 'KE'
811   XLES_TEMP_SAMPLING      = XUNDEF
812   XLES_TEMP_MEAN_START    = XUNDEF
813   XLES_TEMP_MEAN_END      = XUNDEF
814   XLES_TEMP_MEAN_STEP     = 3600.
815   LLES_CART_MASK          = .FALSE.
816   NLES_IINF               = NUNDEF
817   NLES_ISUP               = NUNDEF
818   NLES_JINF               = NUNDEF
819   NLES_JSUP               = NUNDEF
820   LLES_NEB_MASK           = .FALSE.
821   LLES_CORE_MASK          = .FALSE.
822   LLES_MY_MASK            = .FALSE.
823   NLES_MASKS_USER         = NUNDEF
824   LLES_CS_MASK            = .FALSE.
825
826   LLES_PDF               = .FALSE.
827   NPDF                   = 1
828   XTH_PDF_MIN            = 270.
829   XTH_PDF_MAX            = 350.
830   XW_PDF_MIN             = -10.
831   XW_PDF_MAX             = 10.
832   XTHV_PDF_MIN           = 270.
833   XTHV_PDF_MAX           = 350.
834   XRV_PDF_MIN            = 0.
835   XRV_PDF_MAX            = 20.
836   XRC_PDF_MIN            = 0.
837   XRC_PDF_MAX            = 1.
838   XRR_PDF_MIN            = 0.
839   XRR_PDF_MAX            = 1.
840   XRI_PDF_MIN            = 0.
841   XRI_PDF_MAX            = 1.
842   XRS_PDF_MIN            = 0.
843   XRS_PDF_MAX            = 1.
844   XRG_PDF_MIN            = 0.
845   XRG_PDF_MAX            = 1.
846   XRT_PDF_MIN            = 0.
847   XRT_PDF_MAX            = 20.
848   XTHL_PDF_MIN           = 270.
849   XTHL_PDF_MAX           = 350.
850 END IF
851 !
852 !-------------------------------------------------------------------------------
853 !
854 !*      13.   SET DEFAULT VALUES FOR MODD_PARAM_RAD_n :
855 !             ---------------------------------------
856 !
857 XDTRAD        = XTSTEP
858 XDTRAD_CLONLY = XTSTEP
859 LCLEAR_SKY    =.FALSE.
860 NRAD_COLNBR   = 1000
861 NRAD_DIAG     = 0
862 CLW ='RRTM'
863 CAER='SURF'
864 CAOP='CLIM'
865 CEFRADL='MART'
866 CEFRADI='LIOU'
867 COPWSW = 'FOUQ'
868 COPISW = 'EBCU'
869 COPWLW = 'SMSH'
870 COPILW = 'EBCU'
871 XFUDG = 1.
872 LAERO_FT=.FALSE.
873 LFIX_DAT=.FALSE.
874 !
875 !-------------------------------------------------------------------------------
876 !
877 !*      14.   SET DEFAULT VALUES FOR MODD_BLANK :
878 !             -----------------------------------
879 !
880 IF (KMI == 1) THEN
881   XDUMMY1       = 0.
882   XDUMMY2       = 0.
883   XDUMMY3       = 0.
884   XDUMMY4       = 0.
885   XDUMMY5       = 0.
886   XDUMMY6       = 0.
887   XDUMMY7       = 0.
888   XDUMMY8       = 0.
889   XDUMMY=0.
890 !
891   NDUMMY1       = 0
892   NDUMMY2       = 0
893   NDUMMY3       = 0
894   NDUMMY4       = 0
895   NDUMMY5       = 0
896   NDUMMY6       = 0
897   NDUMMY7       = 0
898   NDUMMY8       = 0
899   NDUMMY=0
900 !
901   LDUMMY1       = .TRUE.
902   LDUMMY2       = .TRUE.
903   LDUMMY3       = .TRUE.
904   LDUMMY4       = .TRUE.
905   LDUMMY5       = .TRUE.
906   LDUMMY6       = .TRUE.
907   LDUMMY7       = .TRUE.
908   LDUMMY8       = .TRUE.
909   LDUMMY=.TRUE.
910 !
911   CDUMMY1       = ' '
912   CDUMMY2       = ' '
913   CDUMMY3       = ' '
914   CDUMMY4       = ' '
915   CDUMMY5       = ' '
916   CDUMMY6       = ' '
917   CDUMMY7       = ' '
918   CDUMMY8       = ' '
919   CDUMMY= ' '
920 END IF
921 !
922 !------------------------------------------------------------------------------
923 !
924 !*      15.   SET DEFAULT VALUES FOR MODD_FRC :
925 !             ---------------------------------
926 !
927 IF (KMI == 1) THEN
928   LGEOST_UV_FRC      = .FALSE.
929   LGEOST_TH_FRC      = .FALSE.
930   LTEND_THRV_FRC      = .FALSE.
931   LVERT_MOTION_FRC   = .FALSE.
932   LRELAX_THRV_FRC    = .FALSE.
933   LRELAX_UV_FRC      = .FALSE.
934   XRELAX_TIME_FRC    = 10800.
935   XRELAX_HEIGHT_FRC  = 0.
936   CRELAX_HEIGHT_TYPE = "FIXE"
937   LTRANS             = .FALSE.
938   XUTRANS            = 0.0
939   XVTRANS            = 0.0
940   LPGROUND_FRC       = .FALSE.
941 END IF
942 !
943 !-------------------------------------------------------------------------------
944 !
945 !
946 !*      16.   SET DEFAULT VALUES FOR MODD_PARAM_ICE :
947 !             ---------------------------------------
948 !
949 IF (KMI == 1) THEN
950   LWARM = .TRUE.
951   CPRISTINE_ICE = 'PLAT'
952   LSEDIC  = .FALSE.
953   LCONVHG = .FALSE.
954   CSEDIM  = 'SPLI'
955 END IF
956 !
957 !-------------------------------------------------------------------------------
958 !
959 !
960 !*      17.   SET DEFAULT VALUES FOR MODD_PARAM_KAFR_n :
961 !             --------------------------------------------
962 !
963 XDTCONV       = MAX( 300.0,XTSTEP )
964 NICE          = 1
965 LREFRESH_ALL  = .TRUE.
966 LCHTRANS      = .FALSE.
967 LDOWN         = .TRUE.
968 LSETTADJ      = .FALSE.
969 XTADJD        = 3600.
970 XTADJS        = 10800.
971 LDIAGCONV     = .FALSE.
972 NENSM         = 0
973 !
974 !-------------------------------------------------------------------------------
975 !
976 !
977 !*      18.   SET DEFAULT VALUES FOR MODD_PARAM_MFSHALL_n :
978 !             --------------------------------------------
979 !
980 XIMPL_MF    = 1.
981 CMF_UPDRAFT = 'EDKF'
982 CMF_CLOUD   = 'DIRE'
983 LMIXUV      = .TRUE. 
984 LMF_FLX     = .FALSE.
985 !
986 XALP_PERT   = 0.3  
987 XABUO       = 1.     
988 XBENTR      = 1.   
989 XBDETR      = 0.   
990 XCMF        = 0.065 
991 XENTR_MF    = 0.035
992 XCRAD_MF    = 50.  
993 XENTR_DRY   = 0.55 
994 XDETR_DRY   = 10.  
995 XDETR_LUP   = 1.  
996 XKCF_MF     = 2.75
997 XKRC_MF     = 1.   
998 XTAUSIGMF   = 600.  
999 XPRES_UV    = 0.5  
1000 XFRAC_UP_MAX= 0.33
1001 XALPHA_MF = 2.     
1002 XSIGMA_MF = 20.  
1003 !
1004 XA1    =  2./3.  
1005 XB     =  0.002       
1006 XC     =  0.012     
1007 XBETA1 =  0.9         
1008 !
1009 !-------------------------------------------------------------------------------
1010 !
1011 !*      19.   SET DEFAULT VALUES FOR MODD_PARAM_C2R2 :
1012 !             ----------------------------------------
1013 !
1014 IF (KMI == 1) THEN
1015   XNUC    = 1.0
1016   XALPHAC = 3.0
1017   XNUR    = 2.0
1018   XALPHAR = 1.0
1019 !
1020   LRAIN = .TRUE.
1021   LSEDC = .TRUE.
1022   LACTIT = .FALSE.
1023   LSUPSAT= .FALSE.
1024 !
1025   HPARAM_CCN = 'XXX'
1026   HINI_CCN   = 'XXX'
1027   HTYPE_CCN  = 'X'
1028 !
1029   XCHEN      = 0.0
1030   XKHEN      = 0.0
1031   XMUHEN     = 0.0
1032   XBETAHEN   = 0.0
1033 !
1034   XCONC_CCN   = 0.0
1035   XAERDIFF    = 0.0
1036   XAERHEIGHT  = 2000
1037   XR_MEAN_CCN = 0.0
1038   XLOGSIG_CCN = 0.0
1039   XFSOLUB_CCN = 1.0
1040   XACTEMP_CCN = 280.
1041 ENDIF
1042 !
1043 !-------------------------------------------------------------------------------
1044 !
1045 !*      20.   SET DEFAULT VALUES FOR MODD_CH_MNHC_n
1046 !             -------------------------------------
1047 !
1048 LUSECHEM            = .FALSE.
1049 LUSECHAQ            = .FALSE.
1050 LUSECHIC            = .FALSE.
1051 LCH_INIT_FIELD      = .FALSE.
1052 LCH_SURFACE_FLUX    = .FALSE.
1053 LCH_CONV_SCAV       = .FALSE.
1054 LCH_CONV_LINOX      = .FALSE.
1055 LCH_PH              = .FALSE.
1056 LCH_RET_ICE         = .FALSE.
1057 XCH_PHINIT          = 5.2
1058 XRTMIN_AQ           = 5.e-8
1059 CCHEM_INPUT_FILE    = 'EXSEG1.nam'
1060 CCH_TDISCRETIZATION = 'SPLIT'
1061 NCH_SUBSTEPS        = 1
1062 LCH_TUV_ONLINE      = .FALSE.
1063 CCH_TUV_LOOKUP      = 'PHOTO.TUV39'
1064 CCH_TUV_CLOUDS      = 'NONE'
1065 XCH_TUV_ALBNEW      = -1.
1066 XCH_TUV_DOBNEW      = -1.
1067 XCH_TUV_TUPDATE     = 600.
1068 CCH_VEC_METHOD      = 'MAX'
1069 NCH_VEC_LENGTH      = 1000
1070 XCH_TS1D_TSTEP      = 600.
1071 CCH_TS1D_COMMENT    = 'no comment'
1072 CCH_TS1D_FILENAME   = 'IO1D'
1073 !
1074 !-------------------------------------------------------------------------------
1075 !
1076 !*      21.   SET DEFAULT VALUES FOR MODD_SERIES AND MODD_SERIE_n
1077 !             ---------------------------------------------------
1078 !
1079 IF (KMI == 1) THEN
1080   LSERIES      = .FALSE.
1081   LMASKLANDSEA = .FALSE.
1082   LWMINMAX     = .FALSE.
1083   LSURF        = .FALSE.
1084 ENDIF
1085 !
1086 NIBOXL = 1 !+ JPHEXT
1087 NIBOXH = 1 !+ 2*JPHEXT
1088 NJBOXL = 1 !+ JPHEXT
1089 NJBOXH = 1 !+ 2*JPHEXT
1090 NKCLS  = 1 !+ JPVEXT
1091 NKLOW  = 1 !+ JPVEXT
1092 NKMID  = 1 !+ JPVEXT
1093 NKUP   = 1 !+ JPVEXT
1094 NKCLA  = 1 !+ JPVEXT
1095 NBJSLICE = 1
1096 NJSLICEL(:) = 1 !+ JPHEXT
1097 NJSLICEH(:) = 1 !+ 2*JPHEXT
1098 NFREQSERIES  = INT(XSEGLEN /(100.*XTSTEP) )
1099 NFREQSERIES  = MAX(NFREQSERIES,1)
1100 !
1101 !-------------------------------------------------------------------------------
1102 !
1103 !*      22.   SET DEFAULT VALUES FOR MODD_TURB_CLOUD
1104 !             --------------------------------------
1105 !
1106 IF (KMI == 1) THEN
1107   NMODEL_CLOUD = NUNDEF
1108   CTURBLEN_CLOUD = 'DELT'
1109   XCOEF_AMPL_SAT = 5.
1110   XCEI_MIN = 0.001E-06
1111   XCEI_MAX = 0.01E-06
1112 ENDIF
1113 !-------------------------------------------------------------------------------
1114 !
1115 !*      22.   SET DEFAULT VALUES FOR MODD_MEAN_FIELD
1116 !             --------------------------------------
1117 !
1118 IF (KMI == 1) THEN
1119   LMEAN_FIELD = .FALSE.
1120 ENDIF
1121 !
1122 !-------------------------------------------------------------------------------
1123 !
1124 !*      22.   SET DEFAULT VALUES FOR MODD_AEROSOL
1125 !             -----------------------------------
1126 IF (KMI == 1) THEN ! other values are defined in modd_ch_aerosol
1127 !
1128 ! aerosol lognormal parameterization
1129
1130 LVARSIGI  = .FALSE.   ! switch to active pronostic dispersion for I mode
1131 LVARSIGJ  = .FALSE.   ! switch to active pronostic dispersion for J mode
1132 LHETEROSO4 = .FALSE.  ! switch to active sulfates heteronegeous
1133                       ! production
1134 LSEDIMAERO = .FALSE.  ! switch to active aerosol sedimentation
1135 LAERINIT   = .FALSE.  ! switch to initialize aerosol in arome
1136 CMINERAL      = "NONE"   ! mineral equilibrium scheme
1137 CORGANIC      = "NONE"   ! mineral equilibrium scheme
1138 CNUCLEATION   = "NONE" ! sulfates nucleation scheme
1139 LDEPOS_AER(:) = .FALSE.
1140
1141 ENDIF
1142
1143 !*      23.   SET DEFAULT VALUES FOR MODD_DUST and MODD_SALT
1144 !             ----------------------------------------------
1145 !
1146 IF (KMI == 1) THEN ! other values initialized in modd_dust
1147   LDUST      = .FALSE.
1148   NMODE_DST  = 3
1149   LVARSIG    = .FALSE.
1150   LSEDIMDUST = .FALSE.
1151   LDEPOS_DST(:) = .FALSE.
1152
1153   LSALT      = .FALSE.
1154   NMODE_DST  = 3
1155   LVARSIG_SLT= .FALSE.
1156   LSEDIMSALT = .FALSE.
1157   LDEPOS_SLT(:)     = .FALSE.
1158 ENDIF
1159 !
1160 !-------------------------------------------------------------------------------
1161 !
1162 !
1163 !*      24.   SET DEFAULT VALUES FOR MODD_PASPOL           
1164 !             ----------------------------------
1165 !
1166 ! other values initialized in modd_paspol
1167 !
1168 IF (KMI == 1) THEN 
1169   LPASPOL    = .FALSE.
1170   NRELEASE   = 0
1171   CPPINIT(:) ='1PT'
1172   XPPLAT(:)  = 0.
1173   XPPLON (:) = 0.
1174   XPPMASS(:) = 0.
1175   XPPBOT(:)  = 0.
1176   XPPTOP(:)  = 0.
1177   CPPT1(:)   = "20010921090000"
1178   CPPT2(:)   = "20010921090000"
1179   CPPT3(:)   = "20010921091500"
1180   CPPT4(:)   = "20010921091500"
1181 ENDIF
1182 !
1183 !-------------------------------------------------------------------------------
1184 !
1185 !
1186 !*      25.   SET DEFAULT VALUES FOR MODD_CONDSAMP         
1187 !             ----------------------------------
1188 !
1189 ! other values initialized in modd_condsamp
1190 !
1191 IF (KMI == 1) THEN 
1192   LCONDSAMP  = .FALSE.
1193   NCONDSAMP  = 3
1194   XRADIO(:)  = 900.
1195   XSCAL(:)   = 1.
1196   XHEIGHT_BASE = 100.
1197   XDEPTH_BASE  = 100.
1198   XHEIGHT_TOP  = 100.
1199   XDEPTH_TOP   = 100.
1200 ENDIF  
1201 !-------------------------------------------------------------------------------
1202 !
1203 !
1204 !*      26.   SET DEFAULT VALUES FOR MODD_LATZ_EDFLX      
1205 !             ----------------------------------
1206 !
1207 IF (KMI == 1) THEN 
1208   LUV_FLX=.FALSE.
1209   XUV_FLX1=3.E+14
1210   XUV_FLX2=0.
1211   LTH_FLX=.FALSE.
1212   XTH_FLX=0.75
1213 ENDIF  
1214 #ifdef MNH_FOREFIRE
1215 !-------------------------------------------------------------------------------
1216 !
1217 !*      27.   SET DEFAULT VALUES FOR MODD_FOREFIRE         
1218 !             ----------------------------------
1219 !
1220 ! other values initialized in modd_forefire
1221 !
1222 IF (KMI == 1) THEN 
1223   LFOREFIRE = .FALSE.
1224   LFFCHEM = .FALSE.
1225   COUPLINGRES = 100.
1226   NFFSCALARS = 0
1227 ENDIF  
1228 #endif                 
1229 !-------------------------------------------------------------------------------
1230 !
1231 END SUBROUTINE DEFAULT_DESFM_n