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

palm/trunk/SOURCE/user_init_grid.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! +Generic tunnel added
+! +Example of setting user-defined topography
 ! 
 ! Former revisions:
@@ -55 +56 @@
 !> Execution of user-defined grid initializing actions
 !------------------------------------------------------------------------------!
- SUBROUTINE user_init_grid( nzb_local )
+ SUBROUTINE user_init_grid( topo_3d )
  
 
@@ -68 +69 @@
     IMPLICIT NONE
 
-    INTEGER(iwp), DIMENSION(nysg:nyng,nxlg:nxrg) ::  nzb_local   !< 
+    INTEGER(iwp)                                           ::  k_topo      !< topography top index
+    INTEGER(iwp), DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ::  topo_3d     !< 3D topography field
+
+    REAL(wp) ::  h_topo !< user-defined topography height
 
 !
@@ -78 +82 @@
     SELECT CASE ( TRIM( topography ) )
 
-       CASE ( 'flat', 'single_building', 'single_street_canyon' )
+       CASE ( 'flat', 'single_building', 'single_street_canyon', 'tunnel' )
 !
 !--       Not allowed here since these are the standard cases used in init_grid.
@@ -88 +92 @@
           message_string = 'topography "' // topography // '" not available yet'
           CALL message( 'user_init_grid', 'UI0005', 1, 2, 0, 6, 0 )
+!
+!--       The user is allowed to set surface-mounted as well as non-surface
+!--       mounted topography (e.g. overhanging structures). For both, use 
+!--       3D array topo_3d and set bit 0. The convention is: bit is zero inside 
+!--       topography, bit is 1 for atmospheric grid point. 
+!--       The following example shows how to prescribe sine-like topography 
+!--       along x-direction with amplitude of 10 * dz and wavelength 10 * dy.
+!           DO  i = nxlg, nxrg
+!              h_topo = 10.0_wp * dz * (SIN(3.14_wp*0.5_wp)*i*dx / ( 5.0_wp * dy ) )**2
+! 
+!              k_topo = MINLOC( ABS( zw - h_topo ), 1 ) - 1
+! 
+!              topo_3d(k_topo+1:nzt+1,:,i) =                                     &
+!                                          IBSET( topo_3d(k_topo+1:nzt+1,:,i), 0 ) 
+!           ENDDO 
+! 
+!           CALL exchange_horiz_int( topo_3d, nbgp )
 
        CASE DEFAULT
@@ -100 +121 @@
 
 
+
+
  END SUBROUTINE user_init_grid