Changeset 3379

Timestamp:

Oct 19, 2018 12:38:28 PM (6 years ago)

Author:

knoop

Message:

Added topography flags to subroutine init_pt_anomaly

File:

• palm/trunk/SOURCE/init_pt_anomaly.f90 (modified) (10 diffs)

palm/trunk/SOURCE/init_pt_anomaly.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Added topography flags
+! 
+! 3035 2018-05-24 09:35:20Z schwenkel
 ! Add option to initialize warm air bubble close to surface
 ! 
@@ -70 +73 @@
 !------------------------------------------------------------------------------!
  SUBROUTINE init_pt_anomaly
- 
+
 
     USE arrays_3d,                                                             &
         ONLY:  pt, zu
 
-    USE control_parameters    
-        
+    USE control_parameters
+
     USE grid_variables,                                                        &
         ONLY:  dx, dy
 
     USE indices,                                                               &
-        ONLY:  nbgp, nx, nxl, nxr, ny, nyn, nys, nzb, nzt
-        
+        ONLY:  nbgp, nx, nxl, nxr, ny, nyn, nys, nzb, nzt, wall_flags_0
+
     USE kinds
 
@@ -93 +96 @@
     INTEGER(iwp) ::  k  !< grid index along z
     INTEGER(iwp) ::  kc !< center index along z
-    
+
     REAL(wp)     ::  amount                               !< amount of temperature perturbation
     REAL(wp)     ::  bubble_center_y                      !< center of bubble in y
@@ -99 +102 @@
     REAL(wp)     ::  bubble_sigma_y = 300.0               !< width of bubble in y
     REAL(wp)     ::  bubble_sigma_z = 150.0               !< width of bubble in z
+    REAL(wp)     ::  flag                                 !< flag to mask topography grid points
     REAL(wp)     ::  initial_temperature_difference = 0.4 !< temperature perturbation for bubble in K
     REAL(wp)     ::  radius                               !< radius of pt anomaly
@@ -105 +109 @@
     REAL(wp)     ::  y                                    !< y dimension of pt anomaly
     REAL(wp)     ::  z                                    !< z dimension of pt anomaly
-    
-    
+
+
 !
 !-- Defaults: radius rc, strength z,
@@ -114 +118 @@
     jc =  ic
     kc =  nzt / 2
-    
+
     IF ( INDEX( initializing_actions, 'initialize_ptanom' ) /= 0 )  THEN
 !
@@ -121 +125 @@
           DO  j = nys, nyn
              DO  k = nzb+1, nzt
+!
+!--             Predetermine flag to mask topography
+                flag = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_0(k,j,i), 0 ) )
+
                 x = ( i - ic ) * dx
                 y = ( j - jc ) * dy
@@ -131 +139 @@
                 ENDIF
 
-                pt(k,j,i) = pt(k,j,i) + amount
+                pt(k,j,i) = pt(k,j,i) + amount * flag
 
              ENDDO
           ENDDO
        ENDDO
-        
+
 !
 !-- Initialize warm air bubble close to surface and homogenous elegonated 
@@ -144 +152 @@
 !--    Calculate y-center of model domain
        bubble_center_y = ( ny + 1.0 ) * dy / 2.0
-    
+
 !
 !--    Compute perturbation for potential temperaure
        DO  i = nxl, nxr
           DO  j = nys, nyn
-             DO  k = nzb+1, nzt 
+             DO  k = nzb+1, nzt
+!
+!--             Predetermine flag to mask topography
+                flag = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_0(k,j,i), 0 ) )
+
                 pt(k,j,i) = pt(k,j,i) +                                        &
                                EXP( -0.5 * ( (j* dy  - bubble_center_y) /      &
@@ -155 +167 @@
                                EXP( -0.5 * ( (zu(k)  - bubble_center_z) /      &
                                                        bubble_sigma_z)**2) *   &
-                               initial_temperature_difference
+                               initial_temperature_difference * flag
              ENDDO
           ENDDO