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

Timestamp:

Apr 16, 2019 3:17:43 PM (5 years ago)

Author:

suehring

Message:

Bugfixes in initialization and STG

File:

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

palm/trunk/SOURCE/init_3d_model.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Fix problem with LOD = 2 initialization
+! 
+! 3885 2019-04-11 11:29:34Z kanani
 ! Changes related to global restructuring of location messages and introduction 
 ! of additional debug messages
@@ -1204 +1207 @@
           ENDDO
 !
-!--       Exchange ghost points in case of level-of-detail = 2
-          IF( init_3d%lod_u == 2 )   CALL exchange_horiz( u, nbgp )
-          IF( init_3d%lod_v == 2 )   CALL exchange_horiz( v, nbgp )
-          IF( init_3d%lod_w == 2 )   CALL exchange_horiz( w, nbgp )
+!--       Exchange ghost points and set boundary conditions in case of 
+!--       level-of-detail = 2
+          IF( init_3d%lod_u == 2 )  CALL exchange_horiz( u, nbgp )
+          IF( init_3d%lod_v == 2 )  CALL exchange_horiz( v, nbgp )
+          IF( init_3d%lod_w == 2 )  CALL exchange_horiz( w, nbgp )
           IF( .NOT. neutral  .AND.  init_3d%lod_pt == 2 )                      &
              CALL exchange_horiz( pt, nbgp )
           IF( humidity  .AND.  init_3d%lod_q == 2 )                            &
              CALL exchange_horiz( q, nbgp )
+          
+          IF ( bc_dirichlet_l )  THEN
+             DO  j = nysg, nyng
+                DO  k = nzb, nzt+1
+                   IF( init_3d%lod_u == 2 )  u(k,j,nxlg:nxl)   = u(k,j,nxlu)
+                   IF( init_3d%lod_v == 2 )  v(k,j,nxlg:nxl-1) = v(k,j,nxl)
+                   IF( init_3d%lod_w == 2 )  w(k,j,nxlg:nxl-1) = w(k,j,nxl)
+                   IF( .NOT. neutral  .AND.  init_3d%lod_pt == 2 )                      &
+                      pt(k,j,nxlg:nxl-1) = pt(k,j,nxl)
+                   IF( humidity  .AND.  init_3d%lod_q == 2 )                            &
+                      q(k,j,nxlg:nxl-1)  = q(k,j,nxl) 
+                ENDDO
+             ENDDO
+          ENDIF
+          IF ( bc_dirichlet_r )  THEN
+             DO  j = nysg, nyng
+                DO  k = nzb, nzt+1
+                   IF( init_3d%lod_u == 2 )  u(k,j,nxr+1:nxrg) = u(k,j,nxr)
+                   IF( init_3d%lod_v == 2 )  v(k,j,nxr+1:nxrg) = v(k,j,nxr)
+                   IF( init_3d%lod_w == 2 )  w(k,j,nxr+1:nxrg) = w(k,j,nxr)
+                   IF( .NOT. neutral  .AND.  init_3d%lod_pt == 2 )                      &
+                      pt(k,j,nxr+1:nxrg) = pt(k,j,nxr)
+                   IF( humidity  .AND.  init_3d%lod_q == 2 )                            &
+                      q(k,j,nxr+1:nxrg)  = q(k,j,nxr) 
+                ENDDO
+             ENDDO
+          ENDIF
+          IF ( bc_dirichlet_s )  THEN
+             DO  i = nxlg, nxrg
+                DO  k = nzb, nzt+1
+                   IF( init_3d%lod_u == 2 )  u(k,nysg:nys-1,i) = u(k,nys,i)
+                   IF( init_3d%lod_v == 2 )  v(k,nysg:nys,i)   = v(k,nysv,i)
+                   IF( init_3d%lod_w == 2 )  w(k,nysg:nys-1,i) = w(k,nys,i)
+                   IF( .NOT. neutral  .AND.  init_3d%lod_pt == 2 )                      &
+                      pt(k,nysg:nys-1,i) = pt(k,nys,i)
+                   IF( humidity  .AND.  init_3d%lod_q == 2 )                            &
+                      q(k,nysg:nys-1,i)  = q(k,nys,i)
+                ENDDO
+             ENDDO
+          ENDIF
+          IF ( bc_dirichlet_n )  THEN
+             DO  i = nxlg, nxrg
+                DO  k = nzb, nzt+1
+                   IF( init_3d%lod_u == 2 )  u(k,nyn+1:nyng,i) = u(k,nyn,i)
+                   IF( init_3d%lod_v == 2 )  v(k,nyn+1:nyng,i) = v(k,nyn,i)
+                   IF( init_3d%lod_w == 2 )  w(k,nyn+1:nyng,i) = w(k,nyn,i)
+                   IF( .NOT. neutral  .AND.  init_3d%lod_pt == 2 )                      &
+                      pt(k,nyn+1:nyng,i) = pt(k,nyn,i)
+                   IF( humidity  .AND.  init_3d%lod_q == 2 )                            &
+                      q(k,nyn+1:nyng,i)  = q(k,nyn,i)
+                ENDDO
+             ENDDO
+          ENDIF
 !
 !--       Set geostrophic wind components.