&parallel_nml
 nproma                       = 32
/

&run_nml
 num_lev                      = 65
 lvert_nest                   = .false.  ! use vertical nesting if a nest is active
 dtime                        = 20       ! timestep in seconds
 ldynamics                    = .true.        ! dynamics
 ltransport                   = .true.
 iforcing                     = 3             ! NWP forcing
 lart                         = .false.        ! ICON-ART main switch
 ltestcase                    = .false.       ! false: run with real data
 msg_level                    = 20            ! default: 5, much more: 20
 ltimer                       = .true.        ! set .true. for detailed timer output
 timers_level                 = 1             ! can be increased up to 10 for detailed timer output
 ntracer                      = 5
 check_uuid_gracefully        = .true.
 output                       = "nml"
 ldass_lhn                    = .False.
 debug_check_level            = 1 
/

&grid_nml
 dynamics_parent_grid_id      = 0
 dynamics_grid_filename       = 'icon_grid_0047_R19B07_L.nc'
 radiation_grid_filename      = 'icon_grid_0046_R19B06_LR.nc',
 lredgrid_phys                = .true.
 l_limited_area               = .true.
/

&nwp_phy_nml
 inwp_psp                     = 1   ! PSP scheme 0(turn off) 1(turn one) 2(turn on with non-D wind pert)
 psp_scale                    = 5.0  ! scaling factor for stochastic perturbations scheme
 psp_rnd_seed                 = 31
 psp_apply_ddt_core = .false. ! True if vertical wind perturbations should be applied at the explicit dynamical core timestep
 inwp_gscp                    = 4   ! 1: default, 2: graupel scheme for convection-permitting scales
 !ccn_type_gscp4              =  8   ! CB new namelist parameter: 
 !ccn_type_gscp5              =  8   ! CB new namelist parameter: 
                                    !6: 100 maritime, 7: 500 intermediate, 8: 1700 continental, 9: 3200 polluted continental
 !cloudnue_value              =  0   ! CB new namelist parameter for shape parameter (nu) of CDSD: 
                                    ! -> cloud_cosmo5 (mu=1/3): 0 (REF), 1, 2, 4, 8 
                                    ! -> cloud_nue1mue1 (mu=1): 1 (REF), 2, 4, 8 
 mu_rain                      = 0.5           ! new tuning becoming operational in July 2018
 rain_n0_factor               = 0.1           !       "
 inwp_convection              = 1
 lshallowconv_only            = .true.
 inwp_radiation               = 1
 inwp_cldcover                = 1   ! 0: no cld, 1: new diagnostic, 3: COSMO, 5: grid scale
 inwp_turb                    = 1
 inwp_satad                   = 1
 inwp_sso                     = 1
 inwp_gwd                     = 0
 inwp_surface                 = 1
 icapdcycl                    = 3   ! CAPE modification, improve diurnal cycle over tropical land
 latm_above_top               = .true.
 efdt_min_raylfric            = 7200.
 itype_z0                     = 2   ! 1: default, 2: turn off SSO part of z0
 icpl_aero_conv               = 1
 icpl_aero_gscp               = 0
 icpl_o3_tp                   = 1
 dt_conv                      = 120
 dt_rad                       = 720
 dt_sso                       = 120
 dt_gwd                       = 120
 ldetrain_conv_prec           = .true.  ! .true. should be used in R03B08 EU-nest only
 lrtm_filename                = 'rrtmg_lw.nc'      
 cldopt_filename              = 'ECHAM6_CldOptProps.nc'
/

&turbdiff_nml
 tkhmin                       = 0.5 ! new default since rev. 16527
 tkmmin                       = 0.75  !           
 tkmmin_strat                 = 4
 pat_len                      = 750.
 c_diff                       = 0.2
 rat_sea                      = 0.8      ! operational since May 2017
 ltkesso                      = .true.
 frcsmot                      = 0.2      ! these 2 switches together apply vertical smoothing of the TKE source terms
 imode_frcsmot                = 2  ! in the tropics (only), which reduces the moist bias in the tropical lower troposphere
 itype_sher                   = 2
 ltkeshs                      = .true.
 a_hshr                       = 2.0
 alpha0                       = 0.0123
 alpha0_max                   = 0.0335        ! Charnock parameter tuning
 alpha1                       = 0.75
 icldm_turb                   = 2  ! ** new recommendation for v2.0.15 in conjunction with evaporation fix for grid-scale rain **
 tur_len                      = 300.
 q_crit                       = 2.0
 imode_tkesso                 = 2
/


&radiation_nml
 irad_o3                      = 79
 irad_aero                    = 6
 albedo_type                  = 2        ! 1: default, 2: MODIS
 vmr_co2                      = 390.e-06 ! values representative for 2012
 vmr_ch4                      = 1800.e-09
 vmr_n2o                      = 322.0e-09
 vmr_o2                       = 0.20946
 vmr_cfc11                    = 240.e-12
 vmr_cfc12                    = 532.e-12
/

&nonhydrostatic_nml
 iadv_rhotheta                = 2
 ivctype                      = 2
 itime_scheme                 = 4
 exner_expol                  = 0.333
 vwind_offctr                 = 0.2
 damp_height                  = 12500.
 rayleigh_coeff               = 5.0
 lhdiff_rcf                   = .true.
 divdamp_fac                  = 0.004         ! factor 8 iff divdamp_order=2 !!!
 divdamp_order                = 24       ! use 24 for forecast, 2 for IAU
 divdamp_type                 = 32        ! (3): 2 for assimilation cycle
 divdamp_trans_start          = 12500. ! use 2500. in assimilation cycle
 divdamp_trans_end            = 17500.   ! use 5000. in assimilation cycle
 ndyn_substeps                = 5
 l_open_ubc                   = .false.
 igradp_method                = 3
 l_zdiffu_t                   = .true.
 thslp_zdiffu                 = 0.02
 thhgtd_zdiffu                = 125.
 htop_moist_proc              = 22500.
 hbot_qvsubstep               = 22500   ! revised (to avoid too large CFL)
/

&sleve_nml 
 min_lay_thckn                = 20. 
 max_lay_thckn                = 25000.
 htop_thcknlimit              = 15000.  ! 15 km is the default
 top_height                   = 22000. 
 stretch_fac                  = 0.65
 decay_scale_1                = 4000.
 decay_scale_2                = 2500. 
 decay_exp                    = 1.2 
 flat_height                  = 16000.
/

&dynamics_nml
 iequations                   = 3
 idiv_method                  = 1
 divavg_cntrwgt               = 0.50
 lcoriolis                    = .true.
/

&transport_nml
 ctracer_list                 = '12345'
 ivadv_tracer                 = 3,3,3,3,3
 itype_hlimit                 = 3,4,4,4,4
 ihadv_tracer                 = 52,2,2,2,2
/

&diffusion_nml
 hdiff_order                  = 5
 itype_vn_diffu               = 1
 itype_t_diffu                = 2
 hdiff_efdt_ratio             = 24.     ! 24.0  for R2B6; recommendation for R3B7: 30.0
 hdiff_smag_fac               = 0.025   ! 0.025 for R2B6; recommendation for R3B7: 0.02
 lhdiff_vn                    = .true.
 lhdiff_temp                  = .true.
/

&interpol_nml
 nudge_zone_width             = 10
 nudge_max_coeff              = 0.075
 lsq_high_ord                 = 3
 l_intp_c2l                   = .true.
 l_mono_c2l                   = .true.
 support_baryctr_intp         = .false.
/

&gridref_nml
 grf_intmethod_e              = 6
 grf_intmethod_ct             = 2
 grf_tracfbk                  = 2
 denom_diffu_v                = 150.
/

&nwp_tuning_nml
 itune_albedo                 = 1     ! somewhat reduced albedo (w.r.t. MODIS data) over Sahara
 tune_gkdrag                  = 0.0
 tune_gkwake                  = 0.25
 tune_gfrcrit                 = 0.333
 tune_gust_factor             = 7.25
 tune_dust_abs                = 1.
 tune_zvz0i                   = 1.25    ! new tuning since June 2018
 tune_box_liq_asy             = 3.5  ! new tuning becoming operational in July 2018
 tune_minsnowfrac             = 0.3
 tune_sgsclifac               = 1.0    ! new tuning becoming operational in July 2019
/

&extpar_nml
! extpar_filename              = 'icon_extpar_0044_R19B07_L_20180625_tiles.nc'
 itype_vegetation_cycle       = 1  !operational using 2(ihan)
 itopo                        = 1
 n_iter_smooth_topo           = 1,1
 heightdiff_threshold         = 2250.
 hgtdiff_max_smooth_topo      = 750.,750.,
 read_nc_via_cdi              = .true.
/

&lnd_nml
 ntiles                       = 3
 lmulti_snow                  = .false.
 itype_heatcond               = 3
 idiag_snowfrac               = 20
 itype_snowevap               = 2 
 lsnowtile                    = .true.  !.false.
 lseaice                      = .true.
 llake                        = .true.
 lprog_albsi                  = .true.
 itype_lndtbl                 = 4
 itype_root                   = 2
 itype_evsl                   = 4
 itype_trvg                   = 3 ! operational since 14 March 2018, requires FG cycling of EVAP_PL (internal name plantevap)
 cwimax_ml                    = 5.e-4
 c_soil                       = 1.25
 c_soil_urb                   = 0.5
 sstice_mode                  = 2        ! requires extpar >= 20170202
/

&limarea_nml
 itype_latbc                  = 1
 latbc_filename               = './icbc/M031/latbc_<dddhh>.grib'
 latbc_boundary_grid          = 'icon_grid_0047_R19B07_L_lbc.nc'
 dtime_latbc                  = 10800 !3600
 init_latbc_from_fg           = .true.
 latbc_varnames_map_file      = 'map_file.latbc'
/

&initicon_nml
 init_mode                    = 5 ! 1: operation mode, 5: IAU, 6: 5+sma_inc
 dt_iau                       = 600    ! Window for incr.analysis update
 dt_shift                     = -300  ! Offset for incr.analysis update
 type_iau_wgt                 = 1           ! IAU weighting function (const.)
 iterate_iau                  = .true.      ! iterate IAU cycle with one-sided and centered windows
 zpbl1                        = 500.
 zpbl2                        = 1000.
 lread_ana                    = .true.        ! (T) Read dwdana
 lp2cintp_incr                = .false.    ! parent-to-child-interpolation
 lp2cintp_sfcana              = .false.  ! of analyses (.false. for DET!)
 dwdfg_filename              = 'icbc/M031/fc.grib'
! lread_ana                   = .False.
 dwdana_filename             = 'icbc/M031/an_inc.grib'
 ana_varnames_map_file = 'map_file.ana'
! check_ana(1)%list            = 'P','QV','T','U','V'
! check_ana(2)%list            = 'FRESHSNW','H_SNOW'
 ltile_coldstart              = .false.  ! intermediate mode until tile-based I/O is available
 ltile_init                   = .false.
 use_lakeiceana               = .false. 
 qcana_mode                   = 2           ! (0) no QC-increments,
 qiana_mode                   = 1         ! 0/1: ignore/use QI increments
 qrsgana_mode                 = 1       ! 0/1: ignore/use QR,QS,QG increments
/

&ensemble_pert_nml
 use_ensemble_pert            = .true.
 itype_pert_gen               = 2
 timedep_pert                 = 2 
/


 ! ----------------------------------------------- !
 ! ---  ICON-D2: output fields - regular grid at pressure levels --- !
 ! -----------------------------------------------
&output_nml
 filetype                    = 4
 output_time_unit            = 1      ! 1: seconds
 output_interval             = 'PT3600S'
 output_start                = '#OUTPUT_START#'
 output_end                  = '#OUTPUT_END#'
 steps_per_file       = 1
 taxis_tunit          = 2
 mode                 = 1  ! 1: forecast mode (relative t-axis), 2: climate mode (absolute t-axis)
 include_last         = .true.
 steps_per_file_inclfirst = .false.
 output_filename      = 'idfrf'                      ! file name base
 filename_format      = '<output_filename><ddhhmmss><levtype>'
 remap = 1,
 reg_def_mode                 = 1
 north_pole                   = -170.0,40.0
 reg_lon_def                  = -7.5, 0.02, 5.5
 reg_lat_def                  = -6.3, 0.02, 8.0
 pl_varlist                   = 'geopot','temp','U','V','W','RH','TKE','QV','QC','QR','RHO'
 p_levels  = 95000., 85000., 50000., 25000.,
 output_grid          = .false.
/


 ! ----------------------------------------------- !
 ! ---  ICON-D2: output fields - regular grid at model levels --- !
 ! ----------------------------------------------- 
&output_nml
 filetype                    = 4
 output_time_unit            = 1      ! 1: seconds
 output_interval             = 'PT3600S'
 output_start                = '#OUTPUT_START#'
 output_end                  = '#OUTPUT_END#'
 steps_per_file       = 1
 taxis_tunit          = 2
 mode                 = 1  ! 1: forecast mode (relative t-axis), 2: climate mode (absolute t-axis)
 steps_per_file_inclfirst = .false.
 include_last         = .true.
 output_filename      = 'idfrf'                      ! file name base
 filename_format      = '<output_filename><ddhhmmss><levtype>'
 remap = 1,
 reg_def_mode                 = 1
 north_pole                   = -170.0,40.0
 reg_lon_def                  = -7.5, 0.02, 5.5
 reg_lat_def                  = -6.3, 0.02, 8.0
 ml_varlist                   = 'CAPE_ML', 'CIN_ML','tqc','tqv','tqr','tqi','tqs','tqg','tqh','pres_msl','pres_sfc','T_2m','U_10m','V_10m','TOT_PREC','CLDEPTH','CLCT_MOD','HZEROCL','TMIN_2M','TMAX_2M','dbz_cmax', !,'group:pbl_vars'
 output_grid          = .false.
/


 ! ----------------------------------------------- !
 ! ---  ICON-D2: PSP output fields - regular grid at model levels --- !
 ! ----------------------------------------------- 
&output_nml
 filetype                    = 4
 output_time_unit            = 1      ! 1: seconds
 output_interval             = 'PT3600S'
 output_start                = '#OUTPUT_START#'
 output_end                  = '#OUTPUT_END#'
 steps_per_file       = 1
 taxis_tunit          = 2
 mode                 = 1  ! 1: forecast mode (relative t-axis), 2: climate mode (absolute t-axis)
 steps_per_file_inclfirst = .false.
 include_last         = .true.
 output_filename      = 'idfrf'                      ! file name base
 filename_format      = '<output_filename><ddhhmmss><levtype>_psp'
 remap = 1,
 reg_def_mode                 = 1
 north_pole                   = -170.0,40.0
 reg_lon_def                  = -7.5, 0.02, 5.5
 reg_lat_def                  = -6.3, 0.02, 8.0
! ml_varlist                   = 'group:psp',
 output_grid          = .false.
/