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

Timestamp:

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

Author:

steinfeld

Message:

Update for the plant canopy model

File:

• palm/trunk/SOURCE/modules.f90 (modified) (5 diffs)

palm/trunk/SOURCE/modules.f90

@@ -5 +5 @@
 ! Actual revisions:
 ! -----------------
+! +canopy_heat_flux, cthf, lai,
+! +leaf_surface_concentration, scalar_exchange_coefficient, sec, sls
 ! +hor_index_bounds, hor_index_bounds_previous_run, id_inflow,
 ! inflow_damping_*, mean_inflow_profiles, numprocs_previous_run, nx_on_file,
@@ -135 +137 @@
 
     REAL, DIMENSION(:,:,:), ALLOCATABLE ::                                     &
-          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
+          canopy_heat_flux, cdc, d, diss, lad_s, lad_u, lad_v, lad_w, lai,     &
+          l_wall, sec, sls, 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 ::                             &
@@ -375 +377 @@
              building_height = 50.0, building_length_x = 50.0, &
              building_length_y = 50.0, building_wall_left = 9999999.9, &
-             building_wall_south = 9999999.9, cfl_factor = -1.0, &
+             building_wall_south = 9999999.9, cthf = 0.0, cfl_factor = -1.0, &
              cos_alpha_surface, disturbance_amplitude = 0.25, &
              disturbance_energy_limit = 0.01, &
@@ -396 +398 @@
              f = 0.0, fs = 0.0, g = 9.81, inflow_damping_height = 9999999.9, &
              inflow_damping_width = 9999999.9, kappa = 0.4, km_constant = -1.0,&
-             km_damp_max = -1.0, lad_surface = 0.0, long_filter_factor = 0.0, &
+             km_damp_max = -1.0, lad_surface = 0.0,  &
+             leaf_surface_concentration = 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, &
@@ -412 +415 @@
              rif_max = 1.0, &
              rif_min = -5.0, roughness_length = 0.1, sa_surface = 35.0, &
+             scalar_exchange_coefficient = 0.0, &
              simulated_time = 0.0, simulated_time_at_begin, sin_alpha_surface, &
              skip_time_data_output = 0.0, skip_time_data_output_av = 9999999.9,&