Changeset 139 for palm/trunk/SOURCE


Ignore:
Timestamp:
Nov 29, 2007 9:37:41 AM (14 years ago)
Author:
raasch
Message:

New:
---

Plant canopy model of Watanabe (2004,BLM 112,307-341) added.
It can be switched on by the inipar parameter plant_canopy.
The inipar parameter canopy_mode can be used to prescribe a
plant canopy type. The default case is a homogeneous plant
canopy. Heterogeneous distributions of the leaf area
density and the canopy drag coefficient can be defined in the
new routine user_init_plant_canopy (user_interface).
The inipar parameters lad_surface, lad_vertical_gradient and
lad_vertical_gradient_level can be used in order to
prescribe the vertical profile of leaf area density. The
inipar parameter drag_coefficient determines the canopy
drag coefficient.
Finally, the inipar parameter pch_index determines the
index of the upper boundary of the plant canopy.

Allow new case bc_uv_t = 'dirichlet_0' for channel flow.

For unknown variables (CASE DEFAULT) call new subroutine user_data_output_dvrp

Pressure boundary conditions for vertical walls added to the multigrid solver.
They are applied using new wall flag arrays (wall_flags_..) which are defined
for each grid level. New argument gls added to routine user_init_grid
(user_interface).

Frequence of sorting particles can be controlled with new particles_par
parameter dt_sort_particles. Sorting is moved from the SGS timestep loop in
advec_particles after the end of this loop.

advec_particles, check_parameters, data_output_dvrp, header, init_3d_model, init_grid, init_particles, init_pegrid, modules, package_parin, parin, plant_canopy_model, read_var_list, read_3d_binary, user_interface, write_var_list, write_3d_binary

Changed:


Redefine initial nzb_local as the actual total size of topography (later the
extent of topography in nzb_local is reduced by 1dx at the E topography walls
and by 1dy at the N topography walls to form the basis for nzb_s_inner);
for consistency redefine 'single_building' case.

Vertical profiles now based on nzb_s_inner; they are divided by
ngp_2dh_s_inner (scalars, procucts of scalars) and ngp_2dh (staggered velocity
components and their products, procucts of scalars and velocity components),
respectively.

Allow two instead of one digit to specify isosurface and slicer variables.

Status of 3D-volume NetCDF data file only depends on switch netcdf_64bit_3d (check_open)

prognostic_equations include the respective wall_*flux in the parameter list of
calls of diffusion_s. Same as before, only the values of wall_heatflux(0:4)
can be assigned. At present, wall_humidityflux, wall_qflux, wall_salinityflux,
and wall_scalarflux are kept zero. diffusion_s uses the respective wall_*flux
instead of wall_heatflux. This update serves two purposes:

  • it avoids errors in calculations with humidity/scalar/salinity and prescribed

non-zero wall_heatflux,

  • it prepares PALM for a possible assignment of wall fluxes of

humidity/scalar/salinity in a future release.

buoyancy, check_open, data_output_dvrp, diffusion_s, diffusivities, flow_statistics, header, init_3d_model, init_dvrp, init_grid, modules, prognostic_equations

Errors:


Bugfix: summation of sums_l_l in diffusivities.

Several bugfixes in the ocean part: Initial density rho is calculated
(init_ocean). Error in initializing u_init and v_init removed
(check_parameters). Calculation of density flux now starts from
nzb+1 (production_e).

Bugfix: pleft/pright changed to pnorth/psouth in sendrecv of particle tail
numbers along y, small bugfixes in the SGS part (advec_particles)

Bugfix: model_string needed a default value (combine_plot_fields)

Bugfix: wavenumber calculation for even nx in routines maketri (poisfft)

Bugfix: assignment of fluxes at walls

Bugfix: absolute value of f must be used when calculating the Blackadar mixing length (init_1d_model)

advec_particles, check_parameters, combine_plot_fields, diffusion_s, diffusivities, init_ocean, init_1d_model, poisfft, production_e

Location:
palm/trunk/SOURCE
Files:
28 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/CURRENT_MODIFICATIONS

    r138 r139  
    22---
    33
    4 Plant canopy model of Watanabe (2004,BLM 112,307-341) added.
    5 It can be switched on by the inipar parameter plant_canopy.
    6 The inipar parameter canopy_mode can be used to prescribe a
    7 plant canopy type. The default case is a homogeneous plant
    8 canopy. Heterogeneous distributions of the leaf area
    9 density and the canopy drag coefficient can be defined in the
    10 new routine user_init_plant_canopy (user_interface).
    11 The inipar parameters lad_surface, lad_vertical_gradient and
    12 lad_vertical_gradient_level can be used in order to
    13 prescribe the vertical profile of leaf area density. The
    14 inipar parameter drag_coefficient determines the canopy
    15 drag coefficient.
    16 Finally, the inipar parameter pch_index determines the
    17 index of the upper boundary of the plant canopy. 
    18 
    19 Allow new case bc_uv_t = 'dirichlet_0' for channel flow.
    20 
    21 For unknown variables (CASE DEFAULT) call new subroutine user_data_output_dvrp
    22 
    23 Pressure boundary conditions for vertical walls added to the multigrid solver.
    24 They are applied using new wall flag arrays (wall_flags_..) which are defined
    25 for each grid level. New argument gls added to routine user_init_grid
    26 (user_interface).
    27 
    28 Frequence of sorting particles can be controlled with new particles_par
    29 parameter dt_sort_particles. Sorting is moved from the SGS timestep loop in
    30 advec_particles after the end of this loop.
    31 
    32 advec_particles, check_parameters, data_output_dvrp, header, init_3d_model, init_grid, init_particles, init_pegrid, modules, package_parin, parin, plant_canopy_model, read_var_list, read_3d_binary, user_interface, write_var_list, write_3d_binary
    334
    345
     
    367-------
    378
    38 Redefine initial nzb_local as the actual total size of topography (later the
    39 extent of topography in nzb_local is reduced by 1dx at the E topography walls
    40 and by 1dy at the N topography walls to form the basis for nzb_s_inner);
    41 for consistency redefine 'single_building' case.
    42 
    43 Vertical profiles now based on nzb_s_inner; they are divided by
    44 ngp_2dh_s_inner (scalars, procucts of scalars) and ngp_2dh (staggered velocity
    45 components and their products, procucts of scalars and velocity components),
    46 respectively.
    47 
    48 Allow two instead of one digit to specify isosurface and slicer variables.
    49 
    50 Status of 3D-volume NetCDF data file only depends on switch netcdf_64bit_3d (check_open)
    51 
    52 prognostic_equations include the respective wall_*flux in the parameter list of
    53 calls of diffusion_s. Same as before, only the values of wall_heatflux(0:4)
    54 can be assigned. At present, wall_humidityflux, wall_qflux, wall_salinityflux,
    55 and wall_scalarflux are kept zero. diffusion_s uses the respective wall_*flux
    56 instead of wall_heatflux. This update serves two purposes:
    57 - it avoids errors in calculations with humidity/scalar/salinity and prescribed
    58 non-zero wall_heatflux,
    59 - it prepares PALM for a possible assignment of wall fluxes of
    60 humidity/scalar/salinity in a future release.
    61 
    62 buoyancy, check_open, data_output_dvrp, diffusion_s, diffusivities, flow_statistics, header, init_3d_model, init_dvrp, init_grid, modules, prognostic_equations
    639
    6410
     
    6612------
    6713
    68 Bugfix: summation of sums_l_l in diffusivities.
    6914
    70 Several bugfixes in the ocean part: Initial density rho is calculated
    71 (init_ocean). Error in initializing u_init and v_init removed
    72 (check_parameters). Calculation of density flux now starts from
    73 nzb+1 (production_e).
    74 
    75 Bugfix: pleft/pright changed to pnorth/psouth in sendrecv of particle tail
    76 numbers along y, small bugfixes in the SGS part (advec_particles)
    77 
    78 Bugfix: model_string needed a default value (combine_plot_fields)
    79 
    80 Bugfix: wavenumber calculation for even nx in routines maketri (poisfft)
    81 
    82 Bugfix: assignment of fluxes at walls
    83 
    84 Bugfix: absolute value of f must be used when calculating the Blackadar mixing length (init_1d_model)
    85 
    86 advec_particles, check_parameters, combine_plot_fields, diffusion_s, diffusivities, init_ocean, init_1d_model, poisfft, production_e
  • palm/trunk/SOURCE/advec_particles.f90

    r119 r139  
    44! Actual revisions:
    55! -----------------
     6!
     7! TEST: PRINT statements on unit 9 (commented out)
     8!
     9! Former revisions:
     10! -----------------
     11! $Id$
     12!
     13! 119 2007-10-17 10:27:13Z raasch
    614! Sorting of particles is controlled by dt_sort_particles and moved from
    715! the SGS timestep loop after the end of this loop.
     
    917! numbers along y
    1018! Small bugfixes in the SGS part
    11 ! TEST: PRINT statements on unit 9 (commented out)
    12 !
    13 ! Former revisions:
    14 ! -----------------
    15 ! $Id$
    1619!
    1720! 106 2007-08-16 14:30:26Z raasch
  • palm/trunk/SOURCE/buoyancy.f90

    r132 r139  
    44! Actual revisions:
    55! -----------------
    6 ! Vertical scalar profiles now based on nzb_s_inner and ngp_2dh_s_inner.
     6!
    77!
    88! Former revisions:
    99! -----------------
    1010! $Id$
     11!
     12! 132 2007-11-20 09:46:11Z letzel
     13! Vertical scalar profiles now based on nzb_s_inner and ngp_2dh_s_inner.
    1114!
    1215! 106 2007-08-16 14:30:26Z raasch
  • palm/trunk/SOURCE/check_open.f90

    r120 r139  
    44! Actual revisions:
    55! -----------------
    6 ! Status of 3D-volume NetCDF data file only depends on switch netcdf_64bit_3d
     6!
    77!
    88! Former revisions:
    99! -----------------
    1010! $Id$
     11!
     12! 120 2007-10-17 11:54:43Z raasch
     13! Status of 3D-volume NetCDF data file only depends on switch netcdf_64bit_3d
    1114!
    1215! 105 2007-08-08 07:12:55Z raasch
  • palm/trunk/SOURCE/check_parameters.f90

    r138 r139  
    44! Actual revisions:
    55! -----------------
     6!
     7!
     8! Former revisions:
     9! -----------------
     10! $Id$
     11!
     12! 138 2007-11-28 10:03:58Z letzel
    613! Plant canopy added
    714! Allow new case bc_uv_t = 'dirichlet_0' for channel flow.
    815! Multigrid solver allows topography, checking of dt_sort_particles
    916! Bugfix: initializing u_init and v_init in case of ocean runs
    10 !
    11 ! Former revisions:
    12 ! -----------------
    13 ! $Id$
    1417!
    1518! 109 2007-08-28 15:26:47Z letzel
  • palm/trunk/SOURCE/data_output_dvrp.f90

    r130 r139  
    3232! Actual revisions:
    3333! -----------------
     34!
     35! TEST: different colours for isosurfaces
     36! TEST: write statements
     37!
     38! Former revisions:
     39! -----------------
     40! $Id$
     41!
     42! 130 2007-11-13 14:08:40Z letzel
    3443! allow two instead of one digit to specify isosurface and slicer variables
    3544! for unknown variables (CASE DEFAULT) call new subroutine
    3645! user_data_output_dvrp
    37 ! TEST: different colours for isosurfaces
    38 ! TEST: write statements
    39 !
    40 ! Former revisions:
    41 ! -----------------
    42 ! $Id$
    4346!
    4447! 82 2007-04-16 15:40:52Z raasch
  • palm/trunk/SOURCE/diffusion_s.f90

    r129 r139  
    44! Actual revisions:
    55! -----------------
    6 ! replace wall_heatflux by wall_s_flux that is now included in the parameter
    7 ! list, bugfix for assignment of fluxes at walls
     6!
    87!
    98! Former revisions:
    109! -----------------
    1110! $Id$
     11!
     12! 129 2007-10-30 12:12:24Z letzel
     13! replace wall_heatflux by wall_s_flux that is now included in the parameter
     14! list, bugfix for assignment of fluxes at walls
    1215!
    1316! 20 2007-02-26 00:12:32Z raasch
  • palm/trunk/SOURCE/diffusivities.f90

    r137 r139  
    44! Actual revisions:
    55! -----------------
    6 ! Bugfix for summation of sums_l_l for flow_statistics
    7 ! Vertical scalar profiles now based on nzb_s_inner and ngp_2dh_s_inner.
    86!
    97!
     
    119! -----------------
    1210! $Id$
     11!
     12! 137 2007-11-28 08:50:10Z letzel
     13! Bugfix for summation of sums_l_l for flow_statistics
     14! Vertical scalar profiles now based on nzb_s_inner and ngp_2dh_s_inner.
    1315!
    1416! 97 2007-06-21 08:23:15Z raasch
  • palm/trunk/SOURCE/flow_statistics.f90

    r133 r139  
    44! Actual revisions:
    55! -----------------
     6!
     7!
     8! Former revisions:
     9! -----------------
     10! $Id$
     11!
     12! 133 2007-11-20 10:10:53Z letzel
    613! Vertical profiles now based on nzb_s_inner; they are divided by
    714! ngp_2dh_s_inner (scalars, procucts of scalars) and ngp_2dh (staggered
    815! velocity components and their products, procucts of scalars and velocity
    916! components), respectively.
    10 !
    11 !
    12 ! Former revisions:
    13 ! -----------------
    14 ! $Id$
    1517!
    1618! 106 2007-08-16 14:30:26Z raasch
  • palm/trunk/SOURCE/header.f90

    r138 r139  
    44! Actual revisions:
    55! -----------------
     6!
     7!
     8! Former revisions:
     9! -----------------
     10! $Id$
     11!
     12! 138 2007-11-28 10:03:58Z letzel
    613! Allow new case bc_uv_t = 'dirichlet_0' for channel flow.
    714! Allow two instead of one digit to specify isosurface and slicer variables.
    815! Output of sorting frequency of particles
    9 !
    10 ! Former revisions:
    11 ! -----------------
    12 ! $Id$
    1316!
    1417! 108 2007-08-24 15:10:38Z letzel
  • palm/trunk/SOURCE/init_1d_model.f90

    r135 r139  
    44! Actual revisions:
    55! -----------------
    6 ! Bugfix: absolute value of f must be used when calculating the Blackadar
    7 ! mixing length!
     6!
    87!
    98! Former revisions:
    109! -----------------
    1110! $Id$
     11!
     12! 135 2007-11-22 12:24:23Z raasch
     13! Bugfix: absolute value of f must be used when calculating the Blackadar
     14! mixing length
    1215!
    1316! 82 2007-04-16 15:40:52Z raasch
  • palm/trunk/SOURCE/init_3d_model.f90

    r138 r139  
    77! Actual revisions:
    88! -----------------
     9!
     10!
     11! Former revisions:
     12! -----------------
     13! $Id$
     14!
     15! 138 2007-11-28 10:03:58Z letzel
    916! New counter ngp_2dh_s_inner.
    1017! Allow new case bc_uv_t = 'dirichlet_0' for channel flow.
    1118! Corrected calculation of initial volume flow for 'set_1d-model_profiles' and
    1219! 'set_constant_profiles' in case of buildings in the reference cross-sections.
    13 !
    14 ! Former revisions:
    15 ! -----------------
    16 ! $Id$
    1720!
    1821! 108 2007-08-24 15:10:38Z letzel
  • palm/trunk/SOURCE/init_dvrp.f90

    r130 r139  
    44! Actual revisions:
    55! -----------------
    6 ! allow two instead of one digit to specify isosurface and slicer variables
    7 ! Test output of isosurface on camera file
     6!
    87! TEST: print* statements
    98! ToDo: checking of mode_dvrp for legal values is not correct
     
    1211! -----------------
    1312! $Id$
     13!
     14! 130 2007-11-13 14:08:40Z letzel
     15! allow two instead of one digit to specify isosurface and slicer variables
     16! Test output of isosurface on camera file
    1417!
    1518! 82 2007-04-16 15:40:52Z raasch
  • palm/trunk/SOURCE/init_grid.f90

    r134 r139  
    44! Actual revisions:
    55! -----------------
     6!
     7!
     8! Former revisions:
     9! -----------------
     10! $Id$
     11!
     12! 134 2007-11-21 07:28:38Z letzel
    613! Redefine initial nzb_local as the actual total size of topography (later the
    714! extent of topography in nzb_local is reduced by 1dx at the E topography walls
     
    916! for consistency redefine 'single_building' case.
    1017! Calculation of wall flag arrays
    11 !
    12 ! Former revisions:
    13 ! -----------------
    14 ! $Id$
    1518!
    1619! 94 2007-06-01 15:25:22Z raasch
  • palm/trunk/SOURCE/init_ocean.f90

    r124 r139  
    44! Actual revisions:
    55! -----------------
    6 ! Bugfix: Initial density rho is calculated
     6!
    77!
    88! Former revisions:
    99! ------------------
    1010! $Id$
     11!
     12! 124 2007-10-19 15:47:46Z raasch
     13! Bugfix: Initial density rho is calculated
    1114!
    1215! 97 2007-06-21 08:23:15Z raasch
  • palm/trunk/SOURCE/init_particles.f90

    r117 r139  
    44! Actual revisions:
    55! -----------------
    6 ! Sorting of particles only in case of cloud droplets
     6!
    77!
    88! Former revisions:
    99! -----------------
    1010! $Id$
     11!
     12! 117 2007-10-11 03:27:59Z raasch
     13! Sorting of particles only in case of cloud droplets
    1114!
    1215! 106 2007-08-16 14:30:26Z raasch
  • palm/trunk/SOURCE/init_pegrid.f90

    r114 r139  
    44! Actual revisions:
    55! -----------------
    6 ! Allocation of wall flag arrays for multigrid solver
     6!
    77! TEST OUTPUT (TO BE REMOVED) logging mpi2 ierr values
    88!
     
    1010! -----------------
    1111! $Id$
     12!
     13! 114 2007-10-10 00:03:15Z raasch
     14! Allocation of wall flag arrays for multigrid solver
    1215!
    1316! 108 2007-08-24 15:10:38Z letzel
  • palm/trunk/SOURCE/modules.f90

    r138 r139  
    55! Actual revisions:
    66! -----------------
     7!
     8!
     9! Former revisions:
     10! -----------------
     11! $Id$
     12!
     13! 138 2007-11-28 10:03:58Z letzel
    714! +drag_coefficient, pch_index, lad_surface, lad_vertical_gradient,
    815! lad_vertical_gradient_level, plant_canopy, lad, lad_s, lad_u, lad_v,
    916! lad_w, cdc, lad_vertical_gradient_level_ind, canopy_mode
    10 !
    1117! +dt_sort_particles, ngp_2dh_s_inner, time_sort_particles, flags,
    1218! wall_flags_1..10, wall_humidityflux(0:4), wall_qflux(0:4),
    1319! wall_salinityflux(0:4), wall_scalarflux(0:4)
    14 !
    15 ! Former revisions:
    16 ! -----------------
    17 ! $Id$
    1820!
    1921! 108 2007-08-24 15:10:38Z letzel
  • palm/trunk/SOURCE/package_parin.f90

    r116 r139  
    99! -----------------
    1010! $Id$
     11!
     12! 116 2007-10-11 02:30:27Z raasch
     13! +dt_sort_particles in package_parin
    1114!
    1215! 60 2007-03-11 11:50:04Z raasch
  • palm/trunk/SOURCE/parin.f90

    r138 r139  
    44! Actual revisions:
    55! -----------------
    6 ! +canopy_mode
    7 ! +drag_coefficient, lad_surface, lad_vertical_gradient,
    8 ! +lad_vertical_gradient_level, pch_index, plant_canopy
    9 ! +allocation of leaf area density field
    10 !
     6!
    117!
    128! Former revisions:
    139! -----------------
    1410! $Id$
     11!
     12! 138 2007-11-28 10:03:58Z letzel
     13! +canopy_mode, drag_coefficient, lad_surface, lad_vertical_gradient,
     14! lad_vertical_gradient_level, pch_index, plant_canopy,
     15! +allocation of leaf area density field
    1516!
    1617! 108 2007-08-24 15:10:38Z letzel
  • palm/trunk/SOURCE/plant_canopy_model.f90

    r138 r139  
    44! Actual revisions:
    55! -----------------
    6 ! Initial revision
     6!
    77!
    88! Former revisions:
    99! -----------------
    1010! $Id$
     11!
     12! 138 2007-11-28 10:03:58Z letzel
     13! Initial revision
    1114!
    1215! Description:
  • palm/trunk/SOURCE/poisfft.f90

    r128 r139  
    44! Actual revisions:
    55! -----------------
    6 ! Bugfix: wavenumber calculation for even nx in routines maketri
     6!
    77!
    88! Former revisions:
    99! -----------------
    1010! $Id$
     11!
     12! 128 2007-10-26 13:11:14Z raasch
     13! Bugfix: wavenumber calculation for even nx in routines maketri
    1114!
    1215! 85 2007-05-11 09:35:14Z raasch
  • palm/trunk/SOURCE/poismg.f90

    r114 r139  
    88! Actual revisions:
    99! -----------------
     10!
     11!
     12! Former revisions:
     13! -----------------
     14! $Id$
     15!
     16! 114 2007-10-10 00:03:15Z raasch
    1017! Boundary conditions at walls are implicitly set using flag arrays. Only
    1118! Neumann BC is allowed. Upper walls are still not realized.
    1219! Bottom and top BCs for array f_mg in restrict removed because boundary
    13 ! values are not needed (right hand side of SOR iteration)
    14 !
    15 ! Former revisions:
    16 ! -----------------
    17 ! $Id$
     20! values are not needed (right hand side of SOR iteration).
    1821!
    1922! 75 2007-03-22 09:54:05Z raasch
  • palm/trunk/SOURCE/production_e.f90

    r124 r139  
    44! Actual revisions:
    55! -----------------
    6 ! Bugfix: calculation of density flux in the ocean now starts from nzb+1
     6!
    77!
    88! Former revisions:
    99! -----------------
    1010! $Id$
     11!
     12! 124 2007-10-19 15:47:46Z raasch
     13! Bugfix: calculation of density flux in the ocean now starts from nzb+1
    1114!
    1215! 108 2007-08-24 15:10:38Z letzel
  • palm/trunk/SOURCE/prognostic_equations.f90

    r138 r139  
    44! Actual revisions:
    55! -----------------
    6 ! add call of subroutines that evaluate the canopy drag terms
    7 ! add wall_*flux to parameter list of calls of diffusion_s
     6!
    87!
    98! Former revisions:
    109! -----------------
    1110! $Id$
     11!
     12! 138 2007-11-28 10:03:58Z letzel
     13! add call of subroutines that evaluate the canopy drag terms,
     14! add wall_*flux to parameter list of calls of diffusion_s
    1215!
    1316! 106 2007-08-16 14:30:26Z raasch
  • palm/trunk/SOURCE/read_var_list.f90

    r138 r139  
    44! Actual revisions:
    55! -----------------
    6 ! +canopy_mode
    7 ! +drag_coefficient, lad, lad_surface, lad_vertical_gradient,
    8 ! lad_vertical_gradient_level, lad_vertical_gradient_level_ind,
    9 ! pch_index, plant_canopy
    10 !
    11 ! +time_sort_particles
     6!
    127!
    138! Former revisions:
    149! -----------------
    1510! $Id$
     11!
     12! 138 2007-11-28 10:03:58Z letzel
     13! +canopy_mode, drag_coefficient, lad, lad_surface, lad_vertical_gradient,
     14! lad_vertical_gradient_level, lad_vertical_gradient_level_ind, pch_index,
     15! plant_canopy, time_sort_particles
    1616!
    1717! 102 2007-07-27 09:09:17Z raasch
  • palm/trunk/SOURCE/user_interface.f90

    r138 r139  
    44! Actual revisions:
    55! -----------------
    6 ! new subroutine user_init_plant_canopy
    7 !
    8 ! new subroutine user_data_output_dvrp
    9 ! +argument gls in user_init_grid
     6!
    107!
    118! Former revisions:
    129! -----------------
    1310! $Id$
     11!
     12! 138 2007-11-28 10:03:58Z letzel
     13! new subroutines user_init_plant_canopy, user_data_output_dvrp
     14! +argument gls in user_init_grid
    1415!
    1516! 105 2007-08-08 07:12:55Z raasch
  • palm/trunk/SOURCE/write_var_list.f90

    r138 r139  
    44! Actual revisions:
    55! -----------------
    6 ! +canopy_mode
    7 ! +drag_coefficient, lad, lad_surface, lad_vertical_gradient,
    8 ! lad_vertical_gradient_level, lad_vertical_gradient_level_ind,
    9 ! pch_index, plant_canopy
    10 !
    11 ! +time_sort_particles
     6!
    127!
    138! Former revisions:
    149! -----------------
    1510! $Id$
     11!
     12! 138 2007-11-28 10:03:58Z letzel
     13! +canopy_mode, drag_coefficient, lad, lad_surface, lad_vertical_gradient,
     14! lad_vertical_gradient_level, lad_vertical_gradient_level_ind, pch_index,
     15! plant_canopy, time_sort_particles
    1616!
    1717! 102 2007-07-27 09:09:17Z raasch
Note: See TracChangeset for help on using the changeset viewer.