!---------------------------------- ! ICON standard namelist parameters (NAMELIST_icon_1) !---------------------------------- ! For a complete list see Namelist_overview and Namelist_overview.pdf !----------------------- ! ICON coupling_mode_nml !----------------------- &coupling_mode_nml coupled_mode = .FALSE. / !------------------ ! ICON parallel_nml !------------------ ¶llel_nml nproma = 760 ! SX-9: 1024, Cray: 8, Aurora: 760 p_test_run = .false. num_io_procs = 2 ! 0: no dedicated I/O processor ! num_io_procs_radar = 5 ! 0: no dedicated I/O processor num_restart_procs = 0 itype_comm = 1 iorder_sendrecv = 3 ! (1), 3: best value for CRAY proc0_shift = 1 use_omp_input = .true. l_log_checks = .true. l_test_openmp = .true. num_dist_array_replicas = 4 / !-------------- ! ICON grid_nml !-------------- &grid_nml dynamics_grid_filename = 'icon_grid_9999_R13B07_L.nc', radiation_grid_filename = 'icon_grid_9998_R13B06_LR.nc', dynamics_parent_grid_id = 0,1 lredgrid_phys = .true.,.true. lfeedback = .true. ifeedback_type = 2 l_limited_area = .true. nexlevs_rrg_vnest = 8 start_time = 0, 9999999 end_time = 0, 432000 ! stop nest after 120 h / !------------------ ! ICON initicon_nml !------------------ &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 niter_diffu = 10, niter_divdamp = 25, zpbl1 = 500. zpbl2 = 1000. lread_ana = .true. ! (T) Read dwdana lp2cintp_incr = .false. ! parent-to-child-interpolation lp2cintp_sfcana = .false. ! of analyses (DET:.false.;EPS:.true.) dwdfg_filename = "fg_RB_DOM.grb" dwdana_filename = "an_RB_DOM.grb" ifs2icon_filename = "ifs2icon_RB.nc" ana_varnames_map_file = "/hpc/uhome/extjroqu/bacy_data/const/ana_varnames_map_file.txt" 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 = .true. qcana_mode = 2 ! (0) no QC-increments, ! (1) QC+QV-increments added to QV, ! (2) QC-increments added to QC if already present at gridpoint ! and added to QV if not qiana_mode = 0 ! 0/1: ignore/use QI increments qrsgana_mode = 0 ! 0/1: ignore/use QR,QS,QG increments ! qnxana_2mom_mode = 0 ! 0: diagnose from mass 1: use from an_inc icpl_da_sfcevap = 4 ! use filtered T2m bias and filtered RH increment ... icpl_da_skinc = 2 icpl_da_snowalb = 2 icpl_da_sfcfric = 0 icpl_da_tkhmin = 1 adjust_tso_tsnow = .true. dt_ana = 3600 ! ... to control evaporation relevant parameters / !----------------- ! ICON limarea_nml !----------------- &limarea_nml itype_latbc = 1 nudge_hydro_pres = .true. dtime_latbc = 3600. latbc_path = './' ! latbc_boundary_grid = 'icon_grid_9999_R13B07_L_lbc.nc' latbc_varnames_map_file = "/hpc/uhome/extjroqu/bacy_data/const/ana_varnames_map_file.txt" latbc_filename = 'ilbff' ! latbc_contains_qcqi=.false. init_latbc_from_fg = .true. fac_latbc_presbiascor = 1. / !----------------- ! ICON nudging_nml !----------------- &nudging_nml nudge_type = 0 ! Hauptschalter Nudging Type ! nudge_var = "vn,thermdyn" ! max_nudge_coeff_vn=0.025 ! relaxation time ~1 hours ! max_nudge_coeff_thermdyn=0.0125 ! relaxation time ~2 hours ! nudge_start_height=2000 ! nudge_end_height=80000 ! nudge_scale_height=3500 / !------------- ! ICON run_nml !------------- &run_nml num_lev = 75 lvert_nest = .true. ! use vertical nesting if a nest is active nsteps = 527040 dtime = 30 ! timestep in seconds ldynamics = .true. ! dynamics ltransport = .true. iforcing = 3 ! 3: NWP forcing; 0: no forcing lart = .false. ! ICON-ART main switch ltestcase = .false. ! false: run with real data msg_level = 7 ! default: 5, much more: 20 ltimer = .true. ! set .true. for detailed timer output timers_level = 5 ! can be increased up to 10 for detailed timer output !activate_sync_timers = .false. ntracer = 5 check_uuid_gracefully = .false. output = "nml" ldass_lhn = .FALSE. ! restart_filename = "" ! luse_radarfwo = .FALSE., .FALSE. / !----------------- ! ICON nwp_phy_nml !----------------- &nwp_phy_nml inwp_gscp = 2 ! 1: default, 2: graupel scheme for convection-permitting scales mu_rain = 0.5 ! new tuning becoming operational in July 2018 rain_n0_factor = 0.1 inwp_convection = 1 lshallowconv_only = .true.,.true. lgrayzone_deepconv = .false.,.false. inwp_radiation = 4 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.,.true. efdt_min_raylfric = 7200. itype_z0 = 2 ! 1: default, 2: turn off SSO part of z0 icpl_aero_conv = 1 icpl_aero_gscp = 1 icalc_reff = 0 ! 0 for 1Mom icpl_rad_reff = 0 ! 0 for 1Mom ! ithermo_water = 0 ! 0 for 1Mom icpl_o3_tp = 1 dt_conv = 120 dt_rad = 720 dt_sso = 120 dt_gwd = 120 ldetrain_conv_prec = .true.,.true. ! .true. should be used in R03B08 EU-nest only lrtm_filename = 'rrtmg_lw.nc' ! 'rrtmg_lw.nc' ! not needed for ECRAD (inwp_radiation=4), entries can be removed cldopt_filename = 'rrtm_cldopt.nc' ! 'ECHAM6_CldOptProps.nc' ! not needed for ECRAD (inwp_radiation=4), entries can be removed / !-------------------- ! ICON nwp_tuning_nml !-------------------- &nwp_tuning_nml tune_eiscrit = 7. tune_gkdrag = 0.0 tune_gkwake = 0.25 tune_gfrcrit = 0.333 tune_grcrit = 0.5 ! only for tune_gkdrag != 0.0 tune_gustsso_lim = 20. tune_gust_factor = 7.0 tune_rcucov = 0.075 tune_rhebc_land = 0.825 tune_zvz0i = 0.85 ! new tuning since June 2018 tune_box_liq = 0.04 tune_box_liq_asy = 4.0 ! new tuning becoming operational in July 2018 tune_minsnowfrac = 0.3 tune_sgsclifac = 1.0 ! new tuning becoming operational in July 2019 tune_gfluxlaun = 3.75e-3 icpl_turb_clc = 2 ! tune_minsso = 1 ! tune_blockred = 1.5 itune_gust_diag = 3 ! tune_zceff_min = 0.025 ! lcalib_clcov = .false. max_calibfac_clcl = 2 / !----------------------- ! ICON ensemble_pert_nml !----------------------- &ensemble_pert_nml use_ensemble_pert = .true. ! stdev_sst_pert = 0.0 !range_gkwake = 0.1 ! use default 1.5 with icon-2.6.4-nwp4 range_gkdrag = 0.0 range_gfluxlaun = 0.0 range_capdcfac_et = 0.0 range_tkhmin = 0.1 range_tkmmin = 0.15 range_rlam_heat = 2.5 range_rhebc = 0.0 range_minsnowfrac = 0.0 range_z0_lcc = 0.0 range_rootdp = 0.0 range_rsmin = 0.0 range_laimax = 0.0 range_charnock = 1.0 range_c_soil = 0.0 range_cwimax_ml = 1.0 range_capdcfac_tr = 0.0 range_lowcapefac = 0.0 range_negpblcape = 0.0 range_q_crit = 1.5 range_turlen = 50.0 range_rain_n0fac = 2.0 range_rprcon = 0.0 range_thicklayfac = 0.0 itype_pert_gen = 2 timedep_pert = 1 / !------------------ ! ICON turbdiff_nml !------------------ &turbdiff_nml tkhmin = 0.5 ! new default since rev. 16527 tkmmin = 0.75 ! " tkhmin_strat = 0.75 tkmmin_strat = 4 pat_len = 750. c_diff = 0.2 rlam_heat = 10.0 rat_sea = 0.8 ! operational since May 2017 ltkesso = .true. frcsmot = 0 ! these 2 switches together apply vertical smoothing of the TKE source terms imode_charpar =3 ! in the tropics (only), which reduces the moist bias in the tropical lower troposphere ! use horizontal shear production terms with 1/SQRT(Ri) scaling to prevent unwanted side effects: itype_sher = 2 ltkeshs = .true. a_hshr = 2.0 alpha0 = 0.0123 alpha0_max = 0.0335 ! Charnock parameter tuning alpha1 = 0.125 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 / !------------- ! ICON lnd_nml !------------- &lnd_nml ntiles = 3 lmulti_snow = .false. itype_heatcond = 3 idiag_snowfrac = 20 itype_snowevap = 3 ! operational since late summer 2018, requires cycling of snow_age and hsnow_max 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 itype_canopy = 2 ! operational since May 2020, skin-temperature scheme ! hice_min = 0.05 ! default of nwp sea-ice model ! hice_max = 4.8 ! corresponding to 40% seaice_limit in icon-o / !------------------- ! ICON radiation_nml !------------------- &radiation_nml irad_o3 = 79 irad_aero = 6 izenith = 4 albedo_type = 2 ! 1: default, 2: MODIS direct_albedo = 4 direct_albedo_water = 3 albedo_whitecap = 1 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 ecrad_llw_cloud_scat = .false. / !------------------------ ! ICON nonhydrostatic_nml !------------------------ &nonhydrostatic_nml iadv_rhotheta = 2 ivctype = 2 itime_scheme = 4 exner_expol = 0.6 vwind_offctr = 0.2 damp_height = 12250. 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) / !--------------- ! ICON sleve_nml !--------------- &sleve_nml min_lay_thckn = 20. ! lowest level thickness (between half-levels) max_lay_thckn = 25000. ! maximum layer thickness below htop_thcknlimit 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. / !------------------ ! ICON dynamics_nml !------------------ &dynamics_nml iequations = 3 idiv_method = 1 divavg_cntrwgt = 0.50 lcoriolis = .true. / !------------------- ! ICON transport_nml !------------------- &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 ! iadv_tke = 0 llsq_svd = .true. / !------------------- ! ICON diffusion_nml !------------------- &diffusion_nml hdiff_order = 5 itype_vn_diffu = 1 itype_t_diffu = 2 hdiff_efdt_ratio = 24. hdiff_smag_fac = 0.025 lhdiff_vn = .true. lhdiff_temp = .true. hdiff_multfac = 1.0 hdiff_tv_ratio = 1.0 / !------------------ ! ICON interpol_nml !------------------ &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. / !----------------- ! ICON gridref_nml !----------------- &gridref_nml grf_intmethod_e = 6 grf_intmethod_ct = 2 grf_tracfbk = 2 denom_diffu_v = 150. / !---------------- ! ICON extpar_nml !---------------- &extpar_nml itopo = 1 itype_lwemiss = 2 itype_vegetation_cycle = 2 n_iter_smooth_topo = 1,1 heightdiff_threshold = 2250., 1500. pp_sso = 1 hgtdiff_max_smooth_topo = 750.,750., read_nc_via_cdi = .true. / !------------ ! ICON io_nml !------------ &io_nml ! lflux_avg = .FALSE. itype_pres_msl = 5 ! (1) 4: IFS method with consistency fix dt_checkpoint=864000 ! JL itype_rh = 1 ! (1) 2: mixed phase (water and ice) ! inextra_3d = 2 ! 3D extra variables ! inextra_2d = 10 ! 2D extra variables ! restart_write_mode = "sync" ! write mode of restart-files gust_interval = 3600. !lzaxis_reference = .true. ! (F) use REFERENCE instead of HYBRID lmask_boundary = .true. echotop_meta(1)%time_interval = 900.0 ! Reset der Druckminimumbestimmung auf domain 1 alle 15' echotop_meta(1)%dbzthresh = 18. ! Liste der dbz-thresholds (hier nur einer) fuer domain 1 / !----------------- ! ICON gribout_nml !----------------- &gribout_nml productionStatusOfProcessedData = 2 ! Research products backgroundProcess = 0 ! 0=main,2=assimilation generatingProcessIdentifier = 1,2,3 ! 1=icogl,2=icoeu,3=icode preset = 'deterministic' ! deterministic or ensemble !typeOfProcessedData = 1 ! [WMO Table 1.4] forecast !typeOfGeneratingProcess = 2 ! [WMO Table 4.3] forecast !localDefinitionNumber = 254 ! [DWD] system, 254=det,253=eps localNumberOfExperiment = 1 ! [DWD] Experiment ID generatingcenter = 78 ! [DWD] generatingsubcenter = 255 ! [DWD] tablesVersion = 11 lgribout_24bit = .false. ! (F) write rho,theta,T,p with 24bits lspecialdate_invar = .false. / !-------------------------- ! ICON meteogram_output_nml !-------------------------- &meteogram_output_nml lmeteogram_enabled= .false. !.true. n0_mtgrm = 0 ! meteogram initial time step (0 is first step!) ninc_mtgrm = 300 ! meteogram output interval (in terms of time steps) ldistributed = .false. stationlist_tot = 52.17, 14.12, 'Lindenberg', 51.97, 4.93, 'Cabauw', -10.08, -61.93, 'LBA_Rondonia', 13.50, 2.5 , 'Niamey', 36.61, -97.49, 'ARM_Southern_Great_Plains', -71.32, 156.62, 'ARM_North_Slope_of_Alaska_Barrow', -2.06, 147.43, 'ARM_Tropical_W_Pacific_Manus', -12.43, 130.89, 'ARM_Tropical_W_Pacific_Darwin', 60.00, 80.00, 'Snow Test Russia', ! stationlist_tot = 50.0, 8.6, 'Frankfurt-Flughafen', ! 52.260, 14.313, 'Lindenberg_Obs', ! 52.1873, 13.9894, 'Waldstation_Kehrigk', ! 52.1215, 14.1077, 'Falkenberg', / !----------------------- ! ICON: assimilation_nml !----------------------- &assimilation_nml ! dace_coupling = .false. !# JL ! dace_time_ctrl = 0,3600,1800 ! 60 min. interval ! dace_debug = 1 ! llhn = .false.,.true. ! llhnverif = .false.,.true. nlhn_start = 0 nlhnverif_start = 0 nlhn_end = 20000 nlhnverif_end = 50000 lhn_diag = .true. lhn_hum_adj = .true. lhn_no_ttend = .false. lhn_artif = .true. radardata_file = 'radardata_DOM1.nc','radardata_DOM2.nc', tt_artif_max = 0.009 zlev_artif_max = 3000. lhn_artif_only = .false. lhn_incloud = .true. lhn_filt = .true. lhn_limit = .true. lhn_relax = .true. nlhn_relax = 2 fac_lhn_up = 1.1 fac_lhn_down = 0.7 rqrsgmax = 0.01 lhn_coef = 0.0 start_fadeout = 1.0 lhn_qrs = .true. lhn_logscale = .true. lhn_wweight = .true. nradar = 12 lhn_black = .true. lhn_bright = .true. rttend = 0.0 abs_lhn_lim = 0.0035 blacklist_file = 'blacklist_dx_ilam.nc','blacklist_dx_ilam_nest.nc' height_file = 'radarheight_ilam.nc','radarheight_ilam_nest.nc' / !------------------------------------- ! dummy namelist template !------------------------------------- ! placeholder for additional namelists !---------------------------------------------------------------------------------------- ! OUTPUT: ICON native grid, model levels, all domains: main forecast !---------------------------------------------------------------------------------------- &output_nml filetype = 2 ! output format: 2=GRIB2, 4=NETCDFv2 dom = 1 ! write domains: -1=all, 1=domain 1 only mode = 1 ! 1: forecast output_time_unit = 1 ! 1: seconds output_bounds = 0,15811200,3600 ! start, end, increment taxis_tunit = 1 ! 1: GRIB2 Units = seconds ! 3: GRIB2 Units = hours steps_per_file = 1 include_last = .FALSE. output_filename = '../../output/det/fc' ! file name base filename_format = '_DOM_' ml_varlist = 'pres_msl','asod_s','asodifd_s','athd_s','rh_2m','tot_prec','t_g','t_2m','td_2m','u_10m','v_10m','plantevap','qv','w_so','tqv', 'tqr','lhfl_s','shfl_s'!,'qhfl_s' !output_grid = .TRUE. / !---------------------------------------------------------------------------------------- ! OUTPUT: ICON interpolated to latlon grid, pressure levels, global domain: main forecast !---------------------------------------------------------------------------------------- &output_nml filetype = 2 ! output format: 2=GRIB2, 4=NETCDFv2 dom = 1 ! write domain 1 only mode = 1 ! 1: forecast output_time_unit = 1 ! 1: seconds output_bounds = 0,15811200,3600 ! start, end, increment taxis_tunit = 1 ! 1: GRIB2 Units = seconds ! 3: GRIB2 Units = hours steps_per_file = 1 include_last = .FALSE. output_filename = '../../output/det/fc_latlon' ! file name base filename_format = '_DOM__' pl_varlist = 'temp','u','v','geopot','rh' p_levels = 30000,50000,70000,85000,95000 ml_varlist = 'pres_msl','asod_s','asodifd_s','athd_s','rh_2m','tot_prec','t_g','t_2m','td_2m','u_10m','v_10m','plantevap','qv','w_so','tqv', 'tqr','lhfl_s','shfl_s' output_grid = .false. remap = 1 reg_def_mode = 1 north_pole = 198.0,39.25 reg_lon_def = -28.21, 0.03, 18.43 reg_lat_def = -23.65, 0.03, 21.67 / !---------------------------------------------------------------------------------------- ! OUTPUT: ICON native grid, model levels, all domains: ranged variables !---------------------------------------------------------------------------------------- !&output_nml ! filetype = 2 ! output format: 2=GRIB2, 4=NETCDFv2 ! dom = 1 ! write domains: -1=all, 1=domain 1 only ! mode = 1 ! 1: forecast ! output_time_unit = 1 ! 1: seconds ! output_bounds = 0,15811200,3600 ! start, end, increment ! taxis_tunit = 1 ! 1: GRIB2 Units = seconds ! ! 3: GRIB2 Units = hours ! steps_per_file = 1 ! JL To have one file per timestep ! include_last = .FALSE. ! output_filename = '../../output/det/fc_ranged' ! file name base ! filename_format = '_DOM_' ! JL add jfile ! ml_varlist = 'tot_prec','asodird_s','asodifd_s','asodifu_s' !output_grid = .TRUE. /