Ignore:
Timestamp:
Apr 16, 2019 3:17:43 PM (2 years ago)
Author:
suehring
Message:

Bugfixes in initialization and STG

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/synthetic_turbulence_generator_mod.f90

    r3891 r3900  
    2525! -----------------
    2626! $Id$
     27! Missing re-calculation of perturbation seeds in case of restarts
     28!
     29! 3891 2019-04-12 17:52:01Z suehring
    2730! Bugfix in initialization in case of restart runs.
    2831!
     
    899902                ELSE
    900903                   fw_yz = 0._wp
     904                ENDIF
     905
     906             ENDDO
     907          ENDDO
     908       ENDIF
     909       
     910       IF ( myidy == id_stg_south  .OR.  myidy == id_stg_north )  THEN
     911
     912          IF ( myidy == id_stg_south )  j = -1
     913          IF ( myidy == id_stg_north )  j = nyn+1
     914
     915          DO  i = nxlg, nxrg
     916             DO  k = nzb, nzt+1
     917
     918                IF  ( a11(k) .NE. 0._wp ) THEN
     919                   fu_xz(k,i) = ( u(k,j,i) / mc_factor - u_init(k) ) / a11(k)
     920                ELSE
     921                   fu_xz(k,i) = 0._wp
     922                ENDIF
     923
     924                IF  ( a22(k) .NE. 0._wp ) THEN
     925                   fv_xz(k,i) = ( v(k,j,i) / mc_factor - a21(k) * fu_yz(k,i) - &
     926                               v_init(k) ) / a22(k)
     927                ELSE
     928                   fv_xz(k,i) = 0._wp
     929                ENDIF
     930
     931                IF  ( a33(k) .NE. 0._wp ) THEN
     932                   fw_xz(k,i) = ( w(k,j,i) / mc_factor - a31(k) * fu_yz(k,i) - &
     933                               a32(k) * fv_yz(k,i) ) / a33(k)
     934                ELSE
     935                   fw_xz = 0._wp
    901936                ENDIF
    902937
Note: See TracChangeset for help on using the changeset viewer.