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

Timestamp:

Mar 19, 2008 9:41:30 AM (16 years ago)

Author:

steinfeld

Message:

Update for the plant canopy model

File:

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

palm/trunk/SOURCE/prognostic_equations.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
-! 
+! add call of plant_canopy_model in the prognostic equation for
+! the potential temperature
+! add call of plant_canopy_model in the prognostic equation for
+! the concentration of a passive scalar
 !
 ! Former revisions:
@@ -422 +425 @@
              CALL impact_of_latent_heat( i, j )
           ENDIF
+
+!
+!--       Consideration of heat sources within the plant canopy
+          IF ( plant_canopy .AND. ( cthf /= 0.0 ) ) THEN
+             CALL plant_canopy_model( i, j, 4 )
+          ENDIF
+
           CALL user_actions( i, j, 'pt-tendency' )
 
@@ -612 +622 @@
                 CALL calc_precipitation( i, j )
              ENDIF
+
+!
+!--          Sink or source of scalar concentration due to canopy elements
+             IF ( plant_canopy ) CALL plant_canopy_model( i, j, 5 )
+
              CALL user_actions( i, j, 'q-tendency' )
 
@@ -754 +769 @@
 !
 !--          Additional sink term for flows through plant canopies
-             IF ( plant_canopy )  CALL plant_canopy_model( i, j, 4 )         
+             IF ( plant_canopy )  CALL plant_canopy_model( i, j, 6 )         
 
              CALL user_actions( i, j, 'e-tendency' )
@@ -1034 +1049 @@
              CALL impact_of_latent_heat( i, j )
           ENDIF
+
+!
+!--       Consideration of heat sources within the plant canopy
+          IF ( plant_canopy .AND. ( cthf /= 0.0 ) ) THEN
+             CALL plant_canopy_model( i, j, 4 )
+          ENDIF
+
           CALL user_actions( i, j, 'pt-tendency' )
 
@@ -1142 +1164 @@
                 CALL calc_precipitation( i, j )
              ENDIF
+
+!
+!--          Sink or source of scalar concentration due to canopy elements
+             IF ( plant_canopy ) CALL plant_canopy_model( i, j, 5 )
+
+
              CALL user_actions( i, j, 'q-tendency' )
 
@@ -1219 +1247 @@
 !
 !--          Additional sink term for flows through plant canopies
-             IF ( plant_canopy )  CALL plant_canopy_model( i, j, 4 ) 
+             IF ( plant_canopy )  CALL plant_canopy_model( i, j, 6 ) 
 
              CALL user_actions( i, j, 'e-tendency' )
@@ -1582 +1610 @@
        CALL impact_of_latent_heat
     ENDIF
+
+!
+!-- Consideration of heat sources within the plant canopy
+    IF ( plant_canopy .AND. ( cthf /= 0.0 ) ) THEN
+       CALL plant_canopy_model( 4 )
+    ENDIF
+
+
     CALL user_actions( 'pt-tendency' )
 
@@ -1782 +1818 @@
           CALL calc_precipitation
        ENDIF
+
+!
+!--    Sink or source of scalar concentration due to canopy elements
+       IF ( plant_canopy ) CALL plant_canopy_model( 5 )
+
        CALL user_actions( 'q-tendency' )
 
@@ -1919 +1960 @@
 !
 !--    Additional sink term for flows through plant canopies
-       IF ( plant_canopy )  CALL plant_canopy_model( 4 )
+       IF ( plant_canopy )  CALL plant_canopy_model( 6 )
        CALL user_actions( 'e-tendency' )