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

Timestamp:

Nov 28, 2007 10:03:58 AM (16 years ago)

Author:

letzel

Message:

Plant canopy model of Watanabe (2004,BLM 112,307-341) added.

File:

• palm/trunk/SOURCE/prognostic_equations.f90 (modified) (14 diffs)

palm/trunk/SOURCE/prognostic_equations.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
+! add call of subroutines that evaluate the canopy drag terms
 ! add wall_*flux to parameter list of calls of diffusion_s
 !
@@ -80 +81 @@
     USE diffusion_w_mod
     USE impact_of_latent_heat_mod
+    USE plant_canopy_model_mod
     USE production_e_mod
     USE user_actions_mod
@@ -160 +162 @@
           CALL coriolis( i, j, 1 )
           IF ( sloping_surface )  CALL buoyancy( i, j, pt, pt_reference, 1, 4 )
+
+!
+!--       Drag by plant canopy
+          IF ( plant_canopy )  CALL plant_canopy_model( i, j, 1 )
           CALL user_actions( i, j, 'u-tendency' )
 
@@ -227 +233 @@
           ENDIF
           CALL coriolis( i, j, 2 )
+
+!
+!--       Drag by plant canopy
+          IF ( plant_canopy )  CALL plant_canopy_model( i, j, 2 )     
+
           CALL user_actions( i, j, 'v-tendency' )
 
@@ -303 +314 @@
              ENDIF
           ENDIF
+
+!
+!--       Drag by plant canopy
+          IF ( plant_canopy )  CALL plant_canopy_model( i, j, 3 )
+
           CALL user_actions( i, j, 'w-tendency' )
 
@@ -732 +748 @@
              ENDIF
              CALL production_e( i, j )
+
+!
+!--          Additional sink term for flows through plant canopies
+             IF ( plant_canopy )  CALL plant_canopy_model( i, j, 4 )         
+
              CALL user_actions( i, j, 'e-tendency' )
 
@@ -833 +854 @@
              IF ( sloping_surface )  CALL buoyancy( i, j, pt, pt_reference, 1, &
                                                     4 )
+
+!
+!--          Drag by plant canopy
+             IF ( plant_canopy )  CALL plant_canopy_model( i, j, 1 )
+
              CALL user_actions( i, j, 'u-tendency' )
 
@@ -882 +908 @@
              ENDIF
              CALL coriolis( i, j, 2 )
+
+!
+!--          Drag by plant canopy
+             IF ( plant_canopy )  CALL plant_canopy_model( i, j, 2 )        
+
              CALL user_actions( i, j, 'v-tendency' )
 
@@ -938 +969 @@
              ENDIF
           ENDIF
+
+!
+!--       Drag by plant canopy
+          IF ( plant_canopy )  CALL plant_canopy_model( i, j, 3 )
+
           CALL user_actions( i, j, 'w-tendency' )
 
@@ -1177 +1213 @@
              ENDIF
              CALL production_e( i, j )
+
+!
+!--          Additional sink term for flows through plant canopies
+             IF ( plant_canopy )  CALL plant_canopy_model( i, j, 4 ) 
+
              CALL user_actions( i, j, 'e-tendency' )
 
@@ -1269 +1310 @@
     CALL coriolis( 1 )
     IF ( sloping_surface )  CALL buoyancy( pt, pt_reference, 1, 4 )
+
+!
+!-- Drag by plant canopy
+    IF ( plant_canopy )  CALL plant_canopy_model( 1 )
+
     CALL user_actions( 'u-tendency' )
 
@@ -1340 +1386 @@
     ENDIF
     CALL coriolis( 2 )
+
+!
+!-- Drag by plant canopy
+    IF ( plant_canopy )  CALL plant_canopy_model( 2 )
     CALL user_actions( 'v-tendency' )
 
@@ -1420 +1470 @@
        ENDIF
     ENDIF
+
+!
+!-- Drag by plant canopy
+    IF ( plant_canopy )  CALL plant_canopy_model( 3 )
+
     CALL user_actions( 'w-tendency' )
 
@@ -1858 +1913 @@
        ENDIF
        CALL production_e
+
+!
+!--    Additional sink term for flows through plant canopies
+       IF ( plant_canopy )  CALL plant_canopy_model( 4 )
        CALL user_actions( 'e-tendency' )