Newer
Older
1001
1002
1003
1004
1005
1006
1007
1008
1009
1010
1011
1012
1013
1014
1015
1016
1017
1018
1019
1020
1021
1022
1023
1024
1025
1026
1027
1028
1029
1030
1031
1032
1033
1034
1035
1036
1037
1038
1039
1040
1041
1042
1043
1044
1045
1046
1047
XDUMMY4 = 0.
XDUMMY5 = 0.
XDUMMY6 = 0.
XDUMMY7 = 0.
XDUMMY8 = 0.
XDUMMY=0.
!
NDUMMY1 = 0
NDUMMY2 = 0
NDUMMY3 = 0
NDUMMY4 = 0
NDUMMY5 = 0
NDUMMY6 = 0
NDUMMY7 = 0
NDUMMY8 = 0
NDUMMY=0
!
LDUMMY1 = .TRUE.
LDUMMY2 = .TRUE.
LDUMMY3 = .TRUE.
LDUMMY4 = .TRUE.
LDUMMY5 = .TRUE.
LDUMMY6 = .TRUE.
LDUMMY7 = .TRUE.
LDUMMY8 = .TRUE.
LDUMMY=.TRUE.
!
CDUMMY1 = ' '
CDUMMY2 = ' '
CDUMMY3 = ' '
CDUMMY4 = ' '
CDUMMY5 = ' '
CDUMMY6 = ' '
CDUMMY7 = ' '
CDUMMY8 = ' '
CDUMMY= ' '
END IF
!
!------------------------------------------------------------------------------
!
!* 15. SET DEFAULT VALUES FOR MODD_FRC :
! ---------------------------------
!
IF (KMI == 1) THEN
LGEOST_UV_FRC = .FALSE.
LGEOST_TH_FRC = .FALSE.
LTEND_THRV_FRC = .FALSE.
LVERT_MOTION_FRC = .FALSE.
LRELAX_THRV_FRC = .FALSE.
LRELAX_UV_FRC = .FALSE.
XRELAX_TIME_FRC = 10800.
XRELAX_HEIGHT_FRC = 0.
CRELAX_HEIGHT_TYPE = "FIXE"
LTRANS = .FALSE.
XUTRANS = 0.0
XVTRANS = 0.0
LPGROUND_FRC = .FALSE.
END IF
!
!-------------------------------------------------------------------------------
!
!
!* 16. SET DEFAULT VALUES FOR MODD_PARAM_ICE :
! ---------------------------------------
!
IF (KMI == 1) THEN
Gaelle DELAUTIER
committed
LRED = .FALSE.
LWARM = .TRUE.
CPRISTINE_ICE = 'PLAT'
Gaelle DELAUTIER
committed
LSEDIC = .TRUE.
Gaelle Tanguy
committed
LCONVHG = .FALSE.
CSEDIM = 'SPLI'
Gaelle DELAUTIER
committed
1074
1075
1076
1077
1078
1079
1080
1081
1082
1083
1084
1085
1086
1087
1088
1089
1090
1091
1092
1093
1094
LFEEDBACKT = .TRUE.
LEVLIMIT = .TRUE.
LNULLWETG = .TRUE.
LWETGPOST = .TRUE.
LNULLWETH = .TRUE.
LWETHPOST = .TRUE.
CSNOWRIMING = 'M90 '
CSUBG_RC_RR_ACCR = 'NONE'
CSUBG_RR_EVAP = 'NONE'
CSUBG_PR_PDF = 'SIGM'
XFRACM90 = 0.1
LCRFLIMIT = .TRUE.
NMAXITER = 5
XMRSTEP = 0.00005
XTSTEP_TS = 0.
LADJ_BEFORE = .TRUE.
LADJ_AFTER = .TRUE.
CFRAC_ICE_ADJUST = 'S'
XSPLIT_MAXCFL = 0.8
CFRAC_ICE_SHALLOW_MF = 'S'
LSEDIM_AFTER = .FALSE.
LDEPOSC = .FALSE.
XVDEPOSC= 0.02 ! 2 cm/s
1097
1098
1099
1100
1101
1102
1103
1104
1105
1106
1107
1108
1109
1110
1111
1112
1113
1114
1115
1116
1117
1118
1119
1120
1121
1122
1123
1124
1125
1126
1127
1128
1129
1130
1131
1132
1133
1134
1135
1136
1137
1138
1139
1140
1141
1142
1143
1144
1145
1146
1147
1148
END IF
!
!-------------------------------------------------------------------------------
!
!
!* 17. SET DEFAULT VALUES FOR MODD_PARAM_KAFR_n :
! --------------------------------------------
!
XDTCONV = MAX( 300.0,XTSTEP )
NICE = 1
LREFRESH_ALL = .TRUE.
LCHTRANS = .FALSE.
LDOWN = .TRUE.
LSETTADJ = .FALSE.
XTADJD = 3600.
XTADJS = 10800.
LDIAGCONV = .FALSE.
NENSM = 0
!
!-------------------------------------------------------------------------------
!
!
!* 18. SET DEFAULT VALUES FOR MODD_PARAM_MFSHALL_n :
! --------------------------------------------
!
XIMPL_MF = 1.
CMF_UPDRAFT = 'EDKF'
CMF_CLOUD = 'DIRE'
LMIXUV = .TRUE.
LMF_FLX = .FALSE.
!
XALP_PERT = 0.3
XABUO = 1.
XBENTR = 1.
XBDETR = 0.
XCMF = 0.065
XENTR_MF = 0.035
XCRAD_MF = 50.
XENTR_DRY = 0.55
XDETR_DRY = 10.
XDETR_LUP = 1.
XKCF_MF = 2.75
XKRC_MF = 1.
XTAUSIGMF = 600.
XPRES_UV = 0.5
XFRAC_UP_MAX= 0.33
XALPHA_MF = 2.
XSIGMA_MF = 20.
!
XA1 = 2./3.
XB = 0.002
XC = 0.012
XR = 2.
RODIER Quentin
committed
LGZ = .TRUE.
XGZ = 1.83 ! between 1.83 and 1.33
!
!-------------------------------------------------------------------------------
!
!* 19. SET DEFAULT VALUES FOR MODD_PARAM_C2R2 :
! ----------------------------------------
!
IF (KMI == 1) THEN
XNUC = 1.0
XALPHAC = 3.0
XNUR = 2.0
XALPHAR = 1.0
!
Gaelle TANGUY
committed
LRAIN = .TRUE.
LSEDC = .TRUE.
LACTIT = .FALSE.
LSUPSAT = .FALSE.
LDEPOC = .FALSE.
XVDEPOC = 0.02 ! 2 cm/s
1173
1174
1175
1176
1177
1178
1179
1180
1181
1182
1183
1184
1185
1186
1187
1188
1189
1190
1191
1192
1193
!
HPARAM_CCN = 'XXX'
HINI_CCN = 'XXX'
HTYPE_CCN = 'X'
!
XCHEN = 0.0
XKHEN = 0.0
XMUHEN = 0.0
XBETAHEN = 0.0
!
XCONC_CCN = 0.0
XAERDIFF = 0.0
XAERHEIGHT = 2000
XR_MEAN_CCN = 0.0
XLOGSIG_CCN = 0.0
XFSOLUB_CCN = 1.0
XACTEMP_CCN = 280.
ENDIF
!
!-------------------------------------------------------------------------------
!
!* 19.BIS SET DEFAULT VALUES FOR MODD_PARAM_LIMA :
! ----------------------------------------
!
LPTSPLIT = .FALSE.
L_LFEEDBACKT = .TRUE.
L_NMAXITER = 1
L_XMRSTEP = 0.
L_XTSTEP_TS = 0.
!
IF (KMI == 1) THEN
YNUC = 1.0
YALPHAC = 3.0
YNUR = 2.0
YALPHAR = 1.0
!
OWARM = .TRUE.
LACTI = .TRUE.
ORAIN = .TRUE.
!
OVDEPOC = 0.02 ! 2 cm/s
CINI_CCN = 'AER'
CTYPE_CCN(:) = 'M'
!
YAERDIFF = 0.0
YAERHEIGHT = 2000.
! YR_MEAN_CCN = 0.0 ! In case of 'CCN' initialization
! YLOGSIG_CCN = 0.0
YFSOLUB_CCN = 1.0
YACTEMP_CCN = 280.
!
NMOD_CCN = 1
!
!* AP Scavenging
!
LSCAV = .FALSE.
LAERO_MASS = .FALSE.
!
LCCN_HOM = .TRUE.
CCCN_MODES = 'COPT'
XCCN_CONC(:)=300.
ENDIF
!
IF (KMI == 1) THEN
LHHONI = .FALSE.
LCOLD = .TRUE.
LNUCL = .TRUE.
LSNOW = .TRUE.
LHAIL = .FALSE.
CPRISTINE_ICE_LIMA = 'PLAT'
CHEVRIMED_ICE_LIMA = 'GRAU'
XFACTNUC_DEP = 1.0
XFACTNUC_CON = 1.0
NMOD_IFN = 1
NIND_SPECIE = 1
LMEYERS = .FALSE.
LIFN_HOM = .TRUE.
CIFN_SPECIES = 'PHILLIPS'
CINT_MIXING = 'DM2'
XIFN_CONC(:) = 100.
NMOD_IMM = 0
NPHILLIPS=8
ENDIF
!
!-------------------------------------------------------------------------------
!
1266
1267
1268
1269
1270
1271
1272
1273
1274
1275
1276
1277
1278
1279
1280
1281
1282
1283
1284
1285
1286
1287
1288
!* 20. SET DEFAULT VALUES FOR MODD_CH_MNHC_n
! -------------------------------------
!
LUSECHEM = .FALSE.
LUSECHAQ = .FALSE.
LUSECHIC = .FALSE.
LCH_INIT_FIELD = .FALSE.
LCH_CONV_SCAV = .FALSE.
LCH_CONV_LINOX = .FALSE.
LCH_PH = .FALSE.
LCH_RET_ICE = .FALSE.
XCH_PHINIT = 5.2
XRTMIN_AQ = 5.e-8
CCHEM_INPUT_FILE = 'EXSEG1.nam'
CCH_TDISCRETIZATION = 'SPLIT'
NCH_SUBSTEPS = 1
LCH_TUV_ONLINE = .FALSE.
CCH_TUV_LOOKUP = 'PHOTO.TUV39'
CCH_TUV_CLOUDS = 'NONE'
XCH_TUV_ALBNEW = -1.
XCH_TUV_DOBNEW = -1.
XCH_TUV_TUPDATE = 600.
CCH_VEC_METHOD = 'MAX'
XCH_TS1D_TSTEP = 600.
CCH_TS1D_COMMENT = 'no comment'
CCH_TS1D_FILENAME = 'IO1D'
Gaelle DELAUTIER
committed
CSPEC_PRODLOSS = ''
CSPEC_BUDGET = ''
!
!-------------------------------------------------------------------------------
!
!* 21. SET DEFAULT VALUES FOR MODD_SERIES AND MODD_SERIE_n
! ---------------------------------------------------
!
IF (KMI == 1) THEN
LSERIES = .FALSE.
LMASKLANDSEA = .FALSE.
LWMINMAX = .FALSE.
LSURF = .FALSE.
ENDIF
!
NIBOXL = 1 !+ JPHEXT
NIBOXH = 1 !+ 2*JPHEXT
NJBOXL = 1 !+ JPHEXT
NJBOXH = 1 !+ 2*JPHEXT
NKCLS = 1 !+ JPVEXT
NKLOW = 1 !+ JPVEXT
NKMID = 1 !+ JPVEXT
NKUP = 1 !+ JPVEXT
NKCLA = 1 !+ JPVEXT
NJSLICEL(:) = 1 !+ JPHEXT
NJSLICEH(:) = 1 !+ 2*JPHEXT
1320
1321
1322
1323
1324
1325
1326
1327
1328
1329
1330
1331
1332
1333
1334
1335
1336
1337
1338
1339
1340
1341
1342
1343
1344
1345
1346
1347
1348
1349
1350
1351
1352
1353
1354
1355
1356
1357
1358
1359
1360
1361
1362
1363
1364
1365
1366
1367
1368
1369
1370
1371
1372
1373
1374
1375
1376
1377
1378
1379
1380
1381
1382
1383
1384
1385
1386
1387
1388
1389
1390
1391
1392
1393
1394
1395
1396
1397
1398
1399
1400
1401
1402
1403
1404
1405
1406
1407
1408
1409
1410
1411
1412
1413
1414
1415
1416
1417
1418
1419
1420
NFREQSERIES = INT(XSEGLEN /(100.*XTSTEP) )
NFREQSERIES = MAX(NFREQSERIES,1)
!
!-------------------------------------------------------------------------------
!
!* 22. SET DEFAULT VALUES FOR MODD_TURB_CLOUD
! --------------------------------------
!
IF (KMI == 1) THEN
NMODEL_CLOUD = NUNDEF
CTURBLEN_CLOUD = 'DELT'
XCOEF_AMPL_SAT = 5.
XCEI_MIN = 0.001E-06
XCEI_MAX = 0.01E-06
ENDIF
!-------------------------------------------------------------------------------
!
!* 22. SET DEFAULT VALUES FOR MODD_MEAN_FIELD
! --------------------------------------
!
IF (KMI == 1) THEN
LMEAN_FIELD = .FALSE.
ENDIF
!
!-------------------------------------------------------------------------------
!
!* 22. SET DEFAULT VALUES FOR MODD_AEROSOL
! -----------------------------------
IF (KMI == 1) THEN ! other values are defined in modd_ch_aerosol
!
! aerosol lognormal parameterization
LVARSIGI = .FALSE. ! switch to active pronostic dispersion for I mode
LVARSIGJ = .FALSE. ! switch to active pronostic dispersion for J mode
LHETEROSO4 = .FALSE. ! switch to active sulfates heteronegeous
! production
LSEDIMAERO = .FALSE. ! switch to active aerosol sedimentation
LAERINIT = .FALSE. ! switch to initialize aerosol in arome
CMINERAL = "NONE" ! mineral equilibrium scheme
CORGANIC = "NONE" ! mineral equilibrium scheme
CNUCLEATION = "NONE" ! sulfates nucleation scheme
LDEPOS_AER(:) = .FALSE.
ENDIF
!* 23. SET DEFAULT VALUES FOR MODD_DUST and MODD_SALT
! ----------------------------------------------
!
IF (KMI == 1) THEN ! other values initialized in modd_dust
LDUST = .FALSE.
NMODE_DST = 3
LVARSIG = .FALSE.
LSEDIMDUST = .FALSE.
LDEPOS_DST(:) = .FALSE.
LSALT = .FALSE.
LVARSIG_SLT= .FALSE.
LSEDIMSALT = .FALSE.
LDEPOS_SLT(:) = .FALSE.
ENDIF
!
!-------------------------------------------------------------------------------
!
!
!* 24. SET DEFAULT VALUES FOR MODD_PASPOL
! ----------------------------------
!
! other values initialized in modd_paspol
!
IF (KMI == 1) THEN
LPASPOL = .FALSE.
NRELEASE = 0
CPPINIT(:) ='1PT'
XPPLAT(:) = 0.
XPPLON (:) = 0.
XPPMASS(:) = 0.
XPPBOT(:) = 0.
XPPTOP(:) = 0.
CPPT1(:) = "20010921090000"
CPPT2(:) = "20010921090000"
CPPT3(:) = "20010921091500"
CPPT4(:) = "20010921091500"
ENDIF
!
!-------------------------------------------------------------------------------
!
!
!* 25. SET DEFAULT VALUES FOR MODD_CONDSAMP
! ----------------------------------
!
! other values initialized in modd_condsamp
!
IF (KMI == 1) THEN
LCONDSAMP = .FALSE.
NCONDSAMP = 3
XRADIO(:) = 900.
XSCAL(:) = 1.
XHEIGHT_BASE = 100.
XDEPTH_BASE = 100.
XHEIGHT_TOP = 100.
XDEPTH_TOP = 100.
RODIER Quentin
committed
NFINDTOP = 0
XTHVP = 0.25
LTPLUS = .TRUE.
ENDIF
!-------------------------------------------------------------------------------
!
!
!* 26. SET DEFAULT VALUES FOR MODD_LATZ_EDFLX
! ----------------------------------
!
Gaelle Tanguy
committed
IF (KMI == 1) THEN
LUV_FLX=.FALSE.
XUV_FLX1=3.E+14
XUV_FLX2=0.
LTH_FLX=.FALSE.
XTH_FLX=0.75
ENDIF
#ifdef MNH_FOREFIRE
!-------------------------------------------------------------------------------
!
!* 27. SET DEFAULT VALUES FOR MODD_FOREFIRE
! ----------------------------------
!
! other values initialized in modd_forefire
!
IF (KMI == 1) THEN
LFOREFIRE = .FALSE.
LFFCHEM = .FALSE.
COUPLINGRES = 100.
NFFSCALARS = 0
ENDIF
#endif
!-------------------------------------------------------------------------------
!
!* 28. SET DEFAULT VALUES FOR MODD_BLOWSNOW AND MODD_BLOWSNOW_n
! ----------------------------------------
!
IF (KMI == 1) THEN
LBLOWSNOW = .FALSE.
XALPHA_SNOW = 3.
XRSNOW = 4.
CSNOWSEDIM = 'TABC'
END IF
LSNOWSUBL = .FALSE.
!
!
Gaelle DELAUTIER
committed
!-------------------------------------------------------------------------------
!
!* 29. SET DEFAULT VALUES FOR MODD_VISC
! ----------------------------------
!
! other values initialized in modd_VISC
!
IF (KMI == 1) THEN
LVISC = .FALSE.
LVISC_UVW = .FALSE.
LVISC_TH = .FALSE.
LVISC_SV = .FALSE.
LVISC_R = .FALSE.
Gaelle DELAUTIER
committed
XPRANDTL = 0.
ENDIF
!
!-------------------------------------------------------------------------------
!
!
!* 30. SET DEFAULT VALUES FOR MODD_DRAG
! ----------------------------------
!
! other values initialized in modd_DRAG
!
IF (KMI == 1) THEN
LDRAG = .FALSE.
LMOUNT = .FALSE.
NSTART = 1
XHSTART = 0.
ENDIF
!