Changeset 4675 for palm

Timestamp:

Sep 11, 2020 10:00:26 AM (4 years ago)

Author:

eckhard

Message:

Support for homogeneous (domain-averaged) boundary conditions and soil profile initialization

Location:

palm/trunk/UTIL/inifor/src

Files:

• inifor.f90 (modified) (8 diffs)
• inifor_control.f90 (modified) (4 diffs)
• inifor_defs.f90 (modified) (3 diffs)
• inifor_grid.f90 (modified) (127 diffs)
• inifor_io.f90 (modified) (11 diffs)
• inifor_transform.f90 (modified) (3 diffs)
• inifor_types.f90 (modified) (2 diffs)

palm/trunk/UTIL/inifor/src/inifor.f90

@@ -26 +26 @@
 ! -----------------
 ! $Id$
+! Support for soil profile initialization
+! Improved code formatting
+!
+!
+! 4659 2020-08-31 11:21:17Z eckhard
 ! List all warnings after successful run
 ! Only define netCDF variables in enabled IO groups
@@ -202 +207 @@
 !-- Set up the tables containing the input and output variables and set
 !-- the corresponding netCDF dimensions for each output variable
-    CALL setup_variable_tables( cfg%ic_mode )
+    CALL setup_variable_tables
     CALL setup_io_groups
     CALL log_runtime( 'time', 'init' )
@@ -273 +278 @@
                          CASE( 'init soil' )
    
-                            ALLOCATE( output_arr(0:output_var%grid%nx,            &
-                                                 0:output_var%grid%ny,            &
+                            ALLOCATE( output_arr(0:output_var%grid%nx,         &
+                                                 0:output_var%grid%ny,         &
                                                  SIZE( output_var%grid%depths )) )
    
                          CASE ( 'surface forcing' )
    
-                            ALLOCATE( output_arr(0:output_var%grid%nx,            &
+                            ALLOCATE( output_arr(0:output_var%grid%nx,         &
                                                  0:output_var%grid%ny, 1) )
    
@@ -299 +304 @@
                       CASE ( 'interpolate_3d' )
    
-                         ALLOCATE( output_arr(0:output_var%grid % nx,           &
-                                              0:output_var%grid % ny,           &
+                         ALLOCATE( output_arr(0:output_var%grid % nx,          &
+                                              0:output_var%grid % ny,          &
                                               1:output_var%grid % nz) )
    
                          CALL log_runtime( 'time', 'alloc' )
-                         CALL interpolate_3d(                                     &
-                            input_buffer(output_var%input_id)%array(:,:,:),       &
-                            output_arr(:,:,:),                                    &
-                            output_var%intermediate_grid,                         &
+                         CALL interpolate_3d(                                  &
+                            input_buffer(output_var%input_id)%array(:,:,:),    &
+                            output_arr(:,:,:),                                 &
+                            output_var%intermediate_grid,                      &
                             output_var%grid)
                          CALL log_runtime( 'time', 'comp' )
@@ -315 +320 @@
                       CASE ( 'average profile' )
    
-                         ALLOCATE( output_arr(0:output_var%grid%nx,               &
-                                              0:output_var%grid%ny,               &
-                                              1:output_var%grid%nz) )
+                         ALLOCATE( output_arr(0:0, 0:0, 1:output_var%averaging_grid%nz) )
                          CALL log_runtime( 'time', 'alloc' )
-                         
-                         CALL interp_average_profile(                             &
-                            input_buffer(output_var%input_id)%array(:,:,:),     &
-                            output_arr(0,0,:),                                    &
+                         CALL interp_average_profile(                          &
+                            input_buffer(output_var%input_id)%array(:,:,:),    &
+                            output_arr(0,0,:),                                 &
                             output_var%averaging_grid )
    
-                         IF ( TRIM( output_var%name ) ==                          &
+                         IF ( TRIM( output_var%name ) ==                       &
                               'surface_forcing_surface_pressure' )  THEN
    
@@ -331 +333 @@
                                output_arr(0,0,1) = p0
                             ELSE
-                               CALL get_surface_pressure(                         &
-                                  output_arr(0,0,:), rho_centre,                  &
+                               CALL get_surface_pressure(                      &
+                                  output_arr(0,0,:), rho_centre,               &
                                   output_var%averaging_grid )
                             ENDIF
    
                          ENDIF
+                         CALL log_runtime( 'time', 'comp' )
+
+                      CASE ( 'average levels' )
+
+                         ALLOCATE( output_arr(0:0, 0:0, 1:output_var%averaging_grid%nz) )
+                         CALL log_runtime( 'time', 'alloc' )
+                         CALL average_profile(                                 &
+                            input_buffer(output_var%input_id)%array(:,:,:),    &
+                            output_arr(0,0,:),                                 &
+                            output_var%averaging_grid                          &
+                         )
                          CALL log_runtime( 'time', 'comp' )
    
@@ -393 +406 @@
                          SELECT CASE( TRIM( output_var%name ) )
    
-                         CASE( 'internal_pressure_north',                         &
-                               'internal_pressure_south',                         &
-                               'internal_pressure_east',                          &
+                         CASE( 'internal_pressure_north',                      &
+                               'internal_pressure_south',                      &
+                               'internal_pressure_east',                       &
                                'internal_pressure_west' )
    
-                            CALL average_pressure_perturbation(                   &
-                               input_buffer(output_var%input_id) % array(:,:,:),&
-                               internal_arr(:),                                   &
-                               cosmo_grid, output_var%averaging_grid            &
+                            CALL average_pressure_perturbation(                &
+                               input_buffer(output_var%input_id)%array(:,:,:), &
+                               internal_arr(:),                                &
+                               cosmo_grid, output_var%averaging_grid           &
                             )
    
                          CASE DEFAULT
    
-                            CALL average_profile(                                 &
-                               input_buffer(output_var%input_id) % array(:,:,:),&
-                               internal_arr(:),                                   &
-                               output_var%averaging_grid                        &
+                            CALL average_profile(                              &
+                               input_buffer(output_var%input_id)%array(:,:,:), &
+                               internal_arr(:),                                &
+                               output_var%averaging_grid                       &
                             )
 
                       END SELECT
 
-
-!
-!--                   Output of geostrophic pressure profiles (with --debug
-!--                   option) is currently deactivated, since they are now
-!--                   defined on averaged COSMO levels instead of PALM levels
-!--                   (requires definiton of COSMO levels in netCDF output.)
-                      !IF (.TRUE.)  THEN
-                      !   ALLOCATE( output_arr(1,1,1:output_var%grid % nz) )
-                      !   output_arr(1,1,:) = internal_arr(:)
-                      !ENDIF
                       CALL log_runtime( 'time', 'comp' )
 
@@ -450 +453 @@
                                                     ug_cosmo, vg_cosmo )
 
-                            CALL interpolate_1d( ug_cosmo, ug_palm,             &
+                            CALL interpolate_1d( ug_cosmo, ug_palm,            &
                                                  output_var%grid )
 
-                            CALL interpolate_1d( vg_cosmo, vg_palm,             &
+                            CALL interpolate_1d( vg_cosmo, vg_palm,            &
                                                  output_var%grid )
                          ENDIF

palm/trunk/UTIL/inifor/src/inifor_control.f90

@@ -26 +26 @@
 ! -----------------
 ! $Id$
+! Improve visibility of errors by printing them below warnings listing
+!
+!
+! 4659 2020-08-31 11:21:17Z eckhard
 ! List warnings after successful run or abort
 ! Produce failure exit code (1) on program abort for test automation
@@ -207 +211 @@
  SUBROUTINE report_warnings()
 
-    INTEGER(iwp) ::  warning_idx
+    INTEGER(iwp)        ::  warning_idx
+    CHARACTER (LEN=500) ::  warning = ''
 
     IF (n_wrngs > 0)  THEN
-       message = 'Encountered the following '// TRIM( str( n_wrngs ) ) // " warning(s) during this run:"
-       CALL report( 'report_warnings', message )
+       warning = 'Encountered the following '// TRIM( str( n_wrngs ) ) // " warning(s) during this run:"
+       CALL report( 'report_warnings', warning)
 
        DO warning_idx = 1, n_wrngs
@@ -233 +238 @@
     message = TRIM( message ) // "' successfully."
     IF (n_wrngs > 0)  THEN
-       message = TRIM( message ) // " Some warnings were encountered."
+       message = TRIM( message ) // " Some warnings were encountered, see above."
     ENDIF
     CALL report( 'main loop', message )
@@ -268 +273 @@
     CHARACTER(LEN=*), INTENT(IN) ::  message !< log message
 
+    CALL report_warnings
     CALL report( routine, "ERROR: " // TRIM( message ) // " Stopping." )
-    CALL report_warnings
     CALL close_log
     CALL EXIT(1)

palm/trunk/UTIL/inifor/src/inifor_defs.f90

@@ -21 +21 @@
 ! Current revisions:
 ! -----------------
-! 
-! 
+!
+!
 ! Former revisions:
 ! -----------------
 ! $Id$
+! New command-line options for soil profile initialization and homogeneous
+!    boundary conditions
+! Improved code formatting
+! 
+! 
+! 4659 2020-08-31 11:21:17Z eckhard
 ! Bumped version number
 ! 
@@ -177 +183 @@
 !
 !-- PALM static driver attribute names (PIDS 1.9)
-CHARACTER(SNAME), PARAMETER ::  PIDS_ORIGIN_LON = 'origin_lon'
-CHARACTER(SNAME), PARAMETER ::  PIDS_ORIGIN_LAT = 'origin_lat'
-CHARACTER(SNAME), PARAMETER ::  PIDS_ORIGIN_Z   = 'origin_z'
+ CHARACTER(SNAME), PARAMETER ::  PIDS_ORIGIN_LON = 'origin_lon'
+ CHARACTER(SNAME), PARAMETER ::  PIDS_ORIGIN_LAT = 'origin_lat'
+ CHARACTER(SNAME), PARAMETER ::  PIDS_ORIGIN_Z   = 'origin_z'
 
 ! 
 !-- COSMO netCDF parameters
-INTEGER, PARAMETER          ::  NC_DEPTH_DIM_IDX = 3
-CHARACTER(SNAME), PARAMETER ::  NC_DEPTH_NAME = 'depth_2'
-CHARACTER(SNAME), PARAMETER ::  NC_HHL_NAME = 'HHL'
-CHARACTER(SNAME), PARAMETER ::  NC_RLAT_NAME = 'rlat'
-CHARACTER(SNAME), PARAMETER ::  NC_RLON_NAME = 'rlon'
-CHARACTER(SNAME), PARAMETER ::  NC_ROTATED_POLE_NAME = 'rotated_pole'
-CHARACTER(SNAME), PARAMETER ::  NC_POLE_LATITUDE_NAME = 'grid_north_pole_latitude'
-CHARACTER(SNAME), PARAMETER ::  NC_POLE_LONGITUDE_NAME = 'grid_north_pole_longitude'
+ INTEGER, PARAMETER          ::  NC_DEPTH_DIM_IDX = 3
+ CHARACTER(SNAME), PARAMETER ::  NC_DEPTH_NAME = 'depth_2'
+ CHARACTER(SNAME), PARAMETER ::  NC_HHL_NAME = 'HHL'
+ CHARACTER(SNAME), PARAMETER ::  NC_RLAT_NAME = 'rlat'
+ CHARACTER(SNAME), PARAMETER ::  NC_RLON_NAME = 'rlon'
+ CHARACTER(SNAME), PARAMETER ::  NC_ROTATED_POLE_NAME = 'rotated_pole'
+ CHARACTER(SNAME), PARAMETER ::  NC_POLE_LATITUDE_NAME = 'grid_north_pole_latitude'
+ CHARACTER(SNAME), PARAMETER ::  NC_POLE_LONGITUDE_NAME = 'grid_north_pole_longitude'
 
 !
 !-- INIFOR parameters
+ CHARACTER(LEN=*), PARAMETER ::  CFG_INIT_PROFILE = 'profile'
+ CHARACTER(LEN=*), PARAMETER ::  CFG_INIT_VOLUME = 'volume'
+ CHARACTER(LEN=*), PARAMETER ::  CFG_INIT_SOIL_PROFILE = 'profile'
+ CHARACTER(LEN=*), PARAMETER ::  CFG_INIT_SOIL_VOLUME = 'volume'
+ CHARACTER(LEN=*), PARAMETER ::  CFG_FORCING_HETERO = 'hetero'
+ CHARACTER(LEN=*), PARAMETER ::  CFG_FORCING_HOMO = 'homo'
+ CHARACTER(LEN=*), PARAMETER ::  CFG_FORCING_NUDGING = 'nudging'
  INTEGER(iwp), PARAMETER     ::  FILL_ITERATIONS = 5          !< Number of iterations for extrapolating soil data into COSMO-DE
                                                               !< water cells [-]
@@ -201 +214 @@
     ACHAR( 10 ) // ' Copyright 2017-2020 Deutscher Wetterdienst Offenbach' !< Copyright notice
  CHARACTER(LEN=*), PARAMETER ::  LOG_FILE_NAME = 'inifor.log' !< Name of INIFOR's log file
- CHARACTER(LEN=*), PARAMETER ::  VERSION = '2.0.0'            !< INIFOR version number
+ CHARACTER(LEN=*), PARAMETER ::  VERSION = '2.1.0'            !< INIFOR version number
  
  END MODULE inifor_defs

palm/trunk/UTIL/inifor/src/inifor_grid.f90

@@ -26 +26 @@
 ! -----------------
 ! $Id$
+! Support for profile initialization of soil and for homogeneous (profile)
+!    boundary conditions, including respective command-line options
+! Improved code formatting
+!
+!
+! 4659 2020-08-31 11:21:17Z eckhard
 ! Only define netCDF variables in enabled IO groups
 ! Added cloud and precipitation quantiteis (cloud water, cloud ice, rain, snow
@@ -181 +187 @@
     USE inifor_control
     USE inifor_defs,                                                           &
-        ONLY:  DATE, EARTH_RADIUS, TO_RADIANS, TO_DEGREES, PI,                 &
+        ONLY:  CFG_INIT_PROFILE, CFG_INIT_SOIL_VOLUME, CFG_INIT_SOIL_PROFILE,  &
+               CFG_FORCING_HETERO, CFG_FORCING_HOMO, CFG_FORCING_NUDGING,      &
+               DATE, EARTH_RADIUS, TO_RADIANS, TO_DEGREES, PI,                 &
                SNAME, LNAME, PATH, FORCING_STEP, FILL_ITERATIONS,              &
                BETA, P_SL, T_SL, BETA, RD, RV, G, P_REF, RD_PALM, CP_PALM,     &
@@ -297 +305 @@
     INTEGER(iwp) ::  ndepths !< number of COSMO-DE soil layers
     INTEGER(iwp) ::  start_hour_flow          !< start of flow forcing in number of hours relative to start_date
-    INTEGER(iwp) ::  start_hour_soil          !< start of soil forcing in number of hours relative to start_date, typically equals start_hour_flow
     INTEGER(iwp) ::  start_hour_radiation     !< start of radiation forcing in number of hours relative to start_date, 0 to 2 hours before start_hour_flow to reconstruct hourly averages from one- to three hourly averages of the input data
     INTEGER(iwp) ::  start_hour_precipitation !< start of forcing for precipitaiton forcing in number of hours relative to start_date
@@ -362 +369 @@
     TYPE(grid_definition), TARGET ::  w_south_intermediate            !< intermediate grid for southern w boundary condition
     TYPE(grid_definition), TARGET ::  w_top_intermediate              !< intermediate grid for top w boundary condition
-    TYPE(grid_definition), TARGET ::  north_averaged_scalar_profile   !< grid of the northern geostrophic scalar averaging region
-    TYPE(grid_definition), TARGET ::  south_averaged_scalar_profile   !< grid of the southern geostrophic scalar averaging region
-    TYPE(grid_definition), TARGET ::  west_averaged_scalar_profile    !< grid of the western geostrophic scalar averaging region
-    TYPE(grid_definition), TARGET ::  east_averaged_scalar_profile    !< grid of the eastern geostrophic scalar averaging region
-    TYPE(grid_definition), TARGET ::  averaged_scalar_profile         !< grid of the central geostrophic scalar averaging region
-    TYPE(grid_definition), TARGET ::  averaged_w_profile              !< grid of the central geostrophic w-velocity averaging region
-    TYPE(grid_definition), TARGET ::  averaged_initial_scalar_profile !< averaging grid for initial scalar profiles
-    TYPE(grid_definition), TARGET ::  averaged_initial_w_profile      !< averaging grid for the initial w profile
+    TYPE(grid_definition), TARGET ::  north_geostrophic_scalar_profile!< grid of the northern geostrophic scalar averaging region
+    TYPE(grid_definition), TARGET ::  south_geostrophic_scalar_profile!< grid of the southern geostrophic scalar averaging region
+    TYPE(grid_definition), TARGET ::  west_geostrophic_scalar_profile !< grid of the western geostrophic scalar averaging region
+    TYPE(grid_definition), TARGET ::  east_geostrophic_scalar_profile !< grid of the eastern geostrophic scalar averaging region
+    TYPE(grid_definition), TARGET ::  geostrophic_scalar_profile      !< grid of the central geostrophic scalar averaging region
+    TYPE(grid_definition), TARGET ::  geostrophic_w_profile           !< grid of the central geostrophic w-velocity averaging region
+    TYPE(grid_definition), TARGET ::  averaged_soil_profile           !< averaging grid for initial soil profiles
+    TYPE(grid_definition), TARGET ::  averaged_scalar_profile         !< averaging grid for initial and boundary condition scalar profiles
+    TYPE(grid_definition), TARGET ::  averaged_w_profile              !< averaging grid for initial and boundary condition scalar profiles
+    TYPE(grid_definition), TARGET ::  averaged_scalar_top_point       !< averaging grid for top scalar boundary conditions for homogeneous forcing mode
+    TYPE(grid_definition), TARGET ::  averaged_w_top_point            !< averaging grid for top w boundary condition for homogeneous forcing mode
 
     TYPE(io_group), ALLOCATABLE, TARGET ::  io_group_list(:)  !< List of I/O groups, which group together output variables that share the same input variable
@@ -432 +442 @@
 !-- Parameters for file names
     start_hour_flow = 0
-    start_hour_soil = 0
     start_hour_radiation = 0
     start_hour_precipitation = start_hour_flow 
@@ -464 +473 @@
     cfg%static_driver_file = ''
     cfg%output_file = './palm-4u-input.nc'
-    cfg%ic_mode = 'profile'
-    cfg%bc_mode = 'real'
+    cfg%ic_mode = CFG_INIT_PROFILE
+    cfg%isc_mode = CFG_INIT_SOIL_VOLUME
+    cfg%bc_mode = CFG_FORCING_HETERO
     cfg%averaging_mode = 'level'
 
@@ -485 +495 @@
 
     init_variables_required = .TRUE.
-    boundary_variables_required = TRIM( cfg%bc_mode ) == 'real'
-    ls_forcing_variables_required = TRIM( cfg%bc_mode ) == 'ideal'
+    boundary_variables_required = (                                            &
+       ( TRIM( cfg%bc_mode ) == CFG_FORCING_HETERO )  .OR.                     &
+       ( TRIM( cfg%bc_mode ) == CFG_FORCING_HOMO )                             &
+    )
+    ls_forcing_variables_required = TRIM( cfg%bc_mode ) == CFG_FORCING_NUDGING
     surface_forcing_required = .TRUE.
 
@@ -503 +516 @@
     CALL validate_config( cfg ) 
 
-    CALL report('setup_parameters', "initialization mode: " // TRIM(cfg%ic_mode))
-    CALL report('setup_parameters', "       forcing mode: " // TRIM(cfg%bc_mode))
-    CALL report('setup_parameters', "     averaging mode: " // TRIM(cfg%averaging_mode))
-    CALL report('setup_parameters', "    averaging angle: " // real_to_str(cfg%averaging_angle))
-    CALL report('setup_parameters', "    averaging angle: " // real_to_str(cfg%averaging_angle))
-    CALL report('setup_parameters', "          data path: " // TRIM(cfg%input_path))
-    CALL report('setup_parameters', "           hhl file: " // TRIM(cfg%hhl_file))
-    CALL report('setup_parameters', "       soiltyp file: " // TRIM(cfg%soiltyp_file))
-    CALL report('setup_parameters', "      namelist file: " // TRIM(cfg%namelist_file))
-    CALL report('setup_parameters', "   output data file: " // TRIM(output_file%name))
+    CALL report('setup_parameters', "atmosphere initialization mode: " // TRIM(cfg%ic_mode))
+    CALL report('setup_parameters', "      soil initialization mode: " // TRIM(cfg%isc_mode))
+    CALL report('setup_parameters', "                  forcing mode: " // TRIM(cfg%bc_mode))
+    CALL report('setup_parameters', "                averaging mode: " // TRIM(cfg%averaging_mode))
+    CALL report('setup_parameters', "               averaging angle: " // real_to_str(cfg%averaging_angle))
+    CALL report('setup_parameters', "               averaging angle: " // real_to_str(cfg%averaging_angle))
+    CALL report('setup_parameters', "                     data path: " // TRIM(cfg%input_path))
+    CALL report('setup_parameters', "                      hhl file: " // TRIM(cfg%hhl_file))
+    CALL report('setup_parameters', "                  soiltyp file: " // TRIM(cfg%soiltyp_file))
+    CALL report('setup_parameters', "                 namelist file: " // TRIM(cfg%namelist_file))
+    CALL report('setup_parameters', "              output data file: " // TRIM(output_file%name))
     IF (cfg%process_precipitation )  THEN
         CALL report('setup_parameters', "      precipitation: enabled")
@@ -557 +571 @@
        cfg%input_path, flow_prefix, flow_suffix, flow_files)
     CALL get_input_file_list(                                                  &
-       cfg%start_date, start_hour_soil, end_hour, step_hour,                   &
+       cfg%start_date, start_hour_flow, start_hour_flow, step_hour,            &
        cfg%input_path, soil_prefix, soil_suffix, soil_files)
     CALL get_input_file_list(                                                  &
@@ -1097 +1111 @@
     ENDIF
 
-
-    CALL init_averaging_grid(averaged_initial_scalar_profile, cosmo_grid,   &
+    CALL init_averaging_grid(averaged_soil_profile, cosmo_grid,             &
+            x = 0.5_wp * lx, y = 0.5_wp * ly, z = depths, z0 = z0,          &
+            lonmin = lonmin_palm, lonmax = lonmax_palm,                     &
+            latmin = latmin_palm, latmax = latmax_palm,                     &
+            kind='scalar', name='averaged soil profile')
+
+    CALL init_averaging_grid(averaged_scalar_profile, cosmo_grid,           &
             x = 0.5_wp * lx, y = 0.5_wp * ly, z = z, z0 = z0,               &
             lonmin = lonmin_palm, lonmax = lonmax_palm,                     &
             latmin = latmin_palm, latmax = latmax_palm,                     &
-            kind='scalar', name='averaged initial scalar')
-
-    CALL init_averaging_grid(averaged_initial_w_profile, cosmo_grid,        &
+            kind='scalar', name='averaged scalar profile')
+
+    CALL init_averaging_grid(averaged_w_profile, cosmo_grid,                &
             x = 0.5_wp * lx, y = 0.5_wp * ly, z = zw, z0 = z0,              &
             lonmin = lonmin_palm, lonmax = lonmax_palm,                     &
             latmin = latmin_palm, latmax = latmax_palm,                     &
-            kind='w', name='averaged initial w')
-
-    CALL init_averaging_grid(averaged_scalar_profile, cosmo_grid,           &
+            kind='w', name='averaged w profile')
+
+    CALL init_averaging_grid(averaged_scalar_top_point, cosmo_grid,         &
+            x = 0.5_wp * lx, y = 0.5_wp * ly, z = z_top, z0 = z0,           &
+            lonmin = lonmin_palm, lonmax = lonmax_palm,                     &
+            latmin = latmin_palm, latmax = latmax_palm,                     &
+            kind='scalar', name='averaged scalar top point')
+
+    CALL init_averaging_grid(averaged_w_top_point, cosmo_grid,              &
+            x = 0.5_wp * lx, y = 0.5_wp * ly, z = zw_top, z0 = z0,          &
+            lonmin = lonmin_palm, lonmax = lonmax_palm,                     &
+            latmin = latmin_palm, latmax = latmax_palm,                     &
+            kind='w', name='averaged w top point')
+
+    CALL init_averaging_grid(geostrophic_scalar_profile, cosmo_grid,           &
             x = 0.5_wp * lx, y = 0.5_wp * ly, z = z, z0 = z0,               &
             lonmin = lam_west, lonmax = lam_east,                           &
             latmin = phi_south, latmax = phi_north,                         &
-            kind='scalar', name='centre geostrophic scalar')
-
-    CALL init_averaging_grid(averaged_w_profile, cosmo_grid,                &
+            kind='scalar', name='centre geostrophic scalar profile')
+
+    CALL init_averaging_grid(geostrophic_w_profile, cosmo_grid,                &
             x = 0.5_wp * lx, y = 0.5_wp * ly, z = zw, z0 = z0,              &
             lonmin = lam_west, lonmax = lam_east,                           &
             latmin = phi_south, latmax = phi_north,                         &
-            kind='w', name='centre geostrophic w')
-
-    CALL init_averaging_grid(south_averaged_scalar_profile, cosmo_grid,     &
+            kind='w', name='centre geostrophic w profile')
+
+    CALL init_averaging_grid(south_geostrophic_scalar_profile, cosmo_grid,     &
             x = 0.5_wp * lx, y = 0.5_wp * ly, z = z, z0 = z0,               &
             lonmin = lam_west, lonmax = lam_east,                           &
             latmin = phi_centre - averaging_angle,                          &
             latmax = phi_centre,                                            &
-            kind='scalar', name='south geostrophic scalar')
-
-    CALL init_averaging_grid(north_averaged_scalar_profile, cosmo_grid,     &
+            kind='scalar', name='south geostrophic scalar profile')
+
+    CALL init_averaging_grid(north_geostrophic_scalar_profile, cosmo_grid,     &
             x = 0.5_wp * lx, y = 0.5_wp * ly, z = z, z0 = z0,               &
             lonmin = lam_west, lonmax = lam_east,                           &
             latmin = phi_centre,                                            &
             latmax = phi_centre + averaging_angle,                          &
-            kind='scalar', name='north geostrophic scalar')
-
-    CALL init_averaging_grid(west_averaged_scalar_profile, cosmo_grid,      &
+            kind='scalar', name='north geostrophic scalar profile')
+
+    CALL init_averaging_grid(west_geostrophic_scalar_profile, cosmo_grid,      &
             x = 0.5_wp * lx, y = 0.5_wp * ly, z = z, z0 = z0,               &
             lonmin = lam_centre - averaging_angle,                          &
             lonmax = lam_centre,                                            &
             latmin = phi_south, latmax = phi_north,                         &
-            kind='scalar', name='west geostrophic scalar')
-
-    CALL init_averaging_grid(east_averaged_scalar_profile, cosmo_grid,      &
+            kind='scalar', name='west geostrophic scalar profile')
+
+    CALL init_averaging_grid(east_geostrophic_scalar_profile, cosmo_grid,      &
             x = 0.5_wp * lx, y = 0.5_wp * ly, z = z, z0 = z0,               &
             lonmin = lam_centre,                                            &
             lonmax = lam_centre + averaging_angle,                          &
             latmin = phi_south, latmax = phi_north,                         &
-            kind='scalar', name='east geostrophic scalar')
+            kind='scalar', name='east geostrophic scalar profile')
 
 
@@ -1195 +1226 @@
     ENDIF
 
-    IF (TRIM(cfg%ic_mode) == 'profile')  THEN
+    IF (TRIM(cfg%ic_mode) == CFG_INIT_PROFILE)  THEN
         !TODO: remove this conditional if not needed. 
     ENDIF
@@ -1285 +1316 @@
     setup_volumetric = .TRUE.
     IF (PRESENT(ic_mode))  THEN
-       IF (TRIM(ic_mode) == 'profile')  setup_volumetric = .FALSE.
+       IF (TRIM(ic_mode) == CFG_INIT_PROFILE)  setup_volumetric = .FALSE.
     ENDIF
 
@@ -1374 +1405 @@
 !      
 !--       Allocate neighbour indices and weights
-          IF (TRIM(ic_mode) .NE. 'profile')  THEN
+          IF (TRIM(ic_mode) .NE. CFG_INIT_PROFILE)  THEN
              ALLOCATE( grid%kk(0:nx, 0:ny, 1:nz, 2) )
              grid%kk(:,:,:,:) = -1
@@ -1446 +1477 @@
 !      
 !--       Allocate neighbour indices and weights
-          IF (TRIM(ic_mode) .NE. 'profile')  THEN
+          IF (TRIM(ic_mode) .NE. CFG_INIT_PROFILE)  THEN
              ALLOCATE( grid%kk(0:nx, 0:ny, 1:nz, 2) )
              grid%kk(:,:,:,:) = -1
@@ -1858 +1889 @@
     IF ( ANY( min_dz_stretch_level_end(1:number_stretch_level_start) >      &
               dz_stretch_level_end(1:number_stretch_level_start) ) )  THEN
-    !IF ( ANY( min_dz_stretch_level_end >      &
-    !          dz_stretch_level_end ) ) THEN
           message = 'Each dz_stretch_level_end has to be larger '  // &
                     'than its corresponding value for ' //            &
@@ -2369 +2398 @@
                       invar%to_be_processed = .FALSE.
                       outvar%to_be_processed = .FALSE.
-                      message = "Optional output variable '" // TRIM( outvar%name ) // "' was deactivated" // &
-                         " because the corresponding input variable was not found in file '" // TRIM( filename ) // "'."
+                      message = "Skipping optional output variable '" // TRIM( outvar%name ) // &
+                         ", corresponding input variable not available in file '" // TRIM( filename ) // "'."
                       CALL warn( 'validate_io_groups', message )
                 ENDIF
              ENDIF
 
-             !IF ( .NOT.netcdf_variable_present_in_file( invar%name, filename ) )  THEN
-!--          !      abort due to missing input variable
-             !      message = "The mandatory input variable '" // TRIM( invar%name ) // &
-             !         " was not found in the input file '" // TRIM( filename ) // "'."
-             !      CALL inifor_abort( 'validate_io_groups', message  )
-             !   ENDIF
-             !IF ( .NOT.netcdf_variable_present_in_file( invar%name, filename ) )  THEN
-!--          !      abort due to missing input variable
-             !      message = "The mandatory input variable '" // TRIM( invar%name ) // &
-             !         " was not found in the input file '" // TRIM( filename ) // "'."
-             !      CALL inifor_abort( 'validate_io_groups', message  )
-             !   ENDIF
-             !ENDIF
 !--       outvar loop
           ENDDO
@@ -2493 +2509 @@
 !> Initializes the variable list.
 !------------------------------------------------------------------------------!
- SUBROUTINE setup_variable_tables(ic_mode)
-    CHARACTER(LEN=*), INTENT(IN) ::  ic_mode
-    INTEGER(iwp)                 ::  n_invar = 0  !< number of variables in the input variable table
-    INTEGER(iwp)                 ::  n_outvar = 0 !< number of variables in the output variable table
-    TYPE(nc_var), POINTER        ::  var
+ SUBROUTINE setup_variable_tables
+    INTEGER(iwp)                    ::  n_invar = 0  !< number of variables in the input variable table
+    INTEGER(iwp)                    ::  n_outvar = 0 !< number of variables in the output variable table
+    TYPE(nc_var), POINTER           ::  var
 
     IF (TRIM(cfg%start_date) == '')  THEN
@@ -2626 +2641 @@
        output_file       = output_file,                                     &
        grid              = palm_grid,                                       &
-       intermediate_grid = palm_intermediate                                &
+       intermediate_grid = palm_intermediate,                               &
+       is_profile = (TRIM(cfg%isc_mode) == CFG_INIT_SOIL_PROFILE)           &
     )
 
@@ -2638 +2654 @@
        output_file       = output_file,                                     &
        grid              = palm_grid,                                       &
-       intermediate_grid = palm_intermediate                                &
+       intermediate_grid = palm_intermediate,                               &
+       is_profile = (TRIM(cfg%isc_mode) == CFG_INIT_SOIL_PROFILE)           &
     )
 
@@ -2651 +2668 @@
        grid              = palm_grid,                                       &
        intermediate_grid = palm_intermediate,                               &
-       is_profile = (TRIM(ic_mode) == 'profile')                            &
-    )
-    IF (TRIM(ic_mode) == 'profile')  THEN
-       output_var_table(3)%averaging_grid => averaged_initial_scalar_profile
-    ENDIF
+       is_profile = (TRIM(cfg%ic_mode) == CFG_INIT_PROFILE)                 &
+    )
 
     output_var_table(4) = init_nc_var(                                      &
@@ -2666 +2680 @@
        grid              = scalars_west_grid,                               &
        intermediate_grid = scalars_west_intermediate,                       &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO),                &
        output_file = output_file                                            &
     )
@@ -2678 +2693 @@
        grid              = scalars_east_grid,                               &
        intermediate_grid = scalars_east_intermediate,                       &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO),                &
        output_file = output_file                                            &
     )
@@ -2690 +2706 @@
        grid              = scalars_north_grid,                              &
        intermediate_grid = scalars_north_intermediate,                      &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO),                &
        output_file = output_file                                            &
     )
@@ -2702 +2719 @@
        grid              = scalars_south_grid,                              &
        intermediate_grid = scalars_south_intermediate,                      &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO),                &
        output_file = output_file                                            &
     )
@@ -2714 +2732 @@
        grid              = scalars_top_grid,                                &
        intermediate_grid = scalars_top_intermediate,                        &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO),                &
        output_file = output_file                                            &
     )
@@ -2727 +2746 @@
        grid              = palm_grid,                                       &
        intermediate_grid = palm_intermediate,                               &
-       is_profile = (TRIM(ic_mode) == 'profile')                            &
-    )
-    IF (TRIM(ic_mode) == 'profile')  THEN
-       output_var_table(9)%averaging_grid => averaged_initial_scalar_profile
-    ENDIF
+       is_profile = (TRIM(cfg%ic_mode) == CFG_INIT_PROFILE)                 &
+    )
 
     output_var_table(10) = init_nc_var(                                     &
@@ -2742 +2758 @@
        output_file       = output_file,                                     &
        grid              = scalars_west_grid,                               &
-       intermediate_grid = scalars_west_intermediate                        &
+       intermediate_grid = scalars_west_intermediate,                       &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO)                 &
     )
 
@@ -2754 +2771 @@
        output_file       = output_file,                                     &
        grid              = scalars_east_grid,                               &
-       intermediate_grid = scalars_east_intermediate                        &
+       intermediate_grid = scalars_east_intermediate,                       &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO)                 &
     )
 
@@ -2766 +2784 @@
        output_file       = output_file,                                     &
        grid              = scalars_north_grid,                              &
-       intermediate_grid = scalars_north_intermediate                       &
+       intermediate_grid = scalars_north_intermediate,                      &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO)                 &
     )
 
@@ -2778 +2797 @@
        output_file       = output_file,                                     &
        grid              = scalars_south_grid,                              &
-       intermediate_grid = scalars_south_intermediate                       &
+       intermediate_grid = scalars_south_intermediate,                      &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO)                 &
     )
 
@@ -2790 +2810 @@
        output_file       = output_file,                                     &
        grid              = scalars_top_grid,                                &
-       intermediate_grid = scalars_top_intermediate                         &
+       intermediate_grid = scalars_top_intermediate,                        &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO)                 &
     )
 
@@ -2803 +2824 @@
        grid              = u_initial_grid,                                  &
        intermediate_grid = u_initial_intermediate,                          &
-       is_profile = (TRIM(ic_mode) == 'profile')                            &
-    )
-    IF (TRIM(ic_mode) == 'profile')  THEN
-       output_var_table(15)%averaging_grid => averaged_initial_scalar_profile
-    ENDIF 
+       is_profile = (TRIM(cfg%ic_mode) == CFG_INIT_PROFILE)                 &
+    )
 
     output_var_table(16) = init_nc_var(                                     &
@@ -2818 +2836 @@
        output_file       = output_file,                                     &
        grid              = u_west_grid,                                     &
-       intermediate_grid = u_west_intermediate                              &
+       intermediate_grid = u_west_intermediate,                             &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO)                 &
     )
 
@@ -2830 +2849 @@
        output_file       = output_file,                                     &
        grid              = u_east_grid,                                     &
-       intermediate_grid = u_east_intermediate                              &
+       intermediate_grid = u_east_intermediate,                             &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO)                 &
     )
 
@@ -2842 +2862 @@
        output_file       = output_file,                                     &
        grid              = u_north_grid,                                    &
-       intermediate_grid = u_north_intermediate                             &
+       intermediate_grid = u_north_intermediate,                            &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO)                 &
     )
 
@@ -2854 +2875 @@
        output_file       = output_file,                                     &
        grid              = u_south_grid,                                    &
-       intermediate_grid = u_south_intermediate                             &
+       intermediate_grid = u_south_intermediate,                            &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO)                 &
     )
 
@@ -2866 +2888 @@
        output_file       = output_file,                                     &
        grid              = u_top_grid,                                      &
-       intermediate_grid = u_top_intermediate                               &
+       intermediate_grid = u_top_intermediate,                              &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO)                 &
     )
 
@@ -2879 +2902 @@
        grid              = v_initial_grid,                                  &
        intermediate_grid = v_initial_intermediate,                          &
-       is_profile = (TRIM(ic_mode) == 'profile')                            &
-    )
-    IF (TRIM(ic_mode) == 'profile')  THEN
-       output_var_table(21)%averaging_grid => averaged_initial_scalar_profile
-    ENDIF
+       is_profile = (TRIM(cfg%ic_mode) == CFG_INIT_PROFILE)                 &
+    )
 
     output_var_table(22) = init_nc_var(                                     &
@@ -2894 +2914 @@
        output_file       = output_file,                                     &
        grid              = v_west_grid,                                     &
-       intermediate_grid = v_west_intermediate                              &
+       intermediate_grid = v_west_intermediate,                             &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO)                 &
     )
 
@@ -2906 +2927 @@
        output_file       = output_file,                                     &
        grid              = v_east_grid,                                     &
-       intermediate_grid = v_east_intermediate                              &
+       intermediate_grid = v_east_intermediate,                             &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO)                 &
     )
 
@@ -2918 +2940 @@
        output_file       = output_file,                                     &
        grid              = v_north_grid,                                    &
-       intermediate_grid = v_north_intermediate                             &
+       intermediate_grid = v_north_intermediate,                            &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO)                 &
     )
 
@@ -2930 +2953 @@
        output_file       = output_file,                                     &
        grid              = v_south_grid,                                    &
-       intermediate_grid = v_south_intermediate                             &
+       intermediate_grid = v_south_intermediate,                            &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO)                 &
     )
 
@@ -2942 +2966 @@
        output_file       = output_file,                                     &
        grid              = v_top_grid,                                      &
-       intermediate_grid = v_top_intermediate                               &
+       intermediate_grid = v_top_intermediate,                              &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO)                 &
     )
 
@@ -2955 +2980 @@
        grid              = w_initial_grid,                                  &
        intermediate_grid = w_initial_intermediate,                          &
-       is_profile = (TRIM(ic_mode) == 'profile')                            &
-    )
-    IF (TRIM(ic_mode) == 'profile')  THEN
-       output_var_table(27)%averaging_grid => averaged_initial_w_profile
-    ENDIF
+       is_profile = (TRIM(cfg%ic_mode) == CFG_INIT_PROFILE)                 &
+    )
 
     output_var_table(28) = init_nc_var(                                     &
@@ -2970 +2992 @@
        output_file       = output_file,                                     &
        grid              = w_west_grid,                                     &
-       intermediate_grid = w_west_intermediate                              &
+       intermediate_grid = w_west_intermediate,                             &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO)                 &
     )
 
@@ -2982 +3005 @@
        output_file       = output_file,                                     &
        grid              = w_east_grid,                                     &
-       intermediate_grid = w_east_intermediate                              &
+       intermediate_grid = w_east_intermediate,                             &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO)                 &
     )
 
@@ -2994 +3018 @@
        output_file       = output_file,                                     &
        grid              = w_north_grid,                                    &
-       intermediate_grid = w_north_intermediate                             &
+       intermediate_grid = w_north_intermediate,                            &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO)                 &
     )
 
@@ -3006 +3031 @@
        output_file       = output_file,                                     &
        grid              = w_south_grid,                                    &
-       intermediate_grid = w_south_intermediate                             &
+       intermediate_grid = w_south_intermediate,                            &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO)                 &
     )
 
@@ -3018 +3044 @@
        output_file       = output_file,                                     &
        grid              = w_top_grid,                                      &
-       intermediate_grid = w_top_intermediate                               &
+       intermediate_grid = w_top_intermediate,                              &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO)                 &
     )
 
@@ -3143 +3170 @@
        intermediate_grid = palm_intermediate                                &
     )
-    output_var_table(42)%averaging_grid => averaged_scalar_profile
+    output_var_table(42)%averaging_grid => geostrophic_scalar_profile
 
     output_var_table(43) = init_nc_var(                                     &
@@ -3153 +3180 @@
        input_id          = 1_iwp,                                           &
        output_file       = output_file,                                     &
-       grid              = averaged_scalar_profile,                         &
-       intermediate_grid = averaged_scalar_profile                          &
+       grid              = geostrophic_scalar_profile,                         &
+       intermediate_grid = geostrophic_scalar_profile                          &
     )
 
@@ -3165 +3192 @@
        input_id          = 1_iwp,                                           &
        output_file       = output_file,                                     &
-       grid              = averaged_scalar_profile,                         &
-       intermediate_grid = averaged_scalar_profile                          &
+       grid              = geostrophic_scalar_profile,                         &
+       intermediate_grid = geostrophic_scalar_profile                          &
     )
 
@@ -3177 +3204 @@
        input_id          = 1_iwp,                                           &
        output_file       = output_file,                                     &
-       grid              = averaged_scalar_profile,                         &
-       intermediate_grid = averaged_scalar_profile                          &
+       grid              = geostrophic_scalar_profile,                         &
+       intermediate_grid = geostrophic_scalar_profile                          &
     )
     output_var_table(45)%to_be_processed = ls_forcing_variables_required
@@ -3190 +3217 @@
        input_id          = 1_iwp,                                           &
        output_file       = output_file,                                     &
-       grid              = averaged_scalar_profile,                         &
-       intermediate_grid = averaged_scalar_profile                          &
+       grid              = geostrophic_scalar_profile,                         &
+       intermediate_grid = geostrophic_scalar_profile                          &
     )
     output_var_table(46)%to_be_processed = ls_forcing_variables_required
@@ -3203 +3230 @@
        input_id          = 1_iwp,                                           &
        output_file       = output_file,                                     &
-       grid              = averaged_scalar_profile,                         &
-       intermediate_grid = averaged_scalar_profile                          &
+       grid              = geostrophic_scalar_profile,                         &
+       intermediate_grid = geostrophic_scalar_profile                          &
     )
     output_var_table(47)%to_be_processed = ls_forcing_variables_required
@@ -3216 +3243 @@
        input_id          = 1_iwp,                                           &
        output_file       = output_file,                                     &
-       grid              = averaged_w_profile,                              &
-       intermediate_grid = averaged_w_profile                               &
+       grid              = geostrophic_w_profile,                              &
+       intermediate_grid = geostrophic_w_profile                               &
     )
     output_var_table(48)%to_be_processed = ls_forcing_variables_required
@@ -3230 +3257 @@
        input_id          = 1_iwp,                                           &
        output_file       = output_file,                                     &
-       grid              = averaged_scalar_profile,                         &
-       intermediate_grid = averaged_scalar_profile                          &
+       grid              = geostrophic_scalar_profile,                         &
+       intermediate_grid = geostrophic_scalar_profile                          &
     )
     output_var_table(49)%to_be_processed = ls_forcing_variables_required
@@ -3243 +3270 @@
        input_id          = 1_iwp,                                           &
        output_file       = output_file,                                     &
-       grid              = averaged_scalar_profile,                         &
-       intermediate_grid = averaged_scalar_profile                          &
+       grid              = geostrophic_scalar_profile,                         &
+       intermediate_grid = geostrophic_scalar_profile                          &
     )
     output_var_table(50)%to_be_processed = ls_forcing_variables_required
@@ -3256 +3283 @@
        input_id          = 1_iwp,                                           &
        output_file       = output_file,                                     &
-       grid              = averaged_scalar_profile,                         &
-       intermediate_grid = averaged_scalar_profile                          &
+       grid              = geostrophic_scalar_profile,                         &
+       intermediate_grid = geostrophic_scalar_profile                          &
     )
     output_var_table(51)%to_be_processed = ls_forcing_variables_required
@@ -3269 +3296 @@
        input_id          = 3_iwp,                                           &
        output_file       = output_file,                                     &
-       grid              = averaged_scalar_profile,                         &
-       intermediate_grid = averaged_scalar_profile                          &
+       grid              = geostrophic_scalar_profile,                         &
+       intermediate_grid = geostrophic_scalar_profile                          &
     )
     output_var_table(52)%to_be_processed = ls_forcing_variables_required
@@ -3283 +3310 @@
        input_id          = 3_iwp,                                           &
        output_file       = output_file,                                     &
-       grid              = averaged_scalar_profile,                         &
-       intermediate_grid = averaged_scalar_profile                          &
+       grid              = geostrophic_scalar_profile,                         &
+       intermediate_grid = geostrophic_scalar_profile                          &
     )
     output_var_table(53)%to_be_processed = ls_forcing_variables_required
@@ -3296 +3323 @@
        input_id          = 3_iwp,                                           &
        output_file       = output_file,                                     &
-       grid              = averaged_scalar_profile,                         &
-       intermediate_grid = averaged_scalar_profile                          &
+       grid              = geostrophic_scalar_profile,                         &
+       intermediate_grid = geostrophic_scalar_profile                          &
     )
     output_var_table(54)%to_be_processed = ls_forcing_variables_required
@@ -3309 +3336 @@
        input_id          = 1_iwp,                                           &
        output_file       = output_file,                                     &
-       grid              = averaged_scalar_profile,                         &
-       intermediate_grid = averaged_scalar_profile                          &
+       grid              = geostrophic_scalar_profile,                         &
+       intermediate_grid = geostrophic_scalar_profile                          &
     )
     output_var_table(55)%to_be_processed = ls_forcing_variables_required
@@ -3316 +3343 @@
 
     output_var_table(56) = init_nc_var(                                     &
-       name              = 'internal_density_centre',                              &
+       name              = 'internal_density_centre',                       &
        std_name          = "",                                              &
        long_name         = "",                                              &
@@ -3323 +3350 @@
        input_id          = 4_iwp,                                           &
        output_file       = output_file,                                     &
-       grid              = averaged_scalar_profile,                         &
-       intermediate_grid = averaged_scalar_profile                          &
-    )
-    output_var_table(56)%averaging_grid => averaged_scalar_profile
+       grid              = geostrophic_scalar_profile,                         &
+       intermediate_grid = geostrophic_scalar_profile                          &
+    )
+    output_var_table(56)%averaging_grid => geostrophic_scalar_profile
 
 
     output_var_table(57) = init_nc_var(                                     &
-       name              = 'internal_density_north',                       &
+       name              = 'internal_density_north',                        &
        std_name          = "",                                              &
        long_name         = "",                                              &
@@ -3337 +3364 @@
        input_id          = 4_iwp,                                           &
        output_file       = output_file,                                     &
-       grid              = north_averaged_scalar_profile,                   &
-       intermediate_grid = north_averaged_scalar_profile                    &
-    )
-    output_var_table(57)%averaging_grid => north_averaged_scalar_profile
+       grid              = north_geostrophic_scalar_profile,                   &
+       intermediate_grid = north_geostrophic_scalar_profile                    &
+    )
+    output_var_table(57)%averaging_grid => north_geostrophic_scalar_profile
     output_var_table(57)%to_be_processed = .NOT. cfg%ug_defined_by_user
 
 
     output_var_table(58) = init_nc_var(                                     &
-       name              = 'internal_density_south',                       &
+       name              = 'internal_density_south',                        &
        std_name          = "",                                              &
        long_name         = "",                                              &
@@ -3352 +3379 @@
        input_id          = 4_iwp,                                           &
        output_file       = output_file,                                     &
-       grid              = south_averaged_scalar_profile,                   &
-       intermediate_grid = south_averaged_scalar_profile                    &
-    )
-    output_var_table(58)%averaging_grid => south_averaged_scalar_profile
+       grid              = south_geostrophic_scalar_profile,                   &
+       intermediate_grid = south_geostrophic_scalar_profile                    &
+    )
+    output_var_table(58)%averaging_grid => south_geostrophic_scalar_profile
     output_var_table(58)%to_be_processed = .NOT. cfg%ug_defined_by_user
 
 
     output_var_table(59) = init_nc_var(                                     &
-       name              = 'internal_density_east',                        &
+       name              = 'internal_density_east',                         &
        std_name          = "",                                              &
        long_name         = "",                                              &
@@ -3367 +3394 @@
        input_id          = 4_iwp,                                           &
        output_file       = output_file,                                     &
-       grid              = east_averaged_scalar_profile,                    &
-       intermediate_grid = east_averaged_scalar_profile                     &
-    )
-    output_var_table(59)%averaging_grid => east_averaged_scalar_profile
+       grid              = east_geostrophic_scalar_profile,                    &
+       intermediate_grid = east_geostrophic_scalar_profile                     &
+    )
+    output_var_table(59)%averaging_grid => east_geostrophic_scalar_profile
     output_var_table(59)%to_be_processed = .NOT. cfg%ug_defined_by_user
 
 
     output_var_table(60) = init_nc_var(                                     &
-       name              = 'internal_density_west',                        &
+       name              = 'internal_density_west',                         &
        std_name          = "",                                              &
        long_name         = "",                                              &
@@ -3382 +3409 @@
        input_id          = 4_iwp,                                           &
        output_file       = output_file,                                     &
-       grid              = west_averaged_scalar_profile,                    &
-       intermediate_grid = west_averaged_scalar_profile                     &
-    )
-    output_var_table(60)%averaging_grid => west_averaged_scalar_profile
+       grid              = west_geostrophic_scalar_profile,                    &
+       intermediate_grid = west_geostrophic_scalar_profile                     &
+    )
+    output_var_table(60)%averaging_grid => west_geostrophic_scalar_profile
     output_var_table(60)%to_be_processed = .NOT. cfg%ug_defined_by_user
 
@@ -3396 +3423 @@
        input_id          = 2_iwp,                                           &
        output_file       = output_file,                                     &
-       grid              = north_averaged_scalar_profile,                   &
-       intermediate_grid = north_averaged_scalar_profile                    &
-    )
-    output_var_table(61)%averaging_grid => north_averaged_scalar_profile
+       grid              = north_geostrophic_scalar_profile,                   &
+       intermediate_grid = north_geostrophic_scalar_profile                    &
+    )
+    output_var_table(61)%averaging_grid => north_geostrophic_scalar_profile
     output_var_table(61)%to_be_processed = .NOT. cfg%ug_defined_by_user
 
@@ -3411 +3438 @@
        input_id          = 2_iwp,                                           &
        output_file       = output_file,                                     &
-       grid              = south_averaged_scalar_profile,                   &
-       intermediate_grid = south_averaged_scalar_profile                    &
-    )
-    output_var_table(62)%averaging_grid => south_averaged_scalar_profile
+       grid              = south_geostrophic_scalar_profile,                   &
+       intermediate_grid = south_geostrophic_scalar_profile                    &
+    )
+    output_var_table(62)%averaging_grid => south_geostrophic_scalar_profile
     output_var_table(62)%to_be_processed = .NOT. cfg%ug_defined_by_user
 
@@ -3426 +3453 @@
        input_id          = 2_iwp,                                           &
        output_file       = output_file,                                     &
-       grid              = east_averaged_scalar_profile,                    &
-       intermediate_grid = east_averaged_scalar_profile                     &
-    )
-    output_var_table(63)%averaging_grid => east_averaged_scalar_profile
+       grid              = east_geostrophic_scalar_profile,                    &
+       intermediate_grid = east_geostrophic_scalar_profile                     &
+    )
+    output_var_table(63)%averaging_grid => east_geostrophic_scalar_profile
     output_var_table(63)%to_be_processed = .NOT. cfg%ug_defined_by_user
 
@@ -3441 +3468 @@
        input_id          = 2_iwp,                                           &
        output_file       = output_file,                                     &
-       grid              = west_averaged_scalar_profile,                    &
-       intermediate_grid = west_averaged_scalar_profile                     &
-    )
-    output_var_table(64)%averaging_grid => west_averaged_scalar_profile
+       grid              = west_geostrophic_scalar_profile,                    &
+       intermediate_grid = west_geostrophic_scalar_profile                     &
+    )
+    output_var_table(64)%averaging_grid => west_geostrophic_scalar_profile
     output_var_table(64)%to_be_processed = .NOT. cfg%ug_defined_by_user
 
-    output_var_table(65) = init_nc_var(                                      &
+    output_var_table(65) = init_nc_var(                                     &
        name              = 'init_atmosphere_qc',                            &
        std_name          = "",                                              &
@@ -3457 +3484 @@
        grid              = palm_grid,                                       &
        intermediate_grid = palm_intermediate,                               &
-       is_profile = (TRIM(ic_mode) == 'profile')                            &
-    )
-    IF (TRIM(ic_mode) == 'profile')  THEN
-       output_var_table(65)%averaging_grid => averaged_initial_scalar_profile
-    ENDIF
+       is_profile = (TRIM(cfg%ic_mode) == CFG_INIT_PROFILE)                 &
+    )
     output_var_table(65)%is_optional = .TRUE.
 
@@ -3473 +3497 @@
        output_file       = output_file,                                     &
        grid              = scalars_west_grid,                               &
-       intermediate_grid = scalars_west_intermediate                        &
+       intermediate_grid = scalars_west_intermediate,                       &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO)                 &
     )
     output_var_table(66)%is_optional = .TRUE.
@@ -3486 +3511 @@
        output_file       = output_file,                                     &
        grid              = scalars_east_grid,                               &
-       intermediate_grid = scalars_east_intermediate                        &
+       intermediate_grid = scalars_east_intermediate,                       &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO)                 &
     )
     output_var_table(67)%is_optional = .TRUE.
@@ -3499 +3525 @@
        output_file       = output_file,                                     &
        grid              = scalars_north_grid,                              &
-       intermediate_grid = scalars_north_intermediate                       &
+       intermediate_grid = scalars_north_intermediate,                      &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO)                 &
     )
     output_var_table(68)%is_optional = .TRUE.
@@ -3512 +3539 @@
        output_file       = output_file,                                     &
        grid              = scalars_south_grid,                              &
-       intermediate_grid = scalars_south_intermediate                       &
+       intermediate_grid = scalars_south_intermediate,                      &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO)                 &
     )
     output_var_table(69)%is_optional = .TRUE.
@@ -3525 +3553 @@
        output_file       = output_file,                                     &
        grid              = scalars_top_grid,                                &
-       intermediate_grid = scalars_top_intermediate                         &
+       intermediate_grid = scalars_top_intermediate,                        &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO)                 &
     )
     output_var_table(70)%is_optional = .TRUE.
@@ -3539 +3568 @@
        grid              = palm_grid,                                       &
        intermediate_grid = palm_intermediate,                               &
-       is_profile = (TRIM(ic_mode) == 'profile')                            &
-    )
-    IF (TRIM(ic_mode) == 'profile')  THEN
-       output_var_table(71)%averaging_grid => averaged_initial_scalar_profile
-    ENDIF
+       is_profile = (TRIM(cfg%ic_mode) == CFG_INIT_PROFILE)                 &
+    )
     output_var_table(71)%is_optional = .TRUE.
 
@@ -3555 +3581 @@
        output_file       = output_file,                                     &
        grid              = scalars_west_grid,                               &
-       intermediate_grid = scalars_west_intermediate                        &
+       intermediate_grid = scalars_west_intermediate,                       &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO)                 &
     )
     output_var_table(72)%is_optional = .TRUE.
@@ -3568 +3595 @@
        output_file       = output_file,                                     &
        grid              = scalars_east_grid,                               &
-       intermediate_grid = scalars_east_intermediate                        &
+       intermediate_grid = scalars_east_intermediate,                       &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO)                 &
     )
     output_var_table(73)%is_optional = .TRUE.
@@ -3581 +3609 @@
        output_file       = output_file,                                     &
        grid              = scalars_north_grid,                              &
-       intermediate_grid = scalars_north_intermediate                       &
+       intermediate_grid = scalars_north_intermediate,                      &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO)                 &
     )
     output_var_table(74)%is_optional = .TRUE.
@@ -3594 +3623 @@
        output_file       = output_file,                                     &
        grid              = scalars_south_grid,                              &
-       intermediate_grid = scalars_south_intermediate                       &
+       intermediate_grid = scalars_south_intermediate,                      &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO)                 &
     )
     output_var_table(75)%is_optional = .TRUE.
@@ -3607 +3637 @@
        output_file       = output_file,                                     &
        grid              = scalars_top_grid,                                &
-       intermediate_grid = scalars_top_intermediate                         &
+       intermediate_grid = scalars_top_intermediate,                        &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO)                 &
     )
     output_var_table(76)%is_optional = .TRUE.
@@ -3614 +3645 @@
        name              = 'init_atmosphere_qr',                            &
        std_name          = "",                                              &
-       long_name         = "initial rain water mixture fraction",            &
+       long_name         = "initial rain water mixture fraction",           &
        units             = "kg/kg",                                         &
        kind              = "init scalar",                                   &
@@ -3621 +3652 @@
        grid              = palm_grid,                                       &
        intermediate_grid = palm_intermediate,                               &
-       is_profile = (TRIM(ic_mode) == 'profile')                            &
-    )
-    IF (TRIM(ic_mode) == 'profile')  THEN
-       output_var_table(77)%averaging_grid => averaged_initial_scalar_profile
-    ENDIF
+       is_profile = (TRIM(cfg%ic_mode) == CFG_INIT_PROFILE)                 &
+    )
     output_var_table(77)%is_optional = .TRUE.
 
@@ -3637 +3665 @@
        output_file       = output_file,                                     &
        grid              = scalars_west_grid,                               &
-       intermediate_grid = scalars_west_intermediate                        &
+       intermediate_grid = scalars_west_intermediate,                       &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO)                 &
     )
     output_var_table(78)%is_optional = .TRUE.
@@ -3650 +3679 @@
        output_file       = output_file,                                     &
        grid              = scalars_east_grid,                               &
-       intermediate_grid = scalars_east_intermediate                        &
+       intermediate_grid = scalars_east_intermediate,                       &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO)                 &
     )
     output_var_table(79)%is_optional = .TRUE.
@@ -3663 +3693 @@
        output_file       = output_file,                                     &
        grid              = scalars_north_grid,                              &
-       intermediate_grid = scalars_north_intermediate                       &
+       intermediate_grid = scalars_north_intermediate,                      &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO)                 &
     )
     output_var_table(80)%is_optional = .TRUE.
@@ -3676 +3707 @@
        output_file       = output_file,                                     &
        grid              = scalars_south_grid,                              &
-       intermediate_grid = scalars_south_intermediate                       &
+       intermediate_grid = scalars_south_intermediate,                      &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO)                 &
     )
     output_var_table(81)%is_optional = .TRUE.
@@ -3689 +3721 @@
        output_file       = output_file,                                     &
        grid              = scalars_top_grid,                                &
-       intermediate_grid = scalars_top_intermediate                         &
+       intermediate_grid = scalars_top_intermediate,                        &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO)                 &
     )
     output_var_table(82)%is_optional = .TRUE.
@@ -3696 +3729 @@
        name              = 'init_atmosphere_qs',                            &
        std_name          = "",                                              &
-       long_name         = "initial snow mixture fraction",            &
+       long_name         = "initial snow mixture fraction",                 &
        units             = "kg/kg",                                         &
        kind              = "init scalar",                                   &
@@ -3703 +3736 @@
        grid              = palm_grid,                                       &
        intermediate_grid = palm_intermediate,                               &
-       is_profile = (TRIM(ic_mode) == 'profile')                            &
-    )
-    IF (TRIM(ic_mode) == 'profile')  THEN
-       output_var_table(83)%averaging_grid => averaged_initial_scalar_profile
-    ENDIF
+       is_profile = (TRIM(cfg%ic_mode) == CFG_INIT_PROFILE)                 &
+    )
     output_var_table(83)%is_optional = .TRUE.
 
@@ -3719 +3749 @@
        output_file       = output_file,                                     &
        grid              = scalars_west_grid,                               &
-       intermediate_grid = scalars_west_intermediate                        &
+       intermediate_grid = scalars_west_intermediate,                       &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO)                 &
     )
     output_var_table(84)%is_optional = .TRUE.
@@ -3732 +3763 @@
        output_file       = output_file,                                     &
        grid              = scalars_east_grid,                               &
-       intermediate_grid = scalars_east_intermediate                        &
+       intermediate_grid = scalars_east_intermediate,                       &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO)                 &
     )
     output_var_table(85)%is_optional = .TRUE.
@@ -3745 +3777 @@
        output_file       = output_file,                                     &
        grid              = scalars_north_grid,                              &
-       intermediate_grid = scalars_north_intermediate                       &
+       intermediate_grid = scalars_north_intermediate,                      &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO)                 &
     )
     output_var_table(86)%is_optional = .TRUE.
@@ -3758 +3791 @@
        output_file       = output_file,                                     &
        grid              = scalars_south_grid,                              &
-       intermediate_grid = scalars_south_intermediate                       &
+       intermediate_grid = scalars_south_intermediate,                      &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO)                 &
     )
     output_var_table(87)%is_optional = .TRUE.
@@ -3771 +3805 @@
        output_file       = output_file,                                     &
        grid              = scalars_top_grid,                                &
-       intermediate_grid = scalars_top_intermediate                         &
+       intermediate_grid = scalars_top_intermediate,                        &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO)                 &
     )
     output_var_table(88)%is_optional = .TRUE.
@@ -3778 +3813 @@
        name              = 'init_atmosphere_qg',                            &
        std_name          = "",                                              &
-       long_name         = "initial graupel mixture fraction",            &
+       long_name         = "initial graupel mixture fraction",              &
        units             = "kg/kg",                                         &
        kind              = "init scalar",                                   &
@@ -3785 +3820 @@
        grid              = palm_grid,                                       &
        intermediate_grid = palm_intermediate,                               &
-       is_profile = (TRIM(ic_mode) == 'profile')                            &
-    )
-    IF (TRIM(ic_mode) == 'profile')  THEN
-       output_var_table(89)%averaging_grid => averaged_initial_scalar_profile
-    ENDIF
+       is_profile = (TRIM(cfg%ic_mode) == CFG_INIT_PROFILE)                 &
+    )
     output_var_table(89)%is_optional = .TRUE.
 
@@ -3801 +3833 @@
        output_file       = output_file,                                     &
        grid              = scalars_west_grid,                               &
-       intermediate_grid = scalars_west_intermediate                        &
+       intermediate_grid = scalars_west_intermediate,                       &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO)                 &
     )
     output_var_table(90)%is_optional = .TRUE.
@@ -3814 +3847 @@
        output_file       = output_file,                                     &
        grid              = scalars_east_grid,                               &
-       intermediate_grid = scalars_east_intermediate                        &
+       intermediate_grid = scalars_east_intermediate,                       &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO)                 &
     )
     output_var_table(91)%is_optional = .TRUE.
@@ -3827 +3861 @@
        output_file       = output_file,                                     &
        grid              = scalars_north_grid,                              &
-       intermediate_grid = scalars_north_intermediate                       &
+       intermediate_grid = scalars_north_intermediate,                      &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO)                 &
     )
     output_var_table(92)%is_optional = .TRUE.
@@ -3840 +3875 @@
        output_file       = output_file,                                     &
        grid              = scalars_south_grid,                              &
-       intermediate_grid = scalars_south_intermediate                       &
+       intermediate_grid = scalars_south_intermediate,                      &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO)                 &
     )
     output_var_table(93)%is_optional = .TRUE.
@@ -3853 +3889 @@
        output_file       = output_file,                                     &
        grid              = scalars_top_grid,                                &
-       intermediate_grid = scalars_top_intermediate                         &
+       intermediate_grid = scalars_top_intermediate,                        &
+       is_profile = (TRIM(cfg%bc_mode) == CFG_FORCING_HOMO)                 &
     )
     output_var_table(94)%is_optional = .TRUE.
@@ -3969 +4006 @@
           var%task            = "interpolate_3d"
 
-       CASE( 'init scalar profile', 'init u profile', 'init v profile')
+       CASE( 'init soil profile' )
+          var%nt              = 1
+          var%lod             = 1
+          var%ndim            = 1
+          var%dimids(1)       = output_file%dimids_soil(3)
+          var%dimvarids(1)    = output_file%dimvarids_soil(3)
+          var%to_be_processed = init_variables_required
+          var%is_internal     = .FALSE.
+          var%task            = "average levels"
+          var%averaging_grid  => averaged_soil_profile
+
+       CASE( 'init scalar profile', 'init u profile', 'init v profile' )
           var%nt              = 1
           var%lod             = 1
@@ -3978 +4026 @@
           var%is_internal     = .FALSE.
           var%task            = "average profile"
+          var%averaging_grid  => averaged_scalar_profile
 
        CASE( 'init w profile')
@@ -3988 +4037 @@
           var%is_internal     = .FALSE.
           var%task            = "average profile"
+          var%averaging_grid  => averaged_w_profile
 
        CASE( 'surface forcing' )
-          var%lod             = -1
+          var%lod             = 2
           var%ndim            = 3
           var%dimids(3)       = output_file%dimid_time
@@ -4001 +4051 @@
 
        CASE( 'left scalar', 'right scalar')
-          var%lod             = -1
+          var%lod             = 2
           var%ndim            = 3
           var%dimids(3)       = output_file%dimid_time
@@ -4014 +4064 @@
 
        CASE( 'north scalar', 'south scalar')
-          var%lod             = -1
+          var%lod             = 2
           var%ndim            = 3
           var%dimids(3)       = output_file%dimid_time
@@ -4027 +4077 @@
 
        CASE( 'top scalar', 'top w' )
-          var%lod             = -1
+          var%lod             = 2
           var%ndim            = 3
           var%dimids(3)       = output_file%dimid_time
@@ -4040 +4090 @@
 
        CASE( 'left u', 'right u' )
-          var%lod             = -1
+          var%lod             = 2
           var%ndim            = 3
           var%dimids(3)       = output_file%dimid_time
@@ -4053 +4103 @@
 
        CASE( 'north u', 'south u' )
-          var%lod             = -1
+          var%lod             = 2
           var%ndim            = 3
           var%dimids(3)       = output_file%dimid_time    !t
@@ -4066 +4116 @@
 
        CASE( 'top u' )
-          var%lod             = -1
+          var%lod             = 2
           var%ndim            = 3
           var%dimids(3)       = output_file%dimid_time    !t
@@ -4079 +4129 @@
 
        CASE( 'left v', 'right v' )
-          var%lod             = -1
+          var%lod             = 2
           var%ndim            = 3
           var%dimids(3)       = output_file%dimid_time
@@ -4092 +4142 @@
 
        CASE( 'north v', 'south v' )
-          var%lod             = -1
+          var%lod             = 2
           var%ndim            = 3
           var%dimids(3)       = output_file%dimid_time    !t
@@ -4105 +4155 @@
 
        CASE( 'top v' )
-          var%lod             = -1
+          var%lod             = 2
           var%ndim            = 3
           var%dimids(3)       = output_file%dimid_time    !t
@@ -4118 +4168 @@
 
        CASE( 'left w', 'right w')
-          var%lod             = -1
+          var%lod             = 2
           var%ndim            = 3
           var%dimids(3)       = output_file%dimid_time
@@ -4131 +4181 @@
 
        CASE( 'north w', 'south w' )
-          var%lod             = -1
+          var%lod             = 2
           var%ndim            = 3
           var%dimids(3)       = output_file%dimid_time    !t
@@ -4143 +4193 @@
           var%task            = "interpolate_3d"
 
+       CASE( 'left scalar profile', 'right scalar profile',                    &
+             'north scalar profile', 'south scalar profile',                   &
+             'left u profile', 'right u profile',                              &
+             'north u profile', 'south u profile',                             &
+             'left v profile', 'right v profile',                              &
+             'north v profile', 'south v profile' )
+          var%lod             = 1
+          var%ndim            = 2
+          var%dimids(2)       = output_file%dimid_time    !t
+          var%dimids(1)       = output_file%dimids_scl(3) !z
+          var%dimvarids(2)    = output_file%dimvarid_time
+          var%dimvarids(1)    = output_file%dimvarids_scl(3)
+          var%to_be_processed = boundary_variables_required
+          var%is_internal     = .FALSE.
+          var%task            = "average profile"
+          var%averaging_grid  => averaged_scalar_profile
+
+       CASE( 'top scalar profile', 'top u profile', 'top v profile' )
+          var%lod             = 0
+          var%ndim            = 1
+          var%dimids(1)       = output_file%dimid_time    !t
+          var%dimvarids(1)    = output_file%dimvarid_time
+          var%to_be_processed = boundary_variables_required
+          var%is_internal     = .FALSE.
+          var%task            = "average profile"
+          var%averaging_grid  => averaged_scalar_top_point
+
+       CASE( 'left w profile', 'right w profile',                              &
+             'north w profile', 'south w profile' )
+          var%lod             = 1
+          var%ndim            = 2
+          var%dimids(2)       = output_file%dimid_time    !t
+          var%dimids(1)       = output_file%dimids_vel(3) !z
+          var%dimvarids(2)    = output_file%dimvarid_time
+          var%dimvarids(1)    = output_file%dimvarids_vel(3)
+          var%to_be_processed = boundary_variables_required
+          var%is_internal     = .FALSE.
+          var%task            = "average profile"
+          var%averaging_grid  => averaged_w_profile
+
+       CASE( 'top w profile' )
+          var%lod             = 0
+          var%ndim            = 1
+          var%dimids(1)       = output_file%dimid_time    !t
+          var%dimvarids(1)    = output_file%dimvarid_time
+          var%to_be_processed = boundary_variables_required
+          var%is_internal     = .FALSE.
+          var%task            = "average profile"
+          var%averaging_grid  => averaged_w_top_point
+
        CASE( 'time series' )
           var%lod             = 0
@@ -4153 +4253 @@
 
        CASE( 'constant scalar profile' )
-          var%lod             = -1
+          var%lod             = 1
           var%ndim            = 2
           var%dimids(2)       = output_file%dimid_time    !t
@@ -4164 +4264 @@
 
        CASE( 'large-scale scalar forcing' )
-          var%lod             = -1
+          var%lod             = 1
           var%ndim            = 2
           var%dimids(2)       = output_file%dimid_time    !t
@@ -4175 +4275 @@
 
        CASE( 'geostrophic' )
-          var%lod             = -1
+          var%lod             = 1
           var%ndim            = 2
           var%dimids(2)       = output_file%dimid_time    !t
@@ -4186 +4286 @@
 
        CASE( 'large-scale w forcing' )
-          var%lod             = -1
+          var%lod             = 1
           var%ndim            = 2
           var%dimids(2)       = output_file%dimid_time    !t
@@ -4208 +4308 @@
 
        CASE DEFAULT
-           message = "Variable kind '" // TRIM(kind) // "' not recognized."
+           message = "Variable kind '" // TRIM(out_var_kind) // "' not recognized."
            CALL inifor_abort ('init_nc_var', message)
 
@@ -4350 +4450 @@
              DO  i = 1, nx
 
-                CALL get_basic_state(cosmo_grid%hfl(i,j,:), BETA, P_SL, T_SL,&
-                                     RD, G, basic_state_pressure)
+                CALL get_basic_state( cosmo_grid%hfl(i,j,:), BETA, P_SL, T_SL, &
+                                      RD, G, basic_state_pressure )
 
 !

palm/trunk/UTIL/inifor/src/inifor_io.f90

@@ -26 +26 @@
 ! -----------------
 ! $Id$
+! New command-line options for soil profile initialization and homogeneous
+!    boundary conditions
+!
+!
+! 4659 2020-08-31 11:21:17Z eckhard
 ! New command-line option '--precipitation' enables output of precipitation
 !    surface forcing
@@ -160 +165 @@
     USE inifor_control
     USE inifor_defs,                                                           &
-        ONLY:  DATE, SNAME, PATH, PI, TO_RADIANS, TO_DEGREES, VERSION,         &
+        ONLY:  CFG_INIT_PROFILE, CFG_INIT_VOLUME,                              &
+               CFG_INIT_SOIL_PROFILE, CFG_INIT_SOIL_VOLUME,                    &
+               CFG_FORCING_HETERO, CFG_FORCING_HOMO, CFG_FORCING_NUDGING,      &
+               DATE, SNAME, PATH, PI, TO_RADIANS, TO_DEGREES, VERSION,         &
                NC_DEPTH_DIM_IDX, NC_DEPTH_NAME, NC_HHL_NAME, NC_RLAT_NAME,     &
                NC_RLON_NAME, NC_ROTATED_POLE_NAME, NC_POLE_LATITUDE_NAME,      &
@@ -562 +570 @@
                 cfg%namelist_file = TRIM(arg)
 
-             CASE( '-mode', '-i', '--init-mode' )
+             CASE( '-i', '--init-mode' )
                 CALL get_option_argument( i, arg )
                 cfg%ic_mode = TRIM(arg)
+
+             CASE( '--soil-init-mode' )
+                CALL get_option_argument( i, arg )
+                cfg%isc_mode = TRIM(arg)
 
              CASE( '-f', '--forcing-mode' )
@@ -733 +745 @@
 
     SELECT CASE( TRIM( cfg%ic_mode ) )
-       CASE( 'profile', 'volume' )
+       CASE( CFG_INIT_PROFILE, CFG_INIT_VOLUME )
        CASE DEFAULT
           message = "Initialization mode '" // TRIM( cfg%ic_mode ) //&
                     "' is not supported. " //&
-                    "Please select either 'profile' or 'volume', " //&
-                    "or omit the -i/--init-mode/-mode option entirely, which corresponds "//&
+                    "Please select either '" // CFG_INIT_PROFILE //"' or '" // &
+                    CFG_INIT_VOLUME //"', " //&
+                    "or omit the -i/--init-mode option entirely, which corresponds "//&
                     "to the latter."
           CALL inifor_abort( 'validate_config', message )
     END SELECT
 
+    SELECT CASE( TRIM( cfg%isc_mode ) )
+       CASE( CFG_INIT_SOIL_PROFILE, CFG_INIT_SOIL_VOLUME )
+       CASE DEFAULT
+          message = "Soil initialization mode '" // TRIM( cfg%isc_mode ) //&
+                    "' is not supported. " //&
+                    "Please select either '" // CFG_INIT_SOIL_PROFILE //"' or '" // &
+                    CFG_INIT_SOIL_VOLUME //"', " //&
+                    "or omit the --soil-init-mode option entirely, which corresponds "//&
+                    "to the latter."
+          CALL inifor_abort( 'validate_config', message )
+    END SELECT
+
     SELECT CASE( TRIM(cfg%bc_mode) )
-       CASE( 'real', 'ideal')
+       CASE( CFG_FORCING_HOMO, CFG_FORCING_HETERO, CFG_FORCING_NUDGING )
        CASE DEFAULT
           message = "Forcing mode '" // TRIM( cfg%bc_mode ) //& 
                     "' is not supported. " //&
-                    "Please select either 'real' or 'ideal', " //&
+                    "Please select either '" // CFG_FORCING_NUDGING //          &
+                    "', '" // CFG_FORCING_HOMO // "', or '" //                 &
+                    CFG_FORCING_HETERO // "' " //&
                     "or omit the -f/--forcing-mode option entirely, which corresponds "//&
                     "to the latter."
@@ -1420 +1447 @@
     LOGICAL      ::  var_is_time_dependent
 
-    var_is_time_dependent = (                                               &
-       var%dimids( var%ndim ) == output_file%dimid_time               &
+    var_is_time_dependent = (                                                  &
+       var%dimids( var%ndim ) == output_file%dimid_time                        &
     )
 
@@ -1439 +1466 @@
     SELECT CASE (TRIM(var%kind))
        
-       CASE ( 'init scalar profile', 'init u profile', 'init v profile',       &
-              'init w profile' )
+       CASE( 'init scalar profile', 'init u profile', 'init v profile',        &
+             'init w profile', 'init soil profile' )
 
           CALL check(nf90_put_var( ncid, var%varid, array(1,1,:) ) )
 
-       CASE ( 'init soil', 'init scalar', 'init u', 'init v', 'init w' )
+       CASE( 'init soil', 'init scalar', 'init u', 'init v', 'init w' )
 
           CALL check(nf90_put_var( ncid, var%varid, array(:,:,:) ) )
@@ -1457 +1484 @@
           IF (.NOT. SIZE(array, 2) .EQ. var%dimlen(1))  THEN
              PRINT *, "inifor: update_output: Dimension ", 1, " of variable ", &
-                 TRIM(var%name), " (", var%dimlen(1),                      &
+                 TRIM(var%name), " (", var%dimlen(1),                          &
                  ") does not match the dimension of the output array (",       &
                  SIZE(array, 2), ")."
@@ -1477 +1504 @@
                                    count=count(1:ndim) ) )
           
-       CASE ( 'left u', 'right u', 'left v', 'right v' )
+       CASE( 'left u', 'right u', 'left v', 'right v' )
 
           CALL check(nf90_put_var( ncid, var%varid, array(1,:,:),              &
@@ -1483 +1510 @@
                                    count=count(1:ndim) ) )
 
-       CASE ( 'north u', 'south u', 'north v', 'south v' )
+       CASE( 'north u', 'south u', 'north v', 'south v' )
 
           CALL check(nf90_put_var( ncid, var%varid, array(:,1,:),              &
@@ -1489 +1516 @@
                                    count=count(1:ndim) ) )
 
-       CASE ( 'top u', 'top v' )
+       CASE( 'top u', 'top v' )
 
           CALL check(nf90_put_var( ncid, var%varid, array(:,:,1),              &
@@ -1495 +1522 @@
                                    count=count(1:ndim) ) )
 
-       CASE ( 'time series' )
+       CASE( 'time series' )
 
           CALL check(nf90_put_var( ncid, var%varid, array(1,1,1),              &
                                    start=start(1:ndim+1) ) )
 
-       CASE ( 'constant scalar profile', 'geostrophic' )
+       CASE( 'constant scalar profile', 'geostrophic',                         &
+             'left scalar profile', 'right scalar profile',                    &
+             'north scalar profile', 'south scalar profile',                   &
+             'left u profile', 'right u profile',                              &
+             'north u profile', 'south u profile',                             &
+             'left v profile', 'right v profile',                              &
+             'north v profile', 'south v profile',                             &
+             'top scalar profile', 'top u profile', 'top v profile',           &
+             'top w profile',                                                  &
+             'left w profile', 'right w profile',                              &
+             'north w profile', 'south w profile' )
 
           CALL check(nf90_put_var( ncid, var%varid, array(1,1,:),              &

palm/trunk/UTIL/inifor/src/inifor_transform.f90

@@ -21 +21 @@
 ! Current revisions:
 ! -----------------
-! 
-! 
+!
+!
 ! Former revisions:
 ! -----------------
 ! $Id$
+! Improved code formatting
+! 
+! 
+! 4553 2020-06-03 16:34:15Z eckhard
 ! Improved code readability and documentation
 ! 
@@ -431 +435 @@
 
     TYPE(grid_definition), INTENT(IN)          ::  avg_grid
-    REAL(wp), DIMENSION(:,:,:), INTENT(IN)     ::  source_array 
+    REAL(wp), DIMENSION(:,:,:), INTENT(IN)     ::  source_array
     REAL(wp), DIMENSION(:), INTENT(OUT)        ::  profile_array
 
@@ -445 +449 @@
                  "cosmo_pressure(" // TRIM( str( nlev ) ) //                &
                  "), profile_array(" // TRIM( str( nz ) )  // ")."
-       CALL inifor_abort('average_pressure_perturbation', message)
+       CALL inifor_abort( 'average_profile', message )
     ENDIF
     

palm/trunk/UTIL/inifor/src/inifor_types.f90

@@ -26 +26 @@
 ! -----------------
 ! $Id$
+! Added INIFOR configuration flag for soil profile initialization
+!
+!
+! 4659 2020-08-31 11:21:17Z eckhard
 ! Added flag in support of new command-line option '--precipitation'
 ! Improved code formatting
@@ -133 +137 @@
     CHARACTER(LEN=SNAME) ::  bc_mode              !< destinguishes realistic and idealistic forcing
     CHARACTER(LEN=SNAME) ::  ic_mode              !< destinguishes volume and profile initialization
+    CHARACTER(LEN=SNAME) ::  isc_mode             !< destinguishes volume and profile soil initialization
     CHARACTER(LEN=SNAME) ::  rotation_method      !< selects method for velocity rotation