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

Timestamp:

Feb 18, 2009 5:50:38 PM (16 years ago)

Author:

letzel

Message:

• External pressure gradient (check_parameters, init_3d_model, header, modules, parin, prognostic_equations)
• New topography case 'single_street_canyon'

File:

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

palm/trunk/SOURCE/init_grid.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
-! 
+! new topography case 'single_street_canyon'
 !
 ! Former revisions:
@@ -54 +54 @@
     IMPLICIT NONE
 
-    INTEGER ::  bh, blx, bly, bxl, bxr, byn, bys, gls, i, inc, i_center, j, &
-                j_center, k, l, nxl_l, nxr_l, nyn_l, nys_l, nzb_si, nzt_l, vi
+    INTEGER ::  bh, blx, bly, bxl, bxr, byn, bys, ch, cwx, cwy, cxl, cxr, cyn, &
+                cys, gls, i, inc, i_center, j, j_center, k, l, nxl_l, nxr_l, &
+                nyn_l, nys_l, nzb_si, nzt_l, vi
 
     INTEGER, DIMENSION(:), ALLOCATABLE   ::  vertical_influence
@@ -379 +380 @@
           nzb_local(bys:byn,bxl:bxr) = bh
 
+       CASE ( 'single_street_canyon' )
+!
+!--       Single quasi-2D street canyon of infinite length in x or y direction.
+!--       The canyon is centered in the other direction by default.
+          IF ( canyon_width_x /= 9999999.9 )  THEN
+!
+!--          Street canyon in y direction
+             cwx = NINT( canyon_width_x / dx )
+             IF ( canyon_wall_left == 9999999.9 )  THEN
+                canyon_wall_left = ( nx + 1 - cwx ) / 2 * dx
+             ENDIF
+             cxl = NINT( canyon_wall_left / dx )
+             cxr = cxl + cwx
+
+          ELSEIF ( canyon_width_y /= 9999999.9 )  THEN
+!
+!--          Street canyon in x direction
+             cwy = NINT( canyon_width_y / dy )
+             IF ( canyon_wall_south == 9999999.9 )  THEN
+                canyon_wall_south = ( ny + 1 - cwy ) / 2 * dy
+             ENDIF
+             cys = NINT( canyon_wall_south / dy )
+             cyn = cys + cwy
+
+          ELSE
+             IF ( myid == 0 )  THEN
+                PRINT*, '+++ user_init_grid: no street canyon width given'
+             ENDIF
+             CALL local_stop
+          ENDIF
+
+          ch             = NINT( canyon_height / dz )
+          dp_level_ind_b = ch
+!
+!--       Street canyon size has to meet some requirements
+          IF ( canyon_width_x /= 9999999.9 )  THEN
+             IF ( ( cxl < 1 ) .OR. ( cxr > nx-1 ) .OR. ( cwx < 3 ) .OR.  &
+               ( ch < 3 ) )  THEN
+                IF ( myid == 0 )  THEN
+                   PRINT*, '+++ user_init_grid: inconsistent canyon parameters:'
+                   PRINT*, '                    cxl=', cxl, 'cxr=', cxr,  &
+                                              'cwx=', cwx,  &
+                                              'ch=', ch, 'nx=', nx, 'ny=', ny
+                ENDIF
+                CALL local_stop
+             ENDIF
+          ELSEIF ( canyon_width_y /= 9999999.9 )  THEN
+             IF ( ( cys < 1 ) .OR. ( cyn > ny-1 ) .OR. ( cwy < 3 ) .OR.  &
+               ( ch < 3 ) )  THEN
+                IF ( myid == 0 )  THEN
+                   PRINT*, '+++ user_init_grid: inconsistent canyon parameters:'
+                   PRINT*, '                    cys=', cys, 'cyn=', cyn,  &
+                                              'cwy=', cwy,  &
+                                              'ch=', ch, 'nx=', nx, 'ny=', ny
+                ENDIF
+                CALL local_stop
+             ENDIF
+          ENDIF
+          IF ( canyon_width_x /= 9999999.9 .AND. canyon_width_y /= 9999999.9 )  &
+               THEN
+             IF ( myid == 0 )  THEN
+                PRINT*, '+++ user_init_grid: inconsistent canyon parameters:'
+                PRINT*, '                    street canyon can only be oriented'
+                PRINT*, '                    either in x- or in y-direction'
+             ENDIF
+             CALL local_stop
+          ENDIF
+
+          nzb_local = ch
+          IF ( canyon_width_x /= 9999999.9 )  THEN
+             nzb_local(:,cxl+1:cxr-1) = 0
+          ELSEIF ( canyon_width_y /= 9999999.9 )  THEN
+             nzb_local(cys+1:cyn-1,:) = 0
+          ENDIF
+
        CASE ( 'read_from_file' )
 !