Changeset 2037 for palm/trunk/SOURCE/init_grid.f90

Timestamp:

Oct 26, 2016 11:15:40 AM (8 years ago)

Author:

knoop

Message:

Anelastic approximation implemented

File:

• palm/trunk/SOURCE/init_grid.f90 (modified) (4 diffs)

palm/trunk/SOURCE/init_grid.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! Anelastic approximation implemented
 ! 
 ! Former revisions:
@@ -222 +222 @@
          
     USE grid_variables,                                                        &
-        ONLY:  ddx, ddx2, ddx2_mg, ddy, ddy2, ddy2_mg, dx, dx2, dy, dy2, fwxm, &
+        ONLY:  ddx, ddx2, ddy, ddy2, dx, dx2, dy, dy2, fwxm,                   &
                fwxp, fwym, fwyp, fxm, fxp, fym, fyp, wall_e_x, wall_e_y,       &
                wall_u, wall_v, wall_w_x, wall_w_y, zu_s_inner, zw_w_inner
@@ -293 +293 @@
 
     REAL(wp) ::  dum           !< dummy variable to skip columns while reading topography file    
-    REAL(wp) ::  dx_l          !< grid spacing along x on different multigrid level
-    REAL(wp) ::  dy_l          !< grid spacing along y on different multigrid level
     REAL(wp) ::  dz_stretched  !< stretched vertical grid spacing
 
@@ -431 +429 @@
        ddzu_pres = ddzu
        ddzu_pres(1) = ddzu_pres(2)  ! change for lowest level
-    ENDIF   
-
-!
-!-- In case of multigrid method, compute grid lengths and grid factors for the
-!-- grid levels
-    IF ( psolver(1:9) == 'multigrid' )  THEN
-
-       ALLOCATE( ddx2_mg(maximum_grid_level), ddy2_mg(maximum_grid_level), &
-                 dzu_mg(nzb+1:nzt+1,maximum_grid_level),                   &
-                 dzw_mg(nzb+1:nzt+1,maximum_grid_level),                   &
-                 f1_mg(nzb+1:nzt,maximum_grid_level),                      &
-                 f2_mg(nzb+1:nzt,maximum_grid_level),                      &
-                 f3_mg(nzb+1:nzt,maximum_grid_level) )
-
-       dzu_mg(:,maximum_grid_level) = dzu
-!       
-!--    Next line to ensure an equally spaced grid. 
-       dzu_mg(1,maximum_grid_level) = dzu(2)
-
-       dzw_mg(:,maximum_grid_level) = dzw
-       nzt_l = nzt
-       DO  l = maximum_grid_level-1, 1, -1
-           dzu_mg(nzb+1,l) = 2.0_wp * dzu_mg(nzb+1,l+1)
-           dzw_mg(nzb+1,l) = 2.0_wp * dzw_mg(nzb+1,l+1)
-           nzt_l = nzt_l / 2
-           DO  k = 2, nzt_l+1
-              dzu_mg(k,l) = dzu_mg(2*k-2,l+1) + dzu_mg(2*k-1,l+1)
-              dzw_mg(k,l) = dzw_mg(2*k-2,l+1) + dzw_mg(2*k-1,l+1)
-           ENDDO
-       ENDDO
-
-       nzt_l = nzt
-       dx_l  = dx
-       dy_l  = dy
-       DO  l = maximum_grid_level, 1, -1
-          ddx2_mg(l) = 1.0_wp / dx_l**2
-          ddy2_mg(l) = 1.0_wp / dy_l**2
-          DO  k = nzb+1, nzt_l
-             f2_mg(k,l) = 1.0_wp / ( dzu_mg(k+1,l) * dzw_mg(k,l) )
-             f3_mg(k,l) = 1.0_wp / ( dzu_mg(k,l)   * dzw_mg(k,l) )
-             f1_mg(k,l) = 2.0_wp * ( ddx2_mg(l) + ddy2_mg(l) ) + &
-                          f2_mg(k,l) + f3_mg(k,l)
-          ENDDO
-          nzt_l = nzt_l / 2
-          dx_l  = dx_l * 2.0_wp
-          dy_l  = dy_l * 2.0_wp
-       ENDDO
-
     ENDIF