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

Timestamp:

Jan 18, 2021 11:15:37 AM (3 years ago)

Author:

raasch

Message:

maximum phase velocities are alwasy used for radiation boundary conditions, parameter use_cmax removed

File:

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

palm/trunk/SOURCE/init_3d_model.f90

@@ -24 +24 @@
 ! -----------------
 ! $Id$
+! removed allocation and initialization of arrays required for radiation boundary conditions
+!
+! 4828 2021-01-05 11:21:41Z Giersch
 ! Add possibility to initialize surface flux of passive scalar via static driver
 !
@@ -560 +563 @@
 
 !
-!-- Arrays to store velocity data from t-dt and the phase speeds which are needed for radiation
-!-- boundary conditions.
-    IF ( bc_radiation_l )  THEN
-       ALLOCATE( u_m_l(nzb:nzt+1,nysg:nyng,1:2),                                                   &
-                 v_m_l(nzb:nzt+1,nysg:nyng,0:1),                                                   &
-                 w_m_l(nzb:nzt+1,nysg:nyng,0:1) )
-    ENDIF
-    IF ( bc_radiation_r )  THEN
-       ALLOCATE( u_m_r(nzb:nzt+1,nysg:nyng,nx-1:nx),                                               &
-                 v_m_r(nzb:nzt+1,nysg:nyng,nx-1:nx),                                               &
-                 w_m_r(nzb:nzt+1,nysg:nyng,nx-1:nx) )
-    ENDIF
-    IF ( bc_radiation_l  .OR.  bc_radiation_r )  THEN
-       ALLOCATE( c_u(nzb:nzt+1,nysg:nyng), c_v(nzb:nzt+1,nysg:nyng), c_w(nzb:nzt+1,nysg:nyng) )
-    ENDIF
-    IF ( bc_radiation_s )  THEN
-       ALLOCATE( u_m_s(nzb:nzt+1,0:1,nxlg:nxrg),                                                   &
-                 v_m_s(nzb:nzt+1,1:2,nxlg:nxrg),                                                   &
-                 w_m_s(nzb:nzt+1,0:1,nxlg:nxrg) )
-    ENDIF
-    IF ( bc_radiation_n )  THEN
-       ALLOCATE( u_m_n(nzb:nzt+1,ny-1:ny,nxlg:nxrg),                                               &
-                 v_m_n(nzb:nzt+1,ny-1:ny,nxlg:nxrg),                                               &
-                 w_m_n(nzb:nzt+1,ny-1:ny,nxlg:nxrg) )
-    ENDIF
-    IF ( bc_radiation_s  .OR.  bc_radiation_n )  THEN
-       ALLOCATE( c_u(nzb:nzt+1,nxlg:nxrg), c_v(nzb:nzt+1,nxlg:nxrg), c_w(nzb:nzt+1,nxlg:nxrg) )
-    ENDIF
-    IF ( bc_radiation_l  .OR.  bc_radiation_r  .OR.  bc_radiation_s  .OR.  bc_radiation_n )  THEN
-       ALLOCATE( c_u_m_l(nzb:nzt+1), c_v_m_l(nzb:nzt+1), c_w_m_l(nzb:nzt+1) )
-       ALLOCATE( c_u_m(nzb:nzt+1), c_v_m(nzb:nzt+1), c_w_m(nzb:nzt+1) )
-    ENDIF
-
-!
 !-- Initial assignment of the pointers
     IF ( .NOT. neutral )  THEN
@@ -1344 +1313 @@
        message_string = 'unknown initializing problem'
        CALL message( 'init_3d_model', 'PA0193', 1, 2, 0, 6, 0 )
-    ENDIF
-
-
-    IF (  TRIM( initializing_actions ) /= 'read_restart_data' )  THEN
-!
-!--    Initialize old timelevels needed for radiation boundary conditions
-       IF ( bc_radiation_l )  THEN
-          u_m_l(:,:,:) = u(:,:,1:2)
-          v_m_l(:,:,:) = v(:,:,0:1)
-          w_m_l(:,:,:) = w(:,:,0:1)
-       ENDIF
-       IF ( bc_radiation_r )  THEN
-          u_m_r(:,:,:) = u(:,:,nx-1:nx)
-          v_m_r(:,:,:) = v(:,:,nx-1:nx)
-          w_m_r(:,:,:) = w(:,:,nx-1:nx)
-       ENDIF
-       IF ( bc_radiation_s )  THEN
-          u_m_s(:,:,:) = u(:,0:1,:)
-          v_m_s(:,:,:) = v(:,1:2,:)
-          w_m_s(:,:,:) = w(:,0:1,:)
-       ENDIF
-       IF ( bc_radiation_n )  THEN
-          u_m_n(:,:,:) = u(:,ny-1:ny,:)
-          v_m_n(:,:,:) = v(:,ny-1:ny,:)
-          w_m_n(:,:,:) = w(:,ny-1:ny,:)
-       ENDIF
-
     ENDIF