Changeset 4414 for palm/trunk/SOURCE/advec_ws.f90

Timestamp:

Feb 19, 2020 8:16:04 PM (5 years ago)

Author:

suehring

Message:

Remove deprecated topography arrays; Move basic initialization of numerics into an extra module interface

File:

• palm/trunk/SOURCE/advec_ws.f90 (modified) (5 diffs)

palm/trunk/SOURCE/advec_ws.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Move call for initialization of control flags to ws_init
+! 
+! 4360 2020-01-07 11:25:50Z suehring
 ! Introduction of wall_flags_total_0, which currently sets bits based on static
 ! topography information used in wall_flags_static_0
@@ -379 +382 @@
 
        ENDIF
+!
+!--    Initialize the flag arrays controlling degradation near walls, i.e. 
+!--    to decrease the numerical stencil appropriately. The order of the scheme
+!--    is degraded near solid walls as well as near non-cyclic inflow and outflow
+!--    boundaries. Do this separately for momentum and scalars.
+       IF ( ws_scheme_mom )  CALL ws_init_flags_momentum
+
+       IF ( ws_scheme_sca )  THEN
+          ALLOCATE( advc_flags_s(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
+          advc_flags_s = 0
+          CALL ws_init_flags_scalar( bc_dirichlet_l  .OR.  bc_radiation_l,     &
+                                     bc_dirichlet_n  .OR.  bc_radiation_n,     &
+                                     bc_dirichlet_r  .OR.  bc_radiation_r,     &
+                                     bc_dirichlet_s  .OR.  bc_radiation_s,     &
+                                     advc_flags_s )
+       ENDIF
 
     END SUBROUTINE ws_init
@@ -398 +417 @@
        INTEGER(iwp) ::  k_pp  !< dummy index along z
        INTEGER(iwp) ::  k_ppp !< dummy index along z
-       
+
        LOGICAL      ::  flag_set !< steering variable for advection flags
-   
+
+       ALLOCATE( advc_flags_m(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
        advc_flags_m = 0
-
 !
 !--    Set advc_flags_m to steer the degradation of the advection scheme in advec_ws
@@ -850 +869 @@
        INTEGER(iwp) ::  k_pp  !< dummy index along z
        INTEGER(iwp) ::  k_ppp !< dummy index along z
-       
+
        INTEGER(iwp), INTENT(INOUT), DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ::&
                                                   advc_flag !< flag array to control order of scalar advection
-       
+
        LOGICAL ::  flag_set     !< steering variable for advection flags
        LOGICAL ::  non_cyclic_l !< flag that indicates non-cyclic boundary on the left
@@ -859 +878 @@
        LOGICAL ::  non_cyclic_r !< flag that indicates non-cyclic boundary on the right
        LOGICAL ::  non_cyclic_s !< flag that indicates non-cyclic boundary on the south
-       
+
        LOGICAL, OPTIONAL ::  extensive_degrad !< flag indicating that extensive degradation is required, e.g. for 
                                               !< passive scalars nearby topography along the horizontal directions,