Changeset 4603

Timestamp:

Jul 14, 2020 4:08:30 PM (4 years ago)

Author:

suehring

Message:

Synthetic turbulence generator: Bugfix in initialization from ASCII file - x-length scales at the bottom boundary were not initialized properly

File:

• palm/trunk/SOURCE/synthetic_turbulence_generator_mod.f90 (modified) (6 diffs)

palm/trunk/SOURCE/synthetic_turbulence_generator_mod.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
-! 
+!
+!
 ! Former revisions:
 ! -----------------
 ! $Id$
+! Bugfix in initialization from ASCII file - x-length scales at the bottom boundary were not
+! initialized properly
+! 
+! 4566 2020-06-16 10:11:51Z suehring
 ! - revise parametrization for reynolds-stress components, turbulent length- and time scales
 ! - revise computation of velocity disturbances to be consistent to Xie and Castro (2008)
@@ -813 +817 @@
        READ( 90, * )
 
-       DO  k = nzb+1, nzt+1
+       DO  k = nzb, nzt
           READ( 90, * ) zz, luy, luz, tu(k), lvy, lvz, tv(k), lwy, lwz, tw(k), r11(k), r21(k),     &
                         r22(k), r31(k), r32(k), r33(k), d1, d2, d3, d5
@@ -819 +823 @@
 !
 !--       Convert length scales from meter to number of grid points.
+          IF ( k /= nzb )  THEN
+             nuz(k) = INT( luz * ddzw(k) )
+             nvz(k) = INT( lvz * ddzw(k) )
+             nwz(k) = INT( lwz * ddzw(k) )
+          ELSE
+             nuz(k) = INT( luz * ddzw(k+1) )
+             nvz(k) = INT( lvz * ddzw(k+1) )
+             nwz(k) = INT( lwz * ddzw(k+1) )
+          ENDIF
+
           nuy(k) = INT( luy * ddy )
-          nuz(k) = INT( luz * ddzw(k) )
           nvy(k) = INT( lvy * ddy )
-          nvz(k) = INT( lvz * ddzw(k) )
           nwy(k) = INT( lwy * ddy )
-          nwz(k) = INT( lwz * ddzw(k) )
-!
-!--       Workaround, assume isotropic turbulence
+!
+!--       Set length scales in x-direction. As a workaround assume isotropic turbulence in x- and
+!--       y-direction.
           nwx(k) = nwy(k)
           nvx(k) = nvy(k)
@@ -840 +852 @@
        ENDDO
 !
-!--    Set lenght scales at surface grid point
+!--    Set length scales at the surface and top boundary. At the surface the lengths scales are
+!--    simply overwritten.
+       nwx(nzb) = nwy(nzb+1)
+       nvx(nzb) = nvy(nzb+1)
+       nux(nzb) = nuy(nzb+1)
        nuy(nzb) = nuy(nzb+1)
        nuz(nzb) = nuz(nzb+1)
@@ -848 +864 @@
        nwz(nzb) = nwz(nzb+1)
 
+       nwx(nzt+1) = nwy(nzt)
+       nvx(nzt+1) = nvy(nzt)
+       nux(nzt+1) = nuy(nzt)
+       nuy(nzt+1) = nuy(nzt)
+       nuz(nzt+1) = nuz(nzt)
+       nvy(nzt+1) = nvy(nzt)
+       nvz(nzt+1) = nvz(nzt)
+       nwy(nzt+1) = nwy(nzt)
+       nwz(nzt+1) = nwz(nzt)
+
        CLOSE( 90 )
 !
@@ -887 +913 @@
 !-- Initial filter functions
     CALL stg_setup_filter_function
-
 !
 !-- Allocate velocity seeds for turbulence at xz-layer