Run with land surface model and coupled to the RRTMG radiation model

When using the land surface model (LSM) in combination with the RRTMG radiation scheme, at least two additional NetCDF files are required in the INPUT directory with suffixes _rlw (longwave radiation) and _rsw (shortwave radiation). They have to be placed in the INPUT directory on the host where the simulation will run. The layout of the data files must be equal to the standard files provided by RRTMG (see here). A sounding file for trace gas data can be added optionally, having the suffic _rsnd.

In addition, the PALM configuration file .mrun.config has to be extended as follows

...
#----------------------------------------------------------------------------
# List of input-files
#----------------------------------------------------------------------------
...
rrtmg_lw.nc              in:locopt     d3#:d3f  $base_data/$fname/INPUT    _rlw nc
rrtmg_sw.nc              in:locopt     d3#:d3f  $base_data/$fname/INPUT    _rsw nc
RAD_SND_DATA             in:locopt     d3#:d3f  $base_data/$fname/INPUT    _rsnd nc
...

With the following parameter file example_lsm_rrtmg_p3d a simulation with LSM and RRTMG can be performed.

 &inipar nx = 63, ny = 63, nz = 96,
         dx = 50.0, dy = 50.0, dz = 25.0,
         dz_stretch_level = 1800.0,

         fft_method = 'fftw',

         initializing_actions = 'set_constant_profiles',

         ug_surface = 2.5, vg_surface = 0.0,
         pt_surface = 292.5,
         pt_vertical_gradient       =   0.0,  12.0,   0.5,
         pt_vertical_gradient_level =   0.0, 800.0, 825.0, 

         humidity = .TRUE., 
         cloud_physics = .TRUE.,
         precipitation = .TRUE.,

         q_surface = 0.005,

         surface_pressure = 1000.0,

         phi = 51.96,

         prandtl_layer = .TRUE.,
         bc_pt_b = 'dirichlet', 
         bc_q_b = 'dirichlet', 
         bc_e_b = 'neumann', /


 &d3par  end_time = 18000.0,

         create_disturbances = .T.,
         dt_disturb = 150.0, disturbance_energy_limit = 0.01,

         data_output_2d_on_each_pe = .T.,

         dt_run_control = 1.0,

         data_output = 'ql', 'qv', 'vpt', 'pt', 'shf*_xy', 'qsws*_xy', 
                       'lwp*_xy', 'prr', 'qc', 'qr',
                       't_soil', 't_soil_av', 'pt_av', 
                       'm_soil', 'm_soil_av',
                       'c_liq*_xy', 'c_liq*_xy_av',
                       'c_soil*_xy', 'c_soil*_xy_av',
                       'c_veg*_xy', 'c_veg*_xy_av',
                       'lai*_xy', 'lai*_xy_av',
                       'ghf_eb*_xy', 'ghf_eb*_xy_av', 
                       'shf_eb*_xy', 'shf_eb*_xy_av', 
                       'qsws_eb*_xy', 'qsws_eb*_xy_av', 
                       'qsws_liq_eb*_xy', 'qsws_liq_eb*_xy_av',
                       'qsws_soil_eb*_xy', 'qsws_soil_eb*_xy_av',
                       'qsws_veg_eb*_xy', 'qsws_veg_eb*_xy_av',
                       'rad_net*_xy', 'rad_net*_xy_av',
                       'rad_sw_in', 'rad_sw_out', 'rad_lw_in', 'rad_lw_out', 

        data_output_pr = '#pt', '#qv', 'wpt', '#t_soil', '#m_soil', 'u', 'v', 'rad_sw_in', 'rad_sw_out', 'rad_lw_in', 'rad_lw_out',
         
         dt_data_output = 600.0,

         dt_data_output_av = 600.0,
         averaging_interval = 10.0,
         dt_dopr = 600.0,
/

 &lsm_par veg_type = 2, soil_type = 7,

          conserve_water_content = .T.,
          dewfall = .T.,

          c_surface = 20000.0,
          vegetation_coverage = 0.9,
          z0_eb = 0.1,
          z0h_eb = 0.02,
          lambda_surface_stable = 5.0,
          lambda_surface_unstable = 5.0,
          min_canopy_resistance = 110.0,
          min_soil_resistance = 50.0,
          leaf_area_index = 2.0,

          zs = 0.07, 0.34, 1.47, 2.86,
          soil_moisture = 0.3, 0.3, 0.3, 0.3,
          root_fraction = 0.25, 0.25, 0.25, 0.25,
          soil_temperature = 290.0, 287.0, 285.0, 283.0, 283.0, 
/

 &radiation_par radiation_scheme = 'rrtmg',
                lambda = 0.0, 
                albedo_type = 5,
                constant_albedo = .F.,
                day_init = 287, 
                time_utc_init = 28800.0, 
                dt_radiation = 60.0, 
/