Changeset 1762 for palm/trunk/SOURCE/boundary_conds.f90

Timestamp:

Feb 25, 2016 12:31:13 PM (8 years ago)

Author:

hellstea

Message:

Introduction of nested domain system

File:

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

palm/trunk/SOURCE/boundary_conds.f90

@@ -19 +19 @@
 ! Current revisions:
 ! -----------------
-!
+! Introduction of nested domain feature
 !
 ! Former revisions:
@@ -134 +134 @@
                ibc_pt_b, ibc_pt_t, ibc_q_b, ibc_q_t, ibc_sa_t, ibc_uv_b,       &
                ibc_uv_t, icloud_scheme, inflow_l, inflow_n, inflow_r, inflow_s,&
-               intermediate_timestep_count, large_scale_forcing, ocean,        &
+               intermediate_timestep_count, large_scale_forcing, nest_domain,  &
+               nest_bound_l, nest_bound_s, nudging, ocean,                     &
                outflow_l, outflow_n, outflow_r, outflow_s, passive_scalar,     &
-               precipitation, tsc, use_cmax, &
-               nudging
+               precipitation, tsc, use_cmax
 
     USE grid_variables,                                                        &
@@ -179 +179 @@
 
 !
-!-- Top boundary
+!-- Top boundary. A nested domain ( ibc_uv_t = 3 ) does not require settings.
     IF ( ibc_uv_t == 0 )  THEN
        !$acc kernels present( u_init, u_p, v_init, v_p )
@@ -185 +185 @@
         v_p(nzt+1,:,:) = v_init(nzt+1)
        !$acc end kernels
-    ELSE
+    ELSEIF ( ibc_uv_t == 1 )  THEN
        !$acc kernels present( u_p, v_p )
         u_p(nzt+1,:,:) = u_p(nzt,:,:)
@@ -191 +191 @@
        !$acc end kernels
     ENDIF
-    !$acc kernels present( w_p )
-    w_p(nzt:nzt+1,:,:) = 0.0_wp  ! nzt is not a prognostic level (but cf. pres)
-    !$acc end kernels
+
+    IF ( .NOT. nest_domain )  THEN
+       !$acc kernels present( w_p )
+       w_p(nzt:nzt+1,:,:) = 0.0_wp  ! nzt is not a prognostic level (but cf. pres)
+       !$acc end kernels
+    ENDIF
 
 !
@@ -255 +258 @@
           ENDDO
        ENDDO
-       e_p(nzt+1,:,:) = e_p(nzt,:,:)
+       IF ( .NOT. nest_domain )  THEN
+          e_p(nzt+1,:,:) = e_p(nzt,:,:)
+       ENDIF
        !$acc end kernels
     ENDIF
@@ -325 +330 @@
     ENDIF
 !
-!-- In case of inflow at the south boundary the boundary for v is at nys
-!-- and in case of inflow at the left boundary the boundary for u is at nxl.
-!-- Since in prognostic_equations (cache optimized version) these levels are
-!-- handled as a prognostic level, boundary values have to be restored here.
+!-- In case of inflow or nest boundary at the south boundary the boundary for v
+!-- is at nys and in case of inflow or nest boundary at the left boundary the
+!-- boundary for u is at nxl. Since in prognostic_equations (cache optimized
+!-- version) these levels are handled as a prognostic level, boundary values
+!-- have to be restored here.
 !-- For the SGS-TKE, Neumann boundary conditions are used at the inflow.
     IF ( inflow_s )  THEN
@@ -340 +346 @@
     ELSEIF ( inflow_r )  THEN
        IF ( .NOT. constant_diffusion ) e_p(:,:,nxr+1) = e_p(:,:,nxr)
+    ENDIF
+
+!
+!-- The same restoration for u at i=nxl and v at j=nys as above must be made
+!-- in case of nest boundaries. Note however, that the above ELSEIF-structure is
+!-- not appropriate here as there may be more than one nest boundary on a
+!-- PE-domain. Furthermore Neumann conditions for SGS-TKE are not required here.
+    IF ( nest_bound_s )  THEN
+       v_p(:,nys,:) = v_p(:,nys-1,:)
+    ENDIF
+    IF ( nest_bound_l )  THEN
+       u_p(:,:,nxl) = u_p(:,:,nxl-1)
     ENDIF