Changeset 3182 for palm

Timestamp:

Jul 27, 2018 1:36:03 PM (5 years ago)

Author:

suehring

Message:

New Inifor features: grid stretching, improved command-interface, support start dates in different formats in both YYYYMMDD and YYYYMMDDHH, Ability to manually control input file prefixes (--radiation-prefix, --soil-preifx, --flow-prefix, --soilmoisture-prefix) for compatiblity with DWD forcast naming scheme; GNU-style short and long option; Prepared output of large-scale forcing profiles (no computation yet); Added preprocessor flag netcdf4 to switch output format between netCDF 3 and 4; Updated netCDF variable names and attributes to comply with PIDS v1.9; Inifor bugfixes: Improved compatibility with older Intel Intel compilers by avoiding implicit array allocation; Added origin_lon/_lat values and correct reference time in dynamic driver global attributes; corresponding PALM changes: adjustments to revised Inifor; variables names in dynamic driver adjusted; enable geostrophic forcing also in offline nested mode; variable names in LES-LES and COSMO offline nesting changed; lateral boundary flags for nesting, in- and outflow conditions renamed

Location:

palm/trunk

Files:

• SOURCE/advec_ws.f90 (modified) (19 diffs)
• SOURCE/boundary_conds.f90 (modified) (17 diffs)
• SOURCE/check_parameters.f90 (modified) (12 diffs)
• SOURCE/chemistry_model_mod.f90 (modified) (3 diffs)
• SOURCE/coriolis.f90 (modified) (9 diffs)
• SOURCE/exchange_horiz_2d.f90 (modified) (5 diffs)
• SOURCE/init_3d_model.f90 (modified) (11 diffs)
• SOURCE/init_grid.f90 (modified) (7 diffs)
• SOURCE/init_pegrid.f90 (modified) (10 diffs)
• SOURCE/large_scale_forcing_nudging_mod.f90 (modified) (41 diffs)
• SOURCE/modules.f90 (modified) (5 diffs)
• SOURCE/netcdf_data_input_mod.f90 (modified) (42 diffs)
• SOURCE/palm.f90 (modified) (2 diffs)
• SOURCE/parin.f90 (modified) (12 diffs)
• SOURCE/pmc_interface_mod.f90 (modified) (37 diffs)
• SOURCE/pmc_parent_mod.f90 (modified) (4 diffs)
• SOURCE/poismg_mod.f90 (modified) (17 diffs)
• SOURCE/poismg_noopt_mod.f90 (modified) (18 diffs)
• SOURCE/pres.f90 (modified) (7 diffs)
• SOURCE/prognostic_equations.f90 (modified) (2 diffs)
• SOURCE/read_restart_data_mod.f90 (modified) (13 diffs)
• SOURCE/sor.f90 (modified) (4 diffs)
• SOURCE/synthetic_turbulence_generator_mod.f90 (modified) (18 diffs)
• SOURCE/time_integration.f90 (modified) (8 diffs)
• SOURCE/turbulence_closure_mod.f90 (modified) (8 diffs)
• SOURCE/virtual_flight_mod.f90 (modified) (3 diffs)
• SOURCE/wind_turbine_model_mod.f90 (modified) (1 diff)
• UTIL/inifor/Makefile (modified) (2 diffs)
• UTIL/inifor/Makefile.gnu (modified) (5 diffs)
• UTIL/inifor/Makefile.ifort (modified) (4 diffs)
• UTIL/inifor/README (modified) (2 diffs)
• UTIL/inifor/src/control.f90 (modified) (3 diffs)
• UTIL/inifor/src/defs.f90 (modified) (3 diffs)
• UTIL/inifor/src/get_netcdf_variable.inc (added)
• UTIL/inifor/src/grid.f90 (modified) (93 diffs)
• UTIL/inifor/src/inifor.f90 (modified) (14 diffs)
• UTIL/inifor/src/io.f90 (modified) (21 diffs)
• UTIL/inifor/src/transform.f90 (modified) (16 diffs)
• UTIL/inifor/src/types.f90 (modified) (8 diffs)
• UTIL/inifor/src/util.f90 (modified) (4 diffs)
• UTIL/inifor/tests/test-boundaries.f90 (modified) (4 diffs)
• UTIL/inifor/tests/test-grid.f90 (modified) (4 diffs)
• UTIL/inifor/tests/test-input-files.f90 (modified) (4 diffs)
• UTIL/inifor/tests/test-interpolation.f90 (modified) (6 diffs)
• UTIL/inifor/tests/test-prototype.f90 (modified) (2 diffs)
• UTIL/inifor/tests/test-stretching.f90 (added)
• UTIL/inifor/tests/util.f90 (modified) (4 diffs)

palm/trunk/SOURCE/advec_ws.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-! 
+! Rename variables for boundary-condition flags and for offline nesting mode 
 ! 
 ! Former revisions:
@@ -460 +460 @@
 
        USE control_parameters,                                                 &
-           ONLY:  force_bound_l, force_bound_n, force_bound_r, force_bound_s,  &
-                  inflow_l, inflow_n, inflow_r, inflow_s, momentum_advec,      &
-                  nest_bound_l, nest_bound_n, nest_bound_r, nest_bound_s,      &
-                  outflow_l, outflow_n, outflow_r, outflow_s, scalar_advec
+           ONLY:  bc_dirichlet_l, bc_dirichlet_n, bc_dirichlet_r,              &
+                  bc_dirichlet_s, bc_radiation_l, bc_radiation_n,              &
+                  bc_radiation_r, bc_radiation_s, momentum_advec, scalar_advec
 
        USE indices,                                                            &
@@ -501 +500 @@
                     .OR.  .NOT. BTEST(wall_flags_0(k,j,i+2),0)                 &   
                     .OR.  .NOT. BTEST(wall_flags_0(k,j,i-1),0) )               &
-                      .OR.  ( ( inflow_l .OR. outflow_l .OR. nest_bound_l  .OR.&
-                                force_bound_l )                                &
+                      .OR.  ( ( bc_dirichlet_l .OR. bc_radiation_l )           &
                             .AND.  i == nxl   )                                &
-                      .OR.  ( ( inflow_r .OR. outflow_r .OR. nest_bound_r  .OR.&
-                                force_bound_r )                                &
+                      .OR.  ( ( bc_dirichlet_r .OR. bc_radiation_r )           &
                             .AND.  i == nxr   ) )                              &
                    THEN
@@ -520 +517 @@
                             )                                                  &
                                                     .OR.                       &
-                            ( ( inflow_r .OR. outflow_r .OR. nest_bound_r  .OR.&
-                                force_bound_r )                                &
+                            ( ( bc_dirichlet_r .OR. bc_radiation_r )           &
                               .AND. i == nxr-1 )    .OR.                       &
-                            ( ( inflow_l .OR. outflow_l .OR. nest_bound_l  .OR.&
-                                force_bound_l )                                &
+                            ( ( bc_dirichlet_l .OR. bc_radiation_l )           &
                               .AND. i == nxlu  ) )                             & ! why not nxl+1
                    THEN
@@ -542 +537 @@
                     .OR.  .NOT. BTEST(wall_flags_0(k,j+2,i),0)                 &   
                     .OR.  .NOT. BTEST(wall_flags_0(k,j-1,i),0))                &
-                      .OR.  ( ( inflow_s .OR. outflow_s .OR. nest_bound_s  .OR.&
-                                force_bound_s )                                &
+                      .OR.  ( ( bc_dirichlet_s .OR. bc_radiation_s )           &
                             .AND.  j == nys   )                                &
-                      .OR.  ( ( inflow_n .OR. outflow_n .OR. nest_bound_n  .OR.&
-                                force_bound_n )                                &
+                      .OR.  ( ( bc_dirichlet_n .OR. bc_radiation_n )           &
                             .AND.  j == nyn   ) )                              &
                    THEN
@@ -563 +556 @@
                             )                                                  &
                                                     .OR.                       &
-                            ( ( inflow_s .OR. outflow_s .OR. nest_bound_s  .OR.&
-                                force_bound_s )                                &
+                            ( ( bc_dirichlet_s .OR. bc_radiation_s )           &
                               .AND. j == nysv  )    .OR.                       & ! why not nys+1
-                            ( ( inflow_n .OR. outflow_n .OR. nest_bound_n  .OR.&
-                                force_bound_n )                                &
+                            ( ( bc_dirichlet_n .OR. bc_radiation_n )           &
                               .AND. j == nyn-1 ) )                             &          
                    THEN
@@ -655 +646 @@
 !--                WS1 (9), WS3 (10), WS5 (11)
                    IF ( .NOT. BTEST(wall_flags_0(k,j,i+1),1)  .OR.             &
-                            ( ( inflow_l .OR. outflow_l .OR. nest_bound_l  .OR.&
-                                force_bound_l )                                &
+                            ( ( bc_dirichlet_l .OR. bc_radiation_l )           &
                               .AND. i <= nxlu  )    .OR.                       &
-                            ( ( inflow_r .OR. outflow_r .OR. nest_bound_r  .OR.&
-                                force_bound_r )                                &
+                            ( ( bc_dirichlet_r .OR. bc_radiation_r )           &
                               .AND. i == nxr   ) )                             &
                    THEN
@@ -667 +656 @@
                               .NOT. BTEST(wall_flags_0(k,j,i-1),1) )           &
                                                         .OR.                   &
-                            ( ( inflow_r .OR. outflow_r .OR. nest_bound_r  .OR.&
-                                force_bound_r )                                &
+                            ( ( bc_dirichlet_r .OR. bc_radiation_r )           &
                               .AND. i == nxr-1 )    .OR.                       &
-                            ( ( inflow_l .OR. outflow_l .OR. nest_bound_l  .OR.&
-                                force_bound_l )                                &
+                            ( ( bc_dirichlet_l .OR. bc_radiation_l )           &
                               .AND. i == nxlu+1) )                             &
                    THEN
@@ -691 +678 @@
 !--                WS1 (12), WS3 (13), WS5 (14)
                    IF ( .NOT. BTEST(wall_flags_0(k,j+1,i),1)   .OR.            &
-                            ( ( inflow_s .OR. outflow_s .OR. nest_bound_s  .OR.&
-                                force_bound_s )                                &
+                            ( ( bc_dirichlet_s .OR. bc_radiation_s )           &
                               .AND. j == nys   )    .OR.                       &
-                            ( ( inflow_n .OR. outflow_n .OR. nest_bound_n  .OR.&
-                                force_bound_n )                                &
+                            ( ( bc_dirichlet_n .OR. bc_radiation_n )           &
                               .AND. j == nyn   ) )                             &
                    THEN
@@ -703 +688 @@
                               .NOT. BTEST(wall_flags_0(k,j-1,i),1) )           &
                                                         .OR.                   &
-                            ( ( inflow_s .OR. outflow_s .OR. nest_bound_s  .OR.&
-                                force_bound_s )                                &
+                            ( ( bc_dirichlet_s .OR. bc_radiation_s )           &
                               .AND. j == nysv  )    .OR.                       &
-                            ( ( inflow_n .OR. outflow_n .OR. nest_bound_n  .OR.&
-                                force_bound_n )                                &
+                            ( ( bc_dirichlet_n .OR. bc_radiation_n )           &
                               .AND. j == nyn-1 ) )                             &
                    THEN
@@ -789 +772 @@
 !--                WS1 (18), WS3 (19), WS5 (20)
                    IF ( .NOT. BTEST(wall_flags_0(k,j,i+1),2)  .OR.             &
-                            ( ( inflow_l .OR. outflow_l .OR. nest_bound_l  .OR.&
-                                force_bound_l )                                &
+                            ( ( bc_dirichlet_l .OR. bc_radiation_l )           &
                               .AND. i == nxl   )    .OR.                       &
-                            ( ( inflow_r .OR. outflow_r .OR. nest_bound_r  .OR.&
-                                force_bound_r )                                &
+                            ( ( bc_dirichlet_r .OR. bc_radiation_r )           &
                               .AND. i == nxr   ) )                             &
                   THEN
@@ -803 +784 @@
                               .NOT. BTEST(wall_flags_0(k,j,i-1),2)             &
                                                     .OR.                       &
-                            ( ( inflow_r .OR. outflow_r .OR. nest_bound_r  .OR.&
-                                force_bound_r )                                &
+                            ( ( bc_dirichlet_r .OR. bc_radiation_r )           &
                               .AND. i == nxr-1 )    .OR.                       &
-                            ( ( inflow_l .OR. outflow_l .OR. nest_bound_l  .OR.&
-                                force_bound_l )                                &
+                            ( ( bc_dirichlet_l .OR. bc_radiation_l )           &
                               .AND. i == nxlu  ) )                             &
                    THEN
@@ -827 +806 @@
 !--                WS1 (21), WS3 (22), WS5 (23)
                    IF ( .NOT. BTEST(wall_flags_0(k,j+1,i),2) .OR.              &
-                            ( ( inflow_s .OR. outflow_s .OR. nest_bound_s  .OR.&
-                                force_bound_s )                                &
+                            ( ( bc_dirichlet_s .OR. bc_radiation_s )           &
                               .AND. j <= nysv  )    .OR.                       &
-                            ( ( inflow_n .OR. outflow_n .OR. nest_bound_n  .OR.&
-                                force_bound_n )                                &
+                            ( ( bc_dirichlet_n .OR. bc_radiation_n )           &
                               .AND. j == nyn   ) )                             &
                    THEN
@@ -839 +816 @@
                               .NOT. BTEST(wall_flags_0(k,j-1,i),2) )           &
                                                         .OR.                   &
-                            ( ( inflow_s .OR. outflow_s .OR. nest_bound_s  .OR.&
-                                force_bound_s )                                &
+                            ( (  bc_dirichlet_s .OR. bc_radiation_s )          &
                               .AND. j == nysv+1)    .OR.                       &
-                            ( ( inflow_n .OR. outflow_n .OR. nest_bound_n  .OR.&
-                                force_bound_n )                                &
+                            ( (  bc_dirichlet_n .OR. bc_radiation_n )          &
                               .AND. j == nyn-1 ) )                             &
                    THEN
@@ -924 +899 @@
 !--                WS1 (27), WS3 (28), WS5 (29)
                    IF ( .NOT. BTEST(wall_flags_0(k,j,i+1),3) .OR.              &
-                            ( ( inflow_l .OR. outflow_l .OR. nest_bound_l  .OR.&
-                                force_bound_l )                                &
+                            ( (  bc_dirichlet_l .OR. bc_radiation_l )          &
                               .AND. i == nxl   )    .OR.                       &
-                            ( ( inflow_r .OR. outflow_r .OR. nest_bound_r  .OR.&
-                                force_bound_r )                                &
+                            ( (  bc_dirichlet_r .OR. bc_radiation_r )          &
                               .AND. i == nxr   ) )                             &
                    THEN
@@ -936 +909 @@
                               .NOT. BTEST(wall_flags_0(k,j,i-1),3) )           &
                                                         .OR.                   &
-                            ( ( inflow_r .OR. outflow_r .OR. nest_bound_r  .OR.&
-                                force_bound_r )                                &
+                            ( ( bc_dirichlet_r .OR. bc_radiation_r )           &
                               .AND. i == nxr-1 )    .OR.                       &
-                            ( ( inflow_l .OR. outflow_l .OR. nest_bound_l  .OR.&
-                                force_bound_l )                                &
+                            ( ( bc_dirichlet_l .OR.  bc_radiation_l )          &
                               .AND. i == nxlu  ) )                             &
                    THEN
@@ -960 +931 @@
 !--                WS1 (30), WS3 (31), WS5 (32)
                    IF ( .NOT. BTEST(wall_flags_0(k,j+1,i),3) .OR.              &
-                            ( ( inflow_s .OR. outflow_s .OR. nest_bound_s  .OR.&
-                                force_bound_s )                                &
+                            ( ( bc_dirichlet_s .OR. bc_radiation_s )           &
                               .AND. j == nys   )    .OR.                       &
-                            ( ( inflow_n .OR. outflow_n .OR. nest_bound_n  .OR.&
-                                force_bound_n )                                &
+                            ( ( bc_dirichlet_n .OR. bc_radiation_n )           &
                               .AND. j == nyn   ) )                             &
                    THEN
@@ -972 +941 @@
                               .NOT. BTEST(wall_flags_0(k,j-1,i),3) )           &
                                                         .OR.                   &
-                            ( ( inflow_s .OR. outflow_s .OR. nest_bound_s  .OR.&
-                                force_bound_s )                                &
+                            ( ( bc_dirichlet_s .OR. bc_radiation_s )           &
                               .AND. j == nysv  )    .OR.                       &
-                            ( ( inflow_n .OR. outflow_n .OR. nest_bound_n  .OR.&
-                                force_bound_n )                                &
+                            ( ( bc_dirichlet_n .OR. bc_radiation_n )           &
                               .AND. j == nyn-1 ) )                             &
                    THEN
@@ -1067 +1034 @@
 !--       Set boundary flags at inflow and outflow boundary in case of 
 !--       non-cyclic boundary conditions.
-          IF ( inflow_l      .OR.  outflow_l  .OR.                             &
-               nest_bound_l  .OR.  force_bound_l )  THEN
+          IF ( bc_dirichlet_l  .OR.  bc_radiation_l )  THEN
              advc_flags_1(:,:,nxl-1) = advc_flags_1(:,:,nxl)
              advc_flags_2(:,:,nxl-1) = advc_flags_2(:,:,nxl)
           ENDIF
 
-          IF ( inflow_r      .OR.  outflow_r  .OR.                             &
-               nest_bound_r  .OR.  force_bound_r )  THEN
+          IF ( bc_dirichlet_r  .OR.  bc_radiation_r )  THEN
             advc_flags_1(:,:,nxr+1) = advc_flags_1(:,:,nxr)
             advc_flags_2(:,:,nxr+1) = advc_flags_2(:,:,nxr)
           ENDIF
 
-          IF ( inflow_n      .OR.  outflow_n  .OR.                             &
-               nest_bound_n  .OR.  force_bound_n )  THEN
+          IF ( bc_dirichlet_n  .OR.  bc_radiation_n )  THEN
              advc_flags_1(:,nyn+1,:) = advc_flags_1(:,nyn,:)
              advc_flags_2(:,nyn+1,:) = advc_flags_2(:,nyn,:)
           ENDIF
 
-          IF ( inflow_s      .OR.  outflow_s  .OR.                             &
-               nest_bound_s  .OR.  force_bound_s )  THEN
+          IF ( bc_dirichlet_s  .OR.  bc_radiation_s )  THEN
              advc_flags_1(:,nys-1,:) = advc_flags_1(:,nys,:)
              advc_flags_2(:,nys-1,:) = advc_flags_2(:,nys,:)

palm/trunk/SOURCE/boundary_conds.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! Rename some variables concerning LES-LES as well as offline nesting
 ! 
 ! Former revisions:
@@ -186 +186 @@
     USE arrays_3d,                                                             &
         ONLY:  c_u, c_u_m, c_u_m_l, c_v, c_v_m, c_v_m_l, c_w, c_w_m, c_w_m_l,  &
-               diss, diss_p, dzu, e_p, nc_p, nr_p, pt, pt_p, q, q_p, qc_p, qr_p, s,  & 
+               diss, diss_p, dzu, e_p, nc_p, nr_p, pt, pt_p, q, q_p, qc_p,     &
+               qr_p, s,                                                        & 
                s_p, sa, sa_p, u, ug, u_init, u_m_l, u_m_n, u_m_r, u_m_s, u_p,  &
                v, vg, v_init, v_m_l, v_m_n, v_m_r, v_m_s, v_p,                 &
@@ -195 +196 @@
               
     USE control_parameters,                                                    &
-        ONLY:  air_chemistry, bc_pt_t_val, bc_q_t_val, bc_s_t_val,             &
-               constant_diffusion, cloud_physics, coupling_mode, dt_3d,        &
-               force_bound_l, force_bound_s, forcing, humidity,                &
-               ibc_pt_b, ibc_pt_t, ibc_q_b, ibc_q_t, ibc_s_b, ibc_s_t,         &
-               ibc_sa_t, ibc_uv_b, ibc_uv_t, inflow_l, inflow_n, inflow_r,     &
-               inflow_s, intermediate_timestep_count, kappa,                   &
-               microphysics_morrison, microphysics_seifert, nest_domain,       &
-               nest_bound_l, nest_bound_n, nest_bound_r, nest_bound_s, nudging,&
-               ocean, outflow_l, outflow_n, outflow_r, outflow_s,              &
-               passive_scalar, rans_mode, rans_tke_e, tsc, use_cmax
+        ONLY:  air_chemistry, bc_dirichlet_l, bc_dirichlet_n, bc_dirichlet_r,  &
+               bc_dirichlet_s, bc_radiation_l, bc_radiation_n, bc_radiation_r, &
+               bc_radiation_s, bc_pt_t_val, bc_q_t_val, bc_s_t_val,            &
+               child_domain, constant_diffusion, cloud_physics, coupling_mode, &
+               dt_3d, humidity, ibc_pt_b, ibc_pt_t, ibc_q_b, ibc_q_t, ibc_s_b, &
+               ibc_s_t,ibc_sa_t, ibc_uv_b, ibc_uv_t,                           &
+               intermediate_timestep_count, kappa,                             &
+               microphysics_morrison, microphysics_seifert,                    &
+               nesting_offline, nudging,                                       &
+               ocean, passive_scalar, rans_mode, rans_tke_e, tsc, use_cmax
 
     USE grid_variables,                                                        &
@@ -275 +276 @@
 !
 !-- Vertical nesting: Vertical velocity not zero at the top of the fine grid
-    IF (  .NOT.  nest_domain  .AND.                                            &
+    IF (  .NOT.  child_domain  .AND.                                           &
                  TRIM(coupling_mode) /= 'vnested_fine' )  THEN
        w_p(nzt:nzt+1,:,:) = 0.0_wp  !< nzt is not a prognostic level (but cf. pres)
@@ -409 +410 @@
        ENDIF
 
-       IF ( .NOT. nest_domain )  THEN
+       IF ( .NOT. child_domain )  THEN
           e_p(nzt+1,:,:) = e_p(nzt,:,:)
 !
@@ -421 +422 @@
           IF ( rans_mode_parent  .AND.  .NOT. rans_mode )  THEN
 
-
-
              e_p(nzt+1,:,:) = e_p(nzt,:,:)
-             IF ( nest_bound_l )  e_p(:,:,nxl-1) = e_p(:,:,nxl)
-             IF ( nest_bound_r )  e_p(:,:,nxr+1) = e_p(:,:,nxr)
-             IF ( nest_bound_s )  e_p(:,nys-1,:) = e_p(:,nys,:)
-             IF ( nest_bound_n )  e_p(:,nyn+1,:) = e_p(:,nyn,:)
+             IF ( bc_dirichlet_l )  e_p(:,:,nxl-1) = e_p(:,:,nxl)
+             IF ( bc_dirichlet_r )  e_p(:,:,nxr+1) = e_p(:,:,nxr)
+             IF ( bc_dirichlet_s )  e_p(:,nys-1,:) = e_p(:,nys,:)
+             IF ( bc_dirichlet_n )  e_p(:,nyn+1,:) = e_p(:,nyn,:)
 
           ENDIF
@@ -510 +509 @@
        ENDDO
 
-       IF ( .NOT. nest_domain )  THEN
+       IF ( .NOT. child_domain )  THEN
           diss_p(nzt+1,:,:) = diss_p(nzt,:,:)
        ENDIF
@@ -697 +696 @@
 !-- have to be restored here.
 !-- For the SGS-TKE, Neumann boundary conditions are used at the inflow.
-    IF ( inflow_s )  THEN
+    IF ( bc_dirichlet_s )  THEN
        v_p(:,nys,:) = v_p(:,nys-1,:)
        IF ( .NOT. constant_diffusion ) e_p(:,nys-1,:) = e_p(:,nys,:)
-    ELSEIF ( inflow_n )  THEN
+    ELSEIF ( bc_dirichlet_n )  THEN
        IF ( .NOT. constant_diffusion ) e_p(:,nyn+1,:) = e_p(:,nyn,:)
-    ELSEIF ( inflow_l ) THEN
+    ELSEIF ( bc_dirichlet_l ) THEN
        u_p(:,:,nxl) = u_p(:,:,nxl-1)
        IF ( .NOT. constant_diffusion ) e_p(:,:,nxl-1) = e_p(:,:,nxl)
-    ELSEIF ( inflow_r )  THEN
+    ELSEIF ( bc_dirichlet_r )  THEN
        IF ( .NOT. constant_diffusion ) e_p(:,:,nxr+1) = e_p(:,:,nxr)
     ENDIF
@@ -712 +711 @@
 !-- The same restoration for u at i=nxl and v at j=nys as above must be made
 !-- in case of nest boundaries. This must not be done in case of vertical nesting 
-!-- mode as in that case the lateral boundaries are actually cyclic.
-    IF ( nesting_mode /= 'vertical'  .OR.  forcing )  THEN
-       IF ( nest_bound_s  .OR.  force_bound_s )  THEN
+!-- mode as in that case the lateral boundaries are actually cyclic. 
+!-- @todo: Is this really needed? Boundary values will be overwritten in 
+!--        coupler or by Inifor data.
+    IF ( nesting_mode /= 'vertical'  .OR.  nesting_offline )  THEN
+       IF ( bc_dirichlet_s )  THEN
           v_p(:,nys,:) = v_p(:,nys-1,:)
        ENDIF
-       IF ( nest_bound_l  .OR.  force_bound_l )  THEN
+       IF ( bc_dirichlet_l )  THEN
           u_p(:,:,nxl) = u_p(:,:,nxl-1)
        ENDIF
@@ -724 +725 @@
 !
 !-- Lateral boundary conditions for scalar quantities at the outflow
-    IF ( outflow_s )  THEN
+    IF ( bc_radiation_s )  THEN
        pt_p(:,nys-1,:)     = pt_p(:,nys,:)
        IF ( .NOT. constant_diffusion )  e_p(:,nys-1,:) = e_p(:,nys,:)
@@ -740 +741 @@
        ENDIF
        IF ( passive_scalar )  s_p(:,nys-1,:) = s_p(:,nys,:)
-    ELSEIF ( outflow_n )  THEN
+    ELSEIF ( bc_radiation_n )  THEN
        pt_p(:,nyn+1,:)     = pt_p(:,nyn,:)
        IF ( .NOT. constant_diffusion )  e_p(:,nyn+1,:) = e_p(:,nyn,:)
@@ -756 +757 @@
        ENDIF
        IF ( passive_scalar )  s_p(:,nyn+1,:) = s_p(:,nyn,:)
-    ELSEIF ( outflow_l )  THEN
+    ELSEIF ( bc_radiation_l )  THEN
        pt_p(:,:,nxl-1)     = pt_p(:,:,nxl)
        IF ( .NOT. constant_diffusion )  e_p(:,:,nxl-1) = e_p(:,:,nxl)
@@ -772 +773 @@
        ENDIF
        IF ( passive_scalar )  s_p(:,:,nxl-1) = s_p(:,:,nxl)
-    ELSEIF ( outflow_r )  THEN
+    ELSEIF ( bc_radiation_r )  THEN
        pt_p(:,:,nxr+1)     = pt_p(:,:,nxr)
        IF ( .NOT. constant_diffusion )  e_p(:,:,nxr+1) = e_p(:,:,nxr)
@@ -799 +800 @@
 !-- ensures numerical stability (CFL-condition) or calculated after
 !-- Orlanski(1976) and averaged along the outflow boundary.
-    IF ( outflow_s )  THEN
+    IF ( bc_radiation_s )  THEN
 
        IF ( use_cmax )  THEN
@@ -938 +939 @@
     ENDIF
 
-    IF ( outflow_n )  THEN
+    IF ( bc_radiation_n )  THEN
 
        IF ( use_cmax )  THEN
@@ -1077 +1078 @@
     ENDIF
 
-    IF ( outflow_l )  THEN
+    IF ( bc_radiation_l )  THEN
 
        IF ( use_cmax )  THEN
@@ -1216 +1217 @@
     ENDIF
 
-    IF ( outflow_r )  THEN
+    IF ( bc_radiation_r )  THEN
 
        IF ( use_cmax )  THEN

palm/trunk/SOURCE/check_parameters.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! Rename boundary conditions in offline nesting
 ! 
 ! Former revisions:
@@ -1069 +1069 @@
 !-- Check turbulence closure setup
     CALL tcm_check_parameters
+
 !
 !-- Check approximation
@@ -1437 +1438 @@
        CALL MPI_ALLREDUCE( dt_spinup,   dt_spinup_max,   1, MPI_REAL,          &
                            MPI_MAX, MPI_COMM_WORLD, ierr )
+
        IF ( spinup_time /= spinup_time_max  .OR.  dt_spinup /= dt_spinup_max ) &
        THEN
@@ -1781 +1783 @@
 !
 !-- Overwrite latitude if necessary and compute Coriolis parameter. 
-!-- To do - move initialization of f and fs to coriolis_mod.
+!-- @todo - move initialization of f and fs to coriolis_mod.
     IF ( input_pids_static )  THEN
        latitude  = init_model%latitude
@@ -1976 +1978 @@
 !-- TO_DO: later set bc_p_t to neumann before, in case of nested domain
     ELSEIF ( bc_p_t == 'neumann' .OR. bc_p_t == 'nested'  .OR.                 &
-             bc_p_t == 'forcing'  )  THEN
+             bc_p_t == 'nesting_offline'  )  THEN
        ibc_p_t = 1
     ELSE
@@ -2006 +2008 @@
     ELSEIF ( bc_pt_t == 'initial_gradient' )  THEN
        ibc_pt_t = 2
-    ELSEIF ( bc_pt_t == 'nested'  .OR.  bc_pt_t == 'forcing' )  THEN
+    ELSEIF ( bc_pt_t == 'nested'  .OR.  bc_pt_t == 'nesting_offline' )  THEN
        ibc_pt_t = 3
     ELSE
@@ -2234 +2236 @@
        ELSEIF ( bc_uv_t == 'neumann' )  THEN
           ibc_uv_t = 1
-       ELSEIF ( bc_uv_t == 'nested'  .OR.  bc_uv_t == 'forcing' )  THEN
+       ELSEIF ( bc_uv_t == 'nested'  .OR.  bc_uv_t == 'nesting_offline' )  THEN
           ibc_uv_t = 3
        ELSE
@@ -3926 +3928 @@
           dist_nxl    = MAX( inflow_disturbance_begin, nxl )
           dist_nxr(1) = MIN( inflow_disturbance_end, nxr )
-       ELSEIF ( bc_lr == 'nested'  .OR.  bc_lr == 'forcing' )  THEN
+       ELSEIF ( bc_lr == 'nested'  .OR.  bc_lr == 'nesting_offline' )  THEN
           dist_nxl    = MAX( inflow_disturbance_begin, nxl )
           dist_nxr    = MIN( nx - inflow_disturbance_begin, nxr )
@@ -3936 +3938 @@
           dist_nys    = MAX( inflow_disturbance_begin, nys )
           dist_nyn(1) = MIN( inflow_disturbance_end, nyn )
-       ELSEIF ( bc_ns == 'nested'  .OR.  bc_ns == 'forcing' )  THEN
+       ELSEIF ( bc_ns == 'nested'  .OR.  bc_ns == 'nesting_offline' )  THEN
           dist_nys    = MAX( inflow_disturbance_begin, nys )
           dist_nyn    = MIN( ny - inflow_disturbance_begin, nyn )
@@ -3949 +3951 @@
           dist_nxl    = inflow_disturbance_begin
           dist_nxr(1) = inflow_disturbance_end
-       ELSEIF ( bc_lr == 'nested'  .OR.  bc_lr == 'forcing' )  THEN
+       ELSEIF ( bc_lr == 'nested'  .OR.  bc_lr == 'nesting_offline' )  THEN
           dist_nxr    = nx - inflow_disturbance_begin
           dist_nxl    = inflow_disturbance_begin
@@ -3959 +3961 @@
           dist_nys    = inflow_disturbance_begin
           dist_nyn(1) = inflow_disturbance_end
-       ELSEIF ( bc_ns == 'nested'  .OR.  bc_ns == 'forcing' )  THEN
+       ELSEIF ( bc_ns == 'nested'  .OR.  bc_ns == 'nesting_offline' )  THEN
           dist_nyn    = ny - inflow_disturbance_begin
           dist_nys    = inflow_disturbance_begin
@@ -4403 +4405 @@
        ELSEIF ( bc_t == 'initial_gradient' )  THEN
           ibc_t = 2
-       ELSEIF ( bc_t == 'nested'  .OR.  bc_t == 'forcing' )  THEN
+       ELSEIF ( bc_t == 'nested'  .OR.  bc_t == 'nesting_offline' )  THEN
           ibc_t = 3
        ELSE

palm/trunk/SOURCE/chemistry_model_mod.f90

@@ -22 +22 @@
 ! Current revisions:
 ! -----------------
-! 
+! Rename flags indicating outflow boundary conditions
 ! 
 ! Former revisions:
@@ -291 +291 @@
                                                                                   
     USE control_parameters,                                                    &  
-        ONLY: air_chemistry, outflow_l, outflow_n, outflow_r, outflow_s                 
+        ONLY:  air_chemistry, bc_radiation_l, bc_radiation_n, bc_radiation_r,  &
+               bc_radiation_s                
     USE indices,                                                               &  
-        ONLY: nxl, nxr,  nxlg, nxrg, nyng, nysg, nzt                              
+        ONLY:  nxl, nxr,  nxlg, nxrg, nyng, nysg, nzt                              
                                                                                   
 !    USE prognostic_equations_mod,                                             &
 
     USE arrays_3d,                                                             &     
-        ONLY: dzu                                                
+        ONLY:  dzu                                                
     USE surface_mod,                                                           &
-        ONLY: bc_h                                                           
+        ONLY:  bc_h                                                           
 
     CHARACTER (len=*), INTENT(IN) :: mode
@@ -420 +421 @@
 !--          Lateral boundary conditions for chem species at outflow boundary 
 
-          IF ( outflow_s )  THEN
+          IF ( bc_radiation_s )  THEN
              DO lsp = 1, nspec
                 chem_species(lsp)%conc_p(:,nys-1,:) = chem_species(lsp)%conc_p(:,nys,:)
              ENDDO
-          ELSEIF ( outflow_n )  THEN
+          ELSEIF ( bc_radiation_n )  THEN
              DO lsp = 1, nspec
                 chem_species(lsp)%conc_p(:,nyn+1,:) = chem_species(lsp)%conc_p(:,nyn,:)
              ENDDO
-          ELSEIF ( outflow_l )  THEN
+          ELSEIF ( bc_radiation_l )  THEN
              DO lsp = 1, nspec
                 chem_species(lsp)%conc_p(:,:,nxl-1) = chem_species(lsp)%conc_p(:,:,nxl)
              ENDDO
-          ELSEIF ( outflow_r )  THEN
+          ELSEIF ( bc_radiation_r )  THEN
              DO lsp = 1, nspec
                 chem_species(lsp)%conc_p(:,:,nxr+1) = chem_species(lsp)%conc_p(:,:,nxr)

palm/trunk/SOURCE/coriolis.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! Remove masking of geostrophic wind forcing in offline nesting case
 ! 
 ! Former revisions:
@@ -109 +109 @@
            
        USE control_parameters,                                                 &
-           ONLY:  f, forcing, fs, message_string
+           ONLY:  f, fs, message_string, nesting_offline
            
        USE indices,                                                            &
@@ -123 +123 @@
        INTEGER(iwp) ::  k          !< running index z direction 
 
-       REAL(wp)     ::  flag       !< flag to mask topography
-       REAL(wp)     ::  flag_force !< flag to mask large-scale pressure gradient in case larger-scale forcing is applied
-
-       flag_force = MERGE( 0.0_wp, 1.0_wp, forcing )
+       REAL(wp)     ::  flag           !< flag to mask topography
+
 !
 !--    Compute Coriolis terms for the three velocity components
@@ -144 +142 @@
                       tend(k,j,i) = tend(k,j,i) + f  *    ( 0.25_wp *          &
                                    ( v(k,j,i-1) + v(k,j,i) + v(k,j+1,i-1) +    &
-                                     v(k,j+1,i) ) - vg(k) * flag_force         &
-                                                          ) * flag           &
+                                     v(k,j+1,i) ) - vg(k) ) * flag             &
                                                 - fs *    ( 0.25_wp *          &
                                    ( w(k-1,j,i-1) + w(k-1,j,i) + w(k,j,i-1) +  &
@@ -167 +164 @@
                       tend(k,j,i) = tend(k,j,i) - f *     ( 0.25_wp *          &
                                    ( u(k,j-1,i) + u(k,j,i) + u(k,j-1,i+1) +    &
-                                     u(k,j,i+1) ) - ug(k) * flag_force         &
-                                                          ) * flag
+                                     u(k,j,i+1) ) - ug(k) ) * flag
                    ENDDO
                 ENDDO
@@ -212 +208 @@
            
        USE control_parameters,                                                 &
-           ONLY:  f, forcing, fs, message_string
+           ONLY:  f, fs, message_string, nesting_offline
            
        USE indices,                                                            &
@@ -227 +223 @@
 
        REAL(wp)     ::  flag       !< flag to mask topography
-       REAL(wp)     ::  flag_force !< flag to mask large-scale pressure gradient in case larger-scale forcing is applied
-
-       flag_force = MERGE( 0.0_wp, 1.0_wp, forcing )
+
 !
 !--    Compute Coriolis terms for the three velocity components
@@ -244 +238 @@
                 tend(k,j,i) = tend(k,j,i) + f  *     ( 0.25_wp *               &
                                 ( v(k,j,i-1) + v(k,j,i) + v(k,j+1,i-1) +       &
-                                  v(k,j+1,i) ) - vg(k) * flag_force            &
+                                  v(k,j+1,i) ) - vg(k)                         &
                                                      ) * flag                  &
                                           - fs *     ( 0.25_wp *               &
@@ -261 +255 @@
                 tend(k,j,i) = tend(k,j,i) - f *        ( 0.25_wp *             &
                                 ( u(k,j-1,i) + u(k,j,i) + u(k,j-1,i+1) +       &
-                                  u(k,j,i+1) ) - ug(k)   * flag_force          &
-                                                       ) * flag
+                                  u(k,j,i+1) ) - ug(k) ) * flag
              ENDDO
 

palm/trunk/SOURCE/exchange_horiz_2d.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! Rename variables in offline nesting mode and flags indicating lateral
+! boundary conditions
 ! 
 ! Former revisions:
@@ -84 +85 @@
 
     USE control_parameters,                                                    &
-        ONLY :  bc_lr_cyc, bc_ns_cyc, force_bound_l, force_bound_n,            &
-                force_bound_r, force_bound_s,                                  &
-                inflow_l, inflow_n, inflow_r, inflow_s,                        &
-                nest_bound_l, nest_bound_n, nest_bound_r, nest_bound_s,        &
-                outflow_l, outflow_n, outflow_r, outflow_s
+        ONLY :  bc_dirichlet_l, bc_dirichlet_n, bc_dirichlet_r,                &
+                bc_dirichlet_s, bc_lr_cyc, bc_ns_cyc, bc_radiation_l,          &
+                bc_radiation_n, bc_radiation_r, bc_radiation_s 
                 
     USE cpulog,                                                                &
@@ -184 +183 @@
 !-- Neumann-conditions at inflow/outflow/nested boundaries
     IF ( nesting_mode /= 'vertical' )  THEN
-       IF ( inflow_l .OR. outflow_l .OR. nest_bound_l  .OR.  force_bound_l )   &
-       THEN
+       IF ( bc_dirichlet_l  .OR.  bc_radiation_l )  THEN
           DO  i = nbgp, 1, -1
              ar(:,nxl-i) = ar(:,nxl)
           ENDDO
        ENDIF
-       IF ( inflow_r .OR. outflow_r .OR. nest_bound_r  .OR.  force_bound_r )   &
-       THEN
+       IF ( bc_dirichlet_r  .OR.  bc_radiation_r )  THEN
           DO  i = 1, nbgp
              ar(:,nxr+i) = ar(:,nxr)
           ENDDO
        ENDIF
-       IF ( inflow_s .OR. outflow_s .OR. nest_bound_s  .OR.  force_bound_s )   &
-       THEN
+       IF ( bc_dirichlet_s  .OR.  bc_radiation_s )  THEN
           DO  i = nbgp, 1, -1
              ar(nys-i,:) = ar(nys,:)
           ENDDO
        ENDIF
-       IF ( inflow_n .OR. outflow_n .OR. nest_bound_n  .OR.  force_bound_n )   &
-       THEN
+       IF ( bc_dirichlet_n  .OR.  bc_radiation_n )  THEN
           DO  i = 1, nbgp
              ar(nyn+i,:) = ar(nyn,:)
@@ -227 +222 @@
 
     USE control_parameters,                                                    &
-        ONLY:  bc_lr_cyc, bc_ns_cyc, grid_level, force_bound_l, force_bound_n, &
-               force_bound_r, force_bound_s, nest_bound_l, nest_bound_n,       &
-               nest_bound_r, nest_bound_s
+        ONLY:  bc_dirichlet_l, bc_dirichlet_n, bc_dirichlet_r,                 &
+               bc_dirichlet_s, bc_lr_cyc, bc_ns_cyc, bc_radiation_l,           &
+               bc_radiation_n, bc_radiation_r, bc_radiation_s, grid_level
         
     USE cpulog,                                                                &
@@ -326 +321 @@
 !
 !-- Neumann-conditions at inflow/outflow/nested boundaries
-    IF ( nest_bound_l  .OR.  force_bound_l )  THEN
+    IF ( bc_dirichlet_l )  THEN
        DO  i = nbgp_local, 1, -1
          ar(:,nxl_l-i) = ar(:,nxl_l)
        ENDDO
     ENDIF
-    IF ( nest_bound_r  .OR.  force_bound_r )  THEN
+    IF ( bc_dirichlet_r )  THEN
        DO  i = 1, nbgp_local
           ar(:,nxr_l+i) = ar(:,nxr_l)
        ENDDO
     ENDIF
-    IF ( nest_bound_s  .OR.  force_bound_s )  THEN
+    IF ( bc_dirichlet_s )  THEN
        DO  i = nbgp_local, 1, -1
          ar(nys_l-i,:) = ar(nys_l,:)
        ENDDO
     ENDIF
-    IF ( nest_bound_n  .OR.  force_bound_n )  THEN
+    IF ( bc_dirichlet_n )  THEN
        DO  i = 1, nbgp_local
          ar(nyn_l+i,:) = ar(nyn_l,:)

palm/trunk/SOURCE/init_3d_model.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-! 
+! Revise Inifor initialization
 ! 
 ! Former revisions:
@@ -535 +535 @@
 
     USE netcdf_data_input_mod,                                                 &
-        ONLY:  init_3d, netcdf_data_input_interpolate, netcdf_data_input_init_3d
+        ONLY:  init_3d, netcdf_data_input_init_3d
     
     USE particle_attributes,                                                   &
@@ -611 +611 @@
     REAL(wp)     ::  t_surface !< air temperature at the surface
 
+    REAL(wp), DIMENSION(:), ALLOCATABLE ::  init_l        !< dummy array used for averaging 3D data to obtain inital profiles 
     REAL(wp), DIMENSION(:), ALLOCATABLE ::  p_hydrostatic !< hydrostatic pressure
 
@@ -1006 +1007 @@
 !-- Arrays to store velocity data from t-dt and the phase speeds which
 !-- are needed for radiation boundary conditions
-    IF ( outflow_l )  THEN
+    IF ( bc_radiation_l )  THEN
        ALLOCATE( u_m_l(nzb:nzt+1,nysg:nyng,1:2),                               &
                  v_m_l(nzb:nzt+1,nysg:nyng,0:1),                               &
                  w_m_l(nzb:nzt+1,nysg:nyng,0:1) )
     ENDIF
-    IF ( outflow_r )  THEN
+    IF ( bc_radiation_r )  THEN
        ALLOCATE( u_m_r(nzb:nzt+1,nysg:nyng,nx-1:nx),                           &
                  v_m_r(nzb:nzt+1,nysg:nyng,nx-1:nx),                           &
                  w_m_r(nzb:nzt+1,nysg:nyng,nx-1:nx) )
     ENDIF
-    IF ( outflow_l  .OR.  outflow_r )  THEN
+    IF ( bc_radiation_l  .OR.  bc_radiation_r )  THEN
        ALLOCATE( c_u(nzb:nzt+1,nysg:nyng), c_v(nzb:nzt+1,nysg:nyng),           &
                  c_w(nzb:nzt+1,nysg:nyng) )
     ENDIF
-    IF ( outflow_s )  THEN
+    IF ( bc_radiation_s )  THEN
        ALLOCATE( u_m_s(nzb:nzt+1,0:1,nxlg:nxrg),                               &
                  v_m_s(nzb:nzt+1,1:2,nxlg:nxrg),                               &
                  w_m_s(nzb:nzt+1,0:1,nxlg:nxrg) )
     ENDIF
-    IF ( outflow_n )  THEN
+    IF ( bc_radiation_n )  THEN
        ALLOCATE( u_m_n(nzb:nzt+1,ny-1:ny,nxlg:nxrg),                           &
                  v_m_n(nzb:nzt+1,ny-1:ny,nxlg:nxrg),                           &
                  w_m_n(nzb:nzt+1,ny-1:ny,nxlg:nxrg) )
     ENDIF
-    IF ( outflow_s  .OR.  outflow_n )  THEN
+    IF ( bc_radiation_s  .OR.  bc_radiation_n )  THEN
        ALLOCATE( c_u(nzb:nzt+1,nxlg:nxrg), c_v(nzb:nzt+1,nxlg:nxrg),           &
                  c_w(nzb:nzt+1,nxlg:nxrg) )
     ENDIF
-    IF ( outflow_l  .OR.  outflow_r  .OR.  outflow_s  .OR.  outflow_n )  THEN
+    IF ( bc_radiation_l  .OR.  bc_radiation_r  .OR.  bc_radiation_s  .OR.      &
+         bc_radiation_n )  THEN
        ALLOCATE( c_u_m_l(nzb:nzt+1), c_v_m_l(nzb:nzt+1), c_w_m_l(nzb:nzt+1) )                   
        ALLOCATE( c_u_m(nzb:nzt+1), c_v_m(nzb:nzt+1), c_w_m(nzb:nzt+1) )
@@ -1152 +1154 @@
     sums_l_l           = 0.0_wp
     sums_wsts_bc_l     = 0.0_wp
-
-
-
+    
 !
 !-- Initialize model variables
@@ -1169 +1169 @@
           CALL netcdf_data_input_init_3d
 !
-!--       Please note, at the moment INIFOR assumes only an equidistant vertical
-!--       grid. In case of vertical grid stretching, input of inital data 
-!--       need to be inter- and/or extrapolated. 
-!--       Therefore, check if zu grid on file is identical to numeric zw grid. 
-!--       Please note  
-          IF ( ANY( zu(1:nzt+1) /= init_3d%zu_atmos(1:init_3d%nzu) ) )  THEN
-
-             IF( init_3d%lod_u == 1 )                                          &
-                CALL netcdf_data_input_interpolate(                            &
-                                                 init_3d%u_init(nzb+1:nzt+1),  &
-                                                 zu(nzb+1:nzt+1),              &
-                                                 init_3d%zu_atmos )
-                                                 
-             IF( init_3d%lod_v == 1 )                                          &
-                CALL netcdf_data_input_interpolate(                            &
-                                                 init_3d%v_init(nzb+1:nzt+1),  &
-                                                 zu(nzb+1:nzt+1),              &
-                                                 init_3d%zu_atmos )
-                                                 
-!              CALL netcdf_data_input_interpolate( init_3d%w_init(nzb+1:nzt),    &
-!                                                  zw(nzb+1:nzt),                &
-!                                                  init_3d%zw_atmos )
-
-             IF ( .NOT. neutral  .AND.  init_3d%lod_pt == 1 )                  &
-                CALL netcdf_data_input_interpolate(                            &
-                                             init_3d%pt_init(nzb+1:nzt+1),     &
-                                             zu(nzb+1:nzt+1),                  &
-                                             init_3d%zu_atmos )
-                                             
-             IF ( humidity  .AND.  init_3d%lod_q == 1 )                        &
-                CALL netcdf_data_input_interpolate(                            &
-                                             init_3d%q_init(nzb+1:nzt+1),      &
-                                             zu(nzb+1:nzt+1),                  &
-                                             init_3d%zu_atmos )
-          ENDIF
-!
-!--       In case of LOD=1, initialize 1D profiles and 3D data.  
-          IF( init_3d%lod_u == 1 )  u_init = init_3d%u_init
-          IF( init_3d%lod_v == 1 )  v_init = init_3d%v_init    
-          IF( .NOT. neutral  .AND.  init_3d%lod_pt == 1 )                      &
-             pt_init = init_3d%pt_init
-          IF( humidity  .AND.  init_3d%lod_q == 1 )                            &
-             q_init  = init_3d%q_init
-
-!
-!--       Please note, Inifor provides data from nzb+1 to nzt+1. 
-!--       Initialize pt and q with Neumann condition at nzb.
-          IF( .NOT. neutral )  pt_init(nzb) = pt_init(nzb+1)
-          IF( humidity      )  q_init(nzb)  = q_init(nzb+1)
+!--       Please note, Inifor provides data from nzb+1 to nzt. 
+!--       Bottom and top boundary conditions for Inifor profiles are already 
+!--       set (just after reading), so that this is not necessary here.
+!--       Depending on the provided level-of-detail, initial Inifor data is
+!--       either stored on data type (lod=1), or directly on 3D arrays (lod=2). 
+!--       In order to obtain also initial profiles in case of lod=2 (which 
+!--       is required for e.g. damping), average over 3D data.
+          IF( init_3d%lod_u == 1 )  THEN
+             u_init = init_3d%u_init
+          ELSEIF( init_3d%lod_u == 2 )  THEN 
+             ALLOCATE( init_l(nzb:nzt+1) ) 
+             DO  k = nzb, nzt+1
+                init_l(k) = SUM( u(k,nys:nyn,nxl:nxr) )
+             ENDDO
+             init_l = init_l / REAL( ( nx + 1 ) * ( ny + 1 ), KIND = wp )
+
+#if defined( __parallel )
+             CALL MPI_ALLREDUCE( init_l, u_init, nzt+1-nzb+1,                  &
+                                 MPI_REAL, MPI_SUM, comm2d, ierr )
+#else
+             u_init = init_l
+#endif
+             DEALLOCATE( init_l )
+
+          ENDIF
+           
+          IF( init_3d%lod_v == 1 )  THEN  
+             v_init = init_3d%v_init
+          ELSEIF( init_3d%lod_v == 2 )  THEN 
+             ALLOCATE( init_l(nzb:nzt+1) ) 
+             DO  k = nzb, nzt+1
+                init_l(k) = SUM( v(k,nys:nyn,nxl:nxr) )
+             ENDDO
+             init_l = init_l / REAL( ( nx + 1 ) * ( ny + 1 ), KIND = wp )
+
+#if defined( __parallel )
+             CALL MPI_ALLREDUCE( init_l, v_init, nzt+1-nzb+1,                  &
+                                 MPI_REAL, MPI_SUM, comm2d, ierr )
+#else
+             v_init = init_l
+#endif
+             DEALLOCATE( init_l )
+          ENDIF
+          IF( .NOT. neutral )  THEN
+             IF( init_3d%lod_pt == 1 )  THEN
+                pt_init = init_3d%pt_init
+             ELSEIF( init_3d%lod_pt == 2 )  THEN 
+                ALLOCATE( init_l(nzb:nzt+1) ) 
+                DO  k = nzb, nzt+1
+                   init_l(k) = SUM( pt(k,nys:nyn,nxl:nxr) )
+                ENDDO
+                init_l = init_l / REAL( ( nx + 1 ) * ( ny + 1 ), KIND = wp )
+
+#if defined( __parallel )
+                CALL MPI_ALLREDUCE( init_l, pt_init, nzt+1-nzb+1,               &
+                                    MPI_REAL, MPI_SUM, comm2d, ierr )
+#else
+                pt_init = init_l
+#endif
+                DEALLOCATE( init_l )
+             ENDIF
+          ENDIF
+
+
+          IF( humidity )  THEN
+             IF( init_3d%lod_q == 1 )  THEN
+                q_init = init_3d%q_init
+             ELSEIF( init_3d%lod_q == 2 )  THEN 
+                ALLOCATE( init_l(nzb:nzt+1) ) 
+                DO  k = nzb, nzt+1
+                   init_l(k) = SUM( q(k,nys:nyn,nxl:nxr) )
+                ENDDO
+                init_l = init_l / REAL( ( nx + 1 ) * ( ny + 1 ), KIND = wp )
+
+#if defined( __parallel )
+                CALL MPI_ALLREDUCE( init_l, q_init, nzt+1-nzb+1,               &
+                                    MPI_REAL, MPI_SUM, comm2d, ierr )
+#else
+                q_init = init_l
+#endif
+                DEALLOCATE( init_l )
+             ENDIF
+          ENDIF
+
+!
+!--       Write initial profiles onto 3D arrays. Note, only in case of lod = 1, 
+!--       for lod = 2 data is already on 3D arrays.   
           DO  i = nxlg, nxrg
              DO  j = nysg, nyng
@@ -1224 +1262 @@
                 IF( .NOT. neutral  .AND.  init_3d%lod_pt == 1 )                &
                    pt(:,j,i) = pt_init(:)
-                IF( humidity  .AND.  init_3d%lod_q == 1 )                      &
-                   q(:,j,i)  = q_init(:)
+                IF( humidity  .AND.  init_3d%lod_q == 1 )  q(:,j,i) = q_init(:)
              ENDDO
           ENDDO
 !
-!--       MS: What about the geostrophic wind profiles? Actually these
-!--           are not identical to the initial wind profiles in this case. 
-!--           This need to be further revised.
+!--       Exchange ghost points in case of level-of-detail = 2
+          IF( init_3d%lod_u == 2 )   CALL exchange_horiz( u, nbgp )
+          IF( init_3d%lod_v == 2 )   CALL exchange_horiz( v, nbgp )
+          IF( init_3d%lod_w == 2 )   CALL exchange_horiz( w, nbgp )
+          IF( .NOT. neutral  .AND.  init_3d%lod_pt == 2 )                      &
+             CALL exchange_horiz( pt, nbgp )
+          IF( humidity  .AND.  init_3d%lod_q == 2 )                            &
+             CALL exchange_horiz( q, nbgp )
+!
+!--       Set geostrophic wind components.  
           IF ( init_3d%from_file_ug )  THEN
              ug(:) = init_3d%ug_init(:)
@@ -1238 +1282 @@
              vg(:) = init_3d%vg_init(:)
           ENDIF
-
+          
           ug(nzt+1) = ug(nzt)
           vg(nzt+1) = vg(nzt)
-
 !
 !--       Set inital w to 0
@@ -1815 +1858 @@
 !--       Use these mean profiles at the inflow (provided that Dirichlet
 !--       conditions are used)
-          IF ( inflow_l )  THEN
+          IF ( bc_dirichlet_l )  THEN
              DO  j = nysg, nyng
                 DO  k = nzb, nzt+1
@@ -1963 +2006 @@
 !
 !--    Initialize old timelevels needed for radiation boundary conditions
-       IF ( outflow_l )  THEN
+       IF ( bc_radiation_l )  THEN
           u_m_l(:,:,:) = u(:,:,1:2)
           v_m_l(:,:,:) = v(:,:,0:1)
           w_m_l(:,:,:) = w(:,:,0:1)
        ENDIF
-       IF ( outflow_r )  THEN
+       IF ( bc_radiation_r )  THEN
           u_m_r(:,:,:) = u(:,:,nx-1:nx)
           v_m_r(:,:,:) = v(:,:,nx-1:nx)
           w_m_r(:,:,:) = w(:,:,nx-1:nx)
        ENDIF
-       IF ( outflow_s )  THEN
+       IF ( bc_radiation_s )  THEN
           u_m_s(:,:,:) = u(:,0:1,:)
           v_m_s(:,:,:) = v(:,1:2,:)
           w_m_s(:,:,:) = w(:,0:1,:)
        ENDIF
-       IF ( outflow_n )  THEN
+       IF ( bc_radiation_n )  THEN
           u_m_n(:,:,:) = u(:,ny-1:ny,:)
           v_m_n(:,:,:) = v(:,ny-1:ny,:)
@@ -2301 +2344 @@
 !-- Initialize nudging if required
     IF ( nudging )  CALL nudge_init
-
 !
 !-- Initialize 1D/3D offline-nesting with COSMO model and read data from 
 !-- external file.
-    IF ( large_scale_forcing  .OR.  forcing )  CALL lsf_init
-    
+    IF ( large_scale_forcing  .OR.  nesting_offline )  CALL lsf_init
 !
 !-- Initialize surface forcing corresponding to large-scale forcing. Therein, 

palm/trunk/SOURCE/init_grid.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! Rename variables in mesoscale-offline nesting mode
 ! 
 ! Former revisions:
@@ -342 +342 @@
         
     USE control_parameters,                                                    &
-        ONLY:  bc_lr_cyc, bc_ns_cyc, building_height, building_length_x,       &
+        ONLY:  bc_dirichlet_l, bc_dirichlet_n, bc_dirichlet_r,                 &
+               bc_dirichlet_s, bc_lr_cyc, bc_ns_cyc, bc_radiation_l,           &
+               bc_radiation_n, bc_radiation_r, bc_radiation_s,                 &
+               building_height, building_length_x,                             &
                building_length_y, building_wall_left, building_wall_south,     &
                canyon_height, canyon_wall_left, canyon_wall_south,             &
@@ -350 +353 @@
                dz_stretch_level_end_index, dz_stretch_level_start_index,       &
                dz_stretch_level_start, grid_level,                             &
-               force_bound_l, force_bound_r, force_bound_n, force_bound_s,     &
-               ibc_uv_b, inflow_l, inflow_n, inflow_r, inflow_s,               &
-               masking_method, maximum_grid_level, message_string,             &
-               momentum_advec, nest_domain, nest_bound_l,                      &
-               nest_bound_n, nest_bound_r, nest_bound_s,                       &
-               number_stretch_level_end, number_stretch_level_start, ocean,    &
-               outflow_l, outflow_n, outflow_r, outflow_s, psolver,            & 
-               scalar_advec, topography, topography_grid_convention,           &
-               tunnel_height, tunnel_length, tunnel_width_x, tunnel_width_y,   &
+               ibc_uv_b, masking_method, maximum_grid_level, message_string,   &
+               momentum_advec, number_stretch_level_end,                       &
+               number_stretch_level_start,ocean, psolver, scalar_advec,        &
+               topography, topography_grid_convention, tunnel_height,          &
+               tunnel_length, tunnel_width_x, tunnel_width_y,                  &
                tunnel_wall_depth, use_surface_fluxes, use_top_fluxes,          &
                wall_adjustment_factor
@@ -386 +385 @@
     IMPLICIT NONE
 
-    INTEGER(iwp) ::  i                           !< index variable along x 
-    INTEGER(iwp) ::  j                           !< index variable along y
-    INTEGER(iwp) ::  k                           !< index variable along z
-    INTEGER(iwp) ::  k_top                       !< topography top index on local PE
-    INTEGER(iwp) ::  n                           !< loop variable for stretching
-    INTEGER(iwp) ::  number_dz                   !< number of user-specified dz values       
-    INTEGER(iwp) ::  nzb_local_max               !< vertical grid index of maximum topography height
-    INTEGER(iwp) ::  nzb_local_min               !< vertical grid index of minimum topography height
+    INTEGER(iwp) ::  i             !< index variable along x 
+    INTEGER(iwp) ::  j             !< index variable along y
+    INTEGER(iwp) ::  k             !< index variable along z
+    INTEGER(iwp) ::  k_top         !< topography top index on local PE
+    INTEGER(iwp) ::  l             !< loop variable 
+    INTEGER(iwp) ::  n             !< loop variable for stretching
+    INTEGER(iwp) ::  number_dz     !< number of user-specified dz values       
+    INTEGER(iwp) ::  nzb_local_max !< vertical grid index of maximum topography height
+    INTEGER(iwp) ::  nzb_local_min !< vertical grid index of minimum topography height
                                       
     INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  nzb_local  !< index for topography top at cell-center
@@ -834 +834 @@
     nzb_max = k_top + 1
 #endif
-    IF ( inflow_l  .OR.  outflow_l  .OR.  force_bound_l  .OR.  nest_bound_l  .OR.&
-         inflow_r  .OR.  outflow_r  .OR.  force_bound_r  .OR.  nest_bound_r  .OR.&
-         inflow_n  .OR.  outflow_n  .OR.  force_bound_n  .OR.  nest_bound_n  .OR.&
-         inflow_s  .OR.  outflow_s  .OR.  force_bound_s  .OR.  nest_bound_s )    &
+    IF ( bc_dirichlet_l  .OR.  bc_radiation_l     .OR.                         &
+         bc_dirichlet_r  .OR.  bc_radiation_r     .OR.                         &
+         bc_dirichlet_n  .OR.  bc_radiation_n     .OR.                         &
+         bc_dirichlet_s  .OR.  bc_radiation_s )                                &
          nzb_max = nzt
 !   
 !-- Finally, if topography extents up to the model top, limit nzb_max to nzt.
-    nzb_max = MIN( nzb_max, nzt )
+    nzb_max = MIN( nzb_max, nzt ) 
 !
 !-- Determine minimum index of topography. Usually, this will be nzb. In case
@@ -857 +857 @@
 !-- Initialize boundary conditions via surface type
     CALL init_bc
+
 !
 !-- Allocate and set topography height arrays required for data output
@@ -926 +927 @@
 #if defined( __parallel )
        CALL MPI_ALLREDUCE( MAXVAL( get_topography_top_index( 's' ) ),          &
-                           nzb_local_max, 1, MPI_INTEGER, MPI_MAX, comm2d, ierr )             
+                           nzb_local_max, 1, MPI_INTEGER, MPI_MAX, comm2d, ierr )               
 #else
        nzb_local_max = MAXVAL( get_topography_top_index( 's' ) )

palm/trunk/SOURCE/init_pegrid.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-! 
+! Rename variables and boundary conditions in mesoscale-offline nesting mode
 ! 
 ! Former revisions:
@@ -239 +239 @@
 
     USE control_parameters,                                                    &
-        ONLY:  bc_lr, bc_ns, coupling_mode, coupling_mode_remote,              &
-               coupling_topology, force_bound_l, force_bound_n, force_bound_r, &
-               force_bound_s, gathered_size, grid_level,                       &
-               grid_level_count, inflow_l, inflow_n, inflow_r, inflow_s,       &
-               maximum_grid_level, message_string,                             &
-               mg_switch_to_pe0_level, momentum_advec, nest_bound_l,           &
-               nest_bound_n, nest_bound_r, nest_bound_s, nest_domain, neutral, &
-               psolver, outflow_l, outflow_n, outflow_r, outflow_s,            &
-               outflow_source_plane, recycling_width, scalar_advec,            &
+        ONLY:  bc_dirichlet_l, bc_dirichlet_n, bc_dirichlet_r, bc_dirichlet_s, &
+               bc_lr, bc_ns, bc_radiation_l, bc_radiation_n, bc_radiation_r,   &
+               bc_radiation_s, coupling_mode, coupling_mode_remote,            &
+               coupling_topology, gathered_size, grid_level,                   &
+               grid_level_count, maximum_grid_level, message_string,           &
+               mg_switch_to_pe0_level, momentum_advec, neutral,                &
+               psolver, outflow_source_plane, recycling_width, scalar_advec,   &
                subdomain_size, turbulent_outflow, y_shift
 
@@ -1232 +1230 @@
 !-- Setting of flags for inflow/outflow/nesting conditions.
     IF ( pleft == MPI_PROC_NULL )  THEN
-       IF ( bc_lr == 'dirichlet/radiation' )  THEN
-          inflow_l  = .TRUE.
+       IF ( bc_lr == 'dirichlet/radiation'  .OR.  bc_lr == 'nested'  .OR.      &
+            bc_lr == 'nesting_offline' )  THEN
+          bc_dirichlet_l  = .TRUE.
        ELSEIF ( bc_lr == 'radiation/dirichlet' )  THEN
-          outflow_l = .TRUE.
-       ELSEIF ( bc_lr == 'nested' )  THEN
-          nest_bound_l = .TRUE.
-       ELSEIF ( bc_lr == 'forcing' )  THEN
-          force_bound_l = .TRUE.
+          bc_radiation_l = .TRUE.
        ENDIF
     ENDIF
@@ -1245 +1240 @@
     IF ( pright == MPI_PROC_NULL )  THEN
        IF ( bc_lr == 'dirichlet/radiation' )  THEN
-          outflow_r = .TRUE.
-       ELSEIF ( bc_lr == 'radiation/dirichlet' )  THEN
-          inflow_r  = .TRUE.
-       ELSEIF ( bc_lr == 'nested' )  THEN
-          nest_bound_r = .TRUE.
-       ELSEIF ( bc_lr == 'forcing' )  THEN
-          force_bound_r = .TRUE.
+          bc_radiation_r = .TRUE.
+       ELSEIF ( bc_lr == 'radiation/dirichlet'  .OR.  bc_lr == 'nested'  .OR.  &
+                bc_lr == 'nesting_offline' )  THEN
+          bc_dirichlet_r  = .TRUE.
        ENDIF
     ENDIF
@@ -1257 +1249 @@
     IF ( psouth == MPI_PROC_NULL )  THEN
        IF ( bc_ns == 'dirichlet/radiation' )  THEN
-          outflow_s = .TRUE.
+          bc_radiation_s = .TRUE.
+       ELSEIF ( bc_ns == 'radiation/dirichlet'  .OR.  bc_ns == 'nested'  .OR.  &
+                bc_ns == 'nesting_offline' )  THEN
+          bc_dirichlet_s  = .TRUE.
+       ENDIF
+    ENDIF
+
+    IF ( pnorth == MPI_PROC_NULL )  THEN
+       IF ( bc_ns == 'dirichlet/radiation'  .OR.  bc_ns == 'nested'  .OR.      &
+            bc_ns == 'nesting_offline' )  THEN
+          bc_dirichlet_n  = .TRUE.
        ELSEIF ( bc_ns == 'radiation/dirichlet' )  THEN
-          inflow_s  = .TRUE.
-       ELSEIF ( bc_ns == 'nested' )  THEN
-          nest_bound_s = .TRUE.
-       ELSEIF ( bc_ns == 'forcing' )  THEN
-          force_bound_s = .TRUE.
-       ENDIF
-    ENDIF
-
-    IF ( pnorth == MPI_PROC_NULL )  THEN
-       IF ( bc_ns == 'dirichlet/radiation' )  THEN
-          inflow_n  = .TRUE.
-       ELSEIF ( bc_ns == 'radiation/dirichlet' )  THEN
-          outflow_n = .TRUE.
-       ELSEIF ( bc_ns == 'nested' )  THEN
-          nest_bound_n = .TRUE.
-       ELSEIF ( bc_ns == 'forcing' )  THEN
-          force_bound_n = .TRUE.
+          bc_radiation_n = .TRUE.
        ENDIF
     ENDIF
@@ -1283 +1269 @@
 !-- only at the left lateral boundary.
     IF ( use_syn_turb_gen )  THEN
-       IF ( force_bound_l  .OR.  nest_bound_l  .OR.  inflow_l )  THEN
+       IF ( bc_dirichlet_l )  THEN
           id_stg_left_l = myidx
        ELSE
           id_stg_left_l = 0
        ENDIF
-       IF ( force_bound_r  .OR.  nest_bound_r )  THEN
+       IF ( bc_dirichlet_r )  THEN
           id_stg_right_l = myidx
        ELSE
           id_stg_right_l = 0
        ENDIF
-       IF ( force_bound_s  .OR.  nest_bound_s )  THEN
+       IF ( bc_dirichlet_s )  THEN
           id_stg_south_l = myidy
        ELSE
           id_stg_south_l = 0
        ENDIF
-       IF ( force_bound_n  .OR.  nest_bound_n )  THEN
+       IF ( bc_dirichlet_n )  THEN
           id_stg_north_l = myidy
        ELSE
@@ -1324 +1310 @@
 !
 !-- Broadcast the id of the inflow PE
-    IF ( inflow_l )  THEN
+    IF ( bc_dirichlet_l )  THEN
        id_inflow_l = myidx
     ELSE
@@ -1350 +1336 @@
     IF ( turbulent_outflow )  THEN
 
-       IF ( outflow_r )  THEN
+       IF ( bc_radiation_r )  THEN
           id_outflow_l = myidx
        ELSE
@@ -1375 +1361 @@
 #else
     IF ( bc_lr == 'dirichlet/radiation' )  THEN
-       inflow_l  = .TRUE.
-       outflow_r = .TRUE.
+       bc_dirichlet_l = .TRUE.
+       bc_radiation_r = .TRUE.
     ELSEIF ( bc_lr == 'radiation/dirichlet' )  THEN
-       outflow_l = .TRUE.
-       inflow_r  = .TRUE.
+       bc_radiation_l = .TRUE.
+       bc_dirichlet_r = .TRUE.
     ENDIF
 
     IF ( bc_ns == 'dirichlet/radiation' )  THEN
-       inflow_n  = .TRUE.
-       outflow_s = .TRUE.
+       bc_dirichlet_n = .TRUE.
+       bc_radiation_s = .TRUE.
     ELSEIF ( bc_ns == 'radiation/dirichlet' )  THEN
-       outflow_n = .TRUE.
-       inflow_s  = .TRUE.
+       bc_radiation_n = .TRUE.
+       bc_dirichlet_s = .TRUE.
     ENDIF
 #endif
@@ -1394 +1380 @@
 !-- At the inflow or outflow, u or v, respectively, have to be calculated for
 !-- one more grid point.
-    IF ( inflow_l .OR. outflow_l .OR. nest_bound_l .OR. force_bound_l )  THEN
+    IF ( bc_dirichlet_l  .OR.  bc_radiation_l )  THEN
        nxlu = nxl + 1
     ELSE
        nxlu = nxl
     ENDIF
-    IF ( inflow_s .OR. outflow_s .OR. nest_bound_s .OR. force_bound_s )  THEN
+    IF ( bc_dirichlet_s  .OR.  bc_radiation_s )  THEN
        nysv = nys + 1
     ELSE

palm/trunk/SOURCE/large_scale_forcing_nudging_mod.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-! 
+! * Adjustment to new Inifor version: 
+!   - No vertical interpolation/extrapolation of lateral boundary data required
+!     any more (Inifor can treat grid stretching now
+!   - Revise initialization in case of COSMO forcing 
+! * Rename variables and subroutines for offline nesting
 ! 
 ! Former revisions:
@@ -69 +73 @@
 
     USE arrays_3d,                                                             &
-        ONLY:  dzw, e, heatflux_input_conversion, pt, pt_init, q, q_init, s,   &
-               tend, u, u_init, ug, v, v_init, vg, w, w_subs,                  &
+        ONLY:  dzw, e, diss, heatflux_input_conversion, pt, pt_init, q,        &
+               q_init, s, tend, u, u_init, ug, v, v_init, vg, w, w_subs,       &
                waterflux_input_conversion, zu, zw                  
 
     USE control_parameters,                                                    &
-        ONLY:  bc_lr, bc_ns, bc_pt_b, bc_q_b, constant_diffusion,              &
+        ONLY:  bc_dirichlet_l, bc_dirichlet_n, bc_dirichlet_r, bc_dirichlet_s, &
+               bc_lr, bc_ns, bc_pt_b, bc_q_b, constant_diffusion,              &
                constant_heatflux, constant_waterflux,                          &
-               data_output_pr, dt_3d, end_time, forcing,                       &
-               force_bound_l, force_bound_n, force_bound_r, force_bound_s,     &  
+               data_output_pr, dt_3d, end_time,                                &
                humidity, initializing_actions, intermediate_timestep_count,    &
                ibc_pt_b, ibc_q_b,                                              &
                large_scale_forcing, large_scale_subsidence, lsf_surf, lsf_vert,&
-               lsf_exception, message_string, neutral, nudging, passive_scalar,&
-               pt_surface, ocean, q_surface, surface_heatflux,                 &
-               surface_pressure, surface_waterflux, topography,                &
-               use_subsidence_tendencies
+               lsf_exception, message_string, nesting_offline, neutral,        &
+               nudging, passive_scalar, pt_surface, ocean, q_surface,          &
+               surface_heatflux, surface_pressure, surface_waterflux,          &
+               topography, use_subsidence_tendencies
 
     USE grid_variables
-
-    USE pegrid
 
     USE indices,                                                               &
@@ -96 +98 @@
     USE kinds
 
+    USE netcdf_data_input_mod,                                                 &
+        ONLY:  nest_offl
+
+    USE pegrid
+
     USE surface_mod,                                                           &
         ONLY:  surf_def_h, surf_lsm_h, surf_usm_h
@@ -101 +108 @@
     USE statistics,                                                            &
         ONLY:  hom, statistic_regions, sums_ls_l, weight_substep
-
-    USE netcdf_data_input_mod,                                                 &
-        ONLY:  force, netcdf_data_input_interpolate
 
     INTEGER(iwp) ::  nlsf = 1000                       !< maximum number of profiles in LSF_DATA (large scale forcing)
@@ -140 +144 @@
 !
 !-- Public subroutines
-    PUBLIC ls_forcing_surf, ls_forcing_vert, ls_advec, lsf_init,               &
+    PUBLIC calc_tnudge, ls_forcing_surf, ls_forcing_vert, ls_advec, lsf_init,  &
            lsf_nudging_check_parameters, nudge_init,                           &
            lsf_nudging_check_data_output_pr, lsf_nudging_header,               &
-           calc_tnudge, nudge, nudge_ref, forcing_bc_mass_conservation,        &
-           forcing_bc
+           lsf_nesting_offline, lsf_nesting_offline_mass_conservation,         & 
+           nudge, nudge_ref
+           
 !
 !-- Public variables
     PUBLIC qsws_surf, shf_surf, td_lsa_lpt, td_lsa_q, td_sub_lpt,              &
-           td_sub_q, time_vert, force
+           td_sub_q, time_vert
 
 
@@ -167 +172 @@
 ! Description:
 ! ------------
-!> @todo Missing subroutine description.
-!------------------------------------------------------------------------------!
-    SUBROUTINE forcing_bc_mass_conservation
+!> In this subroutine a constant mass within the model domain is guaranteed. 
+!> Larger-scale models may be based on a compressible equation system, which is
+!> not consistent with PALMs incompressible equation system. In order to avoid
+!> a decrease or increase of mass during the simulation, non-divergent flow
+!> through the lateral and top boundaries is compensated by the vertical wind 
+!> component at the top boundary. 
+!------------------------------------------------------------------------------!
+    SUBROUTINE lsf_nesting_offline_mass_conservation
 
        USE control_parameters,                                                 &
@@ -176 +186 @@
        IMPLICIT NONE
 
-       INTEGER(iwp) ::  i !<
-       INTEGER(iwp) ::  j !<
-       INTEGER(iwp) ::  k !<
-
-       REAL(wp) ::  w_correct !< 
-       REAL(wp), DIMENSION(1:3) ::  volume_flow_l   !<
+       INTEGER(iwp) ::  i !< grid index in x-direction
+       INTEGER(iwp) ::  j !< grid index in y-direction
+       INTEGER(iwp) ::  k !< grid index in z-direction
+
+       REAL(wp) ::  w_correct                       !< vertical velocity increment required to compensate non-divergent flow through the boundaries
+       REAL(wp), DIMENSION(1:3) ::  volume_flow_l   !< local volume flow
 
        volume_flow   = 0.0_wp
@@ -188 +198 @@
        d_area_t = 1.0_wp / ( ( nx + 1 ) * dx * ( ny + 1 ) * dy )
 
-       IF ( force_bound_l )  THEN
+       IF ( bc_dirichlet_l )  THEN
           i = nxl
           DO  j = nys, nyn
@@ -198 +208 @@
           ENDDO
        ENDIF
-       IF ( force_bound_r )  THEN
+       IF ( bc_dirichlet_r )  THEN
           i = nxr+1
           DO  j = nys, nyn
@@ -208 +218 @@
           ENDDO
        ENDIF
-       IF ( force_bound_s )  THEN
+       IF ( bc_dirichlet_s )  THEN
           j = nys
           DO  i = nxl, nxr
@@ -218 +228 @@
           ENDDO
        ENDIF
-       IF ( force_bound_n )  THEN
+       IF ( bc_dirichlet_n )  THEN
           j = nyn+1
           DO  i = nxl, nxr
@@ -255 +265 @@
        ENDDO
 
-write(9,*) "w correction", w_correct
-flush(9)
-
-    END SUBROUTINE forcing_bc_mass_conservation
+    END SUBROUTINE lsf_nesting_offline_mass_conservation
 
 
@@ -264 +271 @@
 ! Description:
 ! ------------
-!> @todo Missing subroutine description.
-!------------------------------------------------------------------------------!
-    SUBROUTINE forcing_bc
+!> Set the lateral and top boundary conditions in case the PALM domain is 
+!> nested offline in a mesoscale model. 
+!------------------------------------------------------------------------------!
+    SUBROUTINE lsf_nesting_offline
 
        USE control_parameters,                                                 &
-           ONLY:  force_bound_l, force_bound_n, force_bound_r, force_bound_s,  &
-                  humidity, neutral, passive_scalar, simulated_time
-
-       USE netcdf_data_input_mod,                                              &
-           ONLY:  force                              
+           ONLY:  bc_dirichlet_l, bc_dirichlet_n, bc_dirichlet_r,              &
+                  bc_dirichlet_s, humidity, neutral, passive_scalar, rans_mode,&
+                  rans_tke_e, time_since_reference_point                      
 
        IMPLICIT NONE
@@ -284 +290 @@
        REAL(wp) ::  ddt_lsf !< inverse value of time resolution of forcing data
        REAL(wp) ::  t_ref   !< time past since last reference step
-
-!
-!--    If required, interpolate and/or extrapolate data vertically. This is 
-!--    required as Inifor outputs only equidistant vertical data.
-       IF ( ANY( zu(1:nzt+1) /= force%zu_atmos(1:force%nzu) ) )  THEN
-          IF ( .NOT. force%interpolated )  THEN
-
-             DO  t = 0, 1
-                IF ( force_bound_l )  THEN
-                   CALL netcdf_data_input_interpolate( force%u_left(t,:,:),    &
-                                                       zu(nzb+1:nzt+1),        &
-                                                       force%zu_atmos )
-                   CALL netcdf_data_input_interpolate( force%v_left(t,:,:),    &
-                                                       zu(nzb+1:nzt+1),        &
-                                                       force%zu_atmos )
-                   CALL netcdf_data_input_interpolate( force%w_left(t,:,:),    &
-                                                       zw(nzb+1:nzt+1),        &
-                                                       force%zw_atmos )
-                   IF ( .NOT. neutral )                                        &
-                      CALL netcdf_data_input_interpolate( force%pt_left(t,:,:),&
-                                                          zu(nzb+1:nzt+1),     &
-                                                          force%zu_atmos )
-                   IF ( humidity )                                             &
-                      CALL netcdf_data_input_interpolate( force%q_left(t,:,:), &
-                                                          zu(nzb+1:nzt+1),     &
-                                                          force%zu_atmos )
-                ENDIF
-                IF ( force_bound_r )  THEN
-                   CALL netcdf_data_input_interpolate( force%u_right(t,:,:),   &
-                                                       zu(nzb+1:nzt+1),        &
-                                                       force%zu_atmos )
-                   CALL netcdf_data_input_interpolate( force%v_right(t,:,:),   &
-                                                       zu(nzb+1:nzt+1),        &
-                                                       force%zu_atmos )
-                   CALL netcdf_data_input_interpolate( force%w_right(t,:,:),   &
-                                                       zw(nzb+1:nzt+1),        &
-                                                       force%zw_atmos )
-                   IF ( .NOT. neutral )                                        &
-                      CALL netcdf_data_input_interpolate( force%pt_right(t,:,:),&
-                                                          zu(nzb+1:nzt+1),     &
-                                                          force%zu_atmos )
-                   IF ( humidity )                                             &
-                      CALL netcdf_data_input_interpolate( force%q_right(t,:,:),&
-                                                          zu(nzb+1:nzt+1),     &
-                                                          force%zu_atmos )
-                ENDIF
-                IF ( force_bound_n )  THEN
-                   CALL netcdf_data_input_interpolate( force%u_north(t,:,:),   &
-                                                       zu(nzb+1:nzt+1),        &
-                                                       force%zu_atmos )
-                   CALL netcdf_data_input_interpolate( force%v_north(t,:,:),   &
-                                                       zu(nzb+1:nzt+1),        &
-                                                       force%zu_atmos )
-                   CALL netcdf_data_input_interpolate( force%w_north(t,:,:),   &
-                                                       zw(nzb+1:nzt+1),        &
-                                                       force%zw_atmos )
-                   IF ( .NOT. neutral )                                        &
-                      CALL netcdf_data_input_interpolate( force%pt_north(t,:,:),&
-                                                          zu(nzb+1:nzt+1),     &
-                                                          force%zu_atmos )
-                   IF ( humidity )                                             &
-                      CALL netcdf_data_input_interpolate( force%q_north(t,:,:),&
-                                                          zu(nzb+1:nzt+1),     &
-                                                          force%zu_atmos )
-                ENDIF
-                IF ( force_bound_s )  THEN
-                   CALL netcdf_data_input_interpolate( force%u_south(t,:,:),   &
-                                                       zu(nzb+1:nzt+1),        &
-                                                       force%zu_atmos )
-                   CALL netcdf_data_input_interpolate( force%v_south(t,:,:),   &
-                                                       zu(nzb+1:nzt+1),        &
-                                                       force%zu_atmos )
-                   CALL netcdf_data_input_interpolate( force%w_south(t,:,:),   &
-                                                       zw(nzb+1:nzt+1),        &
-                                                       force%zw_atmos )
-                   IF ( .NOT. neutral )                                        &
-                      CALL netcdf_data_input_interpolate( force%pt_south(t,:,:),&
-                                                          zu(nzb+1:nzt+1),     &
-                                                          force%zu_atmos )
-                   IF ( humidity )                                             &
-                      CALL netcdf_data_input_interpolate( force%q_south(t,:,:),&
-                                                          zu(nzb+1:nzt+1),     &
-                                                          force%zu_atmos )
-                ENDIF
-             ENDDO
-!
-!--          Note, no interpolation of top boundary. Just use initial value. 
-!--          No physical meaningful extrapolation possible if only one layer is
-!--          given. 
-
-             force%interpolated = .TRUE. 
-          ENDIF
-       ENDIF
-      
+     
 !
 !--    Calculate time interval of forcing data       
-       ddt_lsf = 1.0_wp / ( force%time(force%tind_p) - force%time(force%tind) )
+       ddt_lsf = 1.0_wp / ( nest_offl%time(nest_offl%tind_p) -                 &
+                            nest_offl%time(nest_offl%tind) )
 !
 !--    Calculate reziproke time past since last reference step. Please note, 
-!--    as simulated time is still not updated, the actual time here is 
-!--    simulated time + dt_3d
-       t_ref = simulated_time + dt_3d - force%time(force%tind)
-
-       IF ( force_bound_l )  THEN
-
-          DO  j = nys, nyn
-             DO  k = nzb+1, nzt+1
-                u(k,j,nxlg:nxl) = force%u_left(0,k,j) + ddt_lsf * t_ref *      &
-                         ( force%u_left(1,k,j) - force%u_left(0,k,j) ) *       &
-                           MERGE( 1.0_wp, 0.0_wp,                              &
-                                  BTEST( wall_flags_0(k,j,nxlg:nxl), 1 ) )
-             ENDDO
-          ENDDO
+!--    the time coordinate is still not updated, so that the actual time need
+!--    to be incremented by dt_3d. Moreover, note that the simulation time
+!--    passed since simulation start is time_since_reference_point, not 
+!--    simulated_time!
+       t_ref = time_since_reference_point + dt_3d -                            &
+                                            nest_offl%time(nest_offl%tind)
+                                            
+       IF ( bc_dirichlet_l )  THEN
 
           DO  j = nys, nyn
              DO  k = nzb+1, nzt
-                w(k,j,nxlg:nxl-1) = force%w_left(0,k,j) + ddt_lsf * t_ref *    &
-                         ( force%w_left(1,k,j) - force%w_left(0,k,j) ) *       &
-                           MERGE( 1.0_wp, 0.0_wp,                              &
-                                  BTEST( wall_flags_0(k,j,nxlg:nxl-1), 3 ) )
+                u(k,j,nxlg:nxl) = nest_offl%u_left(0,k,j) + ddt_lsf * t_ref  * &
+                       ( nest_offl%u_left(1,k,j) - nest_offl%u_left(0,k,j) ) * &
+                         MERGE( 1.0_wp, 0.0_wp,                                &
+                                BTEST( wall_flags_0(k,j,nxlg:nxl), 1 ) )
+             ENDDO
+          ENDDO
+
+          DO  j = nys, nyn
+             DO  k = nzb+1, nzt-1
+                w(k,j,nxlg:nxl-1) = nest_offl%w_left(0,k,j) + ddt_lsf * t_ref *&
+                       ( nest_offl%w_left(1,k,j) - nest_offl%w_left(0,k,j) )  *&
+                         MERGE( 1.0_wp, 0.0_wp,                                &
+                                BTEST( wall_flags_0(k,j,nxlg:nxl-1), 3 ) )
              ENDDO
           ENDDO
 
           DO  j = nysv, nyn
-             DO  k = nzb+1, nzt+1
-                v(k,j,nxlg:nxl-1) = force%v_left(0,k,j) + ddt_lsf * t_ref *    &
-                         ( force%v_left(1,k,j) - force%v_left(0,k,j) ) *       &
-                           MERGE( 1.0_wp, 0.0_wp,                              &
-                                  BTEST( wall_flags_0(k,j,nxlg:nxl-1), 2 ) )
+             DO  k = nzb+1, nzt
+                v(k,j,nxlg:nxl-1) = nest_offl%v_left(0,k,j) + ddt_lsf * t_ref *&
+                       ( nest_offl%v_left(1,k,j) - nest_offl%v_left(0,k,j) )  *&
+                         MERGE( 1.0_wp, 0.0_wp,                                &
+                                BTEST( wall_flags_0(k,j,nxlg:nxl-1), 2 ) )
              ENDDO
           ENDDO
@@ -418 +335 @@
           IF ( .NOT. neutral )  THEN
              DO  j = nys, nyn
-                DO  k = nzb+1, nzt+1
-                   pt(k,j,nxlg:nxl-1) = force%pt_left(0,k,j) + ddt_lsf *       &
-                                                                  t_ref   *    &
-                       ( force%pt_left(1,k,j) - force%pt_left(0,k,j) )
+                DO  k = nzb+1, nzt
+                   pt(k,j,nxlg:nxl-1) = nest_offl%pt_left(0,k,j) + ddt_lsf *   &
+                                                                   t_ref   *   &
+                       ( nest_offl%pt_left(1,k,j) - nest_offl%pt_left(0,k,j) )
  
                 ENDDO
@@ -429 +346 @@
           IF ( humidity )  THEN
              DO  j = nys, nyn
-                DO  k = nzb+1, nzt+1
-                   q(k,j,nxlg:nxl-1) = force%q_left(0,k,j) + ddt_lsf   *       &
-                                                                  t_ref   *    &
-                       ( force%q_left(1,k,j) - force%q_left(0,k,j) )
+                DO  k = nzb+1, nzt
+                   q(k,j,nxlg:nxl-1) = nest_offl%q_left(0,k,j) + ddt_lsf   *   &
+                                                                 t_ref     *   &
+                       ( nest_offl%q_left(1,k,j) - nest_offl%q_left(0,k,j) )
  
                 ENDDO
@@ -440 +357 @@
        ENDIF
 
-       IF ( force_bound_r )  THEN
-
-          DO  j = nys, nyn
-             DO  k = nzb+1, nzt+1
-                u(k,j,nxr+1:nxrg) = force%u_right(0,k,j) + ddt_lsf * t_ref *   &
-                        ( force%u_right(1,k,j) - force%u_right(0,k,j) ) *      &
-                           MERGE( 1.0_wp, 0.0_wp,                              &
-                                  BTEST( wall_flags_0(k,j,nxr+1:nxrg), 1 ) )
-
-             ENDDO
-          ENDDO
+       IF ( bc_dirichlet_r )  THEN
+
           DO  j = nys, nyn
              DO  k = nzb+1, nzt
-                w(k,j,nxr+1:nxrg) = force%w_right(0,k,j) + ddt_lsf * t_ref *   &
-                        ( force%w_right(1,k,j) - force%w_right(0,k,j) ) *      &
-                           MERGE( 1.0_wp, 0.0_wp,                              &
-                                  BTEST( wall_flags_0(k,j,nxr+1:nxrg), 3 ) )
+                u(k,j,nxr+1:nxrg) = nest_offl%u_right(0,k,j) + ddt_lsf * t_ref *&
+                      ( nest_offl%u_right(1,k,j) - nest_offl%u_right(0,k,j) )  *&
+                         MERGE( 1.0_wp, 0.0_wp,                                 &
+                                BTEST( wall_flags_0(k,j,nxr+1:nxrg), 1 ) )
+
+             ENDDO
+          ENDDO
+          DO  j = nys, nyn
+             DO  k = nzb+1, nzt-1
+                w(k,j,nxr+1:nxrg) = nest_offl%w_right(0,k,j) + ddt_lsf * t_ref *&
+                      ( nest_offl%w_right(1,k,j) - nest_offl%w_right(0,k,j) )  *&
+                         MERGE( 1.0_wp, 0.0_wp,                                 &
+                                BTEST( wall_flags_0(k,j,nxr+1:nxrg), 3 ) )
              ENDDO
           ENDDO
 
           DO  j = nysv, nyn
-             DO  k = nzb+1, nzt+1
-                v(k,j,nxr+1:nxrg) = force%v_right(0,k,j) + ddt_lsf * t_ref *   &
-                        ( force%v_right(1,k,j) - force%v_right(0,k,j) ) *      &
-                           MERGE( 1.0_wp, 0.0_wp,                              &
-                                  BTEST( wall_flags_0(k,j,nxr+1:nxrg), 2 ) )
+             DO  k = nzb+1, nzt
+                v(k,j,nxr+1:nxrg) = nest_offl%v_right(0,k,j) + ddt_lsf * t_ref *&
+                      ( nest_offl%v_right(1,k,j) - nest_offl%v_right(0,k,j) )  *&
+                         MERGE( 1.0_wp, 0.0_wp,                                 &
+                                BTEST( wall_flags_0(k,j,nxr+1:nxrg), 2 ) )
              ENDDO
           ENDDO
@@ -471 +388 @@
           IF ( .NOT. neutral )  THEN
              DO  j = nys, nyn
-                DO  k = nzb+1, nzt+1
-                   pt(k,j,nxr+1:nxrg) = force%pt_right(0,k,j) + ddt_lsf *      &
-                                                                  t_ref   *    &
-                     ( force%pt_right(1,k,j) - force%pt_right(0,k,j) )
+                DO  k = nzb+1, nzt
+                   pt(k,j,nxr+1:nxrg) = nest_offl%pt_right(0,k,j) + ddt_lsf *  &
+                                                                    t_ref   *  &
+                     ( nest_offl%pt_right(1,k,j) - nest_offl%pt_right(0,k,j) )
  
                 ENDDO
@@ -482 +399 @@
           IF ( humidity )  THEN
              DO  j = nys, nyn
-                DO  k = nzb+1, nzt+1
-                   q(k,j,nxr+1:nxrg) = force%q_right(0,k,j) + ddt_lsf   *      &
-                                                                  t_ref   *    &
-                       ( force%q_right(1,k,j) - force%q_right(0,k,j) )
+                DO  k = nzb+1, nzt
+                   q(k,j,nxr+1:nxrg) = nest_offl%q_right(0,k,j) + ddt_lsf   *  &
+                                                                    t_ref   *  &
+                       ( nest_offl%q_right(1,k,j) - nest_offl%q_right(0,k,j) )
  
                 ENDDO
@@ -493 +410 @@
        ENDIF
 
-       IF ( force_bound_s )  THEN
-
-          DO  i = nxl, nxr
-             DO  k = nzb+1, nzt+1
-                v(k,nysg:nys,i)   = force%v_south(0,k,i) + ddt_lsf * t_ref *   &
-                        ( force%v_south(1,k,i) - force%v_south(0,k,i) ) *      &
-                           MERGE( 1.0_wp, 0.0_wp,                              &
-                                  BTEST( wall_flags_0(k,nysg:nys,i), 2 ) )
-             ENDDO
-          ENDDO
+       IF ( bc_dirichlet_s )  THEN
 
           DO  i = nxl, nxr
              DO  k = nzb+1, nzt
-                w(k,nysg:nys-1,i) = force%w_south(0,k,i) + ddt_lsf * t_ref *   &
-                        ( force%w_south(1,k,i) - force%w_south(0,k,i) ) *      &
-                           MERGE( 1.0_wp, 0.0_wp,                              &
+                v(k,nysg:nys,i)   = nest_offl%v_south(0,k,i) + ddt_lsf * t_ref *&
+                      ( nest_offl%v_south(1,k,i) - nest_offl%v_south(0,k,i) )  *&
+                         MERGE( 1.0_wp, 0.0_wp,                                 &
+                                BTEST( wall_flags_0(k,nysg:nys,i), 2 ) )
+             ENDDO
+          ENDDO
+
+          DO  i = nxl, nxr
+             DO  k = nzb+1, nzt-1
+                w(k,nysg:nys-1,i) = nest_offl%w_south(0,k,i) + ddt_lsf * t_ref  *&
+                        ( nest_offl%w_south(1,k,i) - nest_offl%w_south(0,k,i) ) *&
+                           MERGE( 1.0_wp, 0.0_wp,                                &
                                   BTEST( wall_flags_0(k,nysg:nys-1,i), 3 ) )
              ENDDO
@@ -514 +431 @@
 
           DO  i = nxlu, nxr
-             DO  k = nzb+1, nzt+1
-                u(k,nysg:nys-1,i) = force%u_south(0,k,i) + ddt_lsf * t_ref *   &
-                        ( force%u_south(1,k,i) - force%u_south(0,k,i) ) *      &
-                           MERGE( 1.0_wp, 0.0_wp,                              &
+             DO  k = nzb+1, nzt
+                u(k,nysg:nys-1,i) = nest_offl%u_south(0,k,i) + ddt_lsf * t_ref  *&
+                        ( nest_offl%u_south(1,k,i) - nest_offl%u_south(0,k,i) ) *&
+                           MERGE( 1.0_wp, 0.0_wp,                                &
                                   BTEST( wall_flags_0(k,nysg:nys-1,i), 1 ) )
              ENDDO
@@ -524 +441 @@
           IF ( .NOT. neutral )  THEN
              DO  i = nxl, nxr
-                DO  k = nzb+1, nzt+1
-                   pt(k,nysg:nys-1,i) = force%pt_south(0,k,i) + ddt_lsf *      &
-                                                                  t_ref   *    &
-                     ( force%pt_south(1,k,i) - force%pt_south(0,k,i) )
+                DO  k = nzb+1, nzt
+                   pt(k,nysg:nys-1,i) = nest_offl%pt_south(0,k,i) + ddt_lsf *  &
+                                                                    t_ref   *  &
+                     ( nest_offl%pt_south(1,k,i) - nest_offl%pt_south(0,k,i) )
  
                 ENDDO
@@ -535 +452 @@
           IF ( humidity )  THEN
              DO  i = nxl, nxr
-                DO  k = nzb+1, nzt+1
-                   q(k,nysg:nys-1,i) = force%q_south(0,k,i) + ddt_lsf   *      &
-                                                                  t_ref   *    &
-                       ( force%q_south(1,k,i) - force%q_south(0,k,i) )
+                DO  k = nzb+1, nzt
+                   q(k,nysg:nys-1,i) = nest_offl%q_south(0,k,i) + ddt_lsf   *  &
+                                                                    t_ref   *  &
+                       ( nest_offl%q_south(1,k,i) - nest_offl%q_south(0,k,i) )
  
                 ENDDO
@@ -546 +463 @@
        ENDIF
 
-       IF ( force_bound_n )  THEN
-
-          DO  i = nxl, nxr
-             DO  k = nzb+1, nzt+1
-                v(k,nyn+1:nyng,i)   = force%v_north(0,k,i) + ddt_lsf * t_ref * &
-                        ( force%v_north(1,k,i) - force%v_north(0,k,i) ) *      &
-                           MERGE( 1.0_wp, 0.0_wp,                              &
-                                  BTEST( wall_flags_0(k,nyn+1:nyng,i), 2 ) )
-             ENDDO
-          ENDDO
+       IF ( bc_dirichlet_n )  THEN
+
           DO  i = nxl, nxr
              DO  k = nzb+1, nzt
-                w(k,nyn+1:nyng,i) = force%w_north(0,k,i) + ddt_lsf * t_ref *   &
-                        ( force%w_north(1,k,i) - force%w_north(0,k,i) ) *      &
-                           MERGE( 1.0_wp, 0.0_wp,                              &
+                v(k,nyn+1:nyng,i)   = nest_offl%v_north(0,k,i) + ddt_lsf * t_ref *&
+                        ( nest_offl%v_north(1,k,i) - nest_offl%v_north(0,k,i) )  *&
+                           MERGE( 1.0_wp, 0.0_wp,                                 &
+                                  BTEST( wall_flags_0(k,nyn+1:nyng,i), 2 ) )
+             ENDDO
+          ENDDO
+          DO  i = nxl, nxr
+             DO  k = nzb+1, nzt-1
+                w(k,nyn+1:nyng,i) = nest_offl%w_north(0,k,i) + ddt_lsf * t_ref  *&
+                        ( nest_offl%w_north(1,k,i) - nest_offl%w_north(0,k,i) ) *&
+                           MERGE( 1.0_wp, 0.0_wp,                                &
                                   BTEST( wall_flags_0(k,nyn+1:nyng,i), 3 ) )
              ENDDO
@@ -566 +483 @@
 
           DO  i = nxlu, nxr
-             DO  k = nzb+1, nzt+1
-                u(k,nyn+1:nyng,i) = force%u_north(0,k,i) + ddt_lsf * t_ref *   &
-                        ( force%u_north(1,k,i) - force%u_north(0,k,i) ) *      &
-                           MERGE( 1.0_wp, 0.0_wp,                              &
+             DO  k = nzb+1, nzt
+                u(k,nyn+1:nyng,i) = nest_offl%u_north(0,k,i) + ddt_lsf * t_ref  *&
+                        ( nest_offl%u_north(1,k,i) - nest_offl%u_north(0,k,i) ) *&
+                           MERGE( 1.0_wp, 0.0_wp,                                &
                                   BTEST( wall_flags_0(k,nyn+1:nyng,i), 1 ) )
 
@@ -577 +494 @@
           IF ( .NOT. neutral )  THEN
              DO  i = nxl, nxr
-                DO  k = nzb+1, nzt+1
-                   pt(k,nyn+1:nyng,i) = force%pt_north(0,k,i) + ddt_lsf *      &
-                                                                  t_ref   *    &
-                     ( force%pt_north(1,k,i) - force%pt_north(0,k,i) )
+                DO  k = nzb+1, nzt
+                   pt(k,nyn+1:nyng,i) = nest_offl%pt_north(0,k,i) + ddt_lsf *  &
+                                                                    t_ref   *  &
+                     ( nest_offl%pt_north(1,k,i) - nest_offl%pt_north(0,k,i) )
  
                 ENDDO
@@ -588 +505 @@
           IF ( humidity )  THEN
              DO  i = nxl, nxr
-                DO  k = nzb+1, nzt+1
-                   q(k,nyn+1:nyng,i) = force%q_north(0,k,i) + ddt_lsf   *      &
-                                                                  t_ref   *    &
-                       ( force%q_north(1,k,i) - force%q_north(0,k,i) )
+                DO  k = nzb+1, nzt
+                   q(k,nyn+1:nyng,i) = nest_offl%q_north(0,k,i) + ddt_lsf   *  &
+                                                                    t_ref   *  &
+                       ( nest_offl%q_north(1,k,i) - nest_offl%q_north(0,k,i) )
  
                 ENDDO
@@ -600 +517 @@
 !
 !--    Top boundary. 
-!--    Please note, only map Inifor data on model top in case the numeric is
-!--    identical to the Inifor grid. At the top boundary an extrapolation is 
-!--    not possible. 
        DO  i = nxlu, nxr
           DO  j = nys, nyn
-             u(nzt+1,j,i) = force%u_top(0,j,i) + ddt_lsf * t_ref *             &
-                        ( force%u_top(1,j,i) - force%u_top(0,j,i) ) *          &
-                           MERGE( 1.0_wp, 0.0_wp,                             &
+             u(nzt+1,j,i) = nest_offl%u_top(0,j,i) + ddt_lsf * t_ref *         &
+                        ( nest_offl%u_top(1,j,i) - nest_offl%u_top(0,j,i) ) *  &
+                           MERGE( 1.0_wp, 0.0_wp,                              &
                                   BTEST( wall_flags_0(nzt+1,j,i), 1 ) )
           ENDDO
@@ -614 +528 @@
        DO  i = nxl, nxr
           DO  j = nysv, nyn
-             v(nzt+1,j,i) = force%v_top(0,j,i) + ddt_lsf * t_ref *             &
-                        ( force%v_top(1,j,i) - force%v_top(0,j,i) ) *          &
+             v(nzt+1,j,i) = nest_offl%v_top(0,j,i) + ddt_lsf * t_ref *         &
+                        ( nest_offl%v_top(1,j,i) - nest_offl%v_top(0,j,i) ) *  &
                            MERGE( 1.0_wp, 0.0_wp,                              &
                                   BTEST( wall_flags_0(nzt+1,j,i), 2 ) )
@@ -623 +537 @@
        DO  i = nxl, nxr
           DO  j = nys, nyn
-             w(nzt:nzt+1,j,i) = force%w_top(0,j,i) + ddt_lsf * t_ref *         &
-                        ( force%w_top(1,j,i) - force%w_top(0,j,i) ) *          &
+             w(nzt:nzt+1,j,i) = nest_offl%w_top(0,j,i) + ddt_lsf * t_ref *     &
+                        ( nest_offl%w_top(1,j,i) - nest_offl%w_top(0,j,i) ) *  &
                            MERGE( 1.0_wp, 0.0_wp,                              &
                                   BTEST( wall_flags_0(nzt:nzt+1,j,i), 3 ) )
@@ -634 +548 @@
           DO  i = nxl, nxr
              DO  j = nys, nyn
-                pt(nzt+1,j,i) = force%pt_top(0,j,i) + ddt_lsf * t_ref *        &
-                        ( force%pt_top(1,j,i) - force%pt_top(0,j,i) )
+                pt(nzt+1,j,i) = nest_offl%pt_top(0,j,i) + ddt_lsf * t_ref *    &
+                        ( nest_offl%pt_top(1,j,i) - nest_offl%pt_top(0,j,i) )
              ENDDO
           ENDDO
@@ -643 +557 @@
           DO  i = nxl, nxr
              DO  j = nys, nyn
-                q(nzt+1,j,i) = force%q_top(0,j,i) + ddt_lsf * t_ref *          &
-                        ( force%q_top(1,j,i) - force%q_top(0,j,i) )
+                q(nzt+1,j,i) = nest_offl%q_top(0,j,i) + ddt_lsf * t_ref *      &
+                        ( nest_offl%q_top(1,j,i) - nest_offl%q_top(0,j,i) )
              ENDDO
           ENDDO
@@ -651 +565 @@
 !--    At the edges( left-south, left-north, right-south and right-north) set 
 !--    data on ghost points. 
-       IF ( force_bound_l  .AND.  force_bound_s )  THEN
+       IF ( bc_dirichlet_l  .AND.  bc_dirichlet_s )  THEN
           DO  i = 1, nbgp
              u(:,nys-i,nxlg:nxl)   = u(:,nys,nxlg:nxl)
              w(:,nys-i,nxlg:nxl-1) = w(:,nys,nxlg:nxl-1)
              IF ( .NOT. neutral )  pt(:,nys-i,nxlg:nxl-1) = pt(:,nys,nxlg:nxl-1)
-             IF ( humidity )       q(:,nys-i,nxlg:nxl-1)  = q(:,nys,nxlg:nxl-1)
+             IF ( humidity      )  q(:,nys-i,nxlg:nxl-1)  = q(:,nys,nxlg:nxl-1)
           ENDDO
           DO  i = 1, nbgp+1
@@ -662 +576 @@
           ENDDO
        ENDIF
-       IF ( force_bound_l  .AND.  force_bound_n )  THEN
+       IF ( bc_dirichlet_l  .AND.  bc_dirichlet_n )  THEN
           DO  i = 1, nbgp
              u(:,nyn+i,nxlg:nxl)   = u(:,nyn,nxlg:nxl)
@@ -668 +582 @@
              w(:,nyn+i,nxlg:nxl-1) = w(:,nyn,nxlg:nxl-1)
              IF ( .NOT. neutral )  pt(:,nyn+i,nxlg:nxl-1) = pt(:,nyn,nxlg:nxl-1)
-             IF ( humidity )       q(:,nyn+i,nxlg:nxl-1)  = q(:,nyn,nxlg:nxl-1)
-          ENDDO
-       ENDIF
-       IF ( force_bound_r  .AND.  force_bound_s )  THEN
+             IF ( humidity      )  q(:,nyn+i,nxlg:nxl-1)  = q(:,nyn,nxlg:nxl-1)
+          ENDDO
+       ENDIF
+       IF ( bc_dirichlet_r  .AND.  bc_dirichlet_s )  THEN
           DO  i = 1, nbgp
              u(:,nys-i,nxr+1:nxrg) = u(:,nys,nxr+1:nxrg)
              w(:,nys-i,nxr+1:nxrg) = w(:,nys,nxr+1:nxrg)
              IF ( .NOT. neutral )  pt(:,nys-i,nxr+1:nxrg) = pt(:,nys,nxr+1:nxrg)
-             IF ( humidity )       q(:,nys-i,nxr+1:nxrg)  = q(:,nys,nxr+1:nxrg)
+             IF ( humidity      )  q(:,nys-i,nxr+1:nxrg)  = q(:,nys,nxr+1:nxrg)
           ENDDO
           DO  i = 1, nbgp+1
@@ -682 +596 @@
           ENDDO
        ENDIF
-       IF ( force_bound_r  .AND.  force_bound_n )  THEN
+       IF ( bc_dirichlet_r  .AND.  bc_dirichlet_n )  THEN
           DO  i = 1, nbgp
              u(:,nyn+i,nxr+1:nxrg) = u(:,nyn,nxr+1:nxrg)
@@ -688 +602 @@
              w(:,nyn+i,nxr+1:nxrg) = w(:,nyn,nxr+1:nxrg)
              IF ( .NOT. neutral )  pt(:,nyn+i,nxr+1:nxrg) = pt(:,nyn,nxr+1:nxrg)
-             IF ( humidity )       q(:,nyn+i,nxr+1:nxrg)  = q(:,nyn,nxr+1:nxrg)
-          ENDDO
-       ENDIF
-!
-!--    Moreover, set Neumann boundary condition for subgrid-scale TKE and 
-!--    passive scalar
+             IF ( humidity      )  q(:,nyn+i,nxr+1:nxrg)  = q(:,nyn,nxr+1:nxrg)
+          ENDDO
+       ENDIF
+!
+!--    Moreover, set Neumann boundary condition for subgrid-scale TKE, 
+!--    passive scalar, dissipation, and chemical species if required
+       IF ( rans_mode  .AND.  rans_tke_e )  THEN
+          IF (  bc_dirichlet_l )  diss(:,:,nxl-1) = diss(:,:,nxl)
+          IF (  bc_dirichlet_r )  diss(:,:,nxr+1) = diss(:,:,nxr)
+          IF (  bc_dirichlet_s )  diss(:,nys-1,:) = diss(:,nys,:)
+          IF (  bc_dirichlet_n )  diss(:,nyn+1,:) = diss(:,nyn,:)
+       ENDIF
        IF ( .NOT. constant_diffusion )  THEN
-          IF (  force_bound_l )  e(:,:,nxl-1) = e(:,:,nxl)
-          IF (  force_bound_r )  e(:,:,nxr+1) = e(:,:,nxr)
-          IF (  force_bound_s )  e(:,nys-1,:) = e(:,nys,:)
-          IF (  force_bound_n )  e(:,nyn+1,:) = e(:,nyn,:)
+          IF (  bc_dirichlet_l )  e(:,:,nxl-1) = e(:,:,nxl)
+          IF (  bc_dirichlet_r )  e(:,:,nxr+1) = e(:,:,nxr)
+          IF (  bc_dirichlet_s )  e(:,nys-1,:) = e(:,nys,:)
+          IF (  bc_dirichlet_n )  e(:,nyn+1,:) = e(:,nyn,:)
           e(nzt+1,:,:) = e(nzt,:,:)
        ENDIF
        IF ( passive_scalar )  THEN
-          IF (  force_bound_l )  s(:,:,nxl-1) = s(:,:,nxl)
-          IF (  force_bound_r )  s(:,:,nxr+1) = s(:,:,nxr)
-          IF (  force_bound_s )  s(:,nys-1,:) = s(:,nys,:)
-          IF (  force_bound_n )  s(:,nyn+1,:) = s(:,nyn,:)
-       ENDIF
+          IF (  bc_dirichlet_l )  s(:,:,nxl-1) = s(:,:,nxl)
+          IF (  bc_dirichlet_r )  s(:,:,nxr+1) = s(:,:,nxr)
+          IF (  bc_dirichlet_s )  s(:,nys-1,:) = s(:,nys,:)
+          IF (  bc_dirichlet_n )  s(:,nyn+1,:) = s(:,nyn,:)
+       ENDIF
+
 
 
@@ -720 +641 @@
 !--    treatment of fluxes. 
 !--    For the moment, comment this out! 
-!      surface_pressure = force%surface_pressure(force%tind) +                 &
+!      surface_pressure = nest_offl%surface_pressure(nest_offl%tind) +                 &
 !                                                      ddt_lsf * t_ref *       &
-!                                    ( force%surface_pressure(force%tind_p)    &
-!                                    - force%surface_pressure(force%tind) )
-
-    END SUBROUTINE forcing_bc
+!                                    ( nest_offl%surface_pressure(nest_offl%tind_p)    &
+!                                    - nest_offl%surface_pressure(nest_offl%tind) )
+
+    END SUBROUTINE lsf_nesting_offline
 
 !------------------------------------------------------------------------------!
@@ -1040 +961 @@
        REAL(wp) ::  r_dummy           !<
 
-       IF ( forcing )  THEN
+       IF ( nesting_offline )  THEN
 !
 !--       Allocate arrays for geostrophic wind components. Arrays will 
@@ -1047 +968 @@
 !--       case of cyclic boundary conditions.
           IF ( bc_lr_cyc  .AND.  bc_ns_cyc )  THEN
-             ALLOCATE( force%ug(0:1,nzb:nzt+1) )
-             ALLOCATE( force%vg(0:1,nzb:nzt+1) )
+             ALLOCATE( nest_offl%ug(0:1,nzb:nzt+1) )
+             ALLOCATE( nest_offl%vg(0:1,nzb:nzt+1) )
           ENDIF
 !
 !--       Allocate arrays for reading boundary values. Arrays will incorporate 2 
 !--       time levels in order to interpolate in between. 
-          IF ( force_bound_l )  THEN
-             ALLOCATE( force%u_left(0:1,nzb+1:nzt+1,nys:nyn)  )
-             ALLOCATE( force%v_left(0:1,nzb+1:nzt+1,nysv:nyn) )
-             ALLOCATE( force%w_left(0:1,nzb+1:nzt,nys:nyn)    )
-             IF ( humidity )      ALLOCATE( force%q_left(0:1,nzb+1:nzt+1,nys:nyn)  )
-             IF ( .NOT. neutral ) ALLOCATE( force%pt_left(0:1,nzb+1:nzt+1,nys:nyn) )
-          ENDIF
-          IF ( force_bound_r )  THEN
-             ALLOCATE( force%u_right(0:1,nzb+1:nzt+1,nys:nyn)  )
-             ALLOCATE( force%v_right(0:1,nzb+1:nzt+1,nysv:nyn) )
-             ALLOCATE( force%w_right(0:1,nzb+1:nzt,nys:nyn)    )
-             IF ( humidity )      ALLOCATE( force%q_right(0:1,nzb+1:nzt+1,nys:nyn)  )
-             IF ( .NOT. neutral ) ALLOCATE( force%pt_right(0:1,nzb+1:nzt+1,nys:nyn) )
-          ENDIF
-          IF ( force_bound_n )  THEN
-             ALLOCATE( force%u_north(0:1,nzb+1:nzt+1,nxlu:nxr) )
-             ALLOCATE( force%v_north(0:1,nzb+1:nzt+1,nxl:nxr)  )
-             ALLOCATE( force%w_north(0:1,nzb+1:nzt,nxl:nxr)    )
-             IF ( humidity )      ALLOCATE( force%q_north(0:1,nzb+1:nzt+1,nxl:nxr)  )
-             IF ( .NOT. neutral ) ALLOCATE( force%pt_north(0:1,nzb+1:nzt+1,nxl:nxr) )
-          ENDIF
-          IF ( force_bound_s )  THEN
-             ALLOCATE( force%u_south(0:1,nzb+1:nzt+1,nxlu:nxr) )
-             ALLOCATE( force%v_south(0:1,nzb+1:nzt+1,nxl:nxr)  )
-             ALLOCATE( force%w_south(0:1,nzb+1:nzt,nxl:nxr)    )
-             IF ( humidity )      ALLOCATE( force%q_south(0:1,nzb+1:nzt+1,nxl:nxr)  )
-             IF ( .NOT. neutral ) ALLOCATE( force%pt_south(0:1,nzb+1:nzt+1,nxl:nxr) )
+          IF ( bc_dirichlet_l )  THEN
+             ALLOCATE( nest_offl%u_left(0:1,nzb+1:nzt,nys:nyn)  )
+             ALLOCATE( nest_offl%v_left(0:1,nzb+1:nzt,nysv:nyn) )
+             ALLOCATE( nest_offl%w_left(0:1,nzb+1:nzt-1,nys:nyn) )
+             IF ( humidity )      ALLOCATE( nest_offl%q_left(0:1,nzb+1:nzt,nys:nyn)  )
+             IF ( .NOT. neutral ) ALLOCATE( nest_offl%pt_left(0:1,nzb+1:nzt,nys:nyn) )
+          ENDIF
+          IF ( bc_dirichlet_r )  THEN
+             ALLOCATE( nest_offl%u_right(0:1,nzb+1:nzt,nys:nyn)  )
+             ALLOCATE( nest_offl%v_right(0:1,nzb+1:nzt,nysv:nyn) )
+             ALLOCATE( nest_offl%w_right(0:1,nzb+1:nzt-1,nys:nyn) )
+             IF ( humidity )      ALLOCATE( nest_offl%q_right(0:1,nzb+1:nzt,nys:nyn)  )
+             IF ( .NOT. neutral ) ALLOCATE( nest_offl%pt_right(0:1,nzb+1:nzt,nys:nyn) )
+          ENDIF
+          IF ( bc_dirichlet_n )  THEN
+             ALLOCATE( nest_offl%u_north(0:1,nzb+1:nzt,nxlu:nxr) )
+             ALLOCATE( nest_offl%v_north(0:1,nzb+1:nzt,nxl:nxr)  )
+             ALLOCATE( nest_offl%w_north(0:1,nzb+1:nzt-1,nxl:nxr) )
+             IF ( humidity )      ALLOCATE( nest_offl%q_north(0:1,nzb+1:nzt,nxl:nxr)  )
+             IF ( .NOT. neutral ) ALLOCATE( nest_offl%pt_north(0:1,nzb+1:nzt,nxl:nxr) )
+          ENDIF
+          IF ( bc_dirichlet_s )  THEN
+             ALLOCATE( nest_offl%u_south(0:1,nzb+1:nzt,nxlu:nxr) )
+             ALLOCATE( nest_offl%v_south(0:1,nzb+1:nzt,nxl:nxr)  )
+             ALLOCATE( nest_offl%w_south(0:1,nzb+1:nzt-1,nxl:nxr)    )
+             IF ( humidity )      ALLOCATE( nest_offl%q_south(0:1,nzb+1:nzt,nxl:nxr)  )
+             IF ( .NOT. neutral ) ALLOCATE( nest_offl%pt_south(0:1,nzb+1:nzt,nxl:nxr) )
           ENDIF
           
-          ALLOCATE( force%u_top(0:1,nys:nyn,nxlu:nxr) )
-          ALLOCATE( force%v_top(0:1,nysv:nyn,nxl:nxr) )
-          ALLOCATE( force%w_top(0:1,nys:nyn,nxl:nxr)  )
-          IF ( humidity )      ALLOCATE( force%q_top(0:1,nys:nyn,nxl:nxr)  )
-          IF ( .NOT. neutral ) ALLOCATE( force%pt_top(0:1,nys:nyn,nxl:nxr) )
-
-!
-!--       Initial call of input. Time array, initial 3D data of u, v, w, 
-!--       potential temperature, as well as mixing ratio, will be read. 
-!--       Moreover, data at lateral and top boundary will be read.  
+          ALLOCATE( nest_offl%u_top(0:1,nys:nyn,nxlu:nxr) )
+          ALLOCATE( nest_offl%v_top(0:1,nysv:nyn,nxl:nxr) )
+          ALLOCATE( nest_offl%w_top(0:1,nys:nyn,nxl:nxr)  )
+          IF ( humidity )      ALLOCATE( nest_offl%q_top(0:1,nys:nyn,nxl:nxr)  )
+          IF ( .NOT. neutral ) ALLOCATE( nest_offl%pt_top(0:1,nys:nyn,nxl:nxr) )
+
+!
+!--       Read COSMO data at lateral and top boundaries
           CALL netcdf_data_input_lsf
 !
-!--       Please note, at the moment INIFOR assumes only an equidistant vertical
-!--       grid. In case of vertical grid stretching, input of inital 3D data 
-!--       need to be inter- and/or extrapolated. 
-!--       Therefore, check if zw grid on file is identical to numeric zw grid.   
-          IF ( ANY( zu(1:nzt+1) /= force%zu_atmos(1:force%nzu) ) )  THEN
-!
-!--          Also data at the boundaries need to be inter/extrapolated at both
-!--          time levels
-             DO  t = 0, 1
-                IF ( force_bound_l )  THEN
-                   CALL netcdf_data_input_interpolate( force%u_left(t,:,:),    &
-                                                       zu(1:nzt+1),            &
-                                                       force%zu_atmos )
-                   CALL netcdf_data_input_interpolate( force%v_left(t,:,:),    &
-                                                       zu(1:nzt+1),            &
-                                                       force%zu_atmos )
-                   CALL netcdf_data_input_interpolate( force%w_left(t,:,:),    &
-                                                       zw(1:nzt+1),            &
-                                                       force%zw_atmos )
-                   IF ( .NOT. neutral )                                        &
-                      CALL netcdf_data_input_interpolate( force%pt_left(t,:,:),&
-                                                          zu(1:nzt+1),         &
-                                                          force%zu_atmos )
-                   IF ( humidity )                                             &
-                      CALL netcdf_data_input_interpolate( force%q_left(t,:,:), &
-                                                          zu(1:nzt+1),         &
-                                                          force%zu_atmos )
-                ENDIF
-                IF ( force_bound_r )  THEN
-                   CALL netcdf_data_input_interpolate( force%u_right(t,:,:),   &
-                                                       zu(1:nzt+1),            &
-                                                       force%zu_atmos )
-                   CALL netcdf_data_input_interpolate( force%v_right(t,:,:),   &
-                                                       zu(1:nzt+1),            &
-                                                       force%zu_atmos )
-                   CALL netcdf_data_input_interpolate( force%w_right(t,:,:),   &
-                                                       zw(1:nzt+1),            &
-                                                       force%zw_atmos )
-                   IF ( .NOT. neutral )                                        &
-                      CALL netcdf_data_input_interpolate( force%pt_right(t,:,:),&
-                                                          zu(1:nzt+1),          &
-                                                          force%zu_atmos )
-                   IF ( humidity )                                             &
-                      CALL netcdf_data_input_interpolate( force%q_right(t,:,:),&
-                                                          zu(1:nzt+1),         &
-                                                          force%zu_atmos )
-                ENDIF
-                IF ( force_bound_n )  THEN
-                   CALL netcdf_data_input_interpolate( force%u_north(t,:,:),   &
-                                                       zu(1:nzt+1),            &
-                                                       force%zu_atmos )
-                   CALL netcdf_data_input_interpolate( force%v_north(t,:,:),   &
-                                                       zu(1:nzt+1),            &
-                                                       force%zu_atmos )
-                   CALL netcdf_data_input_interpolate( force%w_north(t,:,:),   &
-                                                       zw(1:nzt+1),            &
-                                                       force%zw_atmos )
-                   IF ( .NOT. neutral )                                        &
-                      CALL netcdf_data_input_interpolate( force%pt_north(t,:,:),&
-                                                          zu(1:nzt+1),          &
-                                                          force%zu_atmos )
-                   IF ( humidity )                                             &
-                      CALL netcdf_data_input_interpolate( force%q_north(t,:,:),&
-                                                          zu(1:nzt+1),         &
-                                                          force%zu_atmos )
-                ENDIF
-                IF ( force_bound_s )  THEN
-                   CALL netcdf_data_input_interpolate( force%u_south(t,:,:),   &
-                                                       zu(1:nzt+1),            &
-                                                       force%zu_atmos )
-                   CALL netcdf_data_input_interpolate( force%v_south(t,:,:),   &
-                                                       zu(1:nzt+1),            &
-                                                       force%zu_atmos )
-                   CALL netcdf_data_input_interpolate( force%w_south(t,:,:),   &
-                                                       zw(1:nzt+1),            &
-                                                       force%zw_atmos )
-                   IF ( .NOT. neutral )                                        &
-                      CALL netcdf_data_input_interpolate( force%pt_south(t,:,:),&
-                                                          zu(1:nzt+1),          &
-                                                          force%zu_atmos )
-                   IF ( humidity )                                             &
-                      CALL netcdf_data_input_interpolate( force%q_south(t,:,:),&
-                                                          zu(1:nzt+1),         &
-                                                          force%zu_atmos )
-                ENDIF
-             ENDDO
-          ENDIF
-
-!
-!--       Exchange ghost points
-          CALL exchange_horiz( u, nbgp )
-          CALL exchange_horiz( v, nbgp )
-          CALL exchange_horiz( w, nbgp )
-          IF ( .NOT. neutral )  CALL exchange_horiz( pt, nbgp )
-          IF ( humidity      )  CALL exchange_horiz( q,  nbgp )
-!
-!--       At lateral boundaries, set also initial boundary conditions
-          IF ( force_bound_l )  THEN
-             u(:,:,nxl)   = u(:,:,nxlu)
-             v(:,:,nxl-1) = v(:,:,nxl)
-             w(:,:,nxl-1) = w(:,:,nxl)
-             IF ( .NOT. neutral )  pt(:,:,nxl-1) = pt(:,:,nxl)
-             IF ( humidity      )  q(:,:,nxl-1)  = q(:,:,nxl)
-          ENDIF
-          IF ( force_bound_r )  THEN
-             u(:,:,nxr+1) = u(:,:,nxr)
-             v(:,:,nxr+1) = v(:,:,nxr)
-             w(:,:,nxr+1) = w(:,:,nxr)
-             IF ( .NOT. neutral )  pt(:,:,nxr+1) = pt(:,:,nxr)
-             IF ( humidity      )  q(:,:,nxr+1)  = q(:,:,nxr)
-          ENDIF
-          IF ( force_bound_s )  THEN
-             u(:,nys-1,:) = u(:,nys,:)
-             v(:,nys,:)   = v(:,nysv,:)
-             w(:,nys-1,:) = w(:,nys,:)
-             IF ( .NOT. neutral )  pt(:,nys-1,:) = pt(:,nys,:)
-             IF ( humidity      )  q(:,nys-1,:)  = q(:,nys,:)
-          ENDIF
-          IF ( force_bound_n )  THEN
-             u(:,nyn+1,:) = u(:,nyn,:)
-             v(:,nyn+1,:) = v(:,nyn,:)
-             w(:,nyn+1,:) = w(:,nyn,:)
-             IF ( .NOT. neutral )  pt(:,nyn+1,:) = pt(:,nyn,:)
-             IF ( humidity      )  q(:,nyn+1,:)  = q(:,nyn,:)
-          ENDIF
-
+!--       Write COSMO data at lateral and top boundaries
+          CALL lsf_nesting_offline
 !
 !--       After 3D data is initialized, ensure mass conservation
-          CALL forcing_bc_mass_conservation
+          CALL lsf_nesting_offline_mass_conservation
 !
 !--       Initialize surface pressure. Please note, time-dependent surface
@@ -1228 +1023 @@
 !--       treatment of fluxes. 
 !--       For the moment, comment this out! 
-!         surface_pressure = force%surface_pressure(0)
+!         surface_pressure = nest_offl%surface_pressure(0)
 
        ELSE

palm/trunk/SOURCE/modules.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-! 
+! Rename offline nesting variables:
+! -inflow_l, inflow_n, inflow_r, inflow_s, 
+!  nest_bound_l, nest_bound_n, nest_bound_r, nest_bound_s, nest_domain, forcing, 
+!  force_bound_l, force_bound_n, force_bound_r, force_bound_s, outflow_l, 
+!  outflow_n, outflow_r, outflow_s
+! +bc_dirichlet_l, bc_dirichlet_n, bc_dirichlet_n, bc_dirichlet_r, 
+!  bc_radiation_l, bc_radiation_n, bc_radiation_n, bc_radiation_r, child_domain 
+!  nesting_offline
 ! 
 ! Former revisions:
@@ -1272 +1279 @@
     LOGICAL ::  agent_time_unlimited = .FALSE.                   !< namelist parameter
     LOGICAL ::  air_chemistry = .FALSE.                          !< chemistry model switch
+    LOGICAL ::  bc_dirichlet_l                                   !< flag indicating dirichlet boundary condition on left model boundary
+    LOGICAL ::  bc_dirichlet_n                                   !< flag indicating dirichlet boundary condition on north model boundary
+    LOGICAL ::  bc_dirichlet_r                                   !< flag indicating dirichlet boundary condition on right model boundary
+    LOGICAL ::  bc_dirichlet_s                                   !< flag indicating dirichlet boundary condition on south model boundary
     LOGICAL ::  bc_lr_cyc =.TRUE.                                !< left-right boundary condition cyclic?
     LOGICAL ::  bc_lr_dirrad = .FALSE.                           !< left-right boundary condition dirichlet/radiation?
@@ -1278 +1289 @@
     LOGICAL ::  bc_ns_dirrad = .FALSE.                           !< north-south boundary condition dirichlet/radiation?
     LOGICAL ::  bc_ns_raddir = .FALSE.                           !< north-south boundary condition radiation/dirichlet?
+    LOGICAL ::  bc_radiation_l = .FALSE.                         !< radiation boundary condition for outflow at left domain boundary
+    LOGICAL ::  bc_radiation_n = .FALSE.                         !< radiation boundary condition for outflow at north domain boundary
+    LOGICAL ::  bc_radiation_r = .FALSE.                         !< radiation boundary condition for outflow at right domain boundary
+    LOGICAL ::  bc_radiation_s = .FALSE.                         !< radiation boundary condition for outflow at south domain boundary
     LOGICAL ::  calc_soil_moisture_during_spinup = .FALSE.       !< namelist parameter
     LOGICAL ::  call_microphysics_at_all_substeps = .FALSE.      !< namelist parameter
     LOGICAL ::  call_psolver_at_all_substeps = .TRUE.            !< namelist parameter
+    LOGICAL ::  child_domain  = .FALSE.                          !< flag indicating that model is nested in a parent domain
     LOGICAL ::  cloud_droplets = .FALSE.                         !< namelist parameter
     LOGICAL ::  cloud_physics = .FALSE.                          !< namelist parameter
@@ -1310 +1326 @@
     LOGICAL ::  first_call_mas = .TRUE.                          !< call mas only once per timestep
     LOGICAL ::  force_print_header = .FALSE.                     !< namelist parameter
-    LOGICAL ::  force_bound_l = .FALSE.                          !< flag indicating domain boundary on left side to set forcing boundary conditions
-    LOGICAL ::  force_bound_n = .FALSE.                          !< flag indicating domain boundary on north side to set forcing boundary conditions
-    LOGICAL ::  force_bound_r = .FALSE.                          !< flag indicating domain boundary on right side to set forcing boundary conditions
-    LOGICAL ::  force_bound_s = .FALSE.                          !< flag indicating domain boundary on south side to set forcing boundary conditions
-    LOGICAL ::  forcing = .FALSE.                                !< flag controlling forcing from large-scale model     
     LOGICAL ::  galilei_transformation = .FALSE.                 !< namelist parameter
     LOGICAL ::  humidity = .FALSE.                               !< namelist parameter
     LOGICAL ::  humidity_remote = .FALSE.                        !< switch for receiving near-surface humidity flux (atmosphere-ocean coupling)
-    LOGICAL ::  inflow_l = .FALSE.                               !< left domain boundary has non-cyclic inflow?
-    LOGICAL ::  inflow_n = .FALSE.                               !< north domain boundary has non-cyclic inflow?
-    LOGICAL ::  inflow_r = .FALSE.                               !< right domain boundary has non-cyclic inflow?
-    LOGICAL ::  inflow_s = .FALSE.                               !< south domain boundary has non-cyclic inflow?
     LOGICAL ::  large_scale_forcing = .FALSE.                    !< namelist parameter
     LOGICAL ::  large_scale_subsidence = .FALSE.                 !< namelist parameter
@@ -1336 +1343 @@
     LOGICAL ::  microphysics_seifert = .FALSE.                   !< use 2-moment Seifert and Beheng scheme
     LOGICAL ::  mg_switch_to_pe0 = .FALSE.                       !< internal multigrid switch for steering the ghost point exchange in case that data has been collected on PE0
-    LOGICAL ::  nest_bound_l = .FALSE.                           !< flag indicating nested domain boundary on left side
-    LOGICAL ::  nest_bound_n = .FALSE.                           !< flag indicating nested domain boundary on north side
-    LOGICAL ::  nest_bound_r = .FALSE.                           !< flag indicating nested domain boundary on right side
-    LOGICAL ::  nest_bound_s = .FALSE.                           !< flag indicating nested domain boundary on south side
-    LOGICAL ::  nest_domain  = .FALSE.                           !< domain is nested into a parent domain?
+    LOGICAL ::  nesting_offline = .FALSE.                        !< flag controlling offline nesting in COSMO model  
     LOGICAL ::  neutral = .FALSE.                                !< namelist parameter
     LOGICAL ::  nudging = .FALSE.                                !< namelist parameter
     LOGICAL ::  ocean = .FALSE.                                  !< namelist parameter
-    LOGICAL ::  outflow_l = .FALSE.                              !< left domain boundary has non-cyclic outflow?
-    LOGICAL ::  outflow_n = .FALSE.                              !< north domain boundary has non-cyclic outflow?
-    LOGICAL ::  outflow_r = .FALSE.                              !< right domain boundary has non-cyclic outflow?
-    LOGICAL ::  outflow_s = .FALSE.                              !< south domain boundary has non-cyclic outflow?
     LOGICAL ::  passive_scalar = .FALSE.                         !< namelist parameter
     LOGICAL ::  plant_canopy = .FALSE.                           !< switch for use of plant canopy model

palm/trunk/SOURCE/netcdf_data_input_mod.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! Adjust input of dynamic driver according to revised Inifor version. 
+! Replace simulated_time by time_since_reference_point.
+! Rename variables in mesoscale-offline nesting mode.
 ! 
 ! Former revisions:
@@ -199 +201 @@
 !-- Define data type for nesting in larger-scale models like COSMO.
 !-- Data type comprises u, v, w, pt, and q at lateral and top boundaries.
-    TYPE force_type
+    TYPE nest_offl_type
 
        CHARACTER(LEN=100), DIMENSION(:), ALLOCATABLE ::  var_names
@@ -206 +208 @@
        INTEGER(iwp) ::  nzu    !< number of vertical levels on scalar grid in dynamic input file
        INTEGER(iwp) ::  nzw    !< number of vertical levels on w grid in dynamic input file
-       INTEGER(iwp) ::  tind   !< time index for reference time in large-scale forcing data
-       INTEGER(iwp) ::  tind_p !< time index for following time in large-scale forcing data
+       INTEGER(iwp) ::  tind   !< time index for reference time in mesoscale-offline nesting
+       INTEGER(iwp) ::  tind_p !< time index for following time in mesoscale-offline nesting
 
        LOGICAL      ::  init         = .FALSE.
-       LOGICAL      ::  interpolated = .FALSE.
        LOGICAL      ::  from_file    = .FALSE.
 
@@ -251 +252 @@
        REAL(wp), DIMENSION(:,:,:), ALLOCATABLE ::  pt_top   !< potentital temperautre at top boundary
 
-    END TYPE force_type
+    END TYPE nest_offl_type
 
     TYPE init_type
@@ -474 +475 @@
     TYPE(dims_xy)    ::  dim_static  !< data structure for x, y-dimension in static input file
 
-    TYPE(force_type) ::  force     !< data structure for data input at lateral and top boundaries (provided by Inifor)
+    TYPE(nest_offl_type) ::  nest_offl  !< data structure for data input at lateral and top boundaries (provided by Inifor)  
 
     TYPE(init_type) ::  init_3d    !< data structure for the initialization of the 3D flow and soil fields
@@ -606 +607 @@
 !-- Public variables
     PUBLIC albedo_pars_f, albedo_type_f, basal_area_density_f, buildings_f,    &
-           building_id_f, building_pars_f, building_type_f, force, init_3d,    &
+           building_id_f, building_pars_f, building_type_f, init_3d,           &
            init_model, input_file_static, input_pids_static,                   &
-           input_pids_dynamic, leaf_area_density_f,                            &
+           input_pids_dynamic, leaf_area_density_f, nest_offl,                 &
            pavement_pars_f, pavement_subsurface_pars_f, pavement_type_f,       &
            root_area_density_lad_f, root_area_density_lsm_f, soil_pars_f,      &
@@ -719 +720 @@
 !--    sun-zenith angles. To avoid this, longitude and latitude in each model
 !--    domain will be set to the values of the root model. Please note, this
-!--    synchronization is required already here.
+!--    synchronization is required already here. 
 #if defined( __parallel )
        CALL MPI_BCAST( init_model%latitude,  1, MPI_REAL, 0,                   &
@@ -726 +727 @@
                        MPI_COMM_WORLD, ierr )
 #endif
-
 
     END SUBROUTINE netcdf_data_input_init
@@ -2057 +2057 @@
 
        USE arrays_3d,                                                          &
-           ONLY:  q, pt, u, v, w
+           ONLY:  q, pt, u, v, w, zu, zw
 
        USE control_parameters,                                                 &
-           ONLY:  bc_lr_cyc, bc_ns_cyc, forcing, humidity, land_surface,       &
-                  message_string, neutral, surface_pressure
+           ONLY:  bc_lr_cyc, bc_ns_cyc, humidity, land_surface, message_string,&
+                  nesting_offline, neutral, surface_pressure
 
        USE indices,                                                            &
@@ -2119 +2119 @@
        CALL get_dimension_length( id_dynamic, init_3d%nzu, 'z'     )
        CALL get_dimension_length( id_dynamic, init_3d%nzw, 'zw'    )
-       CALL get_dimension_length( id_dynamic, init_3d%nzs, 'depth' )
+       CALL get_dimension_length( id_dynamic, init_3d%nzs, 'zsoil' )
 !
 !--    Read also the horizontal dimensions. These are used just used fo
@@ -2142 +2142 @@
        ENDIF
 
-       IF ( init_3d%nzu-1 /= nz )  THEN
+       IF ( init_3d%nzu /= nz )  THEN
           message_string = 'Number of inifor vertical grid points ' //         &
                            'does not match the number of numeric grid '//      &
@@ -2159 +2159 @@
           CALL get_variable( id_dynamic, 'zw', init_3d%zw_atmos )
        ENDIF
-       IF ( check_existence( var_names, 'depth' ) )  THEN
+       IF ( check_existence( var_names, 'zsoil' ) )  THEN
           ALLOCATE( init_3d%z_soil(1:init_3d%nzs) )
-          CALL get_variable( id_dynamic, 'depth', init_3d%z_soil )
-       ENDIF
-!
-!--    Read initial geostrophic wind components at t = 0 (index 1 in file).
-!        IF ( check_existence( var_names, 'tend_ug' ) )  THEN
+          CALL get_variable( id_dynamic, 'zsoil', init_3d%z_soil )
+       ENDIF
+!
+!--    Check for consistency between vertical coordinates in dynamic 
+!--    driver and numeric grid. 
+!--    Please note, depending on compiler options both may be 
+!--    equal up to a certain threshold, and differences between
+!--    the numeric grid and vertical coordinate in the driver can built-
+!--    up to 10E-1-10E-0 m. For this reason, the check is performed not
+!--    for exactly matching values. 
+       IF ( ANY( ABS( zu(1:nzt)   - init_3d%zu_atmos(1:init_3d%nzu) )    &
+                      > 10E-1 )  .OR.                                    &
+            ANY( ABS( zw(1:nzt-1) - init_3d%zw_atmos(1:init_3d%nzw) )    &
+                      > 10E-1 ) )  THEN
+          message_string = 'Vertical grid in dynamic driver does not '// &
+                           'match the numeric grid.'
+          CALL message( 'netcdf_data_input_mod', 'NDI003', 1, 2, 0, 6, 0 )
+       ENDIF
+!
+!--    Read initial geostrophic wind components at 
+!--    t = 0 (index 1 in file). 
        IF ( check_existence( var_names, 'ls_forcing_ug' ) )  THEN
           ALLOCATE( init_3d%ug_init(nzb:nzt+1) )
-!           CALL get_variable_pr( id_dynamic, 'tend_ug', 1,                      &
-!                                 init_3d%ug_init )
-          CALL get_variable_pr( id_dynamic, 'ls_forcing_ug', 1,                &
-                                init_3d%ug_init )
+          init_3d%ug_init = 0.0_wp
+
+          CALL get_variable_pr( id_dynamic, 'ls_forcing_ug', 1,          &
+                                init_3d%ug_init(1:nzt) )
+!
+!--       Set top-boundary condition (Neumann)
+          init_3d%ug_init(nzt+1) = init_3d%ug_init(nzt)
+
           init_3d%from_file_ug = .TRUE.
        ELSE
           init_3d%from_file_ug = .FALSE.
        ENDIF
-!        IF ( check_existence( var_names, 'tend_vg' ) )  THEN
        IF ( check_existence( var_names, 'ls_forcing_vg' ) )  THEN
           ALLOCATE( init_3d%vg_init(nzb:nzt+1) )
-!           CALL get_variable_pr( id_dynamic, 'tend_vg', 1,                      &
-!                                 init_3d%vg_init )
-          CALL get_variable_pr( id_dynamic, 'ls_forcing_vg', 1,                &
-                                init_3d%vg_init )
+          init_3d%vg_init = 0.0_wp
+
+          CALL get_variable_pr( id_dynamic, 'ls_forcing_vg', 1,          &
+                                init_3d%vg_init(1:nzt) )
+!
+!--       Set top-boundary condition (Neumann)
+          init_3d%vg_init(nzt+1) = init_3d%vg_init(nzt)
+
           init_3d%from_file_vg = .TRUE.
        ELSE
@@ -2195 +2218 @@
 !--    into separate loops.  
 !--    Read u-component
-       IF ( check_existence( var_names, 'init_u' ) )  THEN
+       IF ( check_existence( var_names, 'init_atmosphere_u' ) )  THEN
 !
 !--       Read attributes for the fill value and level-of-detail
           CALL get_attribute( id_dynamic, char_fill, init_3d%fill_u,           &
-                              .FALSE., 'init_u' )
+                              .FALSE., 'init_atmosphere_u' )
           CALL get_attribute( id_dynamic, char_lod, init_3d%lod_u,             &
-                              .FALSE., 'init_u' )
+                              .FALSE., 'init_atmosphere_u' )
 !
 !--       level-of-detail 1 - read initialization profile
@@ -2208 +2231 @@
              init_3d%u_init = 0.0_wp
 
-             CALL get_variable( id_dynamic, 'init_u',                          &
-                                init_3d%u_init(nzb+1:nzt+1) )
+             CALL get_variable( id_dynamic, 'init_atmosphere_u',               &
+                                init_3d%u_init(nzb+1:nzt) )
+!
+!--          Set top-boundary condition (Neumann)
+             init_3d%u_init(nzt+1) = init_3d%u_init(nzt)
 !
 !--       level-of-detail 2 - read 3D initialization data
           ELSEIF ( init_3d%lod_u == 2 )  THEN
-
-             CALL get_variable( id_dynamic, 'init_u',                          &
-                                u(nzb+1:nzt+1,nys:nyn,nxlu:nxr),               &
+             CALL get_variable( id_dynamic, 'init_atmosphere_u',               &
+                                u(nzb+1:nzt,nys:nyn,nxlu:nxr),                 &
                                 nxlu, nys+1, nzb+1,                            &
                                 nxr-nxlu+1, nyn-nys+1, init_3d%nzu,            &
                                 dynamic_3d )
+!
+!--          Set value at leftmost model grid point nxl = 0. This is because
+!--          Inifor provides data only from 1:nx-1 since it assumes non-cyclic
+!--          conditions.
+             IF ( nxl == 0 )                                                   &
+                u(nzb+1:nzt,nys:nyn,nxl) = u(nzb+1:nzt,nys:nyn,nxlu)
+!
+!--          Set bottom and top-boundary
+             u(nzb,:,:)   = u(nzb+1,:,:)
+             u(nzt+1,:,:) = u(nzt,:,:)
+             
           ENDIF
           init_3d%from_file_u = .TRUE.
@@ -2224 +2260 @@
 !
 !--    Read v-component
-       IF ( check_existence( var_names, 'init_v' ) )  THEN
+       IF ( check_existence( var_names, 'init_atmosphere_v' ) )  THEN
 !
 !--       Read attributes for the fill value and level-of-detail
           CALL get_attribute( id_dynamic, char_fill, init_3d%fill_v,           &
-                              .FALSE., 'init_v' )
+                              .FALSE., 'init_atmosphere_v' )
           CALL get_attribute( id_dynamic, char_lod, init_3d%lod_v,             &
-                              .FALSE., 'init_v' )
+                              .FALSE., 'init_atmosphere_v' )
 !
 !--       level-of-detail 1 - read initialization profile
@@ -2237 +2273 @@
              init_3d%v_init = 0.0_wp
 
-             CALL get_variable( id_dynamic, 'init_v',                          &
-                                init_3d%v_init(nzb+1:nzt+1) )
-
+             CALL get_variable( id_dynamic, 'init_atmosphere_v',               &
+                                init_3d%v_init(nzb+1:nzt) )
+!
+!--          Set top-boundary condition (Neumann)
+             init_3d%v_init(nzt+1) = init_3d%v_init(nzt)
 !
 !--       level-of-detail 2 - read 3D initialization data
           ELSEIF ( init_3d%lod_v == 2 )  THEN
-
-             CALL get_variable( id_dynamic, 'init_v',                          &
-                                v(nzb+1:nzt+1,nysv:nyn,nxl:nxr),               &
+          
+             CALL get_variable( id_dynamic, 'init_atmosphere_v',               &
+                                v(nzb+1:nzt,nysv:nyn,nxl:nxr),                 &
                                 nxl+1, nysv, nzb+1,                            &
                                 nxr-nxl+1, nyn-nysv+1, init_3d%nzu,            &
                                 dynamic_3d )
+!
+!--          Set value at southmost model grid point nys = 0. This is because
+!--          Inifor provides data only from 1:ny-1 since it assumes non-cyclic
+!--          conditions.
+             IF ( nys == 0 )                                                   &
+                v(nzb+1:nzt,nys,nxl:nxr) = v(nzb+1:nzt,nysv,nxl:nxr)                                
+!
+!--          Set bottom and top-boundary
+             v(nzb,:,:)   = v(nzb+1,:,:)
+             v(nzt+1,:,:) = v(nzt,:,:)
              
           ENDIF
@@ -2255 +2303 @@
 !
 !--    Read w-component
-       IF ( check_existence( var_names, 'init_w' ) )  THEN
+       IF ( check_existence( var_names, 'init_atmosphere_w' ) )  THEN
 !
 !--       Read attributes for the fill value and level-of-detail
           CALL get_attribute( id_dynamic, char_fill, init_3d%fill_w,           &
-                              .FALSE., 'init_w' )
+                              .FALSE., 'init_atmosphere_w' )
           CALL get_attribute( id_dynamic, char_lod, init_3d%lod_w,             &
-                              .FALSE., 'init_w' )
+                              .FALSE., 'init_atmosphere_w' )
 !
 !--       level-of-detail 1 - read initialization profile
@@ -2268 +2316 @@
              init_3d%w_init = 0.0_wp
 
-             CALL get_variable( id_dynamic, 'init_w',                          &
-                                init_3d%w_init(nzb+1:nzt) )
-
+             CALL get_variable( id_dynamic, 'init_atmosphere_w',               &
+                                init_3d%w_init(nzb+1:nzt-1) )
+!
+!--          Set top-boundary condition (Neumann)
+             init_3d%w_init(nzt:nzt+1) = init_3d%w_init(nzt-1)
 !
 !--       level-of-detail 2 - read 3D initialization data
           ELSEIF ( init_3d%lod_w == 2 )  THEN
 
-             CALL get_variable( id_dynamic, 'init_w',                           &
-                                w(nzb+1:nzt,nys:nyn,nxl:nxr),                   &
+             CALL get_variable( id_dynamic, 'init_atmosphere_w',                &
+                                w(nzb+1:nzt-1,nys:nyn,nxl:nxr),                 &
                                 nxl+1, nys+1, nzb+1,                            &
                                 nxr-nxl+1, nyn-nys+1, init_3d%nzw,              &
                                 dynamic_3d )
+!
+!--          Set bottom and top-boundary                                
+             w(nzb,:,:)   = 0.0_wp  
+             w(nzt,:,:)   = w(nzt-1,:,:)
+             w(nzt+1,:,:) = w(nzt-1,:,:)
 
           ENDIF
@@ -2287 +2342 @@
 !--    Read potential temperature
        IF ( .NOT. neutral )  THEN
-          IF ( check_existence( var_names, 'init_pt' ) )  THEN
+          IF ( check_existence( var_names, 'init_atmosphere_pt' ) )  THEN
 !
 !--          Read attributes for the fill value and level-of-detail
              CALL get_attribute( id_dynamic, char_fill, init_3d%fill_pt,       &
-                                 .FALSE., 'init_pt' )
+                                 .FALSE., 'init_atmosphere_pt' )
              CALL get_attribute( id_dynamic, char_lod, init_3d%lod_pt,         &
-                                 .FALSE., 'init_pt' )
+                                 .FALSE., 'init_atmosphere_pt' )
 !
 !--          level-of-detail 1 - read initialization profile
@@ -2299 +2354 @@
                 ALLOCATE( init_3d%pt_init(nzb:nzt+1) )
 
-                CALL get_variable( id_dynamic, 'init_pt',                      &
-                                   init_3d%pt_init(nzb+1:nzt+1) )
-!
-!--             Set Neumann surface boundary condition for initial profil
-                init_3d%pt_init(nzb) = init_3d%pt_init(nzb+1)
+                CALL get_variable( id_dynamic, 'init_atmosphere_pt',           &
+                                   init_3d%pt_init(nzb+1:nzt) )
+!
+!--             Set Neumann top and surface boundary condition for initial 
+!--             profil
+                init_3d%pt_init(nzb)   = init_3d%pt_init(nzb+1)
+                init_3d%pt_init(nzt+1) = init_3d%pt_init(nzt)
 !
 !--          level-of-detail 2 - read 3D initialization data
              ELSEIF ( init_3d%lod_pt == 2 )  THEN
 
-                CALL get_variable( id_dynamic, 'init_pt',                      &
-                                   pt(nzb+1:nzt+1,nys:nyn,nxl:nxr),            &
+                CALL get_variable( id_dynamic, 'init_atmosphere_pt',           &
+                                   pt(nzb+1:nzt,nys:nyn,nxl:nxr),              &
                                    nxl+1, nys+1, nzb+1,                        &
                                    nxr-nxl+1, nyn-nys+1, init_3d%nzu,          &
                                    dynamic_3d )
-
+                                   
+!
+!--             Set bottom and top-boundary
+                pt(nzb,:,:)   = pt(nzb+1,:,:)
+                pt(nzt+1,:,:) = pt(nzt,:,:)             
 
              ENDIF
@@ -2322 +2383 @@
 !--    Read mixing ratio
        IF ( humidity )  THEN
-          IF ( check_existence( var_names, 'init_qv' ) )  THEN
+          IF ( check_existence( var_names, 'init_atmosphere_qv' ) )  THEN
 !
 !--          Read attributes for the fill value and level-of-detail
              CALL get_attribute( id_dynamic, char_fill, init_3d%fill_q,        &
-                                 .FALSE., 'init_qv' )
+                                 .FALSE., 'init_atmosphere_qv' )
              CALL get_attribute( id_dynamic, char_lod, init_3d%lod_q,          &
-                                 .FALSE., 'init_qv' )
+                                 .FALSE., 'init_atmosphere_qv' )
 !
 !--          level-of-detail 1 - read initialization profile
@@ -2334 +2395 @@
                 ALLOCATE( init_3d%q_init(nzb:nzt+1) )
 
-                CALL get_variable( id_dynamic, 'init_qv',                      &
-                                   init_3d%q_init(nzb+1:nzt+1) )
-!
-!--             Set Neumann surface boundary condition for initial profil
-                init_3d%q_init(nzb) = init_3d%q_init(nzb+1)
-
+                CALL get_variable( id_dynamic, 'init_atmosphere_qv',           &
+                                    init_3d%q_init(nzb+1:nzt) )
+!
+!--             Set bottom and top boundary condition (Neumann)
+                init_3d%q_init(nzb)   = init_3d%q_init(nzb+1)
+                init_3d%q_init(nzt+1) = init_3d%q_init(nzt)
 !
 !--          level-of-detail 2 - read 3D initialization data
              ELSEIF ( init_3d%lod_q == 2 )  THEN
              
-                CALL get_variable( id_dynamic, 'init_qv',                      &
-                                   q(nzb+1:nzt+1,nys:nyn,nxl:nxr),             &
+                CALL get_variable( id_dynamic, 'init_atmosphere_qv',           &
+                                   q(nzb+1:nzt,nys:nyn,nxl:nxr),               &
                                    nxl+1, nys+1, nzb+1,                        &
                                    nxr-nxl+1, nyn-nys+1, init_3d%nzu,          &
                                    dynamic_3d )
-
-
-
+                                   
+!
+!--             Set bottom and top-boundary
+                q(nzb,:,:)   = q(nzb+1,:,:)
+                q(nzt+1,:,:) = q(nzt,:,:)
+                
              ENDIF
              init_3d%from_file_q = .TRUE.
@@ -2382 +2446 @@
 
                CALL get_variable( id_dynamic, 'init_soil_m',                   &   
-                                  init_3d%msoil(0:init_3d%nzs-1,nys:nyn,nxl:nxr),&
-                                  nxl, nxr, nys, nyn, 0, init_3d%nzs-1 )
+                                init_3d%msoil(0:init_3d%nzs-1,nys:nyn,nxl:nxr),&
+                                nxl, nxr, nys, nyn, 0, init_3d%nzs-1 )
 
              ENDIF
@@ -2413 +2477 @@
                 
                 CALL get_variable( id_dynamic, 'init_soil_t',                  &   
-                                  init_3d%tsoil(0:init_3d%nzs-1,nys:nyn,nxl:nxr),&
-                                  nxl, nxr, nys, nyn, 0, init_3d%nzs-1 )
+                                init_3d%tsoil(0:init_3d%nzs-1,nys:nyn,nxl:nxr),&
+                                nxl, nxr, nys, nyn, 0, init_3d%nzs-1 )
              ENDIF
              init_3d%from_file_tsoil = .TRUE.
@@ -2439 +2503 @@
           ENDIF
           IF ( .NOT. check_passed )  THEN
-             message_string = 'NetCDF input for u_init must not contain ' //   &
-                              'any _FillValues'
+             message_string = 'NetCDF input for init_atmosphere_u must ' //    &
+                              'not contain any _FillValues'
              CALL message( 'netcdf_data_input_mod', 'NDI004', 2, 2, 0, 6, 0 )
           ENDIF
@@ -2455 +2519 @@
           ENDIF
           IF ( .NOT. check_passed )  THEN
-             message_string = 'NetCDF input for v_init must not contain ' //   &
-                              'any _FillValues'
+             message_string = 'NetCDF input for init_atmosphere_v must ' //    &
+                              'not contain any _FillValues'
              CALL message( 'netcdf_data_input_mod', 'NDI005', 2, 2, 0, 6, 0 )
           ENDIF
@@ -2471 +2535 @@
           ENDIF
           IF ( .NOT. check_passed )  THEN
-             message_string = 'NetCDF input for w_init must not contain ' //   &
-                              'any _FillValues'
+             message_string = 'NetCDF input for init_atmosphere_w must ' //    &
+                              'not contain any _FillValues'
              CALL message( 'netcdf_data_input_mod', 'NDI006', 2, 2, 0, 6, 0 )
           ENDIF
@@ -2487 +2551 @@
           ENDIF
           IF ( .NOT. check_passed )  THEN
-             message_string = 'NetCDF input for pt_init must not contain ' //  &
-                              'any _FillValues'
+             message_string = 'NetCDF input for init_atmosphere_pt must ' //   &
+                              'not contain any _FillValues'
              CALL message( 'netcdf_data_input_mod', 'NDI007', 2, 2, 0, 6, 0 )
           ENDIF
@@ -2503 +2567 @@
           ENDIF
           IF ( .NOT. check_passed )  THEN
-             message_string = 'NetCDF input for q_init must not contain ' //   &
-                              'any _FillValues'
+             message_string = 'NetCDF input for init_atmosphere_q must ' //    &
+                              'not contain any _FillValues'
              CALL message( 'netcdf_data_input_mod', 'NDI008', 2, 2, 0, 6, 0 )
           ENDIF
@@ -2524 +2588 @@
 
        USE control_parameters,                                                 &
-           ONLY:  bc_lr_cyc, bc_ns_cyc, force_bound_l, force_bound_n,          &
-                  force_bound_r, force_bound_s,                                &
-                  forcing, humidity, message_string, neutral, simulated_time
-
+           ONLY:  bc_dirichlet_l, bc_dirichlet_n, bc_dirichlet_r,              &
+                  bc_dirichlet_s, bc_lr_cyc, bc_ns_cyc, humidity,              &
+                  message_string, neutral, nesting_offline,                    &
+                  time_since_reference_point
 
        USE indices,                                                            &
@@ -2533 +2597 @@
 
        IMPLICIT NONE
-
-       LOGICAL      ::  dynamic_3d = .TRUE. !< flag indicating that 3D data is read from dynamic file
        
        INTEGER(iwp) ::  i          !< running index along x-direction
@@ -2545 +2607 @@
        REAL(wp) ::  dum           !< dummy variable to skip columns while reading topography file
 
-       force%from_file = MERGE( .TRUE., .FALSE., input_pids_dynamic )
+       nest_offl%from_file = MERGE( .TRUE., .FALSE., input_pids_dynamic ) 
 !
 !--    Skip input if no forcing from larger-scale models is applied.
-       IF ( .NOT. forcing )  RETURN
+       IF ( .NOT. nesting_offline )  RETURN
 
 !
@@ -2560 +2622 @@
                             TRIM( coupling_char ), id_dynamic )
 !
-!--    Initialize INIFOR forcing.
-       IF ( .NOT. force%init )  THEN
+!--    Initialize INIFOR forcing. 
+       IF ( .NOT. nest_offl%init )  THEN
 !
 !--       At first, inquire all variable names.
@@ -2567 +2629 @@
 !
 !--       Allocate memory to store variable names.
-          ALLOCATE( force%var_names(1:num_vars) )
-          CALL inquire_variable_names( id_dynamic, force%var_names )
+          ALLOCATE( nest_offl%var_names(1:num_vars) )
+          CALL inquire_variable_names( id_dynamic, nest_offl%var_names )
 !
 !--       Read time dimension, allocate memory and finally read time array
-          CALL get_dimension_length( id_dynamic, force%nt, 'time' )
-
-          IF ( check_existence( force%var_names, 'time' ) )  THEN
-             ALLOCATE( force%time(0:force%nt-1) )
-             CALL get_variable( id_dynamic, 'time', force%time )
+          CALL get_dimension_length( id_dynamic, nest_offl%nt, 'time' )
+
+          IF ( check_existence( nest_offl%var_names, 'time' ) )  THEN
+             ALLOCATE( nest_offl%time(0:nest_offl%nt-1) )
+             CALL get_variable( id_dynamic, 'time', nest_offl%time )
           ENDIF
 !
 !--       Read vertical dimension of scalar und w grid
-          CALL get_dimension_length( id_dynamic, force%nzu, 'z' )
-          CALL get_dimension_length( id_dynamic, force%nzw, 'zw' )
-
-          IF ( check_existence( force%var_names, 'z' ) )  THEN
-             ALLOCATE( force%zu_atmos(1:force%nzu) )
-             CALL get_variable( id_dynamic, 'z', force%zu_atmos )
-          ENDIF
-          IF ( check_existence( force%var_names, 'zw' ) )  THEN
-             ALLOCATE( force%zw_atmos(1:force%nzw) )
-             CALL get_variable( id_dynamic, 'zw', force%zw_atmos )
+          CALL get_dimension_length( id_dynamic, nest_offl%nzu, 'z' )
+          CALL get_dimension_length( id_dynamic, nest_offl%nzw, 'zw' )
+
+          IF ( check_existence( nest_offl%var_names, 'z' ) )  THEN
+             ALLOCATE( nest_offl%zu_atmos(1:nest_offl%nzu) )
+             CALL get_variable( id_dynamic, 'z', nest_offl%zu_atmos )
+          ENDIF
+          IF ( check_existence( nest_offl%var_names, 'zw' ) )  THEN
+             ALLOCATE( nest_offl%zw_atmos(1:nest_offl%nzw) )
+             CALL get_variable( id_dynamic, 'zw', nest_offl%zw_atmos )
           ENDIF
 
 !
 !--       Read surface pressure
-          IF ( check_existence( force%var_names,                               &
-                            'surface_forcing_surface_pressure' ) )  THEN
-             ALLOCATE( force%surface_pressure(0:force%nt-1) )
+          IF ( check_existence( nest_offl%var_names,                           &
+                                'surface_forcing_surface_pressure' ) )  THEN
+             ALLOCATE( nest_offl%surface_pressure(0:nest_offl%nt-1) )
              CALL get_variable( id_dynamic,                                    &
                                 'surface_forcing_surface_pressure',            &
-                                force%surface_pressure )
+                                nest_offl%surface_pressure )
           ENDIF
 !
 !--       Set control flag to indicate that initialization is already done
-          force%init = .TRUE.
+          nest_offl%init = .TRUE.
 
        ENDIF
@@ -2610 +2672 @@
 !--    @todo: At the moment time, in INIFOR and simulated time correspond
 !--           to each other. If required, adjust to daytime.
-       force%tind = MINLOC( ABS( force%time - simulated_time ), DIM = 1 )      &
-                    - 1
-       force%tind_p = force%tind + 1       
+       nest_offl%tind = MINLOC( ABS( nest_offl%time -                          &
+                                     time_since_reference_point ), DIM = 1 )   &
+                        - 1
+       nest_offl%tind_p = nest_offl%tind + 1       
 !
 !--    Read geostrophic wind components. In case of forcing, this is only
 !--    required if cyclic boundary conditions are applied.
        IF ( bc_lr_cyc  .AND.  bc_ns_cyc )  THEN
-          DO  t = force%tind, force%tind_p
+          DO  t = nest_offl%tind, nest_offl%tind_p
 !              CALL get_variable_pr( id_dynamic, 'tend_ug', t+1,           &
-!                                    force%ug(t-force%tind,:) )
+!                                    nest_offl%ug(t-nest_offl%tind,:) )
 !              CALL get_variable_pr( id_dynamic, 'tend_vg', t+1,           &
-!                                    force%ug(t-force%tind,:) )
+!                                    nest_offl%ug(t-nest_offl%tind,:) )
              CALL get_variable_pr( id_dynamic, 'ls_forcing_ug', t+1,           &
-                                   force%ug(t-force%tind,:) )
+                                   nest_offl%ug(t-nest_offl%tind,:) )
              CALL get_variable_pr( id_dynamic, 'ls_forcing_vg', t+1,           &
-                                   force%ug(t-force%tind,:) )
+                                   nest_offl%ug(t-nest_offl%tind,:) )
           ENDDO
        ENDIF
@@ -2633 +2696 @@
 !--    For the v-component, the data starts at nysv, while for the other
 !--    quantities the data starts at nys. This is equivalent at the north
-!--    and south domain boundary for the u-component.
-       IF ( force_bound_l )  THEN
+!--    and south domain boundary for the u-component. 
+!--    Further, lateral data is not accessed by parallel IO, indicated by the 
+!--    last passed flag in the subroutine get_variable(). This is because
+!--    not every PE participates in this collective blocking read operation. 
+       IF ( bc_dirichlet_l )  THEN
           CALL get_variable( id_dynamic, 'ls_forcing_left_u',                  &
-                           force%u_left(0:1,nzb+1:nzt+1,nys:nyn),              &
-                           nys+1, nzb+1, force%tind+1,                         &
-                           nyn-nys+1, force%nzu, 2, dynamic_3d )
-          
+                           nest_offl%u_left(0:1,nzb+1:nzt,nys:nyn),            &
+                           nys+1, nzb+1, nest_offl%tind+1,                     &
+                           nyn-nys+1, nest_offl%nzu, 2, .FALSE. )
+      
           CALL get_variable( id_dynamic, 'ls_forcing_left_v',                  &
-                           force%v_left(0:1,nzb+1:nzt+1,nysv:nyn),             &
-                           nysv, nzb+1, force%tind+1,                          &
-                           nyn-nysv+1, force%nzu, 2, dynamic_3d )
+                           nest_offl%v_left(0:1,nzb+1:nzt,nysv:nyn),           &
+                           nysv, nzb+1, nest_offl%tind+1,                      &
+                           nyn-nysv+1, nest_offl%nzu, 2, .FALSE. )
 
           CALL get_variable( id_dynamic, 'ls_forcing_left_w',                  &
-                           force%w_left(0:1,nzb+1:nzt,nys:nyn),                &
-                           nys+1, nzb+1, force%tind+1,                         &
-                           nyn-nys+1, force%nzw, 2, dynamic_3d )
+                           nest_offl%w_left(0:1,nzb+1:nzt-1,nys:nyn),          &
+                           nys+1, nzb+1, nest_offl%tind+1,                     &
+                           nyn-nys+1, nest_offl%nzw, 2, .FALSE. )
 
           IF ( .NOT. neutral )  THEN
              CALL get_variable( id_dynamic, 'ls_forcing_left_pt',              &
-                           force%pt_left(0:1,nzb+1:nzt+1,nys:nyn),             &
-                           nys+1, nzb+1, force%tind+1,                         &
-                           nyn-nys+1, force%nzu, 2, dynamic_3d )
-          ENDIF
+                           nest_offl%pt_left(0:1,nzb+1:nzt,nys:nyn),           &
+                           nys+1, nzb+1, nest_offl%tind+1,                     &
+                           nyn-nys+1, nest_offl%nzu, 2, .FALSE. )
+          ENDIF
+
           IF ( humidity )  THEN
              CALL get_variable( id_dynamic, 'ls_forcing_left_qv',              &
-                           force%q_left(0:1,nzb+1:nzt+1,nys:nyn),              &
-                           nys+1, nzb+1, force%tind+1,                         &
-                           nyn-nys+1, force%nzu, 2, dynamic_3d )
-          ENDIF
-       ENDIF
-
-       IF ( force_bound_r )  THEN
+                           nest_offl%q_left(0:1,nzb+1:nzt,nys:nyn),            &
+                           nys+1, nzb+1, nest_offl%tind+1,                     &
+                           nyn-nys+1, nest_offl%nzu, 2, .FALSE. )
+          ENDIF
+
+       ENDIF
+
+       IF ( bc_dirichlet_r )  THEN
           CALL get_variable( id_dynamic, 'ls_forcing_right_u',                 &
-                           force%u_right(0:1,nzb+1:nzt+1,nys:nyn),             &
-                           nys+1, nzb+1, force%tind+1,                         &
-                           nyn-nys+1, force%nzu, 2, dynamic_3d )
+                           nest_offl%u_right(0:1,nzb+1:nzt,nys:nyn),           &
+                           nys+1, nzb+1, nest_offl%tind+1,                     &
+                           nyn-nys+1, nest_offl%nzu, 2, .FALSE. )
                            
           CALL get_variable( id_dynamic, 'ls_forcing_right_v',                 &
-                           force%v_right(0:1,nzb+1:nzt+1,nysv:nyn),            &
-                           nysv, nzb+1, force%tind+1,                          &
-                           nyn-nysv+1, force%nzu, 2, dynamic_3d )
+                           nest_offl%v_right(0:1,nzb+1:nzt,nysv:nyn),          &
+                           nysv, nzb+1, nest_offl%tind+1,                      &
+                           nyn-nysv+1, nest_offl%nzu, 2, .FALSE. )
                            
           CALL get_variable( id_dynamic, 'ls_forcing_right_w',                 &
-                           force%w_right(0:1,nzb+1:nzt,nys:nyn),               &
-                           nys+1, nzb+1, force%tind+1,                         &
-                           nyn-nys+1, force%nzw, 2, dynamic_3d )
+                           nest_offl%w_right(0:1,nzb+1:nzt-1,nys:nyn),         &
+                           nys+1, nzb+1, nest_offl%tind+1,                     &
+                           nyn-nys+1, nest_offl%nzw, 2, .FALSE. )
                            
           IF ( .NOT. neutral )  THEN
              CALL get_variable( id_dynamic, 'ls_forcing_right_pt',             &
-                           force%pt_right(0:1,nzb+1:nzt+1,nys:nyn),            &
-                           nys+1, nzb+1, force%tind+1,                         &
-                           nyn-nys+1, force%nzu, 2, dynamic_3d )
+                           nest_offl%pt_right(0:1,nzb+1:nzt,nys:nyn),          &
+                           nys+1, nzb+1, nest_offl%tind+1,                     &
+                           nyn-nys+1, nest_offl%nzu, 2, .FALSE. )
           ENDIF
           IF ( humidity )  THEN
              CALL get_variable( id_dynamic, 'ls_forcing_right_qv',             &
-                           force%q_right(0:1,nzb+1:nzt+1,nys:nyn),             &
-                           nys+1, nzb+1, force%tind+1,                         &
-                           nyn-nys+1, force%nzu, 2, dynamic_3d )
-          ENDIF
-       ENDIF
-
-       IF ( force_bound_n )  THEN
+                           nest_offl%q_right(0:1,nzb+1:nzt,nys:nyn),           &
+                           nys+1, nzb+1, nest_offl%tind+1,                     &
+                           nyn-nys+1, nest_offl%nzu, 2, .FALSE. )
+          ENDIF
+       ENDIF
+
+       IF ( bc_dirichlet_n )  THEN
        
           CALL get_variable( id_dynamic, 'ls_forcing_north_u',                 &
-                           force%u_north(0:1,nzb+1:nzt+1,nxlu:nxr),            &
-                           nxlu, nzb+1, force%tind+1,                          &
-                           nxr-nxlu+1, force%nzu, 2, dynamic_3d )
-
+                           nest_offl%u_north(0:1,nzb+1:nzt,nxlu:nxr),          &
+                           nxlu, nzb+1, nest_offl%tind+1,                      &
+                           nxr-nxlu+1, nest_offl%nzu, 2, .FALSE. )
+                           
           CALL get_variable( id_dynamic, 'ls_forcing_north_v',                 &
-                           force%v_north(0:1,nzb+1:nzt+1,nxl:nxr),             &
-                           nxl+1, nzb+1, force%tind+1,                         &
-                           nxr-nxl+1, force%nzu, 2, dynamic_3d )
+                           nest_offl%v_north(0:1,nzb+1:nzt,nxl:nxr),           &
+                           nxl+1, nzb+1, nest_offl%tind+1,                     &
+                           nxr-nxl+1, nest_offl%nzu, 2, .FALSE. )
                            
           CALL get_variable( id_dynamic, 'ls_forcing_north_w',                 &
-                           force%w_north(0:1,nzb+1:nzt,nxl:nxr),               &
-                           nxl+1, nzb+1, force%tind+1,                         &
-                           nxr-nxl+1, force%nzw, 2, dynamic_3d )
+                           nest_offl%w_north(0:1,nzb+1:nzt-1,nxl:nxr),         &
+                           nxl+1, nzb+1, nest_offl%tind+1,                     &
+                           nxr-nxl+1, nest_offl%nzw, 2, .FALSE. )
                            
           IF ( .NOT. neutral )  THEN
              CALL get_variable( id_dynamic, 'ls_forcing_north_pt',             &
-                           force%pt_north(0:1,nzb+1:nzt+1,nxl:nxr),            &
-                           nxl+1, nzb+1, force%tind+1,                         &
-                           nxr-nxl+1, force%nzu, 2, dynamic_3d )
+                           nest_offl%pt_north(0:1,nzb+1:nzt,nxl:nxr),          &
+                           nxl+1, nzb+1, nest_offl%tind+1,                     &
+                           nxr-nxl+1, nest_offl%nzu, 2, .FALSE. )
           ENDIF
           IF ( humidity )  THEN
              CALL get_variable( id_dynamic, 'ls_forcing_north_qv',             &
-                           force%q_north(0:1,nzb+1:nzt+1,nxl:nxr),             &
-                           nxl+1, nzb+1, force%tind+1,                         &
-                           nxr-nxl+1, force%nzu, 2, dynamic_3d )
-          ENDIF
-       ENDIF
-
-       IF ( force_bound_s )  THEN
+                           nest_offl%q_north(0:1,nzb+1:nzt,nxl:nxr),           &
+                           nxl+1, nzb+1, nest_offl%tind+1,                     &
+                           nxr-nxl+1, nest_offl%nzu, 2, .FALSE. )
+          ENDIF
+       ENDIF
+
+       IF ( bc_dirichlet_s )  THEN
           CALL get_variable( id_dynamic, 'ls_forcing_south_u',                 &
-                           force%u_south(0:1,nzb+1:nzt+1,nxlu:nxr),            &
-                           nxlu, nzb+1, force%tind+1,                          &
-                           nxr-nxlu+1, force%nzu, 2, dynamic_3d )
+                           nest_offl%u_south(0:1,nzb+1:nzt,nxlu:nxr),          &
+                           nxlu, nzb+1, nest_offl%tind+1,                      &
+                           nxr-nxlu+1, nest_offl%nzu, 2, .FALSE. )
 
           CALL get_variable( id_dynamic, 'ls_forcing_south_v',                 &
-                           force%v_south(0:1,nzb+1:nzt+1,nxl:nxr),             &
-                           nxl+1, nzb+1, force%tind+1,                         &
-                           nxr-nxl+1, force%nzu, 2, dynamic_3d )
+                           nest_offl%v_south(0:1,nzb+1:nzt,nxl:nxr),           &
+                           nxl+1, nzb+1, nest_offl%tind+1,                     &
+                           nxr-nxl+1, nest_offl%nzu, 2, .FALSE. )
                            
           CALL get_variable( id_dynamic, 'ls_forcing_south_w',                 &
-                           force%w_south(0:1,nzb+1:nzt,nxl:nxr),               &
-                           nxl+1, nzb+1, force%tind+1,                         &
-                           nxr-nxl+1, force%nzw, 2, dynamic_3d )
+                           nest_offl%w_south(0:1,nzb+1:nzt-1,nxl:nxr),         &
+                           nxl+1, nzb+1, nest_offl%tind+1,                     &
+                           nxr-nxl+1, nest_offl%nzw, 2, .FALSE. )
                            
           IF ( .NOT. neutral )  THEN
              CALL get_variable( id_dynamic, 'ls_forcing_south_pt',             &
-                           force%pt_south(0:1,nzb+1:nzt+1,nxl:nxr),            &
-                           nxl+1, nzb+1, force%tind+1,                         &
-                           nxr-nxl+1, force%nzu, 2, dynamic_3d )
+                           nest_offl%pt_south(0:1,nzb+1:nzt,nxl:nxr),          &
+                           nxl+1, nzb+1, nest_offl%tind+1,                     &
+                           nxr-nxl+1, nest_offl%nzu, 2, .FALSE. )
           ENDIF
           IF ( humidity )  THEN
              CALL get_variable( id_dynamic, 'ls_forcing_south_qv',             &
-                           force%q_south(0:1,nzb+1:nzt+1,nxl:nxr),             &
-                           nxl+1, nzb+1, force%tind+1,                         &
-                           nxr-nxl+1, force%nzu, 2, dynamic_3d )
-          ENDIF
-       ENDIF
+                           nest_offl%q_south(0:1,nzb+1:nzt,nxl:nxr),           &
+                           nxl+1, nzb+1, nest_offl%tind+1,                     &
+                           nxr-nxl+1, nest_offl%nzu, 2, .FALSE. )
+          ENDIF
+       ENDIF
+
 !
 !--    Top boundary
        CALL get_variable( id_dynamic, 'ls_forcing_top_u',                      &
-                             force%u_top(0:1,nys:nyn,nxlu:nxr),                &
-                             nxlu, nys+1, force%tind+1,                        &
-                             nxr-nxlu+1, nyn-nys+1, 2, dynamic_3d )
+                             nest_offl%u_top(0:1,nys:nyn,nxlu:nxr),            &
+                             nxlu, nys+1, nest_offl%tind+1,                    &
+                             nxr-nxlu+1, nyn-nys+1, 2, .TRUE. )
 
        CALL get_variable( id_dynamic, 'ls_forcing_top_v',                      &
-                             force%v_top(0:1,nysv:nyn,nxl:nxr),                &
-                             nxl+1, nysv, force%tind+1,                        &
-                             nxr-nxl+1, nyn-nysv+1, 2, dynamic_3d )
+                             nest_offl%v_top(0:1,nysv:nyn,nxl:nxr),            &
+                             nxl+1, nysv, nest_offl%tind+1,                    &
+                             nxr-nxl+1, nyn-nysv+1, 2, .TRUE. )
                              
        CALL get_variable( id_dynamic, 'ls_forcing_top_w',                      &
-                             force%w_top(0:1,nys:nyn,nxl:nxr),                 &
-                             nxl+1, nys+1, force%tind+1,                       &
-                             nxr-nxl+1, nyn-nys+1, 2, dynamic_3d )
+                             nest_offl%w_top(0:1,nys:nyn,nxl:nxr),             &
+                             nxl+1, nys+1, nest_offl%tind+1,                   &
+                             nxr-nxl+1, nyn-nys+1, 2, .TRUE. )
                              
        IF ( .NOT. neutral )  THEN
           CALL get_variable( id_dynamic, 'ls_forcing_top_pt',                  &
-                                force%pt_top(0:1,nys:nyn,nxl:nxr),             &
-                                nxl+1, nys+1, force%tind+1,                    &
-                                nxr-nxl+1, nyn-nys+1, 2, dynamic_3d )
+                                nest_offl%pt_top(0:1,nys:nyn,nxl:nxr),         &
+                                nxl+1, nys+1, nest_offl%tind+1,                &
+                                nxr-nxl+1, nyn-nys+1, 2, .TRUE. )
        ENDIF
        IF ( humidity )  THEN
           CALL get_variable( id_dynamic, 'ls_forcing_top_qv',                  &
-                                force%q_top(0:1,nys:nyn,nxl:nxr),              &
-                                nxl+1, nys+1, force%tind+1,                    &
-                                nxr-nxl+1, nyn-nys+1, 2, dynamic_3d )
+                                nest_offl%q_top(0:1,nys:nyn,nxl:nxr),          &
+                                nxl+1, nys+1, nest_offl%tind+1,                &
+                                nxr-nxl+1, nyn-nys+1, 2, .TRUE. )
        ENDIF
 
@@ -2792 +2861 @@
        CALL cpu_log( log_point_s(86), 'NetCDF input forcing', 'stop' )
 
-!
-!--    Finally, after data input set control flag indicating that vertical
-!--    inter- and/or extrapolation is required.
-!--    Please note, inter/extrapolation of INIFOR data is only a workaroud,
-!--    as long as INIFOR delivers vertically equidistant data.
-       force%interpolated = .FALSE.
-
     END SUBROUTINE netcdf_data_input_lsf
 
@@ -2810 +2872 @@
 
        USE control_parameters,                                                 &
-           ONLY:  initializing_actions, forcing, message_string
+           ONLY:  initializing_actions, message_string, nesting_offline 
 
        IMPLICIT NONE
@@ -2816 +2878 @@
 !
 !--    In case of forcing, check whether dynamic input file is present
-       IF ( .NOT. input_pids_dynamic  .AND.  forcing )  THEN
-          message_string = 'forcing = .TRUE. requires dynamic input file ' //  &
+       IF ( .NOT. input_pids_dynamic  .AND.  nesting_offline  )  THEN
+          message_string = 'nesting_offline = .TRUE. requires dynamic '  //    &
+                            'input file ' //                                   &
                             TRIM( input_file_dynamic ) // TRIM( coupling_char )
           CALL message( 'netcdf_data_input_mod', 'NDI009', 1, 2, 0, 6, 0 )
@@ -4514 +4577 @@
 !------------------------------------------------------------------------------!
     SUBROUTINE get_variable_3d_real_dynamic( id, variable_name, var,           &
-                            i1s, i2s, i3s, count_1, count_2, count_3, dynamic)
+                                             i1s, i2s, i3s,                    &
+                                             count_1, count_2, count_3,        &
+                                             par_access )
                                 
        USE indices
@@ -4523 +4588 @@
        CHARACTER(LEN=*)              ::  variable_name   !< variable name
 
-       LOGICAL                       ::  dynamic         !< additional flag just used to select correct overloaded routine from interface block 
+       LOGICAL                       ::  par_access      !< additional flag indicating whether parallel read operations should be performed or not 
        
        INTEGER(iwp)                  ::  count_1         !< number of elements to be read along 1st dimension (with respect to file)
@@ -4550 +4615 @@
 #if defined( __netcdf )
 !
-!--    Inquire variable id
+!--    Inquire variable id.
        nc_stat = NF90_INQ_VARID( id, TRIM( variable_name ), id_var )
 !
 !--    Check for collective read-operation and set respective NetCDF flags if
 !--    required. 
-       IF ( collective_read )  THEN
+!--    Please note, in contrast to the other input routines where each PEs 
+!--    reads its subdomain data, dynamic input data not by all PEs, only
+!--    by those which encompass lateral model boundaries. Hence, collective
+!--    read operations are only enabled for top-boundary data.
+       IF ( collective_read  .AND.  par_access )  THEN
           nc_stat = NF90_VAR_PAR_ACCESS (id, id_var, NF90_COLLECTIVE)
        ENDIF    

palm/trunk/SOURCE/palm.f90

@@ -238 +238 @@
     USE control_parameters,                                                    &
         ONLY:  air_chemistry,                                                  &
-               cloud_physics, constant_diffusion, coupling_char, coupling_mode,&
+               cloud_physics, constant_diffusion, child_domain, coupling_char, &
+               coupling_mode,                                                  &
                do2d_at_begin, do3d_at_begin, humidity, initializing_actions,   &
                io_blocks, io_group, large_scale_forcing,                       &
                message_string, microphysics_morrison, microphysics_seifert,    &
-               nest_domain, neutral, nudging, passive_scalar, runnr,           &
+               neutral, nudging, passive_scalar, runnr,                        &
                simulated_time, simulated_time_chr, spinup,                     &
                time_since_reference_point,                                     &
@@ -428 +429 @@
 !
 !--    Exchange_horiz is needed after the nest initialization
-          IF ( nest_domain )  THEN
+          IF ( child_domain )  THEN
              CALL exchange_horiz( u, nbgp )
              CALL exchange_horiz( v, nbgp )

palm/trunk/SOURCE/parin.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! Rename variables and boundary conditions in mesoscale-offline nesting mode
 ! 
 ! Former revisions:
@@ -529 +529 @@
              dz_stretch_factor, dz_stretch_level, dz_stretch_level_start,      &
              dz_stretch_level_end, end_time_1d, ensemble_member_nr, e_init,    &
-             e_min, fft_method, flux_input_mode, flux_output_mode, forcing,    &
+             e_min, fft_method, flux_input_mode, flux_output_mode,             &
              galilei_transformation, humidity,                                 &
              inflow_damping_height, inflow_damping_width,                      &
@@ -538 +538 @@
              loop_optimization, lsf_exception, masking_method, mg_cycles,      &
              mg_switch_to_pe0_level, mixing_length_1d, momentum_advec,         &
-             most_method, na_init, nc_const, netcdf_precision, neutral, ngsrb, &
+             most_method, na_init, nc_const, nesting_offline,                  &
+             netcdf_precision, neutral, ngsrb,                                 &
              nsor, nsor_ini, nudging, nx, ny, nz, ocean, omega, omega_sor,     &
              outflow_source_plane, passive_scalar,                             &
@@ -601 +602 @@
              dz_stretch_factor, dz_stretch_level, dz_stretch_level_start,      &
              dz_stretch_level_end, end_time_1d, ensemble_member_nr, e_init,    &
-             e_min, fft_method, flux_input_mode, flux_output_mode, forcing,    &
+             e_min, fft_method, flux_input_mode, flux_output_mode,             &
              galilei_transformation, humidity,                                 &
              inflow_damping_height, inflow_damping_width,                      &
@@ -610 +611 @@
              loop_optimization, lsf_exception, masking_method, mg_cycles,      &
              mg_switch_to_pe0_level, mixing_length_1d, momentum_advec,         &
-             most_method, na_init, nc_const, netcdf_precision, neutral, ngsrb, &
+             most_method, na_init, nc_const, nesting_offline,                  &
+             netcdf_precision, neutral, ngsrb,                                 &
              nsor, nsor_ini, nudging, nx, ny, nz, ocean, omega, omega_sor,     &
              outflow_source_plane, passive_scalar,                             &
@@ -924 +926 @@
                    bc_ns   = 'cyclic'
                 ENDIF
-                IF ( nest_domain )  THEN
+                IF ( child_domain )  THEN
                    bc_uv_t  = 'nested'
                    bc_pt_t  = 'nested'
@@ -936 +938 @@
 !--          nested domains.
              ELSE 
-                IF ( nest_domain )  THEN
+                IF ( child_domain )  THEN
                    bc_lr    = 'nested'
                    bc_ns    = 'nested'
@@ -948 +950 @@
              ENDIF
           ENDIF
-
-          IF ( forcing )  THEN
-             bc_lr    = 'forcing'
-             bc_ns    = 'forcing'
-             bc_uv_t  = 'forcing'
-             bc_pt_t  = 'forcing'
-             bc_q_t   = 'forcing'
-             bc_s_t   = 'forcing'  ! scalar boundary condition is not clear
-             bc_cs_t  = 'forcing'  ! same for chemical species
+!
+!--       Set boundary conditions also in case the model is offline-nested in
+!--       larger-scale models.
+          IF ( nesting_offline )  THEN
+             bc_lr    = 'nesting_offline'
+             bc_ns    = 'nesting_offline'
+             bc_uv_t  = 'nesting_offline'
+             bc_pt_t  = 'nesting_offline'
+             bc_q_t   = 'nesting_offline'
+             bc_s_t   = 'nesting_offline'  ! scalar boundary condition is not clear
+             bc_cs_t  = 'nesting_offline'  ! same for chemical species
              bc_p_t   = 'neumann'
           ENDIF
@@ -968 +972 @@
 !--       are set properly. An exception is made in case of restart runs and
 !--       if user decides to do everything by its own.
-          IF ( nest_domain  .AND.  .NOT. (                                     &
+          IF ( child_domain  .AND.  .NOT. (                                    &
                TRIM( initializing_actions ) == 'read_restart_data'      .OR.   &
                TRIM( initializing_actions ) == 'set_constant_profiles'  .OR.   &
@@ -979 +983 @@
 
              initializing_actions = 'set_constant_profiles'
-          ENDIF
-            
+          ENDIF            
 !
 !--       Check validity of lateral boundary conditions. This has to be done
@@ -987 +990 @@
           IF ( bc_lr /= 'cyclic'  .AND.  bc_lr /= 'dirichlet/radiation'  .AND. &
                bc_lr /= 'radiation/dirichlet'  .AND.  bc_lr /= 'nested'  .AND. &
-               bc_lr /= 'forcing' )  THEN
+               bc_lr /= 'nesting_offline' )  THEN
              message_string = 'unknown boundary condition: bc_lr = "' // &
                               TRIM( bc_lr ) // '"'
@@ -994 +997 @@
           IF ( bc_ns /= 'cyclic'  .AND.  bc_ns /= 'dirichlet/radiation'  .AND. &
                bc_ns /= 'radiation/dirichlet'  .AND.  bc_ns /= 'nested'  .AND. &
-               bc_ns /= 'forcing' )  THEN
+               bc_ns /= 'nesting_offline' )  THEN
              message_string = 'unknown boundary condition: bc_ns = "' // &
                               TRIM( bc_ns ) // '"'

palm/trunk/SOURCE/pmc_interface_mod.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-! 
+! Variable names for nest_bound_x replaced by bc_dirichlet_x.
+! Remove commented prints into debug files.
 ! 
 ! Former revisions:
@@ -306 +307 @@
    USE arrays_3d,                                                              &
         ONLY:  diss, diss_2, dzu, dzw, e, e_p, e_2, nc, nc_2, nc_p, nr, nr_2,  &
-               pt, pt_2, q, q_2, qc, qc_2, qr, qr_2, s, s_2,                       &
+               pt, pt_2, q, q_2, qc, qc_2, qr, qr_2, s, s_2,                   &
                u, u_p, u_2, v, v_p, v_2, w, w_p, w_2, zu, zw
 #endif
 
     USE control_parameters,                                                    &
-        ONLY:  air_chemistry, cloud_physics,                                   &
+        ONLY:  air_chemistry, bc_dirichlet_l, bc_dirichlet_n, bc_dirichlet_r,  &
+               bc_dirichlet_s, cloud_physics, child_domain,                    &
                constant_diffusion, constant_flux_layer,                        &
                coupling_char, dt_3d, dz, humidity, message_string,             &
                microphysics_morrison, microphysics_seifert,                    &
-               nest_bound_l, nest_bound_r, nest_bound_s, nest_bound_n,         &
-               nest_domain, neutral, passive_scalar, rans_mode, rans_tke_e,    &
+               neutral, passive_scalar, rans_mode, rans_tke_e,                 &
                roughness_length, simulated_time, topography, volume_flow
 
@@ -743 +744 @@
 !-- course the root domain (cpl_id = 1) is not nested)
     IF ( cpl_id >= 2 )  THEN
-       nest_domain = .TRUE.
+       child_domain = .TRUE.
        WRITE( coupling_char, '(A2,I2.2)') '_N', cpl_id
     ENDIF
@@ -1192 +1193 @@
         childgrid(m)%uz_coord   = zmax_coarse(2)
         childgrid(m)%uz_coord_b = zmax_coarse(1)
-
-!         WRITE(9,*)                 'edge coordinates for child id ',child_id,m
-!         WRITE(9,*)                 'Number of Boundray cells lpm  ',nbgp_lpm
-!         WRITE(9,'(a,3i7,2f10.2)') ' model size                    ', nx_cl, ny_cl, nz_cl, dx_cl, dy_cl
-!          WRITE(9,'(a,5f10.2)')     ' model edge                    ', childgrid(m)%lx_coord,  &
-!                                childgrid(m)%rx_coord, childgrid(m)%sy_coord,                  &
-!                                childgrid(m)%ny_coord,childgrid(m)%uz_coord
-!          WRITE(9,'(a,4f10.2)')     ' model edge with Boundary      ', childgrid(m)%lx_coord_b,& 
-!                                childgrid(m)%rx_coord_b, childgrid(m)%sy_coord_b,              &
-!                                childgrid(m)%ny_coord_b
 
      END SUBROUTINE set_child_edge_coords
@@ -1714 +1705 @@
 !--    interpolation routines.
        nzt_topo_nestbc_l = nzb
-       IF ( nest_bound_l )  THEN
+       IF ( bc_dirichlet_l )  THEN
           DO  i = nxl-1, nxl
              DO  j = nys, nyn
@@ -1740 +1731 @@
       
        nzt_topo_nestbc_r = nzb
-       IF ( nest_bound_r )  THEN
+       IF ( bc_dirichlet_r )  THEN
           i = nxr + 1
           DO  j = nys, nyn
@@ -1765 +1756 @@
 
        nzt_topo_nestbc_s = nzb
-       IF ( nest_bound_s )  THEN
+       IF ( bc_dirichlet_s )  THEN
           DO  j = nys-1, nys
              DO  i = nxl, nxr
@@ -1791 +1782 @@
 
        nzt_topo_nestbc_n = nzb
-       IF ( nest_bound_n )  THEN
+       IF ( bc_dirichlet_n )  THEN
           j = nyn + 1
           DO  i = nxl, nxr
@@ -1862 +1853 @@
 !--    logc_kbounds_* need to be allocated and initialized here.
 !--    Left boundary
-       IF ( nest_bound_l )  THEN
+       IF ( bc_dirichlet_l )  THEN
 
           ALLOCATE( logc_u_l(1:2,nzb:nzt_topo_nestbc_l,nys:nyn) )
@@ -1923 +1914 @@
 !
 !--    Right boundary
-       IF ( nest_bound_r )  THEN
+       IF ( bc_dirichlet_r )  THEN
            
           ALLOCATE( logc_u_r(1:2,nzb:nzt_topo_nestbc_r,nys:nyn) )
@@ -1985 +1976 @@
 !
 !--    South boundary
-       IF ( nest_bound_s )  THEN
+       IF ( bc_dirichlet_s )  THEN
 
           ALLOCATE( logc_u_s(1:2,nzb:nzt_topo_nestbc_s,nxl:nxr) )
@@ -2044 +2035 @@
 !
 !--    North boundary
-       IF ( nest_bound_n )  THEN
+       IF ( bc_dirichlet_n )  THEN
 
           ALLOCATE( logc_u_n(1:2,nzb:nzt_topo_nestbc_n,nxl:nxr) )
@@ -2116 +2107 @@
 !        
 !--       Left boundary
-          IF ( nest_bound_l )  THEN
+          IF ( bc_dirichlet_l )  THEN
              logc_kbounds_u_l(1:2,nys:nyn) = 0
              logc_kbounds_v_l(1:2,nys:nyn) = 0             
@@ -2255 +2246 @@
              ENDDO
 
-          ENDIF   !  IF ( nest_bound_l )
+          ENDIF   !  IF ( bc_dirichlet_l )
 !        
 !--       Right boundary
-          IF ( nest_bound_r )  THEN
+          IF ( bc_dirichlet_r )  THEN
              logc_kbounds_u_r(1:2,nys:nyn) = 0
              logc_kbounds_v_r(1:2,nys:nyn) = 0             
@@ -2395 +2386 @@
              ENDDO
              
-          ENDIF   !  IF ( nest_bound_r )
+          ENDIF   !  IF ( bc_dirichlet_r )
 !        
 !--       South boundary
-          IF ( nest_bound_s )  THEN
+          IF ( bc_dirichlet_s )  THEN
              logc_kbounds_u_s(1:2,nxl:nxr) = 0
              logc_kbounds_v_s(1:2,nxl:nxr) = 0
@@ -2537 +2528 @@
              ENDDO
 
-          ENDIF   !  IF (nest_bound_s )
+          ENDIF   !  IF (bc_dirichlet_s )
 !        
 !--       North boundary
-          IF ( nest_bound_n )  THEN
+          IF ( bc_dirichlet_n )  THEN
              logc_kbounds_u_n(1:2,nxl:nxr) = 0             
              logc_kbounds_v_n(1:2,nxl:nxr) = 0
@@ -2678 +2669 @@
              ENDDO
 
-          ENDIF   !  IF ( nest_bound_n )
+          ENDIF   !  IF ( bc_dirichlet_n )
 
        ENDIF   !  IF ( topography /= 'flat' )
@@ -2786 +2777 @@
 
        END SELECT
-          
-       !write(9,"('pmci_define_loglaw_correction_parameters: ', 6(i3,2x))")    &
-       !     direction, ij, k, wall_index, inc, lc
 
     END SUBROUTINE pmci_define_loglaw_correction_parameters
@@ -3325 +3313 @@
 !
        IF ( .NOT. rans_mode  .AND.  .NOT. rans_mode_parent )  THEN 
-          IF ( nest_bound_l )  THEN
+          IF ( bc_dirichlet_l )  THEN
              ALLOCATE( tkefactor_l(nzb:nzt+1,nysg:nyng) )
              tkefactor_l = 0.0_wp
@@ -3345 +3333 @@
           ENDIF
 
-          IF ( nest_bound_r )  THEN
+          IF ( bc_dirichlet_r )  THEN
              ALLOCATE( tkefactor_r(nzb:nzt+1,nysg:nyng) )
              tkefactor_r = 0.0_wp
@@ -3365 +3353 @@
           ENDIF
 
-          IF ( nest_bound_s )  THEN
+          IF ( bc_dirichlet_s )  THEN
              ALLOCATE( tkefactor_s(nzb:nzt+1,nxlg:nxrg) )
              tkefactor_s = 0.0_wp
@@ -3386 +3374 @@
           ENDIF
 
-          IF ( nest_bound_n )  THEN
+          IF ( bc_dirichlet_n )  THEN
              ALLOCATE( tkefactor_n(nzb:nzt+1,nxlg:nxrg) )
              tkefactor_n = 0.0_wp
@@ -3428 +3416 @@
 !--    RANS mode
        ELSE
-          IF ( nest_bound_l )  THEN
+          IF ( bc_dirichlet_l )  THEN
              ALLOCATE( tkefactor_l(nzb:nzt+1,nysg:nyng) )
              tkefactor_l = 1.0_wp
           ENDIF
-          IF ( nest_bound_r )  THEN
+          IF ( bc_dirichlet_r )  THEN
              ALLOCATE( tkefactor_r(nzb:nzt+1,nysg:nyng) )
              tkefactor_r = 1.0_wp
           ENDIF
-          IF ( nest_bound_s )  THEN
+          IF ( bc_dirichlet_s )  THEN
              ALLOCATE( tkefactor_s(nzb:nzt+1,nxlg:nxrg) )
              tkefactor_s = 1.0_wp
           ENDIF
-          IF ( nest_bound_n )  THEN
+          IF ( bc_dirichlet_n )  THEN
              ALLOCATE( tkefactor_n(nzb:nzt+1,nxlg:nxrg) )
              tkefactor_n = 1.0_wp
@@ -3980 +3968 @@
        je = nyn
        IF ( nesting_mode /= 'vertical' )  THEN
-          IF ( nest_bound_l )  THEN
+          IF ( bc_dirichlet_l )  THEN
              ib = nxl - 1
 !
@@ -3988 +3976 @@
              ENDIF
           ENDIF
-          IF ( nest_bound_s )  THEN
+          IF ( bc_dirichlet_s )  THEN
              jb = nys - 1
 !
@@ -3996 +3984 @@
              ENDIF
           ENDIF
-          IF ( nest_bound_r )  THEN
+          IF ( bc_dirichlet_r )  THEN
              ie = nxr + 1
           ENDIF
-          IF ( nest_bound_n )  THEN
+          IF ( bc_dirichlet_n )  THEN
              je = nyn + 1
           ENDIF
@@ -4201 +4189 @@
     volume_flow_l(1) = 0.0_wp
 
-    IF ( nest_bound_l )  THEN
+    IF ( bc_dirichlet_l )  THEN
        i = 0
        innor = dy
@@ -4213 +4201 @@
     ENDIF
 
-    IF ( nest_bound_r )  THEN
+    IF ( bc_dirichlet_r )  THEN
        i = nx + 1
        innor = -dy
@@ -4238 +4226 @@
     volume_flow_l(2) = 0.0_wp
 
-    IF ( nest_bound_s )  THEN
+    IF ( bc_dirichlet_s )  THEN
        j = 0
        innor = dx
@@ -4250 +4238 @@
     ENDIF
 
-    IF ( nest_bound_n )  THEN
+    IF ( bc_dirichlet_n )  THEN
        j = ny + 1
        innor = -dx
@@ -4549 +4537 @@
 !
 !--       Left border pe:
-          IF ( nest_bound_l )  THEN
+          IF ( bc_dirichlet_l )  THEN
              
              CALL pmci_interp_tril_lr( u,  uc,  icu, jco, kco, r1xu, r2xu,     &
@@ -4656 +4644 @@
 !
 !--       Right border pe
-          IF ( nest_bound_r )  THEN
+          IF ( bc_dirichlet_r )  THEN
              
              CALL pmci_interp_tril_lr( u,  uc,  icu, jco, kco, r1xu, r2xu,     &
@@ -4769 +4757 @@
 !
 !--       South border pe
-          IF ( nest_bound_s )  THEN
+          IF ( bc_dirichlet_s )  THEN
 
              CALL pmci_interp_tril_sn( u,  uc,  icu, jco, kco, r1xu, r2xu,     &
@@ -4880 +4868 @@
 !
 !--       North border pe
-          IF ( nest_bound_n )  THEN
+          IF ( bc_dirichlet_n )  THEN
              
              CALL pmci_interp_tril_sn( u,  uc,  icu, jco, kco, r1xu, r2xu,     &
@@ -5748 +5736 @@
 !--    comprehensive explanation for this is still pending.
        IF ( nesting_mode == 'vertical' )  THEN
-          IF ( nest_bound_l )  THEN
+          IF ( bc_dirichlet_l )  THEN
              iclp = icl + nhll
           ENDIF
-          IF ( nest_bound_r ) THEN
+          IF ( bc_dirichlet_r ) THEN
              icrm = icr - nhlr
           ENDIF
-          IF ( nest_bound_s )  THEN
+          IF ( bc_dirichlet_s )  THEN
              jcsp = jcs + nhls
           ENDIF
-          IF ( nest_bound_n )  THEN
+          IF ( bc_dirichlet_n )  THEN
              jcnm = jcn - nhln
           ENDIF
        ELSE
-          IF ( nest_bound_l )  THEN
+          IF ( bc_dirichlet_l )  THEN
              IF ( var == 'u' )  THEN
                 iclp = icl + nhll + 1 + 1
@@ -5768 +5756 @@
              ENDIF
           ENDIF
-          IF ( nest_bound_r )  THEN
+          IF ( bc_dirichlet_r )  THEN
              icrm = icr - nhlr - 1
           ENDIF
 
-          IF ( nest_bound_s )  THEN
+          IF ( bc_dirichlet_s )  THEN
              IF ( var == 'v' )  THEN
                 jcsp = jcs + nhls + 1 + 1
@@ -5779 +5767 @@
              ENDIF
           ENDIF
-          IF ( nest_bound_n )  THEN
+          IF ( bc_dirichlet_n )  THEN
              jcnm = jcn - nhln - 1
           ENDIF

palm/trunk/SOURCE/pmc_parent_mod.f90

@@ -203 +203 @@
     INTEGER(iwp) ::  istat     !<
 
-
     DO  i = 1, SIZE( pmc_parent_for_child )-1
 
@@ -219 +218 @@
        CALL MPI_COMM_REMOTE_SIZE( children(childid)%inter_comm,                &
                                   children(childid)%inter_npes, istat )
-
 !
 !--    Intra communicator is used for MPI_GET
@@ -228 +226 @@
 
        ALLOCATE( children(childid)%pes(children(childid)%inter_npes))
-
 !
 !--    Allocate array of TYPE arraydef for all child PEs to store information
@@ -865 +862 @@
     DO
        CALL pmc_bcast( myname%couple_index, 0, comm=m_to_child_comm(childid) )
+
        IF ( myname%couple_index == -1 )  EXIT
+
        CALL pmc_bcast( myname%parentdesc,   0, comm=m_to_child_comm(childid) )
        CALL pmc_bcast( myname%nameonparent, 0, comm=m_to_child_comm(childid) )

palm/trunk/SOURCE/poismg_mod.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! Rename variables in mesoscale-offline nesting mode
 ! 
 ! Former revisions:
@@ -26 +26 @@
 ! $Id$
 ! Set lateral boundary conditions for divergence
-! 
+!
 ! 2937 2018-03-27 14:58:33Z suehring
 ! Corrected "Former revisions" section
@@ -107 +107 @@
  
     USE control_parameters,                                                    &
-        ONLY:  grid_level, force_bound_l, force_bound_n, force_bound_r,        &
-               force_bound_s, forcing, inflow_l, inflow_n, inflow_r, inflow_s, &
-               nest_bound_l, nest_bound_n, nest_bound_r, nest_bound_s,         &
-               outflow_l, outflow_n, outflow_r, outflow_s
+        ONLY:  bc_dirichlet_l, bc_dirichlet_n, bc_dirichlet_r,                 &
+               bc_dirichlet_s, bc_radiation_l, bc_radiation_n, bc_radiation_r, &
+               bc_radiation_s, grid_level, nesting_offline
 
     USE cpulog,                                                                &
@@ -206 +205 @@
 !--    Set lateral boundary conditions in non-cyclic case
        IF ( .NOT. bc_lr_cyc )  THEN
-          IF ( inflow_l .OR. outflow_l .OR. nest_bound_l .OR. force_bound_l )  &
+          IF ( bc_dirichlet_l  .OR.  bc_radiation_l )                          &
              d(:,:,nxl-1) = d(:,:,nxl)
-          IF ( inflow_r .OR. outflow_r .OR. nest_bound_r .OR. force_bound_r )  &
+          IF ( bc_dirichlet_r  .OR.  bc_radiation_r )                          &
              d(:,:,nxr+1) = d(:,:,nxr)
        ENDIF
        IF ( .NOT. bc_ns_cyc )  THEN
-          IF ( inflow_n .OR. outflow_n .OR. nest_bound_n .OR. force_bound_n )  &
+          IF ( bc_dirichlet_n  .OR.  bc_radiation_n )                          &
              d(:,nyn+1,:) = d(:,nyn,:)
-          IF ( inflow_s .OR. outflow_s .OR. nest_bound_s .OR. force_bound_s )  &
+          IF ( bc_dirichlet_s  .OR.  bc_radiation_s )                          &
              d(:,nys-1,:) = d(:,nys,:)
        ENDIF
@@ -375 +374 @@
 
        IF ( .NOT. bc_lr_cyc )  THEN
-          IF ( inflow_l  .OR.  outflow_l  .OR.  nest_bound_l  .OR.             &
-               force_bound_l )  THEN
+          IF ( bc_dirichlet_l  .OR.  bc_radiation_l )  THEN
              r(:,:,nxl_mg(l)-1) = r(:,:,nxl_mg(l))
           ENDIF
-          IF ( inflow_r  .OR.  outflow_r  .OR.  nest_bound_r  .OR.             &
-               force_bound_r )  THEN
+          IF ( bc_dirichlet_r  .OR.  bc_radiation_r )  THEN
              r(:,:,nxr_mg(l)+1) = r(:,:,nxr_mg(l))
           ENDIF
@@ -386 +383 @@
 
        IF ( .NOT. bc_ns_cyc )  THEN
-          IF ( inflow_n  .OR.  outflow_n  .OR.  nest_bound_n  .OR.             &
-               force_bound_n )  THEN
+          IF ( bc_dirichlet_n  .OR.  bc_radiation_n )  THEN
              r(:,nyn_mg(l)+1,:) = r(:,nyn_mg(l),:)
           ENDIF
-          IF ( inflow_s  .OR.  outflow_s  .OR.  nest_bound_s  .OR.             &
-               force_bound_s )  THEN
+          IF ( bc_dirichlet_s  .OR.  bc_radiation_s )  THEN
              r(:,nys_mg(l)-1,:) = r(:,nys_mg(l),:)
           ENDIF
@@ -510 +505 @@
 !--    Horizontal boundary conditions
        IF ( .NOT. bc_lr_cyc )  THEN
-          IF ( inflow_l  .OR.  outflow_l  .OR.  nest_bound_l  .OR.             &
-               force_bound_l )  THEN
+          IF ( bc_dirichlet_l  .OR.  bc_radiation_l )  THEN
              f_mg(:,:,nxl_mg(l)-1) = f_mg(:,:,nxl_mg(l))
           ENDIF
-          IF ( inflow_r  .OR.  outflow_r  .OR.  nest_bound_r  .OR.             &
-               force_bound_r )  THEN
+          IF ( bc_dirichlet_r  .OR.  bc_radiation_r )  THEN
              f_mg(:,:,nxr_mg(l)+1) = f_mg(:,:,nxr_mg(l))
           ENDIF
@@ -521 +514 @@
 
        IF ( .NOT. bc_ns_cyc )  THEN
-          IF ( inflow_n  .OR.  outflow_n  .OR.  nest_bound_n  .OR.             &
-               force_bound_n )  THEN
+          IF ( bc_dirichlet_n  .OR.  bc_radiation_n )  THEN
              f_mg(:,nyn_mg(l)+1,:) = f_mg(:,nyn_mg(l),:)
           ENDIF
-          IF ( inflow_s  .OR.  outflow_s  .OR.  nest_bound_s  .OR.             &
-               force_bound_s )  THEN
+          IF ( bc_dirichlet_s  .OR.  bc_radiation_s )  THEN
              f_mg(:,nys_mg(l)-1,:) = f_mg(:,nys_mg(l),:)
           ENDIF
@@ -681 +672 @@
 
        IF ( .NOT. bc_lr_cyc )  THEN
-          IF ( inflow_l  .OR.  outflow_l  .OR.  nest_bound_l  .OR.             &
-               force_bound_l )  THEN
+          IF ( bc_dirichlet_l  .OR.  bc_radiation_l )  THEN
              temp(:,:,nxl_mg(l)-1) = temp(:,:,nxl_mg(l))
           ENDIF
-          IF ( inflow_r  .OR.  outflow_r  .OR.  nest_bound_r  .OR.             &
-               force_bound_r )  THEN
+          IF ( bc_dirichlet_r  .OR.  bc_radiation_r )  THEN
              temp(:,:,nxr_mg(l)+1) = temp(:,:,nxr_mg(l))
           ENDIF
@@ -692 +681 @@
 
        IF ( .NOT. bc_ns_cyc )  THEN
-          IF ( inflow_n  .OR.  outflow_n  .OR.  nest_bound_n  .OR.             &
-               force_bound_n )  THEN
+          IF ( bc_dirichlet_n  .OR.  bc_radiation_n )  THEN
              temp(:,nyn_mg(l)+1,:) = temp(:,nyn_mg(l),:)
           ENDIF
-          IF ( inflow_s  .OR.  outflow_s  .OR.  nest_bound_s  .OR.             &
-               force_bound_s )  THEN
+          IF ( bc_dirichlet_s  .OR.  bc_radiation_s )  THEN
              temp(:,nys_mg(l)-1,:) = temp(:,nys_mg(l),:)
           ENDIF
@@ -991 +978 @@
 
              IF ( .NOT. bc_lr_cyc )  THEN
-                IF ( inflow_l  .OR.  outflow_l  .OR.  nest_bound_l  .OR.       &
-                     force_bound_l )  THEN
+                IF ( bc_dirichlet_l  .OR.  bc_radiation_l )  THEN
                    p_mg(:,:,nxl_mg(l)-1) = p_mg(:,:,nxl_mg(l))
                 ENDIF
-                IF ( inflow_r  .OR.  outflow_r  .OR.  nest_bound_r  .OR.       &
-                     force_bound_r )  THEN
+                IF ( bc_dirichlet_r  .OR.  bc_radiation_r )  THEN
                    p_mg(:,:,nxr_mg(l)+1) = p_mg(:,:,nxr_mg(l))
                 ENDIF
@@ -1002 +987 @@
 
              IF ( .NOT. bc_ns_cyc )  THEN
-                IF ( inflow_n  .OR.  outflow_n  .OR.  nest_bound_n  .OR.       &
-                     force_bound_n )  THEN
+                IF ( bc_dirichlet_n  .OR.  bc_radiation_n )  THEN
                    p_mg(:,nyn_mg(l)+1,:) = p_mg(:,nyn_mg(l),:)
                 ENDIF
-                IF ( inflow_s  .OR.  outflow_s  .OR.  nest_bound_s  .OR.       &
-                     force_bound_s )  THEN
+                IF ( bc_dirichlet_s  .OR.  bc_radiation_s )  THEN
                    p_mg(:,nys_mg(l)-1,:) = p_mg(:,nys_mg(l),:)
                 ENDIF
@@ -1423 +1406 @@
        USE control_parameters,                                                 &
            ONLY:  bc_lr_dirrad, bc_lr_raddir, bc_ns_dirrad, bc_ns_raddir,      &
-                  gamma_mg, grid_level_count, ibc_p_b, ibc_p_t,                &
+                  child_domain, gamma_mg, grid_level_count, ibc_p_b, ibc_p_t,  &
                   maximum_grid_level,  mg_switch_to_pe0_level,                 &
-                  mg_switch_to_pe0, nest_domain, ngsrb 
+                  mg_switch_to_pe0, ngsrb 
 
        USE indices,                                                            &
@@ -1564 +1547 @@
 !--          because then they have the total domain.
              IF ( bc_lr_dirrad )  THEN
-                inflow_l  = .TRUE.
-                inflow_r  = .FALSE.
-                outflow_l = .FALSE.
-                outflow_r = .TRUE.
+                bc_dirichlet_l  = .TRUE.
+                bc_dirichlet_r  = .FALSE.
+                bc_radiation_l = .FALSE.
+                bc_radiation_r = .TRUE.
              ELSEIF ( bc_lr_raddir )  THEN
-                inflow_l  = .FALSE.
-                inflow_r  = .TRUE.
-                outflow_l = .TRUE.
-                outflow_r = .FALSE.
-             ELSEIF ( nest_domain )  THEN
-                nest_bound_l = .TRUE.
-                nest_bound_r = .TRUE.
-             ELSEIF ( forcing )  THEN
-                force_bound_l = .TRUE.
-                force_bound_r = .TRUE.
+                bc_dirichlet_l  = .FALSE.
+                bc_dirichlet_r  = .TRUE.
+                bc_radiation_l = .TRUE.
+                bc_radiation_r = .FALSE.
+             ELSEIF ( child_domain  .OR.  nesting_offline )  THEN
+                bc_dirichlet_l = .TRUE.
+                bc_dirichlet_r = .TRUE.
              ENDIF
 
              IF ( bc_ns_dirrad )  THEN
-                inflow_n  = .TRUE.
-                inflow_s  = .FALSE.
-                outflow_n = .FALSE.
-                outflow_s = .TRUE.
+                bc_dirichlet_n  = .TRUE.
+                bc_dirichlet_s  = .FALSE.
+                bc_radiation_n = .FALSE.
+                bc_radiation_s = .TRUE.
              ELSEIF ( bc_ns_raddir )  THEN
-                inflow_n  = .FALSE.
-                inflow_s  = .TRUE.
-                outflow_n = .TRUE.
-                outflow_s = .FALSE.
-             ELSEIF ( nest_domain )  THEN
-                nest_bound_s = .TRUE.
-                nest_bound_n = .TRUE.
-             ELSEIF ( forcing )  THEN
-                force_bound_s = .TRUE.
-                force_bound_n = .TRUE.
+                bc_dirichlet_n  = .FALSE.
+                bc_dirichlet_s  = .TRUE.
+                bc_radiation_n = .TRUE.
+                bc_radiation_s = .FALSE.
+             ELSEIF ( child_domain  .OR.  nesting_offline)  THEN
+                bc_dirichlet_s = .TRUE.
+                bc_dirichlet_n = .TRUE.
              ENDIF
 
@@ -1659 +1636 @@
 !--          For non-cyclic lateral boundary conditions and in case of nesting, 
 !--          restore the in-/outflow conditions. 
-             inflow_l  = .FALSE.;  inflow_r  = .FALSE.
-             inflow_n  = .FALSE.;  inflow_s  = .FALSE.
-             outflow_l = .FALSE.;  outflow_r = .FALSE.
-             outflow_n = .FALSE.;  outflow_s = .FALSE.
-!
-!--          In case of nesting or forcing, restore lateral boundary conditions
-             IF ( nest_domain )  THEN
-                nest_bound_l = .FALSE.
-                nest_bound_r = .FALSE.
-                nest_bound_s = .FALSE.
-                nest_bound_n = .FALSE.      
-             ENDIF
-             IF ( forcing )  THEN
-                force_bound_l = .FALSE.
-                force_bound_r = .FALSE.
-                force_bound_s = .FALSE.
-                force_bound_n = .FALSE.      
-             ENDIF
+             bc_dirichlet_l = .FALSE.;  bc_dirichlet_r = .FALSE.
+             bc_dirichlet_n = .FALSE.;  bc_dirichlet_s = .FALSE.
+             bc_radiation_l = .FALSE.;  bc_radiation_r = .FALSE.
+             bc_radiation_n = .FALSE.;  bc_radiation_s = .FALSE.
 
              IF ( pleft == MPI_PROC_NULL )  THEN
-                IF ( bc_lr_dirrad )  THEN
-                   inflow_l  = .TRUE.
+                IF ( bc_lr_dirrad  .OR.  child_domain  .OR.  nesting_offline ) &
+                THEN
+                   bc_dirichlet_l = .TRUE.
                 ELSEIF ( bc_lr_raddir )  THEN
-                   outflow_l = .TRUE.
-                ELSEIF ( nest_domain )  THEN
-                   nest_bound_l = .TRUE.
-                ELSEIF ( forcing )  THEN
-                   force_bound_l = .TRUE.
+                   bc_radiation_l = .TRUE.
                 ENDIF
              ENDIF
@@ -1692 +1652 @@
              IF ( pright == MPI_PROC_NULL )  THEN
                 IF ( bc_lr_dirrad )  THEN
-                   outflow_r = .TRUE.
-                ELSEIF ( bc_lr_raddir )  THEN
-                   inflow_r  = .TRUE.
-                ELSEIF ( nest_domain )  THEN
-                   nest_bound_r = .TRUE.
-                ELSEIF ( forcing )  THEN
-                   force_bound_r = .TRUE.
+                   bc_radiation_r = .TRUE.
+                ELSEIF ( bc_lr_raddir  .OR.  child_domain  .OR.                &
+                         nesting_offline )  THEN
+                   bc_dirichlet_r = .TRUE.
                 ENDIF
              ENDIF
@@ -1704 +1661 @@
              IF ( psouth == MPI_PROC_NULL )  THEN
                 IF ( bc_ns_dirrad )  THEN
-                   outflow_s = .TRUE.
+                   bc_radiation_s = .TRUE.
+                ELSEIF ( bc_ns_raddir  .OR.  child_domain  .OR.                &
+                         nesting_offline )  THEN
+                   bc_dirichlet_s = .TRUE.
+                ENDIF
+             ENDIF
+
+             IF ( pnorth == MPI_PROC_NULL )  THEN
+                IF ( bc_ns_dirrad  .OR.  child_domain  .OR.  nesting_offline ) &
+                THEN
+                   bc_dirichlet_n = .TRUE.
                 ELSEIF ( bc_ns_raddir )  THEN
-                   inflow_s  = .TRUE.
-                ELSEIF ( nest_domain )  THEN
-                   nest_bound_s = .TRUE.
-                ELSEIF ( forcing )  THEN
-                   force_bound_s = .TRUE.
-                ENDIF
-             ENDIF
-
-             IF ( pnorth == MPI_PROC_NULL )  THEN
-                IF ( bc_ns_dirrad )  THEN
-                   inflow_n  = .TRUE.
-                ELSEIF ( bc_ns_raddir )  THEN
-                   outflow_n = .TRUE.
-                ELSEIF ( nest_domain )  THEN
-                   nest_bound_n = .TRUE.
-                ELSEIF ( forcing )  THEN
-                   force_bound_n = .TRUE.
+                   bc_radiation_n = .TRUE.
                 ENDIF
              ENDIF

palm/trunk/SOURCE/poismg_noopt_mod.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! Rename variables in mesoscale-offline nesting mode
 ! 
 ! Former revisions:
@@ -165 +165 @@
  
     USE control_parameters,                                                    &
-        ONLY:  grid_level, force_bound_l, force_bound_n, force_bound_r,        &
-               force_bound_s, forcing, inflow_l, inflow_n, inflow_r, inflow_s, &
-               nest_bound_l, nest_bound_n, nest_bound_r, nest_bound_s,         &
-               outflow_l, outflow_n, outflow_r, outflow_s
+        ONLY:  bc_dirichlet_l, bc_dirichlet_n, bc_dirichlet_r,                 &
+               bc_dirichlet_s, bc_radiation_l, bc_radiation_n, bc_radiation_r, &
+               bc_radiation_s, child_domain, grid_level, nesting_offline
 
     USE cpulog,                                                                &
@@ -261 +260 @@
 !--    Set lateral boundary conditions in non-cyclic case
        IF ( .NOT. bc_lr_cyc )  THEN
-          IF ( inflow_l .OR. outflow_l .OR. nest_bound_l .OR. force_bound_l )  &
+          IF ( bc_dirichlet_l   .OR.  bc_radiation_l )                         &
              d(:,:,nxl-1) = d(:,:,nxl)
-          IF ( inflow_r .OR. outflow_r .OR. nest_bound_r .OR. force_bound_r )  &
+          IF ( bc_dirichlet_r   .OR.  bc_radiation_r )                         &
              d(:,:,nxr+1) = d(:,:,nxr)
        ENDIF
        IF ( .NOT. bc_ns_cyc )  THEN
-          IF ( inflow_n .OR. outflow_n .OR. nest_bound_n .OR. force_bound_n )  &
+          IF ( bc_dirichlet_n   .OR.  bc_radiation_n )                         &
              d(:,nyn+1,:) = d(:,nyn,:)
-          IF ( inflow_s .OR. outflow_s .OR. nest_bound_s .OR. force_bound_s )  &
+          IF ( bc_dirichlet_s   .OR.  bc_radiation_s )                         &
              d(:,nys-1,:) = d(:,nys,:)
        ENDIF
@@ -442 +441 @@
 
        IF ( .NOT. bc_lr_cyc )  THEN
-          IF ( inflow_l .OR. outflow_l .OR. nest_bound_l .OR. force_bound_l )  THEN
+          IF ( bc_dirichlet_l  .OR.  bc_radiation_l )  THEN
              r(:,:,nxl_mg(l)-1) = r(:,:,nxl_mg(l))
           ENDIF
-          IF ( inflow_r .OR. outflow_r .OR. nest_bound_r .OR. force_bound_r )  THEN
+          IF ( bc_dirichlet_r  .OR.  bc_radiation_r )  THEN
              r(:,:,nxr_mg(l)+1) = r(:,:,nxr_mg(l))