Changeset 707 for palm/trunk/SOURCE/init_pegrid.f90

Timestamp:

Mar 29, 2011 11:39:40 AM (13 years ago)

Author:

raasch

Message:

New:
---

In case of multigrid method, on coarse grid levels, gathered data are
identically processed on all PEs (before, on PE0 only), so that the subsequent
scattering of data is not neccessary any more. (modules, init_pegrid, poismg)

Changed:

---

Calculation of weighted average of p is now handled in the same way
regardless of the number of ghost layers (advection scheme). (pres)

multigrid and sor method are using p_loc for iterative
advancements of pressure. p_sub removed. (init_3d_model, modules, poismg, pres, sor)

bc_lr and bc_ns replaced by bc_lr_dirrad, bc_lr_raddir, bc_ns_dirrad, bc_ns_raddir
for speed optimization. (calc_spectra, check_parameters, exchange_horiz,
exchange_horiz_2d, header, init_3d_model, init_grid, init_pegrid, modules,
poismg, pres, sor, time_integration, timestep)

grid_level directly used as index for MPI data type arrays. (exchange_horiz,
poismg)

initial assignments of zero to array p for iterative solvers only (init_3d_model)

Errors:

---

localsum calculation modified for proper OpenMP reduction. (pres)

Bugfix: bottom (nzb) and top (nzt+1) boundary conditions set in routines
resid and restrict. They were missed before, which may have led to
unpredictable results. (poismg)

File:

• palm/trunk/SOURCE/init_pegrid.f90 (modified) (8 diffs)

palm/trunk/SOURCE/init_pegrid.f90

@@ -4 +4 @@
 ! Current revisions:
 ! -----------------
+! bc_lr/ns replaced by bc_lr/ns_cyc/dirrad/raddir
 ! 
 ! ATTENTION: nnz_x undefined problem still has to be solved!!!!!!!!
@@ -186 +187 @@
 !
 !-- If necessary, set horizontal boundary conditions to non-cyclic
-    IF ( bc_lr /= 'cyclic' )  cyclic(1) = .FALSE.
-    IF ( bc_ns /= 'cyclic' )  cyclic(2) = .FALSE.
+    IF ( .NOT. bc_lr_cyc )  cyclic(1) = .FALSE.
+    IF ( .NOT. bc_ns_cyc )  cyclic(2) = .FALSE.
 
 !
@@ -331 +332 @@
     nzt  = MIN( nz, nzta )
     nnz  = nza
+
+!
+!-- Set switches to define if the PE is situated at the border of the virtual
+!-- processor grid
+    IF ( nxl == 0 )   left_border_pe  = .TRUE.
+    IF ( nxr == nx )  right_border_pe = .TRUE.
+    IF ( nys == 0 )   south_border_pe = .TRUE.
+    IF ( nyn == ny )  north_border_pe = .TRUE.
 
 !
@@ -1066 +1075 @@
 !-- horizontal boundary conditions.
     IF ( pleft == MPI_PROC_NULL )  THEN
-       IF ( bc_lr == 'dirichlet/radiation' )  THEN
+       IF ( bc_lr_dirrad )  THEN
           inflow_l  = .TRUE.
-       ELSEIF ( bc_lr == 'radiation/dirichlet' )  THEN
+       ELSEIF ( bc_lr_raddir )  THEN
           outflow_l = .TRUE.
        ENDIF
@@ -1074 +1083 @@
 
     IF ( pright == MPI_PROC_NULL )  THEN
-       IF ( bc_lr == 'dirichlet/radiation' )  THEN
+       IF ( bc_lr_dirrad )  THEN
           outflow_r = .TRUE.
-       ELSEIF ( bc_lr == 'radiation/dirichlet' )  THEN
+       ELSEIF ( bc_lr_raddir )  THEN
           inflow_r  = .TRUE.
        ENDIF
@@ -1082 +1091 @@
 
     IF ( psouth == MPI_PROC_NULL )  THEN
-       IF ( bc_ns == 'dirichlet/radiation' )  THEN
+       IF ( bc_ns_dirrad )  THEN
           outflow_s = .TRUE.
-       ELSEIF ( bc_ns == 'radiation/dirichlet' )  THEN
+       ELSEIF ( bc_ns_raddir )  THEN
           inflow_s  = .TRUE.
        ENDIF
@@ -1090 +1099 @@
 
     IF ( pnorth == MPI_PROC_NULL )  THEN
-       IF ( bc_ns == 'dirichlet/radiation' )  THEN
+       IF ( bc_ns_dirrad )  THEN
           inflow_n  = .TRUE.
-       ELSEIF ( bc_ns == 'radiation/dirichlet' )  THEN
+       ELSEIF ( bc_ns_raddir )  THEN
           outflow_n = .TRUE.
        ENDIF
@@ -1122 +1131 @@
 
 #else
-    IF ( bc_lr == 'dirichlet/radiation' )  THEN
+    IF ( bc_lr_dirrad )  THEN
        inflow_l  = .TRUE.
        outflow_r = .TRUE.
-    ELSEIF ( bc_lr == 'radiation/dirichlet' )  THEN
+    ELSEIF ( bc_lr_raddir )  THEN
        outflow_l = .TRUE.
        inflow_r  = .TRUE.
     ENDIF
 
-    IF ( bc_ns == 'dirichlet/radiation' )  THEN
+    IF ( bc_ns_dirrad )  THEN
        inflow_n  = .TRUE.
        outflow_s = .TRUE.
-    ELSEIF ( bc_ns == 'radiation/dirichlet' )  THEN
+    ELSEIF ( bc_ns_raddir )  THEN
        outflow_n = .TRUE.
        inflow_s  = .TRUE.