&run_nml
    ltestcase                   = .false.
    dtime                       =  300
    output                      = 'nml'
    msg_level                   = 15
    num_lev                     = 90
    lvert_nest                  = .false.
    ldynamics                   = .true.
    ltransport                  = .true.
    ntracer                     = 5
    iforcing                    = 3
/



&nwp_phy_nml
    lupatmo_phy = .FALSE.
    inwp_gscp=1,
    inwp_convection=1,
    inwp_radiation=1,
    inwp_cldcover=1,
    inwp_turb=1,
    inwp_satad=1,
    inwp_sso=1,
    inwp_gwd=1,
    inwp_surface=1,
    icapdcycl=3,
    latm_above_top=.false.
    efdt_min_raylfric=3600.,
    itype_z0=2,
    icpl_aero_conv=1,
    icpl_aero_gscp=1,
    lrtm_filename='rrtmg_lw.nc',
    !cldopt_filename='rrtm_cldopt.nc', ! simply other file name, probably same file content
    dt_rad=2160,
    dt_conv=720,
    dt_sso=1440,
    dt_gwd=1440,
    ldetrain_conv_prec=.false.
    mu_rain=0.5,
    rain_n0_factor=0.1,
/

&radiation_nml
    irad_o3=79,
    irad_aero=6,
    albedo_type=2, ! this requires t_2m climatology in external parameter file
    vmr_co2=390.e-06,
    vmr_ch4=1800.e-09,
    vmr_n2o=322.0e-09,
    vmr_o2=0.20946,
    vmr_cfc11=240.e-12,
    vmr_cfc12=532.e-12,
/

&lnd_nml
    !ntiles=3,! needs more layers of t_g in initial conditions
    nlev_snow=3,
    lmulti_snow=.false.,
    itype_heatcond=3,
    idiag_snowfrac=20,
    lsnowtile=.true.,
    lseaice=.true.,
    llake=.true.,
    lprog_albsi=.true.,
    itype_lndtbl=4,
    itype_root=2,
    itype_evsl=4,
    itype_trvg=3,
    cwimax_ml=5.e-4,
    c_soil=1.25,
    c_soil_urb=0.5,
    sstice_mode=2, !use sst climatology in external parameter file to update ssts
    itype_snowevap=3,
/



&grid_nml
    dynamics_parent_grid_id     = 0
    lredgrid_phys               = .true.
/

&extpar_nml
    itopo                       = 1
/

&initicon_nml
    init_mode                   = 1,
    lconsistency_checks         = .false.
    ana_varnames_map_file       = 'ana_varnames_map_file.txt'
/

! settings from operational setup for vertical coordinate
&sleve_nml
    min_lay_thckn               = 20.
    max_lay_thckn               = 400.
    htop_thcknlimit             = 14000.
    top_height                  = 75000.
    stretch_fac                 = 0.9
    decay_scale_1               = 4000.
    decay_scale_2               = 2500.
    decay_exp                   = 1.2
    flat_height                 = 16000.
/

&io_nml
    lnetcdf_flt64_output        = .true.
/


&nonhydrostatic_nml
    iadv_rhotheta=2,
    ivctype=2,
    itime_scheme=4,
    exner_expol=0.333,
    vwind_offctr=0.2,
    damp_height=44000.,
    rayleigh_coeff=0.5,
    lhdiff_rcf=.true.,
    divdamp_order=24,
    divdamp_type=32,
    divdamp_fac=0.004,
    divdamp_trans_start=12500,
    divdamp_trans_end=17500,
    l_open_ubc=.false.,
    igradp_method=3,
    l_zdiffu_t=.true.,
    thslp_zdiffu=0.02,
    thhgtd_zdiffu=125.,
    htop_moist_proc=22500.,
    hbot_qvsubstep=16000.,
/


&parallel_nml
    nproma                      = 16
/


    ! output for analysis of simulation (surface fields)
    &output_nml
        output_interval             = 'PT21600S'
        steps_per_file              = 1
        filename_format             = "surface_fields_<jfile>"
        output_start                = '#OUTPUT_START#'
        output_end                  = '#OUTPUT_END#'
        include_last                = .true.
        remap                       = 1         ! remap to regular lat/lon grid
        reg_lat_def                 = -20.,1.0,90.         ! start, increment, end latitude in degrees
        reg_lon_def                 = -180.,1.0,179.0        ! start, increment, end longitude in degrees
        ml_varlist                  = 't_2m','tmin_2m','tmax_2m','pres_sfc','pres_msl','u_10m','v_10m','tot_prec','tqc','tqi','tqv','snowfrac'
    /
    

    ! output for analysis of simulation (on pressure levels)
    &output_nml
        output_interval             = 'PT21600S'
        steps_per_file              = 1
        filename_format             = "fields_3d_<jfile>"
        output_start                = '#OUTPUT_START#'
        output_end                  = '#OUTPUT_END#'
        include_last                = .true.
        remap                       = 1         ! remap to regular lat/lon grid
        reg_lat_def                 = -20.,1.0,90.         ! start, increment, end latitude in degrees
        reg_lon_def                 = -180.,1.0,179.0        ! start, increment, end longitude in degrees
        pl_varlist                  = 'u','v','omega','geopot','temp'
        p_levels         =    100000,85000,75000,65000,60000, 55000,50000,47500,45000,42500,40000,37500, 35000,32500,30000,28000,26500, 25000,23500,22000, 21000,20000,19000,18000, 17000,16000, 15000,14000,13000,12000, 11000,10000,8500,7000,5500,4000,3000, 2500,2000,1500,1000, 700,500,300,200,150,100,70,50,30,20,10
    /