Changeset 240 for palm/trunk/SOURCE/header.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/header.f90 (modified) (6 diffs)

palm/trunk/SOURCE/header.f90

@@ -5 +5 @@
 ! -----------------
 ! Output of cluster_size
+! +dp_external, dp_level_b, dp_smooth, dpdxy
+! +canyon_height, canyon_width_x, canyon_width_y, canyon_wall_left,
+! canyon_wall_south
 !
 ! Former revisions:
@@ -99 +102 @@
     CHARACTER (LEN=85) ::  roben, runten
 
-    INTEGER ::  av, bh, blx, bly, bxl, bxr, byn, bys, i, ihost, io, j, l, ll
+    INTEGER ::  av, bh, blx, bly, bxl, bxr, byn, bys, ch, cwx, cwy, cxl, cxr, &
+         cyn, cys, i, ihost, io, j, l, ll
     REAL    ::  cpuseconds_per_simulated_second
 
@@ -285 +289 @@
     ENDIF
     IF ( passive_scalar )  WRITE ( io, 134 )
-    IF ( conserve_volume_flow )  WRITE ( io, 150 )
+    IF ( conserve_volume_flow )  THEN
+       WRITE ( io, 150 )
+    ELSEIF ( dp_external )  THEN
+       IF ( dp_smooth )  THEN
+          WRITE ( io, 151 ) dpdxy, dp_level_b, ', vertically smoothed.'
+       ELSE
+          WRITE ( io, 151 ) dpdxy, dp_level_b, '.'
+       ENDIF
+    ENDIF
     WRITE ( io, 99 )
 
@@ -381 +393 @@
           WRITE ( io, 271 )  building_length_x, building_length_y, &
                              building_height, bxl, bxr, bys, byn
+
+       CASE ( 'single_street_canyon' )
+          ch  = NINT( canyon_height / dz )
+          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
+             WRITE ( io, 272 )  'y', canyon_height, ch, 'u', cxl, cxr
+
+          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
+             WRITE ( io, 272 )  'x', canyon_height, ch, 'v', cys, cyn
+          ENDIF
 
     END SELECT
@@ -1370 +1407 @@
 150 FORMAT (' --> Volume flow at the right and north boundary will be ', &
                   'conserved')
+151 FORMAT (' --> External pressure gradient directly prescribed by the user:'/, &
+              2(1X,E12.5),'Pa/m', &
+             ' in x/y direction starting from dp_level_b =', F6.3, 'm', &
+             A /)
 200 FORMAT (//' Run time and time step information:'/ &
              ' ----------------------------------'/)
@@ -1407 +1448 @@
               ' Horizontal index bounds (l/r/s/n): ',I4,' / ',I4,' / ',I4, &
                 ' / ',I4)
+272 FORMAT (  ' Single quasi-2D street canyon of infinite length in ',A, &
+              ' direction' / &
+              ' Canyon height: ', F6.2, 'm, ch = ', I4, '.'      / &
+              ' Canyon position (',A,'-walls): cxl = ', I4,', cxr = ', I4, '.')
 280 FORMAT (//' Vegetation canopy (drag) model:'/ &
               ' ------------------------------'// &