Changeset 3347 for palm/trunk/SOURCE/pres.f90

Timestamp:

Oct 15, 2018 2:21:08 PM (6 years ago)

Author:

suehring

Message:

Offline nesting revised and separated from large_scale_forcing_mod; Time-dependent synthetic turbulence generator; bugfixes in USM/LSM radiation model and init_pegrid

File:

• palm/trunk/SOURCE/pres.f90 (modified) (4 diffs)

palm/trunk/SOURCE/pres.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Bugfixes in offline nesting.
+! Add comment. 
+! 
+! 3341 2018-10-15 10:31:27Z suehring
 ! Rename variables for boundary flags and nesting
 ! 
@@ -169 +173 @@
                dt_3d, gathered_size, ibc_p_b, ibc_p_t,                         &
                intermediate_timestep_count, intermediate_timestep_count_max,   &
-               mg_switch_to_pe0_level, psolver, subdomain_size,                &
+               mg_switch_to_pe0_level, nesting_offline,                        &
+               psolver, subdomain_size,                                        &
                topography, volume_flow, volume_flow_area, volume_flow_initial
 
@@ -386 +391 @@
     ENDIF
 
-    IF ( ibc_p_b == 1  .AND.  ibc_p_t == 1  .AND.                               &
+    IF ( ibc_p_b == 1  .AND.  ibc_p_t == 1  .AND.  .NOT. nesting_offline  .AND. &
          .NOT. child_domain_nvn  .AND. intermediate_timestep_count /= 0 )       &
     THEN
@@ -719 +724 @@
        volume_flow_l(2) = 0.0_wp
     ENDIF
-
+!
+!-- Add pressure gradients to the velocity components. Note, the loops are
+!-- running over the entire model domain, even though, in case of non-cyclic 
+!-- boundaries u- and v-component are not prognostic at i=0 and j=0, 
+!-- respectiveley. However, in case of Dirichlet boundary conditions for the 
+!-- velocities, zero-gradient conditions for the pressure are set, so that
+!-- no modification is imposed at the boundaries. 
     !$OMP PARALLEL PRIVATE (i,j,k)
     !$OMP DO