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. For this example run, simply copy the two files from here to your INPUT directory and rename them as described above.

Before you proceed, make sure you have installed the RRTMG library as external library and you have re-compiled PALM afterwards. An installation manual can be found here.

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, 
/

This setup requires about 1.5h of time on a quad-core processor with 4x2.66Ghz.