Changeset 106 for palm/trunk/SOURCE/init_3d_model.f90

Timestamp:

Aug 16, 2007 2:30:26 PM (17 years ago)

Author:

raasch

Message:

preliminary update of bugfixes and extensions for non-cyclic BCs

File:

• palm/trunk/SOURCE/init_3d_model.f90 (modified) (3 diffs)

palm/trunk/SOURCE/init_3d_model.f90

@@ -7 +7 @@
 ! Actual revisions:
 ! -----------------
-! Flux initialization in case of coupled runs, +momentum fluxes at top boundary
+! Flux initialization in case of coupled runs, +momentum fluxes at top boundary,
+! +arrays for phase speed c_u, c_v, c_w, indices for u|v|w_m_l|r changed
 !
 ! Former revisions:
@@ -233 +234 @@
 
 !
-!-- Arrays to store velocity data from t-dt needed for radiation boundary
-!-- conditions
+!-- Arrays to store velocity data from t-dt and the phase speeds which
+!-- are needed for radiation boundary conditions
     IF ( outflow_l )  THEN
-       ALLOCATE( u_m_l(nzb:nzt+1,nys-1:nyn+1,-1:1), &
-                 v_m_l(nzb:nzt+1,nys-1:nyn+1,-1:1), &
-                 w_m_l(nzb:nzt+1,nys-1:nyn+1,-1:1) )
+       ALLOCATE( u_m_l(nzb:nzt+1,nys-1:nyn+1,1:2), &
+                 v_m_l(nzb:nzt+1,nys-1:nyn+1,0:1), &
+                 w_m_l(nzb:nzt+1,nys-1:nyn+1,0:1) )
     ENDIF
     IF ( outflow_r )  THEN
-       ALLOCATE( u_m_r(nzb:nzt+1,nys-1:nyn+1,nx-1:nx+1), &
-                 v_m_r(nzb:nzt+1,nys-1:nyn+1,nx-1:nx+1), &
-                 w_m_r(nzb:nzt+1,nys-1:nyn+1,nx-1:nx+1) )
+       ALLOCATE( u_m_r(nzb:nzt+1,nys-1:nyn+1,nx-1:nx), &
+                 v_m_r(nzb:nzt+1,nys-1:nyn+1,nx-1:nx), &
+                 w_m_r(nzb:nzt+1,nys-1:nyn+1,nx-1:nx) )
+    ENDIF
+    IF ( outflow_l  .OR.  outflow_r )  THEN
+       ALLOCATE( c_u(nzb:nzt+1,nys-1:nyn+1), c_v(nzb:nzt+1,nys-1:nyn+1), &
+                 c_w(nzb:nzt+1,nys-1:nyn+1) )
     ENDIF
     IF ( outflow_s )  THEN
-       ALLOCATE( u_m_s(nzb:nzt+1,-1:1,nxl-1:nxr+1), &
-                 v_m_s(nzb:nzt+1,-1:1,nxl-1:nxr+1), &
-                 w_m_s(nzb:nzt+1,-1:1,nxl-1:nxr+1) )
+       ALLOCATE( u_m_s(nzb:nzt+1,0:1,nxl-1:nxr+1), &
+                 v_m_s(nzb:nzt+1,1:2,nxl-1:nxr+1), &
+                 w_m_s(nzb:nzt+1,0:1,nxl-1:nxr+1) )
     ENDIF
     IF ( outflow_n )  THEN
-       ALLOCATE( u_m_n(nzb:nzt+1,ny-1:ny+1,nxl-1:nxr+1), &
-                 v_m_n(nzb:nzt+1,ny-1:ny+1,nxl-1:nxr+1), &
-                 w_m_n(nzb:nzt+1,ny-1:ny+1,nxl-1:nxr+1) )
+       ALLOCATE( u_m_n(nzb:nzt+1,ny-1:ny,nxl-1:nxr+1), &
+                 v_m_n(nzb:nzt+1,ny-1:ny,nxl-1:nxr+1), &
+                 w_m_n(nzb:nzt+1,ny-1:ny,nxl-1:nxr+1) )
+    ENDIF
+    IF ( outflow_s  .OR.  outflow_n )  THEN
+       ALLOCATE( c_u(nzb:nzt+1,nxl-1:nxr+1), c_v(nzb:nzt+1,nxl-1:nxr+1), &
+                 c_w(nzb:nzt+1,nxl-1:nxr+1) )
     ENDIF
 
@@ -808 +817 @@
 !--    Initialize old timelevels needed for radiation boundary conditions
        IF ( outflow_l )  THEN
-          u_m_l(:,:,:) = u(:,:,-1:1)
-          v_m_l(:,:,:) = v(:,:,-1:1)
-          w_m_l(:,:,:) = w(:,:,-1:1)
+          u_m_l(:,:,:) = u(:,:,1:2)
+          v_m_l(:,:,:) = v(:,:,0:1)
+          w_m_l(:,:,:) = w(:,:,0:1)
        ENDIF
        IF ( outflow_r )  THEN
-          u_m_r(:,:,:) = u(:,:,nx-1:nx+1)
-          v_m_r(:,:,:) = v(:,:,nx-1:nx+1)
-          w_m_r(:,:,:) = w(:,:,nx-1:nx+1)
+          u_m_r(:,:,:) = u(:,:,nx-1:nx)
+          v_m_r(:,:,:) = v(:,:,nx-1:nx)
+          w_m_r(:,:,:) = w(:,:,nx-1:nx)
        ENDIF
        IF ( outflow_s )  THEN
-          u_m_s(:,:,:) = u(:,-1:1,:)
-          v_m_s(:,:,:) = v(:,-1:1,:)
-          w_m_s(:,:,:) = w(:,-1:1,:)
+          u_m_s(:,:,:) = u(:,0:1,:)
+          v_m_s(:,:,:) = v(:,1:2,:)
+          w_m_s(:,:,:) = w(:,0:1,:)
        ENDIF
        IF ( outflow_n )  THEN
-          u_m_n(:,:,:) = u(:,ny-1:ny+1,:)
-          v_m_n(:,:,:) = v(:,ny-1:ny+1,:)
-          w_m_n(:,:,:) = w(:,ny-1:ny+1,:)
+          u_m_n(:,:,:) = u(:,ny-1:ny,:)
+          v_m_n(:,:,:) = v(:,ny-1:ny,:)
+          w_m_n(:,:,:) = w(:,ny-1:ny,:)
        ENDIF