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

Timestamp:

Feb 18, 2009 5:50:38 PM (15 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/prognostic_equations.f90 (modified) (7 diffs)

palm/trunk/SOURCE/prognostic_equations.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
-! 
+! external pressure gradient
 !
 ! Former revisions:
@@ -173 +173 @@
 !--       Drag by plant canopy
           IF ( plant_canopy )  CALL plant_canopy_model( i, j, 1 )
+
+!
+!--       External pressure gradient
+          IF ( dp_external )  THEN
+             DO  k = dp_level_ind_b+1, nzt
+                tend(k,j,i) = tend(k,j,i) - dpdxy(1) * dp_smooth_factor(k)
+             ENDDO
+          ENDIF
+
           CALL user_actions( i, j, 'u-tendency' )
 
@@ -244 +253 @@
 !--       Drag by plant canopy
           IF ( plant_canopy )  CALL plant_canopy_model( i, j, 2 )     
+
+!
+!--       External pressure gradient
+          IF ( dp_external )  THEN
+             DO  k = dp_level_ind_b+1, nzt
+                tend(k,j,i) = tend(k,j,i) - dpdxy(2) * dp_smooth_factor(k)
+             ENDDO
+          ENDIF
 
           CALL user_actions( i, j, 'v-tendency' )
@@ -878 +895 @@
              IF ( plant_canopy )  CALL plant_canopy_model( i, j, 1 )
 
+!
+!--          External pressure gradient
+             IF ( dp_external )  THEN
+                DO  k = dp_level_ind_b+1, nzt
+                   tend(k,j,i) = tend(k,j,i) - dpdxy(1) * dp_smooth_factor(k)
+                ENDDO
+             ENDIF
+
              CALL user_actions( i, j, 'u-tendency' )
 
@@ -931 +956 @@
 !--          Drag by plant canopy
              IF ( plant_canopy )  CALL plant_canopy_model( i, j, 2 )        
+
+!
+!--          External pressure gradient
+             IF ( dp_external )  THEN
+                DO  k = dp_level_ind_b+1, nzt
+                   tend(k,j,i) = tend(k,j,i) - dpdxy(2) * dp_smooth_factor(k)
+                ENDDO
+             ENDIF
 
              CALL user_actions( i, j, 'v-tendency' )
@@ -1347 +1380 @@
     IF ( plant_canopy )  CALL plant_canopy_model( 1 )
 
+!
+!-- External pressure gradient
+    IF ( dp_external )  THEN
+       DO  i = nxlu, nxr
+          DO  j = nys, nyn
+             DO  k = dp_level_ind_b+1, nzt
+                tend(k,j,i) = tend(k,j,i) - dpdxy(1) * dp_smooth_factor(k)
+             ENDDO
+          ENDDO
+       ENDDO
+    ENDIF
+
     CALL user_actions( 'u-tendency' )
 
@@ -1422 +1467 @@
 !-- Drag by plant canopy
     IF ( plant_canopy )  CALL plant_canopy_model( 2 )
+
+!
+!-- External pressure gradient
+    IF ( dp_external )  THEN
+       DO  i = nxl, nxr
+          DO  j = nysv, nyn
+             DO  k = dp_level_ind_b+1, nzt
+                tend(k,j,i) = tend(k,j,i) - dpdxy(2) * dp_smooth_factor(k)
+             ENDDO
+          ENDDO
+       ENDDO
+    ENDIF
+
     CALL user_actions( 'v-tendency' )