Changeset 3012 for palm/trunk/SOURCE

Timestamp:

May 8, 2018 11:01:29 AM (7 years ago)

Author:

Giersch

Message:

To avoid jumps while plotting w-profiles w level nzt+1 is set to w level nzt after velocity modifications through the pressure solver were carried out

File:

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

palm/trunk/SOURCE/pres.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-! 
+! To avoid jumps while plotting profiles w level nzt+1 is set to w level nzt 
+! after velocity modifications through the pressure solver are carried out
 ! 
 ! Former revisions:
@@ -154 +155 @@
 
     USE control_parameters,                                                    &
-        ONLY:  bc_lr_cyc, bc_ns_cyc, conserve_volume_flow, dt_3d,              &
-               gathered_size, ibc_p_b, ibc_p_t, intermediate_timestep_count,   &
-               intermediate_timestep_count_max, mg_switch_to_pe0_level,        &
-               nest_domain, outflow_l, outflow_n, outflow_r,                   &
-               outflow_s, psolver, subdomain_size, topography, volume_flow,    &
-               volume_flow_area, volume_flow_initial
+        ONLY:  bc_lr_cyc, bc_ns_cyc, conserve_volume_flow, coupling_mode,      &
+               dt_3d, gathered_size, ibc_p_b, ibc_p_t,                         &
+               intermediate_timestep_count, intermediate_timestep_count_max,   &
+               mg_switch_to_pe0_level, nest_domain, outflow_l, outflow_n,      &
+               outflow_r, outflow_s, psolver, subdomain_size, topography,      &
+               volume_flow, volume_flow_area, volume_flow_initial
 
     USE cpulog,                                                                &
@@ -745 +746 @@
 
 !
+!-- The vertical velocity is not set to zero at nzt + 1 for nested domains
+!-- Instead it is set to the values of nzt (see routine vnest_boundary_conds 
+!-- or pmci_interp_tril_t) BEFORE calling the pressure solver. To avoid jumps
+!-- while plotting profiles w at the top has to be set to the values in the 
+!-- height nzt after above modifications. Hint: w level nzt+1 does not impact
+!-- results. 
+    IF (nest_domain .OR. coupling_mode == 'vnested_fine') THEN
+       w(nzt+1,:,:) = w(nzt,:,:)
+    ENDIF
+
+!
 !-- Sum up the volume flow through the right and north boundary
     IF ( conserve_volume_flow  .AND.  bc_lr_cyc  .AND.  bc_ns_cyc  .AND.       &