Changeset 138 for palm/trunk/SOURCE/check_parameters.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/check_parameters.f90 (modified) (4 diffs)

palm/trunk/SOURCE/check_parameters.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
+! Plant canopy added
 ! Allow new case bc_uv_t = 'dirichlet_0' for channel flow.
 ! Multigrid solver allows topography, checking of dt_sort_particles
@@ -592 +593 @@
     ENDIF
 
+    IF ( plant_canopy .AND. ( drag_coefficient == 0.0 ) ) THEN
+       IF ( myid == 0 )  PRINT*, '+++ check_parameters: plant_canopy = TRUE', &
+                                 'requires a non-zero drag_coefficient'
+       CALL local_stop
+    ENDIF 
+
 !
 !-- In case of no model continuation run, check initialising parameters and
@@ -605 +612 @@
        IF ( ocean )           sa_init = sa_surface
        IF ( passive_scalar )  q_init  = s_surface
+       IF ( plant_canopy )    lad = 0.0
 
 !
@@ -934 +942 @@
        ENDIF
 
-    ENDIF
-
+!
+!--    If required compute the profile of leaf area density used in the plant canopy model 
+       IF ( plant_canopy ) THEN
+       
+          i = 1
+          gradient = 0.0
+
+          IF ( .NOT. ocean ) THEN
+ 
+             lad_vertical_gradient_level_ind(1) = 0
+             DO k = 1, pch_index
+                IF ( lad_vertical_gradient_level(i) < zu(k) .AND.  &
+                     lad_vertical_gradient_level(i) >= 0.0 ) THEN
+                   gradient = lad_vertical_gradient(i)
+                   lad_vertical_gradient_level_ind(i) = k - 1
+                   i = i + 1
+                   IF ( i > 10 ) THEN
+                      IF ( myid == 0 ) THEN
+                         PRINT*, '+++ user_init_3d_model: upper bound 10 of array',  &
+                                 ' "lad_vertical_gradient_level" exceeded'
+                      ENDIF
+                      CALL local_stop
+                   ENDIF
+                ENDIF
+                IF ( gradient /= 0.0 ) THEN
+                   IF ( k /= 1 ) THEN
+                      lad(k) = lad(k-1) + dzu(k) * gradient
+                   ELSE
+                      lad(k) = lad_surface + 0.5 * dzu(k) *gradient
+                   ENDIF
+                ELSE
+                   lad(k) = lad(k-1)
+                ENDIF
+             ENDDO
+
+          ENDIF
+
+!
+!--       In case of no given leaf area density gradients, choose a vanishing gradient
+          IF ( lad_vertical_gradient_level(1) == -9999999.9 ) THEN
+             lad_vertical_gradient_level(1) = 0.0
+          ENDIF
+
+       ENDIF
+         
+    ENDIF
+             
 !
 !-- Compute Coriolis parameter