Newer
Older
!MNH_LIC Copyright 1994-2014 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.
! _Source: /home/cvsroot/MNH-VX-Y-Z/src/MNH/default_desfmn.f90,v _
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
!-----------------------------------------------------------------
! ###########################
MODULE MODI_DEFAULT_DESFM_n
! ###########################
!
INTERFACE
!
SUBROUTINE DEFAULT_DESFM_n(KMI)
INTEGER, INTENT(IN) :: KMI ! Model index
END SUBROUTINE DEFAULT_DESFM_n
!
END INTERFACE
!
END MODULE MODI_DEFAULT_DESFM_n
!
!
!
! ###############################
SUBROUTINE DEFAULT_DESFM_n(KMI)
! ###############################
!
!!**** *DEFAULT_DESFM_n * - set default values for descriptive variables of
!! model KMI
!!
!! PURPOSE
!! -------
! The purpose of this routine is to set default values for the variables
! in descriptor files by filling the corresponding variables which
! are stored in modules.
!
!
!!** METHOD
!! ------
!! Each variable in modules, which can be initialized by reading its
!! value in the descriptor file is set to a default value.
!! When this routine is used during INIT, the modules of the first model
!! are used to temporarily store the variables associated with a nested
!! model.
!! When this routine is used during SPAWNING, the modules of a second
!! model must be initialized.
!! Default values for variables common to all models are set only
!! at the first call of DEFAULT_DESFM_n (i.e. when KMI=1)
!!
!!
!! EXTERNAL
!! --------
!! NONE
!!
!! IMPLICIT ARGUMENTS
!! ------------------
!! Module MODD_PARAMETERS : JPHEXT,JPVEXT
!!
!! Module MODD_CONF : CCONF,L2D,L1D,LFLAT,NMODEL,NVERB
!!
!! Module MODD_DYN : XSEGLEN,XASSELIN,LCORIO,LNUMDIFF
!! XALKTOP,XALZBOT
!!
!! Module MODD_BAKOUT
!!
!! Module MODD_NESTING : NDAD(m),NDTRATIO(m),XWAY(m)
!!
!! Module MODD_CONF_n : LUSERV,LUSERC,LUSERR,LUSERI,LUSERS
!! LUSERG,LUSERH,CSEG,CEXP
!!
!! Module MODD_LUNIT_n : CINIFILE,CCPLFILE
!!
!!
!! Module MODD_DYN_n : XTSTEP,CPRESOPT,NITR,XRELAX,LHO_RELAX
!! LVE_RELAX,XRIMKMAX,NRIMX,NRIMY
!!
!! Module MODD_ADV_n : CUVW_ADV_SCHEME,CMET_ADV_SCHEME,CSV_ADV_SCHEME,NLITER
!!
!! Module MODD_PARAM_n : CTURB,CRAD,CDCONV,CSCONV
!!
Gaelle Tanguy
committed
!! Module MODD_LBC_n : CLBCX, CLBCY,NLBLX,NLBLY,XCPHASE,XCPHASE_PBL,XPOND
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
!!
!! Module MODD_TURB_n : XIMPL,CTURBLEN,CTURBDIM,LTURB_FLX,LTURB_DIAG,LSUBG_COND
!! LTGT_FLX
!!
!!
!! Module MODD_PARAM_RAD_n:
!! XDTRAD,XDTRAD_CLONLY,LCLEAR_SKY,NRAD_COLNBR, NRAD_DIAG
!!
!! Module MODD_BUDGET : CBUTYPE,NBUMOD,XBULEN,NBUKL, NBUKH,LBU_KCP,XBUWRI
!! NBUIL, NBUIH,NBUJL, NBUJH,LBU_ICP,LBU_JCP,NBUMASK
!! LBU_RU,LBU_RV,LBU_RW,LBU_RTH,LBU_RTKE,LBU_RRV,LBU_RRC,LBU_RRR
!! LBU_RRI,LBU_RRS,LBU_RRG,LBU_RRH,LBU_RSVx
!! NADVXU, NADVYU, NADVZU, NCURVU, NCORU, NDIFU, NRELU, NHTURBU,
!! NVTURBU, NPRESU
!! NADVXV, NADVYV, NADVZV, NCURVV, NCORV, NDIFV, NRELV, NHTURBV,
!! NVTURBV, NPRESV
!! NADVXW, NADVYW, NADVZW, NCURVW, NCORW, NGRAVW, NDIFW, NRELW,
!! NHTURBW, NVTURBW, NPRESW
!! NADVXTH, NADVYTH, NADVZTH, NPREFTH, NDIFTH, NRELTH, NHTURBTH,
!! NVTURBTH, NREVATH, NCONDTH
!! NADVXTKE, NADVYTKE, NADVZTKE, NDIFTKE, NDPTKE, NTPTKE, NDISSTKE,
!! NTRTKE
!! NADVXRV, NADVYRV, NADVZRV, NDIFRV, NRELRV, NHTURBRV, NVTURBRV,
!! NREVARV, NCONDRV
!! NADVXRC, NADVYRC, NADVZRC, NDIFRC, NHTURBRC, NVTURBRC, NACCRRC,
!! NAUTORC, NCONDRC
!! NADVXRR, NADVYRR, NADVZRR, NDIFRR, NACCRRR, NAUTORR, NREVARR,
!! NSEDIRR
!! NADVXRI, NADVYRI, NADVZRI, NDIFRI
!! NADVXRS, NADVYRS, NADVZRS, NDIFRS
!! NADVXRG, NADVYRG, NADVZRG, NDIFRG
!! NADVXRH, NADVYRH, NADVZRH, NDIFRH
!! NADVXSVx, NADVYSVx, NADVZSVx, NDIFSVx, NHTURBSVx, NVTURBSVx
!!
!! Module MODD_BLANK :
!!
!! XDUMMYi, NDUMMYi, LDUMMYi, CDUMMYi
!!
!! Module MODD_FRC :
!!
!! LGEOST_UV_FRC,LGEOST_TH_FRC,LTEND_THRV_FRC
!! LVERT_MOTION_FRC,LRELAX_THRV_FRC,LRELAX_UV_FRC,XRELAX_TIME_FRC
!! XRELAX_HEIGHT_FRC,CRELAX_HEIGHT_TYPE,LTRANS,XUTRANS,XVTRANS,
!! LPGROUND_FRC
!!
!! Module MODD_PARAM_ICE :
!!
!! LWARM,CPRISTINE_ICE
!!
!! Module MODD_PARAM_KAFR_n :
!!
!! XDTCONV,LREFRESH_ALL,LDOWN,NICE,LCHTRANS
!!
!! Module MODD_PARAM_MFSHALL_n :
!!
!! CMF_UPDRAFT,LMIXUV,CMF_CLOUD,XIMPL_MF,LMF_FLX
!!
!!
!!
!!
!! REFERENCE
!! ---------
!! Book2 of the documentation (routine DEFAULT_DESFM_n)
!!
!!
!! AUTHOR
!! ------
!! V. Ducrocq * Meteo France *
!!
!! MODIFICATIONS
!! -------------
!! Original 02/06/94
!! Modifications 17/10/94 (Stein) For LCORIO
!! Modifications 06/12/94 (Stein) remove LBOUSS+add LABSLAYER, LNUMDIFF
!! ,LSTEADYLS
!! Modifications 06/12/94 (Stein) remove LABSLAYER, add LHO_RELAX,
!! LVE_RELAX, NRIMX, NRIMY, XRIMKMAX
!! Modifications 09/01/95 (Lafore) add LSTEADY_DMASS
!! Modifications 09/01/95 (Stein) add the turbulence scheme namelist
!! Modifications 09/01/95 (Stein) add the 1D switch
!! Modifications 10/03/95 (Mallet) add the coupling files
!! 29/06/95 ( Stein, Nicolau, Hereil) add the budgets
!! Modifications 25/09/95 ( Stein )add the LES tools
!! Modifications 25/10/95 ( Stein )add the radiations
!! Modifications 23/10/95 (Vila, lafore) new scalar advection scheme
!! Modifications 24/02/96 (Stein) change the default value for CCPLFILE
!! Modifications 12/02/96 (Lafore) transformation to DEFAULT_DESFM_n for
!! spawning
!! Modifications 25/04/96 (Suhre) add the blank module
!! Modifications 29/07/96 (Pinty&Suhre) add module MODD_FRC
!! Modifications 11/04/96 (Pinty) add the rain-ice scheme and modify
!! the splitted arrays in MODD_PARAM_RAD_n
!! Modifications 11/01/97 (Pinty) add the deep convection scheme
!! Modifications 24/11/96 (Masson) add LREFRESH_ALL in deep convection
!! Modifications 12/02/96 (Lafore) transformation to DEFAULT_DESFM_n for spawning
!! Modifications 22/07/96 (Lafore) gridnesting implementation
!! Modifications 29/07/96 (Lafore) add the module MODD_FMOUT (renamed MODD_BAKOUT)
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
!! Modifications 23/06/97 (Stein) add the equation system name
!! Modifications 10/07/97 (Masson) add MODD_PARAM_GROUNDn : CROUGH
!! Modifications 28/07/97 (Masson) remove LREFRESH_ALL and LSTEADY_DMASS
!! Modifications 08/10/97 (Stein) switch (_n=1) to initialize the
!! parameters common to all models
!! Modifications 24/01/98 (Bechtold) add LREFRESH_ALL, LCHTRANS,
!! LTEND_THRV_FR and LSST_FRC
!! Modifications 18/07/99 (Stein) add LRAD_DIAG
!! Modification 15/03/99 (Masson) use of XUNDEF
!! Modification 11/12/00 (Tomasini) Add CSEA_FLUX to MODD_PARAMn
!! Modification 22/01/01 (Gazen) delete NSV and add LHORELAX_SVC2R2
!! LHORELAX_SVCHEM,LHORELAX_SVLG
!! Modification 15/03/02 (Solmon) radiation scheme: remove NSPOT and add
!! default for aerosol and cloud rad. prop. control
!! Modification 22/05/02 (Jabouille) put chimical default here
!! Modification 01/2004 (Masson) removes surface (externalization)
!! 09/04 (M. Tomasini) New namelist to modify the
!! Cloud mixing length
!! 07/05 (P.Tulet) New namelists for dust and aerosol
!! Modification 01/2007 (Malardel, Pergaud) Add MODD_PARAM_MFSHALL_n
!! Modification 10/2009 (Aumond) Add user multimasks for LES
!! Modification 10/2009 (Aumond) Add MEAN_FIELD
!! Modification 12/04/07 (Leriche) add LUSECHAQ for aqueous chemistry
!! Modification 30/05/07 (Leriche) add LCH_PH and XCH_PHINIT for pH
!! Modification 25/04/08 (Leriche) add XRTMIN_AQ LWC threshold for aq. chemistry
!! 16/07/10 add LHORELAX_SVIC
!! 16/09/10 add LUSECHIC
!! 13/01/11 add LCH_RET_ICE
!! 01/07/11 (F.Couvreux) Add CONDSAMP
!! 01/07/11 (B.Aouizerats) Add CAOP
Gaelle Tanguy
committed
!! 07/2013 (C.Lac) add WENO, LCHECK
!! 07/2013 (Bosseur & Filippi) adds Forefire
Gaelle Tanguy
committed
!! 08/2015 (Redelsperger & Pianezze) add XPOND coefficient for LBC
!! Modification 24/03/16 (Leriche) remove LCH_SURFACE_FLUX
!! put NCH_VEC_LENGTH = 50 instead of 1000
!!
!! 04/2016 (C.LAC) negative contribution to the budget splitted between advection, turbulence and microphysics for KHKO/C2R2
!! Modification 24/03/16 (Leriche) remove LCH_SURFACE_FLUX
!! put NCH_VEC_LENGTH = 50 instead of 1000
Gaelle TANGUY
committed
!! 10/2016 (C.Lac) VSIGQSAT change from 0 to 0.02 for coherence with AROME
!! 10/2016 (C.Lac) Add droplet deposition
!! 10/2016 (R.Honnert and S.Riette) : Improvement of EDKF and adaptation to the grey zone
Gaelle DELAUTIER
committed
!! 10/2016 (F Brosse) add prod/loss terms computation for chemistry
!! 07/2017 (V. Masson) adds time step for output files writing.
!! 09/2017 Q.Rodier add LTEND_UV_FRC
!! 02/2018 Q.Libois ECRAD
! Philippe Wautelet: 05/2016-04/2018: new data structures and calls for I/O
Gaelle DELAUTIER
committed
!! 01/2018 (S. Riette) new budgets and variables for ICE3/ICE4
!! 01/2018 (J.Colin) add VISC and DRAG
!! 07/2017 (V. Vionnet) add blowing snow variables
RODIER Quentin
committed
!! 01/2019 (R. Honnert) add reduction of the mass-flux surface closure with the resolution
RODIER Quentin
committed
!! Bielli S. 02/2019 Sea salt : significant sea wave height influences salt emission; 5 salt modes
RODIER Quentin
committed
!! 05/2019 F.Brient add tracer emission from the top of the boundary-layer
!!
!-------------------------------------------------------------------------------
!
!* 0. DECLARATIONS
! ------------
USE MODD_PARAMETERS
USE MODD_CONF ! For INIT only DEFAULT_DESFM1
USE MODD_CONFZ
USE MODD_DYN
USE MODD_NESTING
USE MODD_BAKOUT
USE MODD_SERIES
USE MODD_CONF_n ! modules used to set the default values is only
USE MODD_LUNIT_n ! the one corresponding to model 1. These memory
USE MODD_DIM_n ! addresses will then be filled by the values read in
USE MODD_DYN_n ! the DESFM corresponding to model n which may have
USE MODD_ADV_n ! missing values. This is why we affect default values.
USE MODD_PARAM_n ! For SPAWNING DEFAULT_DESFM2 is also used
USE MODD_LBC_n
USE MODD_OUT_n
USE MODD_TURB_n
USE MODD_BUDGET
USE MODD_LES
USE MODD_PARAM_RAD_n
#ifdef MNH_ECRAD
USE MODD_PARAM_ECRAD_n
#endif
USE MODD_BLANK
USE MODD_FRC
USE MODD_PARAM_ICE
USE MODD_PARAM_C2R2
USE MODD_TURB_CLOUD
USE MODD_PARAM_KAFR_n
USE MODD_PARAM_MFSHALL_n
USE MODD_CH_MNHC_n
USE MODD_SERIES_n
USE MODD_NUDGING_n
USE MODD_CH_AEROSOL
USE MODD_DUST
USE MODD_SALT
USE MODD_PASPOL
USE MODD_CONDSAMP
USE MODD_MEAN_FIELD
USE MODD_DRAGTREE
!
!
USE MODD_PARAM_LIMA, ONLY : LCOLD, LNUCL, LSEDI, LHHONI, LSNOW, LHAIL, LMEYERS,&
Gaelle DELAUTIER
committed
NMOD_IFN, XIFN_CONC, LIFN_HOM, CIFN_SPECIES, &
CINT_MIXING, NMOD_IMM, NIND_SPECIE, &
CPRISTINE_ICE_LIMA, CHEVRIMED_ICE_LIMA, &
XFACTNUC_DEP, XFACTNUC_CON, &
OWARM=>LWARM, LACTI, ORAIN=>LRAIN, OSEDC=>LSEDC, &
Gaelle DELAUTIER
committed
OACTIT=>LACTIT, LBOUND, NMOD_CCN, XCCN_CONC, &
LCCN_HOM, CCCN_MODES, &
YALPHAR=>XALPHAR, YNUR=>XNUR, &
YALPHAC=>XALPHAC, YNUC=>XNUC, CINI_CCN=>HINI_CCN, &
CTYPE_CCN=>HTYPE_CCN, YFSOLUB_CCN=>XFSOLUB_CCN, &
YACTEMP_CCN=>XACTEMP_CCN, YAERDIFF=>XAERDIFF, &
YAERHEIGHT=>XAERHEIGHT, &
LSCAV, LAERO_MASS, NPHILLIPS, &
ODEPOC=>LDEPOC, OVDEPOC=>XVDEPOC, OACTTKE=>LACTTKE,&
LPTSPLIT, L_LFEEDBACKT=>LFEEDBACKT, L_NMAXITER=>NMAXITER, &
L_XMRSTEP=>XMRSTEP, L_XTSTEP_TS=>XTSTEP_TS
USE MODD_LATZ_EDFLX
USE MODD_2D_FRC
USE MODD_BLOWSNOW
USE MODD_BLOWSNOW_n
Gaelle DELAUTIER
committed
USE MODD_DRAG_n
#ifdef MNH_FOREFIRE
USE MODD_FOREFIRE
#endif
!
IMPLICIT NONE
!
!* 0.1 declarations of arguments
!
INTEGER, INTENT(IN) :: KMI ! Model index
!
!* 0.2 declaration of local variables
!
INTEGER :: JM ! loop index
!
!-------------------------------------------------------------------------------
!
!* 1. SET DEFAULT VALUES FOR MODD_LUNIT_n :
! ----------------------------------
!
WAUTELET Philippe
committed
! CINIFILE='INIFILE'
CINIFILEPGD='' !Necessary to keep this line to prevent problems with spawning
CCPLFILE(:)=' '
!
!-------------------------------------------------------------------------------
!
!* 2. SET DEFAULT VALUES FOR MODD_CONF AND MODD_CONF_n :
! ------------------------------------------------
!
IF (KMI == 1) THEN
CCONF ='START'
LTHINSHELL = .FALSE.
L2D = .FALSE.
L1D = .FALSE.
LFLAT = .FALSE.
NMODEL = 1
CEQNSYS = 'DUR'
NVERB = 5
CEXP = 'EXP01'
CSEG = 'SEG01'
LFORCING = .FALSE.
L2D_ADV_FRC= .FALSE.
L2D_REL_FRC= .FALSE.
XRELAX_HEIGHT_BOT = 0.
XRELAX_HEIGHT_TOP = 30000.
XRELAX_TIME = 864000.
LPACK = .TRUE.
NHALO = 1
#ifdef MNH_SX5
CSPLIT ='YSPLITTING' ! NEC vectoriel architecture , low number of PROC
#else
CSPLIT ='BSPLITTING' ! Scalaire architecture , high number of PROC
#endif
NZ_PROC = 0 !JUAN Z_SPLITTING :: number of proc in Z splitting
NZ_SPLITTING = 10 !JUAN Z_SPLITTING :: for debug NZ=1=flat_inv; NZ=10=flat_invz; NZ=1+2 the two
LLG = .FALSE.
LINIT_LG = .FALSE.
CINIT_LG = 'FMOUT'
LNOMIXLG = .FALSE.
END IF
!
CCLOUD = 'NONE'
LUSERV = .TRUE.
LUSERC = .FALSE.
LUSERR = .FALSE.
LUSERI = .FALSE.
LUSERS = .FALSE.
LUSERG = .FALSE.
LUSERH = .FALSE.
!NSV = 0
!NSV_USER = 0
LUSECI = .FALSE.
!
!-------------------------------------------------------------------------------
!
!* 3. SET DEFAULT VALUES FOR MODD_DYN AND MODD_DYN_n :
! -----------------------------------------------
!
IF (KMI == 1) THEN
XSEGLEN = 43200.
XASSELIN = 0.2
XASSELIN_SV = 0.02
LCORIO = .TRUE.
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
427
428
429
430
431
432
433
434
435
436
LNUMDIFTH = .FALSE.
LNUMDIFSV = .FALSE.
XALZBOT = 4000.
XALKTOP = 0.01
XALKGRD = 0.01
XALZBAS = 0.01
END IF
!
XTSTEP = 60.
CPRESOPT = 'CRESI'
NITR = 4
LITRADJ = .TRUE.
LRES = .FALSE.
XRES = 1.E-07
XRELAX = 1.
LVE_RELAX = .FALSE.
LVE_RELAX_GRD = .FALSE.
XRIMKMAX = 0.01 / XTSTEP
XT4DIFU = 1800.
XT4DIFTH = 1800.
XT4DIFSV = 1800.
!
IF (KMI == 1) THEN ! for model 1 we have a Large scale information
NRIMX = JPRIMMAX ! for U,V,W,TH,Rv used for the hor. relaxation
NRIMY = JPRIMMAX
ELSE
NRIMX = 0 ! for inner models we use only surfacic fields to
NRIMY = 0 ! give the lbc and no hor. relaxation is used
END IF
!
LHORELAX_UVWTH = .FALSE.
LHORELAX_RV = .FALSE.
LHORELAX_RC = .FALSE. ! for all these fields, no large scale is usally available
LHORELAX_RR = .FALSE. ! for model 1 and for inner models, we only use surfacic
LHORELAX_RS = .FALSE. ! fiels ( no hor. relax. )
LHORELAX_RI = .FALSE.
LHORELAX_RG = .FALSE.
LHORELAX_RH = .FALSE.
LHORELAX_TKE = .FALSE.
LHORELAX_SV(:) = .FALSE.
LHORELAX_SVC2R2 = .FALSE.
LHORELAX_SVC1R3 = .FALSE.
LHORELAX_SVELEC = .FALSE.
LHORELAX_SVLG = .FALSE.
LHORELAX_SVCHEM = .FALSE.
LHORELAX_SVCHIC = .FALSE.
LHORELAX_SVDST = .FALSE.
LHORELAX_SVSLT = .FALSE.
LHORELAX_SVPP = .FALSE.
LHORELAX_SVCS = .FALSE.
LHORELAX_SVAER = .FALSE.
#ifdef MNH_FOREFIRE
LHORELAX_SVFF = .FALSE.
#endif
LHORELAX_SVSNW = .FALSE.
!
!
!-------------------------------------------------------------------------------
!
!* 4. SET DEFAULT VALUES FOR MODD_NESTING :
! -----------------------------------
!
IF (KMI == 1) THEN
NDAD(1)=1
DO JM=2,JPMODELMAX
NDAD(JM) = JM - 1
END DO
NDTRATIO(:) = 1
XWAY(:) = 2. ! two-way interactive gridnesting
XWAY(1) = 0. ! except for model 1
END IF
!
!-------------------------------------------------------------------------------
!
!* 5. SET DEFAULT VALUES FOR MODD_ADV_n :
! ----------------------------------
!
CUVW_ADV_SCHEME = 'CEN4TH'
CMET_ADV_SCHEME = 'PPM_01'
CSV_ADV_SCHEME = 'PPM_01'
CTEMP_SCHEME = 'RKC4'
NWENO_ORDER = 3
NSPLIT = 1
LSPLIT_CFL = .TRUE.
Gaelle TANGUY
committed
LSPLIT_WENO = .TRUE.
!
!-------------------------------------------------------------------------------
!
!* 6. SET DEFAULT VALUES FOR MODD_PARAM_n :
! -----------------------------------
!
CTURB = 'NONE'
CRAD = 'NONE'
CDCONV = 'NONE'
CSCONV = 'NONE'
CELEC = 'NONE'
CACTCCN = 'NONE'
!
!-------------------------------------------------------------------------------
!
!* 7. SET DEFAULT VALUES FOR MODD_LBC_n :
! ---------------------------------
!
CLBCX(1) ='CYCL'
CLBCX(2) ='CYCL'
CLBCY(1) ='CYCL'
CLBCY(2) ='CYCL'
NLBLX(:) = 1
NLBLY(:) = 1
XCPHASE = 20.
Gaelle Tanguy
committed
XCPHASE_PBL = 0.
XCARPKMAX = XUNDEF
!
!-------------------------------------------------------------------------------
!
!* 8. SET DEFAULT VALUES FOR MODD_NUDGING_n :
! ---------------------------------
!
LNUDGING = .FALSE.
XTNUDGING = 21600.
!
!-------------------------------------------------------------------------------
!
!* 9. SET DEFAULT VALUES FOR MODD_BAKOUT and MODD_OUT_n :
! ------------------------------------------------
!
!
!
!-------------------------------------------------------------------------------
!
!* 10. SET DEFAULT VALUES FOR MODD_TURB_n :
! ----------------------------------
!
XIMPL = 1.
XKEMIN = 0.01
Gaelle DELAUTIER
committed
XCEDIS = 0.84
CTURBLEN = 'BL89'
CTURBDIM = '1DIM'
LTURB_FLX =.FALSE.
LTURB_DIAG=.FALSE.
LSUBG_COND=.FALSE.
CSUBG_AUCV='NONE'
LSIGMAS =.TRUE.
LSIG_CONV =.FALSE.
LRMC01 =.FALSE.
CTOM ='NONE'
Gaelle TANGUY
committed
VSIGQSAT = 0.02
!-------------------------------------------------------------------------------
!
!* 10b. SET DEFAULT VALUES FOR MODD_DRAGTREE :
! ----------------------------------
!
LDRAGTREE = .FALSE.
Gaelle TANGUY
committed
LDEPOTREE = .FALSE.
XVDEPOTREE = 0.02 ! 2 cm/s
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
!
!-------------------------------------------------------------------------------
!
!* 11. SET DEFAULT VALUES FOR MODD_BUDGET :
! ------------------------------------
!
! 11.1 General budget variables
!
IF (KMI == 1) THEN
CBUTYPE = 'NONE'
NBUMOD = 1
XBULEN = XSEGLEN
XBUWRI = XSEGLEN
NBUKL = 1
NBUKH = 0
LBU_KCP = .TRUE.
!
! 11.2 Variables for the cartesian box
!
NBUIL = 1
NBUIH = 0
NBUJL = 1
NBUJH = 0
LBU_ICP = .TRUE.
LBU_JCP = .TRUE.
!
! 11.3 Variables for the mask
!
NBUMASK = 1
!
! 11.4 Variables for budget and processes choice
LBU_RU = .FALSE.
NASSEU = 0
NNESTU = 0
NFRCU = 0
NNUDU = 0
NCURVU = 0
NCORU = 0
NDIFU = 0
NRELU = 0
NMAFLU = 0
NDRAGU = 0
NVTURBU = 0
NHTURBU = 0
NPRESU = 0
!
! Budget of RV
LBU_RV = .FALSE.
NASSEV = 0
NNESTV = 0
NFRCV = 0
NNUDV = 0
NCURVV = 0
NCORV = 0
NDIFV = 0
NRELV = 0
NMAFLV = 0
NDRAGV = 0
NVTURBV = 0
NHTURBV = 0
NPRESV = 0
!
! Budget of RW
LBU_RW = .FALSE.
NASSEW = 0
NNESTW = 0
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
NFRCW = 0
NNUDW = 0
NCURVW = 0
NCORW = 0
NGRAVW = 0
NDIFW = 0
NRELW = 0
NVTURBW = 0
NHTURBW = 0
NPRESW = 0
!
! Budget of RTH
LBU_RTH = .FALSE.
NASSETH = 0
NNESTTH = 0
NADVTH = 0
NFRCTH = 0
NNUDTH = 0
NPREFTH = 0
NDIFTH = 0
NRELTH = 0
NRADTH = 0
NDCONVTH = 0
NMAFLTH = 0
NVTURBTH = 0
NHTURBTH = 0
NDISSHTH = 0
NNEGATH = 0
NNEADVTH = 0
NNETURTH = 0
NNECONTH = 0
NREVATH = 0
NCONDTH = 0
NHENUTH = 0
NHONTH = 0
NSFRTH = 0
NDEPSTH = 0
NDEPGTH = 0
NREVATH = 0
NRIMTH = 0
NACCTH = 0
NCFRZTH = 0
NWETGTH = 0
NDRYGTH = 0
Gaelle DELAUTIER
committed
NWETHTH = 0
NDRYHTH = 0
NGMLTTH = 0
NIMLTTH = 0
NBERFITH = 0
NCDEPITH = 0
Gaelle DELAUTIER
committed
NHMLTTH = 0
NCORRTH = 0
NADJUTH = 0
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
!
! Budget of RTKE
LBU_RTKE = .FALSE.
NASSETKE = 0
NADVTKE = 0
NFRCTKE = 0
NDIFTKE = 0
NRELTKE = 0
NDRAGTKE = 0
NDPTKE = 0
NTPTKE = 0
NDISSTKE = 0
NTRTKE = 0
!
! Budget of RRV
LBU_RRV = .FALSE.
NASSERV = 0
NNESTRV = 0
NADVRV = 0
NFRCRV = 0
NNUDRV = 0
NDIFRV = 0
NRELRV = 0
NDCONVRV = 0
NMAFLRV = 0
NVTURBRV = 0
NHTURBRV = 0
NNEGARV = 0
NNETURRV = 0
NNECONRV = 0
NNEADVRV = 0
NREVARV = 0
NCONDRV = 0
NHENURV = 0
NDEPSRV = 0
NDEPGRV = 0
NREVARV = 0
NCDEPIRV = 0
Gaelle DELAUTIER
committed
NCORRRV = 0
NADJURV = 0
!
! Budget of RRC
LBU_RRC = .FALSE.
NASSERC = 0
NNESTRC = 0
NADVRC = 0
NFRCRC = 0
NDIFRC = 0
NRELRC = 0
NDCONVRC = 0
NVTURBRC = 0
NHTURBRC = 0
NNEGARC = 0
NNETURRC = 0
NNECONRC = 0
NNEADVRC = 0
NACCRRC = 0
NAUTORC = 0
NCONDRC = 0
NAUTORC = 0
NACCRRC = 0
NHONRC = 0
NRIMRC = 0
Gaelle DELAUTIER
committed
NCMELRC = 0
NWETGRC = 0
NDRYGRC = 0
NIMLTRC = 0
NBERFIRC = 0
NCDEPIRC = 0
NWETHRC = 0
Gaelle DELAUTIER
committed
NDRYHRC = 0
NCORRRC = 0
NADJURC = 0
Gaelle TANGUY
committed
NDEPORC = 0
NDEPOTRRC = 0
!
! Budget of RRR
LBU_RRR = .FALSE.
NASSERR = 0
NNESTRR = 0
NADVRR = 0
NFRCRR = 0
NDIFRR = 0
NRELRR = 0
NNEGARR = 0
NACCRRR = 0
NAUTORR = 0
NREVARR = 0
NSEDIRR = 0
NSFRRR = 0
NACCRR = 0
Gaelle DELAUTIER
committed
NCMELRR = 0
NCFRZRR = 0
NWETGRR = 0
NDRYGRR = 0
NGMLTRR = 0
NWETHRR = 0
Gaelle DELAUTIER
committed
NDRYHRR = 0
Gaelle DELAUTIER
committed
NCORRRR = 0
!
! Budget of RRI
LBU_RRI = .FALSE.
NASSERI = 0
NNESTRI = 0
NADVRI = 0
NFRCRI = 0
NDIFRI = 0
NRELRI = 0
NDCONVRI = 0
NVTURBRI = 0
NHTURBRI = 0
NNEGARI = 0
NSEDIRI = 0
NHENURI = 0
NHONRI = 0
NAGGSRI = 0
NAUTSRI = 0
NCFRZRI = 0
NWETGRI = 0
NDRYGRI = 0
NIMLTRI = 0
NBERFIRI = 0
NCDEPIRI = 0
NWETHRI = 0
Gaelle DELAUTIER
committed
NDRYHRI = 0
NCORRRI = 0
NADJURI = 0
!
! Budget of RRS
LBU_RRS = .FALSE.
NASSERS = 0
NNESTRS = 0
NADVRS = 0
NFRCRS = 0
NDIFRS = 0
NRELRS = 0
NNEGARS = 0
NSEDIRS = 0
NDEPSRS = 0
NAGGSRS = 0
NAUTSRS = 0
NRIMRS = 0
NACCRS = 0
NCMELRS = 0
NWETGRS = 0
NDRYGRS = 0
NWETHRS = 0
Gaelle DELAUTIER
committed
NDRYHRS = 0
NCORRRS = 0
!
! Budget of RRG
LBU_RRG = .FALSE.
NASSERG = 0
NNESTRG = 0
NADVRG = 0
NFRCRG = 0
NDIFRG = 0
NRELRG = 0
NNEGARG = 0
NSEDIRG = 0
NSFRRG = 0
NDEPGRG = 0
NRIMRG = 0
NACCRG = 0
NCMELRG = 0
NCFRZRG = 0
NWETGRG = 0
NDRYGRG = 0
NGMLTRG = 0
NWETHRG = 0
Gaelle DELAUTIER
committed
NDRYHRG = 0
NCORRRG = 0
NHGCVRG = 0
NGHCVRG = 0
!
! Budget of RRH
LBU_RRH = .FALSE.
NASSERH = 0
NNESTRH = 0
NADVRH = 0
NFRCRH = 0
NDIFRH = 0
NRELRH = 0
NNEGARH = 0
NWETGRH = 0
NWETHRH = 0
Gaelle DELAUTIER
committed
NDRYHRH = 0
Gaelle DELAUTIER
committed
NCORRRH = 0
NHGCVRH = 0
NGHCVRH = 0
!
! Budget of RSVx
LBU_RSV = .FALSE.
NASSESV = 0
NNESTSV = 0
NADVSV = 0
NFRCSV = 0
NDIFSV = 0
NRELSV = 0
NDCONVSV = 0
NMAFLSV = 0
NVTURBSV = 0
NHTURBSV = 0
NCHEMSV = 0
Gaelle Tanguy
committed
NNEGASV = 0
Gaelle TANGUY
committed
NDEPOSV = 0
NDEPOTRSV = 0
878
879
880
881
882
883
884
885
886
887
888
889
890
891
892
893
894
895
896
897
898
899
900
901
902
903
904
905
906
907
908
909
910
911
912
913
914
915
916
917
918
919
920
921
922
923
924
925
926
927
928
929
930
931
932
933
934
935
936
937
938
939
940
941
942
943
944
945
946
947
948
949
950
951
952
953
954
955
956
957
958
959
960
961
962
963
964
965
966
967
!
!
END IF
!
!-------------------------------------------------------------------------------
!
!* 12. SET DEFAULT VALUES FOR MODD_LES :
! ---------------------------------
!
IF (KMI == 1) THEN
LLES_MEAN = .FALSE.
LLES_RESOLVED = .FALSE.
LLES_SUBGRID = .FALSE.
LLES_UPDRAFT = .FALSE.
LLES_DOWNDRAFT = .FALSE.
LLES_SPECTRA = .FALSE.
!
NLES_LEVELS = NUNDEF
XLES_ALTITUDES = XUNDEF
NSPECTRA_LEVELS = NUNDEF
XSPECTRA_ALTITUDES = XUNDEF
NLES_TEMP_SERIE_I = NUNDEF
NLES_TEMP_SERIE_J = NUNDEF
NLES_TEMP_SERIE_Z = NUNDEF
CLES_NORM_TYPE = 'NONE'
CBL_HEIGHT_DEF = 'KE'
XLES_TEMP_SAMPLING = XUNDEF
XLES_TEMP_MEAN_START = XUNDEF
XLES_TEMP_MEAN_END = XUNDEF
XLES_TEMP_MEAN_STEP = 3600.
LLES_CART_MASK = .FALSE.
NLES_IINF = NUNDEF
NLES_ISUP = NUNDEF
NLES_JINF = NUNDEF
NLES_JSUP = NUNDEF
LLES_NEB_MASK = .FALSE.
LLES_CORE_MASK = .FALSE.
LLES_MY_MASK = .FALSE.
NLES_MASKS_USER = NUNDEF
LLES_CS_MASK = .FALSE.
LLES_PDF = .FALSE.
NPDF = 1
XTH_PDF_MIN = 270.
XTH_PDF_MAX = 350.
XW_PDF_MIN = -10.
XW_PDF_MAX = 10.
XTHV_PDF_MIN = 270.
XTHV_PDF_MAX = 350.
XRV_PDF_MIN = 0.
XRV_PDF_MAX = 20.
XRC_PDF_MIN = 0.
XRC_PDF_MAX = 1.
XRR_PDF_MIN = 0.
XRR_PDF_MAX = 1.
XRI_PDF_MIN = 0.
XRI_PDF_MAX = 1.
XRS_PDF_MIN = 0.
XRS_PDF_MAX = 1.
XRG_PDF_MIN = 0.
XRG_PDF_MAX = 1.
XRT_PDF_MIN = 0.
XRT_PDF_MAX = 20.
XTHL_PDF_MIN = 270.
XTHL_PDF_MAX = 350.
END IF
!
!-------------------------------------------------------------------------------
!
!* 13. SET DEFAULT VALUES FOR MODD_PARAM_RAD_n :
! ---------------------------------------
!
XDTRAD = XTSTEP
XDTRAD_CLONLY = XTSTEP
LCLEAR_SKY =.FALSE.
NRAD_COLNBR = 1000
NRAD_DIAG = 0
CLW ='RRTM'
CAER='SURF'
CAOP='CLIM'
CEFRADL='MART'
CEFRADI='LIOU'
COPWSW = 'FOUQ'
COPISW = 'EBCU'
COPWLW = 'SMSH'
COPILW = 'EBCU'
XFUDG = 1.
LAERO_FT=.FALSE.
LFIX_DAT=.FALSE.
!
#ifdef MNH_ECRAD
!* 13bis. SET DEFAULT VALUES FOR MODD_PARAM_ECRAD_n :
! ---------------------------------------
!
NSWSOLVER = 0 ! 0: 'McICA 1: 'SPARTACUS' 2: 'SPARTACUS' + 3D effect
NLWSOLVER = 0 ! 0: 'McICA 1: 'SPARTACUS' 2: 'SPARTACUS' + 3D effect
! LEFF3D = .TRUE.
! LSIDEM = .TRUE.
NREG = 3 ! Number of cloudy regions (3=TripleClouds)
! LLWCSCA = .TRUE. ! LW cloud scattering
! LLWASCA = .TRUE. ! LW aerosols scattering
NLWSCATTERING = 2
NAERMACC = 0
! CGAS = 'RRTMG-IFS' ! Gas optics model
NOVLP = 1 ! overlap assumption ; 0= 'Max-Ran' ; 1= 'Exp-Ran'; 2 = 'Exp-Exp'
NLIQOPT = 3 ! 1: 'Monochromatic', 2: 'HuStamnes', 3: 'SOCRATES', 4: 'Slingo'
NICEOPT = 3 ! 1: 'Monochromatic', 2: 'Fu-PSRAD', 3: 'Fu-IFS', 4: 'Baran', 5: 'Baran2016', 6: 'Baran2017'
! LSW_ML_E = .FALSE.
! LLW_ML_E = .FALSE.
! LPSRAD = .FALSE.
!
NRADLP = 1 ! 0: ERA-15, 1: Zhang and Rossow, 2: Martin (1994) et Woods (2000)
NRADIP = 1 ! 0: 40 mum, 1: Liou and Ou (1994), 2: Liou and Ou (1994) improved, 3: Sun and Rikus (1999)
XCLOUD_FRAC_STD = 1.0_JPRB ! change to 0.75 for more realistic distribution
#endif
!-------------------------------------------------------------------------------
!
!* 14. SET DEFAULT VALUES FOR MODD_BLANK :
! -----------------------------------
!
IF (KMI == 1) THEN
XDUMMY1 = 0.
XDUMMY2 = 0.