Changeset 3937 for palm/trunk/SOURCE/nesting_offl_mod.f90

Timestamp:

Apr 29, 2019 3:09:07 PM (5 years ago)

Author:

suehring

Message:

Move initialization of STG behind initialization of offline nesting; Bugfix in STG in case of very early restart; calculation of scaling parameters used for parametrization of synthetic turbulence profiles improved; in offline nesting, set boundary value at upper-left and upper-south grid points for u- and v-component, respectively

File:

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

palm/trunk/SOURCE/nesting_offl_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Set boundary conditon on upper-left and upper-south grid point for the u- and
+! v-component, respectively.
+! 
+! 3891 2019-04-12 17:52:01Z suehring
 ! Bugfix, do not overwrite lateral and top boundary data in case of restart 
 ! runs. 
@@ -633 +637 @@
           ENDDO
        ENDDO
+!
+!--    For left boundary set boundary condition for u-component also at top 
+!--    grid point. 
+!--    Note, this has no effect on the numeric solution, only for data output.
+       IF ( bc_dirichlet_l )  u(nzt+1,:,nxl) = u(nzt+1,:,nxlu)
 
        DO  i = nxl, nxr
@@ -644 +653 @@
           ENDDO
        ENDDO
+!
+!--    For south boundary set boundary condition for v-component also at top 
+!--    grid point. 
+!--    Note, this has no effect on the numeric solution, only for data output.
+       IF ( bc_dirichlet_s )  v(nzt+1,nys,:) = v(nzt+1,nysv,:)
 
        DO  i = nxl, nxr
@@ -961 +975 @@
        topo_max     = topo_max_l
 #endif
-
        zi_ribulk = MAX( zi_ribulk, topo_max )
-
+       
     END SUBROUTINE calc_zi