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

palm/trunk/SOURCE/modules.f90

@@ -5 +5 @@
 ! Actual revisions:
 ! -----------------
+! +drag_coefficient, pch_index, lad_surface, lad_vertical_gradient, 
+! lad_vertical_gradient_level, plant_canopy, lad, lad_s, lad_u, lad_v, 
+! lad_w, cdc, lad_vertical_gradient_level_ind, canopy_mode
+!
 ! +dt_sort_particles, ngp_2dh_s_inner, time_sort_particles, flags, 
 ! wall_flags_1..10, wall_humidityflux(0:4), wall_qflux(0:4), 
@@ -107 +111 @@
 
     REAL, DIMENSION(:), ALLOCATABLE ::                                         &
-          ddzu, dd2zu, dzu, ddzw, dzw, hyp, km_damp_x, km_damp_y, l_grid,      &
+          ddzu, dd2zu, dzu, ddzw, dzw, hyp, km_damp_x, km_damp_y, lad, l_grid, &
           pt_init, q_init, rdf, sa_init, ug, u_init, u_nzb_p1_for_vfc, vg,     &
           v_init, v_nzb_p1_for_vfc, zu, zw
@@ -126 +130 @@
 
     REAL, DIMENSION(:,:,:), ALLOCATABLE ::                                     &
-          d, diss, l_wall, tend, u_m_l, u_m_n, u_m_r, u_m_s, v_m_l, v_m_n,     &
-          v_m_r, v_m_s, w_m_l, w_m_n, w_m_r, w_m_s
+          cdc, d, diss, lad_s, lad_u, lad_v, lad_w, l_wall, tend,              &
+          u_m_l, u_m_n, u_m_r, u_m_s, v_m_l, v_m_n, v_m_r, v_m_s, w_m_l,       &
+          w_m_n, w_m_r, w_m_s
 
     REAL, DIMENSION(:,:,:), ALLOCATABLE, TARGET ::                             &
@@ -259 +264 @@
                             bc_sa_t = 'neumann', &
                             bc_uv_b = 'dirichlet', bc_uv_t = 'dirichlet', &
+                            canopy_mode = 'block', &
                             coupling_mode = 'uncoupled', &
                             coupling_mode_remote = 'uncoupled', &
@@ -307 +313 @@
                 do2d_xz_time_count(0:1), do2d_yz_time_count(0:1), &
                 do3d_no(0:1) = 0, do3d_time_count(0:1), &
+                lad_vertical_gradient_level_ind(10) = -9999, &
+                pch_index = 0, &
                 pt_vertical_gradient_level_ind(10) = -9999, &
                 q_vertical_gradient_level_ind(10) = -9999, &
@@ -340 +348 @@
                 outflow_l = .FALSE., outflow_n = .FALSE., outflow_r = .FALSE., &
                 outflow_s = .FALSE., passive_scalar = .FALSE., &
+                plant_canopy = .FALSE., &
                 prandtl_layer = .TRUE., precipitation = .FALSE., &
                 profil_output = .FALSE., radiation = .FALSE., &
@@ -365 +374 @@
              disturbance_level_b = -9999999.9, &
              disturbance_level_t = -9999999.9, &
+             drag_coefficient = 0.0, &
              dt = -1.0, dt_averaging_input = 0.0, &
              dt_averaging_input_pr = 9999999.9, dt_coupling = 9999999.9, &
@@ -379 +389 @@
              end_time = 0.0, &
              f = 0.0, fs = 0.0, g = 9.81, kappa = 0.4, km_constant = -1.0, &
-             km_damp_max = -1.0, long_filter_factor = 0.0, &
+             km_damp_max = -1.0, lad_surface = 0.0, long_filter_factor = 0.0, &
              maximum_cpu_time_allowed = 0.0, molecular_viscosity = 1.461E-5, &
              old_dt = 1.0E-10, omega = 7.29212E-5, omega_sor = 1.8, &
@@ -419 +429 @@
 
     REAL ::  do2d_xy_last_time(0:1) = -1.0, do2d_xz_last_time(0:1) = -1.0, &
-             do2d_yz_last_time(0:1) = -1.0, pt_vertical_gradient(10) = 0.0, &
+             do2d_yz_last_time(0:1) = -1.0, &
+             lad_vertical_gradient(10) = 0.0, &
+             lad_vertical_gradient_level(10) = -9999999.9, &
+             pt_vertical_gradient(10) = 0.0, &
              pt_vertical_gradient_level(10) = -9999999.9, &
              q_vertical_gradient(10) = 0.0, &