Changeset 3900 for palm/trunk/SOURCE/synthetic_turbulence_generator_mod.f90

Timestamp:

Apr 16, 2019 3:17:43 PM (5 years ago)

Author:

suehring

Message:

Bugfixes in initialization and STG

File:

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

palm/trunk/SOURCE/synthetic_turbulence_generator_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Missing re-calculation of perturbation seeds in case of restarts
+! 
+! 3891 2019-04-12 17:52:01Z suehring
 ! Bugfix in initialization in case of restart runs. 
 ! 
@@ -899 +902 @@
                 ELSE
                    fw_yz = 0._wp
+                ENDIF
+
+             ENDDO
+          ENDDO
+       ENDIF
+       
+       IF ( myidy == id_stg_south  .OR.  myidy == id_stg_north )  THEN
+
+          IF ( myidy == id_stg_south )  j = -1
+          IF ( myidy == id_stg_north )  j = nyn+1
+
+          DO  i = nxlg, nxrg
+             DO  k = nzb, nzt+1
+
+                IF  ( a11(k) .NE. 0._wp ) THEN
+                   fu_xz(k,i) = ( u(k,j,i) / mc_factor - u_init(k) ) / a11(k)
+                ELSE
+                   fu_xz(k,i) = 0._wp
+                ENDIF
+
+                IF  ( a22(k) .NE. 0._wp ) THEN
+                   fv_xz(k,i) = ( v(k,j,i) / mc_factor - a21(k) * fu_yz(k,i) - &
+                               v_init(k) ) / a22(k)
+                ELSE
+                   fv_xz(k,i) = 0._wp
+                ENDIF
+
+                IF  ( a33(k) .NE. 0._wp ) THEN
+                   fw_xz(k,i) = ( w(k,j,i) / mc_factor - a31(k) * fu_yz(k,i) - &
+                               a32(k) * fv_yz(k,i) ) / a33(k)
+                ELSE
+                   fw_xz = 0._wp
                 ENDIF