Changeset 2232 for palm/trunk/SOURCE/header.f90

Timestamp:

May 30, 2017 5:47:52 PM (7 years ago)

Author:

suehring

Message:

Adjustments according new topography and surface-modelling concept implemented

File:

• palm/trunk/SOURCE/header.f90 (modified) (8 diffs)

palm/trunk/SOURCE/header.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! Adjustments to new topography and surface concept
+! Generic tunnel setup added
 ! 
 ! Former revisions:
@@ -293 +294 @@
 
     USE arrays_3d,                                                             &
-        ONLY:  pt_init, qsws, q_init, s_init, sa_init, shf, ug, vg, w_subs, zu,&
-               zw
+        ONLY:  pt_init, q_init, s_init, sa_init, ug, vg, w_subs, zu, zw
         
     USE control_parameters
@@ -320 +320 @@
  
     USE land_surface_model_mod,                                                &
-        ONLY: land_surface, lsm_header
+        ONLY: lsm_header
 
     USE microphysics_mod,                                                      &
@@ -362 +362 @@
     USE spectra_mod,                                                           &
         ONLY:  calculate_spectra, spectra_header
+
+    USE surface_mod,                                                           &
+        ONLY:  surf_def_h
 
     IMPLICIT NONE
@@ -933 +936 @@
           ENDIF
 
+       CASE ( 'tunnel' )
+          IF ( tunnel_width_x /= 9999999.9_wp )  THEN
+!
+!--          Tunnel axis in y direction
+             IF ( tunnel_length == 9999999.9_wp  .OR.                          &
+                  tunnel_length >= ( nx + 1 ) * dx )  THEN
+                WRITE ( io, 273 )  'y', tunnel_height, tunnel_wall_depth,      &
+                                        tunnel_width_x
+             ELSE
+                WRITE ( io, 274 )  'y', tunnel_height, tunnel_wall_depth,      &
+                                        tunnel_width_x, tunnel_length
+             ENDIF
+
+          ELSEIF ( tunnel_width_y /= 9999999.9_wp )  THEN
+!
+!--          Tunnel axis in x direction
+             IF ( tunnel_length == 9999999.9_wp  .OR.                          &
+                  tunnel_length >= ( ny + 1 ) * dy )  THEN
+                WRITE ( io, 273 )  'x', tunnel_height, tunnel_wall_depth,      &
+                                        tunnel_width_y
+             ELSE
+                WRITE ( io, 274 )  'x', tunnel_height, tunnel_wall_depth,      &
+                                        tunnel_width_y, tunnel_length
+             ENDIF
+          ENDIF
+
     END SELECT
 
@@ -1065 +1094 @@
        IF ( constant_heatflux )  THEN
           IF ( large_scale_forcing .AND. lsf_surf )  THEN
-             WRITE ( io, 306 )  shf(0,0)
+             IF ( surf_def_h(0)%ns >= 1 )  WRITE ( io, 306 )  surf_def_h(0)%shf(1)
           ELSE
              WRITE ( io, 306 )  surface_heatflux
@@ -1073 +1102 @@
        IF ( humidity  .AND.  constant_waterflux )  THEN
           IF ( large_scale_forcing .AND. lsf_surf )  THEN
-             WRITE ( io, 311 ) qsws(0,0)
+             WRITE ( io, 311 ) surf_def_h(0)%qsws(1)
           ELSE
              WRITE ( io, 311 ) surface_waterflux
@@ -2024 +2053 @@
               ' Canyon height: ', F6.2, 'm, ch = ', I4, '.'      / &
               ' Canyon position (',A,'-walls): cxl = ', I4,', cxr = ', I4, '.')
+273 FORMAT (  ' Tunnel of infinite length in ',A, &
+              ' direction' / &
+              ' Tunnel height: ', F6.2, / &
+              ' Tunnel-wall depth: ', F6.2      / &
+              ' Tunnel width: ', F6.2 )
+274 FORMAT (  ' Tunnel in ', A, ' direction.' / &
+              ' Tunnel height: ', F6.2, / &    
+              ' Tunnel-wall depth: ', F6.2      / &
+              ' Tunnel width: ', F6.2, / &
+              ' Tunnel length: ', F6.2 )
 278 FORMAT (' Topography grid definition convention:'/ &
             ' cell edge (staggered grid points'/  &