!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 _ !----------------------------------------------------------------- ! ########################### 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_FMOUT : XOUT !! !! 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 !! !! Module MODD_LBC_n : CLBCX, CLBCY,NLBLX,NLBLY,XCPHASE,XCPHASE_PBL,XPOND !! !! 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 !! 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 !! 07/2013 (C.Lac) add WENO, LCHECK !! 07/2013 (Bosseur & Filippi) adds Forefire !! 08/2015 (Redelsperger & Pianezze) add XPOND coefficient for LBC !------------------------------------------------------------------------------- ! !* 0. DECLARATIONS ! ------------ ! USE MODD_IBM_PARAM_n !===>FA USE MODD_EMIS_PARAM_n !===>FA ! USE MODD_PARAMETERS USE MODD_CONF ! For INIT only DEFAULT_DESFM1 USE MODD_CONFZ USE MODD_DYN USE MODD_NESTING USE MODD_FMOUT 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 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_LATZ_EDFLX USE MODD_2D_FRC #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 : ! ---------------------------------- ! CINIFILE='INIFILE' CINIFILEPGD='' 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 = 3 #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. LCHECK = .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. LNUMDIFU = .FALSE. 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 ! !------------------------------------------------------------------------------- ! !* 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 = 'WENO_K' CMET_ADV_SCHEME = 'PPM_01' CSV_ADV_SCHEME = 'PPM_01' CTEMP_SCHEME = 'RK53' NWENO_ORDER = 3 NSPLIT = 1 LSPLIT_CFL = .TRUE. XSPLIT_CFL = 1.0 LCFL_WRIT = .FALSE. ! !------------------------------------------------------------------------------- ! !* 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. XCPHASE_PBL = 0. XCARPKMAX = XUNDEF XPOND = 1.0 ! !------------------------------------------------------------------------------- ! !* 8. SET DEFAULT VALUES FOR MODD_NUDGING_n : ! --------------------------------- ! LNUDGING = .FALSE. XTNUDGING = 21600. ! !------------------------------------------------------------------------------- ! !* 9. SET DEFAULT VALUES FOR MODD_FMOUT and MODD_OUT_n : ! ------------------------------------------------ ! IF (KMI == 1) XFMOUT (:,:) = XUNDEF ! ! !------------------------------------------------------------------------------- ! !* 10. SET DEFAULT VALUES FOR MODD_TURB_n : ! ---------------------------------- ! XIMPL = 1. XKEMIN = 0.01 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' VSIGQSAT = 0. !------------------------------------------------------------------------------- ! !* 10b. SET DEFAULT VALUES FOR MODD_DRAGTREE : ! ---------------------------------- ! LDRAGTREE = .FALSE. ! !------------------------------------------------------------------------------- ! !* 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 NADVU = 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 NADVV = 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 NADVW = 0 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 NGMLTTH = 0 NIMLTTH = 0 NBERFITH = 0 NCDEPITH = 0 NWETHTH = 0 NHMLTTH = 0 ! ! 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 ! ! 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 NWETGRC = 0 NDRYGRC = 0 NIMLTRC = 0 NBERFIRC = 0 NCDEPIRC = 0 NWETHRC = 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 NCFRZRR = 0 NWETGRR = 0 NDRYGRR = 0 NGMLTRR = 0 NWETHRR = 0 NHMLTRR = 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 ! ! 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 ! ! 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 NCOHGRG = 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 NCOHGRH = 0 NHMLTRH = 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 NNEGASV = 0 ! ! 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. ! !------------------------------------------------------------------------------- ! !* 14. SET DEFAULT VALUES FOR MODD_BLANK : ! ----------------------------------- ! IF (KMI == 1) THEN XDUMMY1 = 0. XDUMMY2 = 0. XDUMMY3 = 0. 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 LWARM = .TRUE. CPRISTINE_ICE = 'PLAT' LSEDIC = .FALSE. LCONVHG = .FALSE. CSEDIM = 'SPLI' 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 XBETA1 = 0.9 ! !------------------------------------------------------------------------------- ! !* 19. SET DEFAULT VALUES FOR MODD_PARAM_C2R2 : ! ---------------------------------------- ! IF (KMI == 1) THEN XNUC = 1.0 XALPHAC = 3.0 XNUR = 2.0 XALPHAR = 1.0 ! LRAIN = .TRUE. LSEDC = .TRUE. LACTIT = .FALSE. LSUPSAT= .FALSE. ! 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 ! !------------------------------------------------------------------------------- ! !* 20. SET DEFAULT VALUES FOR MODD_CH_MNHC_n ! ------------------------------------- ! LUSECHEM = .FALSE. LUSECHAQ = .FALSE. LUSECHIC = .FALSE. LCH_INIT_FIELD = .FALSE. LCH_SURFACE_FLUX = .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' NCH_VEC_LENGTH = 1000 XCH_TS1D_TSTEP = 600. CCH_TS1D_COMMENT = 'no comment' CCH_TS1D_FILENAME = 'IO1D' ! !------------------------------------------------------------------------------- ! !* 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 NBJSLICE = 1 NJSLICEL(:) = 1 !+ JPHEXT NJSLICEH(:) = 1 !+ 2*JPHEXT 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. NMODE_SLT = 3 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. ENDIF !------------------------------------------------------------------------------- ! ! !* 26. SET DEFAULT VALUES FOR MODD_LATZ_EDFLX ! ---------------------------------- ! 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 !------------------------------------------------------------------------------- !------------------------------------------------------------------------------- !===>FA ! ! !* 28. SET DEFAULT VALUES FOR MODD_EMIS_PARAMn ! -------------------------------------- IF (KMI == 1) THEN LEMIS = .FALSE. NEMIS_FILE = 0 NEMIS_HOUR = 0 XEMIS_RESI = 1e-3 ENDIF !------------------------------------------------------------------------------- !------------------------------------------------------------------------------- !===>FA ! ! !* 28. SET DEFAULT VALUES FOR MODD_IBM_PARAMn ! -------------------------------------- ! IF (KMI == 1) THEN LIBM = .FALSE. LIBM_RLS = .FALSE. LRECYCL = .FALSE. CIBM_TYPE = "NONE" XIBM_EPSI = 1.E-9 XIBM_IEPS = 1.E+9 XIBM_SMOOTH = 1.E-9 XIBM_VIT = 1. XXRECYCLB1 = 0. XXRECYCLE1 = 0. XYRECYCLB1 = 0. XYRECYCLE1 = 0. XXRECYCLB2 = 0. XXRECYCLE2 = 0. XYRECYCLB2 = 0. XYRECYCLE2 = 0. XDRECYCL = 0. XARECYCL = 0. XIBM_RUG = 0.01 XIBM_SFTH = 0. XIBM_RES = 1.E-6 XIBM_RES2 = 1.E-6 NIBM_ITR = 3 NIBM_SMOOTH = 1 XIBM_IDEA_DELT = 1. XIBM_TKE =+1. XIBM_VISC_U=-1. XIBM_VISC_E=-1. XIBM_VISC_T=-1. XIBM_VISC_S=-1. XIBM_VISC_R=-1. XIBM_LATORI = 60. XIBM_LONORI = 0. NIBM_LAYER_P = 2 NIBM_LAYER_Q = 2 NIBM_LAYER_R = 2 NIBM_LAYER_S = 2 NIBM_LAYER_T = 2 NIBM_LAYER_E = 2 NIBM_LAYER_V = 2 XIBM_RADIUS_P = 2. XIBM_RADIUS_Q = 2. XIBM_RADIUS_R = 2. XIBM_RADIUS_S = 2. XIBM_RADIUS_T = 2. XIBM_RADIUS_E = 2. XIBM_RADIUS_V = 2. XIBM_POWERS_P = 1. XIBM_POWERS_Q = 1. XIBM_POWERS_R = 1. XIBM_POWERS_S = 1. XIBM_POWERS_T = 1. XIBM_POWERS_E = 1. XIBM_POWERS_V = 1. CIBM_MODE_INTE3_P = 'LAI' CIBM_MODE_INTE3_Q = 'LAI' CIBM_MODE_INTE3_R = 'LAI' CIBM_MODE_INTE3_S = 'LAI' CIBM_MODE_INTE3_T = 'LAI' CIBM_MODE_INTE3_E = 'LAI' CIBM_MODE_INTE3_V = 'LAI' CIBM_MODE_INTE1_P = 'CL2' CIBM_MODE_INTE1_Q = 'CL2' CIBM_MODE_INTE1_R = 'CL2' CIBM_MODE_INTE1_S = 'CL2' CIBM_MODE_INTE1_T = 'CL2' CIBM_MODE_INTE1_E = 'CL2' CIBM_MODE_INTE1NV = 'CL2' CIBM_MODE_INTE1TV = 'CL2' CIBM_MODE_INTE1CV = 'CL2' CIBM_MODE_BOUND_P = 'ASY' CIBM_MODE_BOUND_Q = 'ASY' CIBM_MODE_BOUND_R = 'ASY' CIBM_MODE_BOUND_S = 'ASY' CIBM_MODE_BOUND_T = 'ASY' CIBM_MODE_BOUND_E = 'ASY' CIBM_MODE_BOUNT_V = 'ASY' CIBM_MODE_BOUNN_V = 'ASY' CIBM_MODE_BOUNC_V = 'ASY' XIBM_FORC_BOUND_P = 0. XIBM_FORC_BOUND_Q = 0. XIBM_FORC_BOUND_R = 0. XIBM_FORC_BOUND_S = 0. XIBM_FORC_BOUND_T = 0. XIBM_FORC_BOUND_E = 0. XIBM_FORC_BOUNN_V = 0. XIBM_FORC_BOUNT_V = 0. XIBM_FORC_BOUNC_V = 0. CIBM_TYPE_BOUND_P = 'NEU' CIBM_TYPE_BOUND_Q = 'NEU' CIBM_TYPE_BOUND_R = 'NEU' CIBM_TYPE_BOUND_S = 'NEU' CIBM_TYPE_BOUND_T = 'NEU' CIBM_TYPE_BOUND_E = 'NEU' CIBM_TYPE_BOUNT_V = 'NEU' CIBM_TYPE_BOUNN_V = 'DIR' CIBM_TYPE_BOUNC_V = 'DIR' CIBM_FORC_BOUND_P = 'CST' CIBM_FORC_BOUND_Q = 'CST' CIBM_FORC_BOUND_R = 'CST' CIBM_FORC_BOUND_S = 'CST' CIBM_FORC_BOUND_T = 'CST' CIBM_FORC_BOUND_E = 'CST' CIBM_FORC_BOUNN_V = 'CST' CIBM_FORC_BOUNT_V = 'CST' CIBM_FORC_BOUNC_V = 'CST' CIBM_FORC_BOUNR_V = 'CST' LIBM_SOLAR = .FALSE. XIBM_SOLAR_MAX = 20.0 XIBM_SOLAR_EXP = 3.00 XIBM_SOLAR_ERR = 0.10 XIBM_SOLAR_REF = 348. XIBM_SOLAR_LBF = 0.01 XIBM_SOLAR_LBG = 1.00 XIBM_SOLAR_DZG = 1.00 XIBM_SOLAR_DZR = 1.00 XIBM_SOLAR_CPG = 1e+6 XIBM_SOLAR_EMI = 1.00 ENDIF ! END SUBROUTINE DEFAULT_DESFM_n